大慶水庫冰表面溫度特征及其在冰工程實踐應(yīng)用方法探討

李潤玲，白乙拉，榮 慧，賈 青，李志軍

（1.大連理工大學(xué) 海岸和近海工程國家重點實驗室，遼寧 大連 116024；2.渤海大學(xué) 數(shù)理學(xué)院，遼寧 錦州 121000；3.黑龍江大學(xué) 水利電力學(xué)院，黑龍江 哈爾濱 150086）

大慶水庫冰表面溫度特征及其在冰工程實踐應(yīng)用方法探討

李潤玲1，白乙拉2，榮 慧2，賈 青3，李志軍1

（1.大連理工大學(xué) 海岸和近海工程國家重點實驗室，遼寧 大連 116024；2.渤海大學(xué) 數(shù)理學(xué)院，遼寧 錦州 121000；3.黑龍江大學(xué) 水利電力學(xué)院，黑龍江 哈爾濱 150086）

根據(jù)冰面熱平衡方程并以大慶紅旗泡水庫冰面觀測的氣溫、太陽輻射、風(fēng)速、相對濕度、冰厚數(shù)據(jù)以及水庫附近安達(dá)氣象站記錄的云量、濕度和海平面氣壓數(shù)據(jù)，求解了2009/2010年度冬季的冰表面溫度，并分析了冰表面溫度特征及其同太陽輻射之間的關(guān)系、同氣溫、風(fēng)速之間的響應(yīng)關(guān)系。在此基礎(chǔ)上，就傳統(tǒng)考慮的抗冰工程設(shè)計所需的設(shè)計冰表面溫度進(jìn)行了評估；就日益增加的安全用冰工程設(shè)計所需冰表面溫度取值也給予定量討論；對氣候變化下的冰內(nèi)和冰下水體生態(tài)適應(yīng)問題所需的冰面溫度日內(nèi)變化速率和開河期冰凌爆破所需冰表面溫度的取值也進(jìn)行探索。盡管受國際氣象資料觀測規(guī)定的限制，歷史資料的應(yīng)用存在利用接近觀測時刻資料的事實，但仍然期望這些分析和討論能夠推動科學(xué)地、更廣闊地應(yīng)用冰表面溫度。

淡水冰；表面溫度；氣溫；冰工程；冰生態(tài)

1 研究背景

隨著氣候變暖，地表冰發(fā)生著諸多形式的變化，但國人依賴于冰的活動卻有增無減。除了眾所周知的冰凌對固定式或浮式結(jié)構(gòu)物（如渤海結(jié)冰區(qū)的海上平臺［1］、沿岸港口［2］、電廠［3］；冰區(qū)航行的船舶［4］）的作用外，突發(fā)性結(jié)冰給渤海、北方冬季水產(chǎn)養(yǎng)殖業(yè)也帶來了重大經(jīng)濟(jì)損失［5］。另外，越來越多的國人喜好冰上娛樂活動和冰上交通，它們也增加了不穩(wěn)定冰期時的安全風(fēng)險［6］。這意味著在冰凌對水工和海洋結(jié)構(gòu)物破壞的防御問題之外，又增加了冬季冰面交通、垂釣和其它娛樂的冰上活動安全問題；以及渤海和北方冰下水體中的水產(chǎn)養(yǎng)殖安全越冬問題。這些將冰應(yīng)用領(lǐng)域由傳統(tǒng)的冰防范拓展到了安全利用冰以及冰內(nèi)及冰下水內(nèi)生態(tài)等新領(lǐng)域。冰內(nèi)溫度取決于冰表面溫度和底面溫度，一般低氣溫和低冰溫發(fā)生在沒有太陽輻射和冰面放熱的凌晨［7］；而高氣溫、高冰溫發(fā)生在凈輻射最高的午后。另外，極端氣象事件的頻繁發(fā)生，導(dǎo)致突變的風(fēng)、氣溫的不規(guī)律變化，也成為不可忽視的因子。對于以上各種情況，利用冰面熱平衡方程計算包含輻射、風(fēng)速、氣溫、冰厚在內(nèi)的冰表面溫度是最可靠的。然后再針對工程應(yīng)用的不同特殊條件，探索利用冰面熱平衡方程的計算結(jié)果建立起簡化的冰表面、冰內(nèi)極端溫度表達(dá)方式。本文依照上述思想，結(jié)合在黑龍江大慶紅旗泡水庫2009—2010年的現(xiàn)場觀測氣溫、輻射、風(fēng)速、冰厚，以及距水庫較近的安達(dá)氣象站有關(guān)云量和海平面氣壓數(shù)據(jù)，模擬計算出2009.12.1—2010.2.10水庫冰表面溫度時間系列。在此基礎(chǔ)上，探討冰表面溫度變化特征和最低、最高冰表面溫度的取值方法，和如何依此作為確定抗冰工程、安全用冰工程、冰下環(huán)境和生態(tài)、冰凌爆破的科學(xué)依據(jù)。

2 表面熱平衡方程及其物理量的計算方法

有關(guān)靜態(tài)水體（湖泊、水庫、水塘、養(yǎng)殖池）的冰生消過程有許多研究，其中冰內(nèi)溫度就是研究的一部分內(nèi)容。國際上有些模式簡化了湖水混合層熱力學(xué)的過程，并將其與日常氣象觀測數(shù)據(jù)和氣候預(yù)報模式結(jié)合，模擬出湖水全年的熱力學(xué)過程及對氣候的影響［8-9］。2012年楊宇等以雪/冰表面溫度和氣溫變化之間的時滯為辨識參量；以計算的雪/冰表面溫度和觀測皮溫之差值為目標(biāo)函數(shù)，建立了最優(yōu)辨識模型，獲得芬蘭湖冰雪/冰表面溫度和氣溫之間時滯的最優(yōu)變化過程［10］。這套計算方法后被引用到計算中國東北冰的表面溫度。中國北方干燥，冬季降雪較少。冰表面熱平衡方程的一般表達(dá)式是：

式中：α為冰表面反照率；QS為入射短波輻射；I0為穿過表層進(jìn)入冰內(nèi)的短波輻射；Tsfc為冰表面溫度；Qd、Qb分別為入射長波輻射和反射長波輻射；Qh為感熱通量；Ql為潛熱通量；Fc為由冰內(nèi)傳遞到冰面的熱通量；Fm為冰面融化的熱通量。所有通量均是以指向冰表面方向為正。計算表面熱通量和表面溫度需要的氣象參數(shù)有風(fēng)速（Va）、氣溫（Ta）、相對濕度（r）、云量（CN）。圍繞式（1），需要通過式（2）—式（11）計算所需物理參量。

反照率計算形式為：

式中：F+為反輻射；F-為總輻射。

選取Efimova的入射長波輻射參數(shù)化方案和Jocobs給出的云量對長波輻射的影響，作為入射長波輻射（Qd）參數(shù)化方案；并根據(jù)Stefan-Boltzmann輻射定律，反射長波輻射（Qb）是表面溫度的函數(shù)，則：

式中：σs為Stefan-Boltzmann常數(shù)，5.68×10-8；e1為表面水蒸氣壓力；ε為表面輻射率；CN為總云量。

式中：r為相對濕度；e為飽和水氣壓。

利用觀測的氣溫數(shù)據(jù)可以計算出飽和水氣壓：

式中Ta為氣溫。

湍流通量包括感熱通量和潛熱通量，它們分別由式（7）和式（8）計算：

式中：ρa(bǔ)為空氣密度；ca為空氣比熱；Rl為蒸發(fā)焓；CH和CE分別為感熱交換系數(shù)和潛熱交換系數(shù)；Va為風(fēng)速；qa為空氣的比濕；qsfc為冰表面比濕。

計算比濕的公式為：

式中：p0為海平面氣壓，分別用Ta、Tsfc代入式（6）和式（9）進(jìn)而計算出qa、qsfc。

冰表面熱傳導(dǎo)通量為：

式中：ki為冰的熱傳導(dǎo)系數(shù)；tf為冰點；h為冰厚度，根據(jù)觀測的冰厚數(shù)據(jù)可以計算出熱傳導(dǎo)通量。由于冬季結(jié)冰期的氣溫很低，認(rèn)為不會發(fā)生冰表面融化。因此，假設(shè)：m

3 基于大慶水庫冬季數(shù)據(jù)的冰表面溫度

求解冰表面熱平衡方程需要氣象和冰觀測數(shù)據(jù)。2008年12月26日至2009年4月9日和2009年11月26日至2010年4月17日，在紅旗泡水庫庫區(qū)距離南站約150m的冰面上安裝自動氣象站［11］，并開展冰物理現(xiàn)場調(diào)查。本文選用的數(shù)據(jù)有：2009/2010年度冰上氣溫、風(fēng)速、太陽總輻射和太陽反輻射；近冰表面及冰層內(nèi)部溫度、冰層厚度。另外，從中國氣象數(shù)據(jù)共享網(wǎng)下載離紅旗泡水庫較近的安達(dá)氣象站云量、海平面氣壓、相對濕度數(shù)據(jù)。

由于中國氣象數(shù)據(jù)共享網(wǎng)上提供的云量、海平面氣壓、相對濕度的時間間隔為6 h；現(xiàn)場觀測氣溫、風(fēng)速、太陽總輻射和太陽反輻射的時間間隔為1 m in；現(xiàn)場冰層厚度記錄時間間隔為15m in，為此需要建立一個對應(yīng)30m in，處在半點的數(shù)據(jù)集。對中國氣象數(shù)據(jù)共享網(wǎng)的數(shù)據(jù)，利用三次樣條插值的方法，把每日00∶00、06∶00、12∶00、18∶00四個時間點的數(shù)據(jù)，插值成每30min的數(shù)據(jù)。把現(xiàn)場觀測數(shù)據(jù)，挑選出對應(yīng)的每30m in數(shù)據(jù)。從2009年12月1日到2010年2月10日，共有2 814組數(shù)據(jù)。求解大慶水庫冰表層溫度需要的其他參數(shù)，根據(jù)文獻(xiàn)收集并匯總于表1。

選取正割迭代法對冰表面熱平衡方程求解，得到對應(yīng)的2 814組數(shù)據(jù)所在時刻的冰表面溫度。其中圖1選取了2009年12月1日至2009年12月31日1 249組輻射、氣溫、冰表面溫度，給出它們的時間變化曲線。

由圖1可以發(fā)現(xiàn)：隨著氣溫的波動，冰表面溫度也隨之波動；冰表面溫度、氣溫和凈輻射之間存在明顯的周期性變化，變化周期為一日。其中凈輻射的周期規(guī)律較為突出，其次是冰表面溫度和氣溫。三者之間的峰值發(fā)生時間的一般規(guī)律是凈輻射峰值出現(xiàn)時間早于冰表面溫度；冰表面溫度早于氣溫。也就是說，太陽輻射是冰表面溫度的控制因子，它的強(qiáng)烈輻射引起冰表面溫度的“滯后”響應(yīng)。氣溫同樣控制冰表面溫度，但受凈輻射和風(fēng)的影響，冰表面溫度對輻射和風(fēng)的“響應(yīng)”敏捷于氣溫。因此，冰表面溫度稍早于氣溫，兩者的時間差相對較小，而輻射的峰值相對于冰表面溫度較早出現(xiàn)。同樣，圖1還可以分析谷值冰表面溫度和谷值氣溫的特征。所獲得的谷值冰表面溫度和谷值氣溫之間的時差。而凌晨谷值氣溫發(fā)生的時間集中在凌晨7∶30，見圖2。兩者相差1.5 h；而太陽輻射的平均峰值發(fā)生在12∶00，每日平均最高氣溫則發(fā)生在14∶00，滯后2.0 h。在16∶00至次日7∶00期間，太陽輻射基本為零。

表1 計算大慶紅旗泡水庫冰表面溫度的相關(guān)參數(shù)取值匯總表

圖1 2009年12月冰面觀測氣溫、凈輻射和計算冰表面溫度歷時過程

4 冰表面溫度與冰溫關(guān)系及其在冰工程實踐應(yīng)用的方法商榷

從冰工程角度出發(fā)，冰基礎(chǔ)科學(xué)參數(shù)有冰厚度和冰強(qiáng)度。而冰強(qiáng)度的主要控制因子是冰溫度，因此提供冰工程活動科學(xué)參數(shù)，就需要給出冰的溫度［12-13］。

在特定溫度環(huán)境中，冰內(nèi)的熱質(zhì)趨于平衡。如果冰內(nèi)部不存在熱源，則冰內(nèi)溫度呈線性分布，但在冰表面升溫和降溫中，冰內(nèi)熱傳導(dǎo)引起冰內(nèi)溫度出現(xiàn)“熱中間層”或者“冷中間層”［14］。冰工程問題，往往考慮的是極限環(huán)境，因此假設(shè)冰溫度在冰內(nèi)線性分布，然后采用分層處理［12］或者直接處理［7］。

分層處理由Timco和Frederking（1990）提出，并應(yīng)用在冰層壓縮強(qiáng)度評估上［12］。其處理思想在渤海引用［7］，具體步驟是：通過氣溫、水溫、冰厚和冰期將平整冰層分為若干小層，確定每小層的冰溫度、鹽度、密度，計算該層的孔隙率，然后計算對應(yīng)力學(xué)參數(shù)。取各小層力學(xué)參數(shù)的平均值作為平整冰層的力學(xué)指標(biāo)。文獻(xiàn)［12］指出把平整冰層分為9層獲得的南極平整冰層壓縮強(qiáng)度具有較高精度，并且計算的強(qiáng)度值收斂。假設(shè)平整冰層溫度剖面每一小層的平均溫度均由表面至底面呈線性分布，第n層的冰溫θn（n=1代表最表層）由式（12）計算：

圖2 最低氣溫凌晨出現(xiàn)時刻統(tǒng)計

式中：θB為平整冰層底面溫度，等于水溫；θs為冰層夜間表面溫度。

4.1 在抗冰工程中的應(yīng)用在抗冰工程中，采用冰作用力上限，即設(shè)計中需要考慮冰強(qiáng)度高的冰凌上限作用力［13］。根據(jù)圖1和式（1）可知，冰表面溫度最低情況出現(xiàn)在太陽輻射為零、風(fēng)速最大的環(huán)境條件。為了覆蓋在冬季出現(xiàn)的這些條件，分析使用了2009年12月1日至2010年2月10日的氣象站規(guī)定觀測時間的觀測數(shù)據(jù)和計算數(shù)據(jù)，共計258組。這些數(shù)據(jù)既覆蓋了整個冬季的觀測時間，又不失去分析問題的能力。為了獲得最低溫度，圖3選擇了凈輻射接近0的18∶00，0∶00和6∶00的數(shù)據(jù)并獲得冰表面溫度和氣溫之間的關(guān)系。從圖3可以發(fā)現(xiàn)輻射影響并不存在，但風(fēng)速的作用仍然存在，即圖3中數(shù)據(jù)點的分散原因。而且風(fēng)速越大，表面溫度越接近氣溫。為此對這批輻射接近0的數(shù)據(jù)（192組）進(jìn)行氣溫、風(fēng)速對冰表面溫度影響的統(tǒng)計分析，得到式（13）。

根據(jù)冬季大慶水庫記錄數(shù)據(jù)出現(xiàn)的最大風(fēng)速為10m/s，將其代入式（13），得到圖3的虛線。該虛線即為抗冰工程設(shè)計時需要的冰表面設(shè)計溫度。它所對應(yīng)的冰內(nèi)溫度應(yīng)該也較低，相應(yīng)的冰強(qiáng)度高。在文獻(xiàn)［7］中，采用了抗冰結(jié)構(gòu)物設(shè)計冰溫取凌晨氣溫等于冰表面溫度的估算，這種作法相當(dāng)于認(rèn)為在無輻射且風(fēng)速最大條件下，最低的表層氣溫等于氣溫的情況。但事實上，根據(jù)圖3，在0～-13℃的時候，這個作法偏于冒險；在-13～-23℃時較好；而在低于-23℃時偏于保守。也就是說，對于冬季氣溫較高的結(jié)冰區(qū)域，這個估算做法具有風(fēng)險。

文獻(xiàn)［13］指出對于橋墩上的冰作用力，主要是春季流凌擠壓作用力，它與每年3～4月份開河時的冰溫有關(guān)，但是每年開河時的冰溫分布有差異，需要選擇數(shù)年觀測數(shù)據(jù)中最低冰溫為依據(jù)。事實上，為了抗冰結(jié)構(gòu)物設(shè)計，要根據(jù)結(jié)構(gòu)物重現(xiàn)期，分析歷史上的開河日，在統(tǒng)計開河日的基礎(chǔ)上，統(tǒng)計這段時間的氣溫，用統(tǒng)計方法確定相同重現(xiàn)期的氣溫，再進(jìn)而確定相同重現(xiàn)期的冰溫，用冰力學(xué)參數(shù)同冰溫的關(guān)系，確定冰力學(xué)參數(shù)［7］。

圖3 2009/2010年度冬季夜間冰表面溫和氣溫之間關(guān)系

4.2 在安全用冰工程中的應(yīng)用冰工程的另一種情況是利用冰，如冰上交通，冰上娛樂等，它們將冰作為承載“地基”。因此，只有選擇冰溫高、強(qiáng)度低的下限冰承載力，才能保證它使用期間的安全［15］。這時應(yīng)該取凈輻射最大，風(fēng)速接近0m/s時的冰表面溫度作為設(shè)計冰溫依據(jù)，然后再評估冰的強(qiáng)度。因此，利用2009年12月1日至2010年2月10日中午12點的觀測和計算數(shù)據(jù)繪制圖4。圖中實心點是夜間無輻射時的計算冰表面溫度，它受到風(fēng)速的影響；圖中空心點是中午12點時的數(shù)據(jù)，它既受到風(fēng)速的影響，又受到輻射的影響。根據(jù)大慶觀測數(shù)據(jù)，將風(fēng)速低于2m/s的視為無風(fēng)，實測最大凈輻射為180W/m2，分析得到氣溫低于-17℃的實線段。另外，計算冰表面溫度中，有些受到感熱影響，計算值大于冰點溫度，即高于淡水冰的冰點溫度0℃，這屬于物理上不合理的數(shù)據(jù)。這種情況下，設(shè)計溫度取冰點溫度0℃。因此，在安全用冰工程中的設(shè)計冰表面溫度由圖4的二段實線構(gòu)成。

事實上，冰上活動多數(shù)處于穩(wěn)定盛冰期，而防凌處于不穩(wěn)定的流冰期。一般冰上活動時段的氣溫低、冰溫低，但是冰上活動往往是臨時性的，沒有重現(xiàn)期的概念。在黑龍江的冰道承載力估算中，只使用當(dāng)日和前三日的氣溫資料評估［16］。這些評估也沒有將輻射考慮在內(nèi)，2011年中國學(xué)者才提出考慮靜態(tài)和動態(tài)水體冰生消中的太陽輻射影響［17］。它雖然沒有分解出對冰表面溫度的貢獻(xiàn)，但考慮了對冰厚度的貢獻(xiàn)。另外，關(guān)于河流冰下流速對冰厚度的貢獻(xiàn)也有半經(jīng)驗公式［18］，這些為推動冰道承載力的精確評估提供開端。

4.3 在結(jié)冰環(huán)境下生態(tài)中的應(yīng)用對于冰內(nèi)和冰下水環(huán)境和生態(tài)而言，既需要考慮極端低冰溫下，也需要考慮極端高冰溫下的冰內(nèi)和水體溫度、溶解氧、太陽穿透冰層的輻照度。在極地惡劣環(huán)境中的生態(tài)問題已有努力［19］，在中國北方水產(chǎn)養(yǎng)殖方面也有經(jīng)驗報道［20-21］。因此，一些學(xué)者啟動了低溫環(huán)境模擬冰下水環(huán)境條件的日變化幅度對控制水生生物的適應(yīng)性能力和響應(yīng)方面的科學(xué)試驗。

由于冰內(nèi)溫度在二個代表性極端環(huán)境的上下限區(qū)域內(nèi)（圖4的實線和虛線）變動，結(jié)冰后，冰內(nèi)和冰下水溫出現(xiàn)24 h周期變化，這種變化影響到水生動物的適應(yīng)性，因此不僅冰溫是影響因子，而且冰溫的變化速率更是重要。圖5在圖4的基礎(chǔ)上，增加依據(jù)這二條極端冰表層溫度上下限曲線所計算的冰內(nèi)溫差變化率。

圖4 2009/2010年度冬季冰表面溫度同氣溫之間關(guān)系

圖5 2009/2010年度冬季冰表面最低、最高溫度和冰內(nèi)溫差日變化率同氣溫之間關(guān)系

依此希望圖5成為指導(dǎo)冰區(qū)和冰下生態(tài)環(huán)境研究的思路之一，但其根本思想是作為生態(tài)適應(yīng)方面的溫度和溫度變化率效應(yīng)的基礎(chǔ)。在這種思路下，溫度調(diào)查是基礎(chǔ)。這也可以同中國學(xué)者近幾年已經(jīng)積累的結(jié)冰環(huán)境中的冰下水溫和冰下水質(zhì)環(huán)境成果［22-23］，開展對接嘗試。

另一個與冰溫度變化率有關(guān)的應(yīng)用是春季冰對水庫護(hù)坡的熱膨脹作用力問題［24］。它如同抗冰問題，是發(fā)生在冰層升溫階段的現(xiàn)象，而生態(tài)問題是發(fā)生在整個冰期，特別是穩(wěn)定的盛冰期。

4.4 在冰凌爆破中的應(yīng)用冰凌爆破屬于黃河防凌的傳統(tǒng)措施之一，直到21世紀(jì)初，冰凌爆破再次“升溫”，其中除了發(fā)展新型火力外，還將爆破效果同冰力學(xué)性質(zhì)聯(lián)系起來［25］，因此爆破效果也同冰溫聯(lián)系起來。

對于冰凌爆破，它類同于冰管理，與實施爆破當(dāng)時冰性質(zhì)有關(guān)。但實踐中，有幾條事實考慮，使得冰表面溫度得以簡化。首先，盡管實施爆破可以不受時間限制，但實際操作在白天，也就是存在太陽輻射對冰表面溫度影響。在白天時間，無論有風(fēng)和無風(fēng)都可以實施，因此白天最大風(fēng)速和最小風(fēng)速對應(yīng)的冰表面溫度決定了變化范圍。另外，冰凌爆破是開河期間出現(xiàn)險情時的應(yīng)急措施，這就規(guī)定了當(dāng)時的氣溫不會很低。圖6中的實線對應(yīng)午間無風(fēng)時的最強(qiáng)輻射影響下的冰表面溫度；虛線則對應(yīng)風(fēng)速最大且無輻射時的冰表面溫度。三角形的最低溫度直線對應(yīng)氣溫-10℃，該溫度是發(fā)生開河時的最低氣溫。它們所構(gòu)成的三角形內(nèi)部范圍是冰凌爆破的關(guān)心區(qū)。

冰凌爆破的冰溫評價則建議直接依據(jù)冰表面熱平衡方程計算的三角形區(qū)域內(nèi)的冰表面溫度結(jié)果。由于冰底面溫度和冰厚度也是影響爆破效果的因子，同時建議未來考慮冰下水流的貢獻(xiàn)。

圖6 2009/2010年度冬季冰表面白天最低、最高溫度同氣溫之間關(guān)系以及冰凌爆破應(yīng)用區(qū)

5 討論與結(jié)論

根據(jù)連續(xù)記錄數(shù)據(jù)分析得到一些基本規(guī)律，但在工程實踐中目前存在資料使用的限制。歷史氣象數(shù)據(jù)記錄采用國際準(zhǔn)則，即北京時間0∶00，6∶00，12∶00，18∶00四次觀測。即使是采用連續(xù)記錄數(shù)據(jù)，一般公開的仍然是這四次數(shù)據(jù)。這些時刻不對應(yīng)分析得到的輻射和氣溫極值時刻，如果需要考慮重現(xiàn)期，就必須使用歷史數(shù)據(jù)，無法回避這一事實，只好選擇與四次觀測時間相近的數(shù)據(jù)統(tǒng)計分析。因此，在未來，需要思考二件工作，一是增加自動記錄設(shè)備投入；二是建立修正歷史數(shù)據(jù)的統(tǒng)計關(guān)系。

由于淡水冰封凍后徹底失去氣-水交換，盡管淡水結(jié)冰區(qū)最低氣溫主要發(fā)生在凌晨7∶30，使用凌晨6∶00觀測氣溫作為冰面溫度存在局限性。因此，針對渤海工程海冰提出的利用凌晨氣溫作為冰表面溫度對于結(jié)冰淡水水庫、水池，河流不完全適宜。盡量利用連續(xù)觀測資料建立起不同時間之間氣溫分布統(tǒng)計關(guān)系，然后再利用歷史6∶00資料并給予修正。

抗冰工程需要對應(yīng)時間的最低冰表面溫度，即沒有輻射和風(fēng)速最大條件下對應(yīng)冰表面溫度狀況，重點關(guān)心凌汛期。另外，如果抗冰結(jié)構(gòu)物設(shè)計需要重現(xiàn)期，凌汛期的氣溫也需要考慮到重現(xiàn)期。對于春季的冰熱膨脹引起的冰對結(jié)構(gòu)物作用，必須使用溫度變化率。安全用冰工程需要考慮應(yīng)用期間的最高冰表面溫度，它對應(yīng)輻射最大和風(fēng)速最小條件下的冰表面溫度，重點關(guān)心盛冰期；冰內(nèi)和冰下生態(tài)需要的是最高和最低溫差以及溫差變化率，重點關(guān)注盛冰期；冰凌爆破則需要白天有風(fēng)和無風(fēng)環(huán)境條件下的冰表面溫度，只需要關(guān)注流凌期。

首先利用冰表面熱平衡方程計算數(shù)據(jù)，然后簡化給出針對不同涉冰領(lǐng)域的冰表面溫度同氣溫的關(guān)系。這些關(guān)系可望成為未來將氣候變化下的氣溫變化引入到實踐應(yīng)用中的基礎(chǔ)。

［1］ 周振濤，苗文東.談海洋石油防抗海冰［J］.中國海事，2011（2）：23-25.

［2］ 付博新，宋向群，郭子堅，等.海冰對港口作業(yè)的影響及應(yīng)對措施［J］.水道港口，2007，28（6）：444-447.

［3］ 賀益英.寒區(qū)工業(yè)取水口防冰工程新方法的研究［J］.水利學(xué)報，2003（1）：7-11.

［4］ 劉強(qiáng)，孫健，王鳳武 .黃、渤海海冰對船舶通航安全的影響及應(yīng)對措施［J］.世界海運，2012，35（6）：30-33.

［5］ 孫劭，蘇潔，史培軍.2010年渤海海冰災(zāi)害特征分析［J］.自然災(zāi)害學(xué)報，2011，20（6）：87-93.

［6］ 姜建春.冰河事故救援處置對策研究［J］.武警學(xué)院學(xué)報，2014，30（8）：33-35.

［7］ 李志軍，Sodhi Devinder S，盧鵬.渤海海冰工程設(shè)計參數(shù)分布［J］.工程力學(xué)，2006，23（6）：167-172.

［8］ Mironov D，Heise E，Kourzeneva E，etal.Implementation of the lake parameterization scheme Flake into the nu?mericalweather predictionmodel COSMO［J］.Boreal EnvironmentResearch，2010，15（2）：218-230.

［9］ Kourzeneva E.External data for lake parameterization in numericalweather prediction and climatemodeling［J］. Boreal EnvironmentResearch，2010，15（2）：165-177.

［10］ Yang Yu，Li Zhijun，Shi Liqiong，etal.An iterative interpolation algorithm formodeling snow/ice surface and in?ner temperature［J］.International Journal of Advancements in Computing Technology，2012，4（12）：195-204.

［11］ 李志軍，楊宇，彭旭明，等.黑龍江紅旗泡水庫冰生長過程現(xiàn)場觀測數(shù)據(jù)的剖析［J］.西安理工大學(xué)學(xué)報，2009，25（3）：270-274.

［12］ Timco GW，F(xiàn)rederking，RMW.Compressive strength of sea ice sheets［J］.Cold Regions Science and Technolo?gy，1990，17（3）：227-240.

［13］ 于天來，雷俊卿，單思迪 .計算春季流冰對橋墩動壓力時冰溫取值的試驗研究［J］.中外公路，2010，30（4）：220-223.

［14］ 李志軍，隋吉學(xué)，嚴(yán)德成，等.遼東灣平整固定冰冰溫及其它物理特性的測定與分析［J］.海洋學(xué)報：中文版，1989，11（4）：525-533.

［15］ 徐伯孟，于生清，盧興良.嫩江冰上交通［J］.冰川凍土，1992，14（4）：375-380.

［16］ 鄭榮，賈俊明，劉臣君，等.黑龍江冰道承載力測算［J］.黑河科技，1996（4）：27-28.

［17］ 練繼建，趙新.靜動水冰厚生長消融全過程的輻射冰凍度-日法預(yù)測研究［J］.水利學(xué)報，2011，42（11）：1261-1267.

［18］ 李志軍，孫萬光，許士國，等 .短期水文氣象資料估算哈爾濱至同江冰厚度［J］.水科學(xué)進(jìn)展，2009，20（3）：428-433.

［19］ 戴芳芳，王自磐，李志軍，等.南極冬季威德爾海海冰物理結(jié)構(gòu)與葉綠素a垂直分布特征［J］.海洋學(xué)報：中文版，2008，30（4）：104-113.

［20］ 邰永芬，孔令杰.北方高寒地區(qū)魚類安全越冬要點［J］.黑龍江水產(chǎn)，2014（5）：38-40.

［21］ 徐廣遠(yuǎn)，鄶明，張恩鵬，等.冬季海冰災(zāi)害期間海參養(yǎng)殖管理方法［J］.中國水產(chǎn)，2010（11）：45-46.

［22］ 脫友才，劉志國，鄧云，等.豐滿水庫水溫的原型觀測及分析［J］.水科學(xué)進(jìn)展，2014，25（5）：731-738.

［23］ 張巖，李暢游，Shen Hongtao，等.烏梁素海湖泊冰生長過程中總氮的遷移規(guī)律［J］.水科學(xué)進(jìn)展，2013，24（5）：728-735.

［24］ 賈青，李志軍，梁景江，等 .利用現(xiàn)場觀測和力學(xué)試驗確定設(shè)計冰力學(xué)參數(shù)［J］.黑龍江水利科技，2014，42（10）：1-3.

［25］ 史興隆，佟錚，王呼和.河道冰體的可爆性研究［J］.內(nèi)蒙古工業(yè)大學(xué)學(xué)報，2014，33（2）：88-94.

LI Runling1，BAI Yila2，RONG Hui2，JIA Qing3，LI Zhijun1

（1.State Key Laboratory of Coastaland Offshore Engineering，Dalian University of Technology，Dalian 116024，China；2.DepartmentofMathematics，BohaiUniversity，Jinzhou 121000，China；3.College ofWater Conservancy and Electric Power，Heilongjiang University，Harbin 150086，China）

Discussions on the features of icesurface temperature in Daqing Reservoir and the methods of their practical utilizations in differentice engineering

The ice surface temperature in the winter of 2009/2010 was solved by using ice surface heat bal?ance equation，measured data of air temperature，solar radiation，wind speed，relative hum idity and ice thick?ness at Daqing reservoir，and recorded cloud cover，humidity and sea-level pressure at Anda Weather Sta?tion.The relationships between ice surface temperature and solar radiation，and air temperature，wind speed were analyzed.Based on these relationships，the ice surface temperature needed for traditional ice-resistance engineering was evaluated；the icesurface temperature used on safe ice-utilization engineering was quantita?tively discussed；daily variation of ice surface temperature in ecological adaptation in ice and under ice-cov?ering water in the climate changes，and the range of ice surface temperature used on ice blasting during iceflood periods were studied.In spit of the fact using the history data near the four observation times in practices because of international rules on weather observations，these analyzed relations are expected to pro?mote scientific and more practical utilizations on ice surface temperature in ice concerned engineering.

freshwater ice；surface temperature；air temperature；ice engineering；ice ecology

P338；X43

Adoi：10.13244/j.cnki.jiwhr.2015.04.006

11672-3031（2015）04-0277-08

（責(zé)任編輯：李福田）

2015-02-28

國家自然科學(xué)基金面上項目（41376186）；大連理工大學(xué)海岸和近海工程國家重點實驗室開放基金項目（LP1217）；黑龍江省教育廳科學(xué)技術(shù)研究項目（12531509）

李潤玲（1979-），女，河北張家口人，博士生，從事冰下水體環(huán)境和生態(tài)研究。E-mail：lirunlings@mail.d lut.edu.cn

李志軍（1960-），男，河北崇禮人，教授，博士，從事冰物理和力學(xué)性質(zhì)研究。E-mail：lizhijun@dlut.edu.cn