李婷　陳笑娟　張靜　孫玉龍　魏軍　李丹

摘要：基于災(zāi)害風(fēng)險(xiǎn)評(píng)估理論，利用極值概率分布模式，結(jié)合致災(zāi)因子和承災(zāi)體內(nèi)在屬性，構(gòu)建大風(fēng)風(fēng)速-日光溫室損失率脆弱性曲線(xiàn)，對(duì)河北省大風(fēng)災(zāi)害日光溫室受損風(fēng)險(xiǎn)進(jìn)行定量評(píng)估。研究表明，隨重現(xiàn)期的增加，風(fēng)災(zāi)危險(xiǎn)性及日光溫室災(zāi)損風(fēng)險(xiǎn)的強(qiáng)度和范圍也增加。空間上，地勢(shì)低平的東南部風(fēng)災(zāi)危險(xiǎn)性最高，主要位于滄州東部、衡水北部和邯鄲東部，日極大風(fēng)速超過(guò)32.6 m/s?；撅L(fēng)壓與三大地貌的分布呈現(xiàn)一致，高值出現(xiàn)在西北高原區(qū)和太行以東-燕山以南的平原區(qū)，太行-燕山山區(qū)較小。脆弱性曲線(xiàn)表明日光溫室開(kāi)始出現(xiàn)損失的臨界風(fēng)速值為20.0 m/s。溫室受損風(fēng)險(xiǎn)受致災(zāi)因子和承災(zāi)體共同作用，總體呈現(xiàn)西北低東南高，高風(fēng)險(xiǎn)區(qū)分布在唐山-秦皇島交界、唐山西部、廊坊及中南部地區(qū)，地均損失超過(guò)50萬(wàn)元/km2。中南部的溫室面積小于東北部，但其致災(zāi)因子危險(xiǎn)性較高，損失風(fēng)險(xiǎn)也較高。

關(guān)鍵詞：日光溫室;風(fēng)險(xiǎn);重現(xiàn)期;脆弱性曲線(xiàn);風(fēng)災(zāi)損失率

中圖分類(lèi)號(hào)：P429? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2020）04-0063-07

Abstract： Based on the disaster risk assessment theory and climate probability distribution model， the wind disaster risk to greenhouse in Hebei province was evaluated. The vulnerability curves of wind speed—solar greenhouse loss rate were constructed， considering disaster-causing factor and bearing body. The results showed that， with the increase of return periods， the strength and incidence of hazard factor and loss risk was greater. In terms of space， the southeastern area with low terrain had the highest hazard factor risk， mainly located in the eastern part of Cangzhou， as well as the northern Hengshui and eastern Handan. The extreme wind velocity was above 32.6 m/s. The basic wind pressure was consistent with the landforms distribution. The basic wind pressure was consistent with the distribution of the three landforms， with the high value occurring in the northwest plateau area and the plain area to the east of Taihang and the south of Yanshan Mountain areas. The vulnerability curve indicated that the threshold of wind speed which might cause loss to greenhouse was 20.0 m/s. The risk of greenhouse loss was affected by the combined effects of disaster-causing factor and bearing body. The risk in the northwest was generally lower than the southeast. The high-risk area was distributed in the junction of Tangshan-Qinhuangdao， western Tangshan， Langfang and south-central regions. The average loss was more than 500 000 yuan/km2. The greenhouse area in the south-central part was smaller than that in the northeast， but the loss was higher owing to its higher hazard factor risk.

Key words： greenhouse; risk; return period; vulnerability curver; wind disaster loss rate

日光溫室是中國(guó)北方地區(qū)特有的設(shè)施農(nóng)業(yè)結(jié)構(gòu)形式，具有效益高的特點(diǎn)[1-3]。其增溫保濕效應(yīng)促進(jìn)蔬果均衡生產(chǎn)，解決了冬季食品供應(yīng)難題，推動(dòng)農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整[4]。河北省是中國(guó)設(shè)施農(nóng)業(yè)生產(chǎn)大省[5]，設(shè)施蔬菜已成為農(nóng)業(yè)中最具競(jìng)爭(zhēng)力的優(yōu)勢(shì)產(chǎn)業(yè)，種植面積超過(guò)1.5萬(wàn)hm2。但是日光溫室結(jié)構(gòu)相對(duì)簡(jiǎn)單，對(duì)外界氣象條件的依賴(lài)性較高，災(zāi)害性天氣可對(duì)溫室作物的產(chǎn)量和品質(zhì)造成影響[6]。據(jù)統(tǒng)計(jì)，1984—2015年河北省平均每年發(fā)生大風(fēng)災(zāi)害40多次，最多年份可達(dá)119縣次，風(fēng)災(zāi)是影響溫室生產(chǎn)的頻發(fā)性災(zāi)害[7-9]，若防御不當(dāng)，風(fēng)災(zāi)可刮破棚膜，破壞保溫性，造成作物減產(chǎn)，重度風(fēng)災(zāi)會(huì)直接損毀溫室結(jié)構(gòu)，造成嚴(yán)重經(jīng)濟(jì)損失[2，9]。因此，研究日光溫室大風(fēng)災(zāi)害風(fēng)險(xiǎn)評(píng)估方法，建立科學(xué)的風(fēng)險(xiǎn)評(píng)估模型，是實(shí)施農(nóng)業(yè)風(fēng)險(xiǎn)管理的前提與基礎(chǔ)[10，11]，對(duì)于規(guī)避風(fēng)險(xiǎn)、穩(wěn)定發(fā)展有重要的實(shí)際意義。

以往對(duì)農(nóng)業(yè)風(fēng)險(xiǎn)評(píng)估研究對(duì)象集中于小麥、玉米、水稻等大田作物[12-14]，并且大多是在氣象數(shù)據(jù)和作物減產(chǎn)率的基礎(chǔ)上來(lái)建立災(zāi)害風(fēng)險(xiǎn)指標(biāo)體系，進(jìn)行等級(jí)劃分和風(fēng)險(xiǎn)評(píng)估。且針對(duì)日光溫室，小氣候試驗(yàn)及建造結(jié)構(gòu)研究較多[15，16]，風(fēng)災(zāi)方面，李楠等[2]通過(guò)風(fēng)洞試驗(yàn)，得到山東省日光溫室各區(qū)域發(fā)生風(fēng)災(zāi)的臨界風(fēng)速，并確定風(fēng)災(zāi)等級(jí)指標(biāo)。黃川容等[17]結(jié)合災(zāi)情數(shù)據(jù)資料，計(jì)算了不同風(fēng)力對(duì)北京溫室可造成災(zāi)害的概率，確定了日光溫室風(fēng)災(zāi)等級(jí)劃分的量化標(biāo)準(zhǔn)。根據(jù)區(qū)域?yàn)?zāi)害系統(tǒng)和風(fēng)險(xiǎn)評(píng)估理論[18]，風(fēng)險(xiǎn)應(yīng)是致災(zāi)因子、承災(zāi)體暴露度和脆弱性三者綜合作用的結(jié)果，但是從這三方面展開(kāi)設(shè)施農(nóng)業(yè)大風(fēng)災(zāi)害的研究基礎(chǔ)尚且薄弱。此外，前人的多數(shù)成果是在指標(biāo)權(quán)重確定方面，此方法仍是基于主觀(guān)判斷的定性評(píng)估方法，對(duì)指標(biāo)資料的數(shù)據(jù)需求量較大，并且該方法的合理性存在爭(zhēng)議[19]。脆弱性曲線(xiàn)法是建立不同致災(zāi)強(qiáng)度與承災(zāi)體損失的關(guān)系，進(jìn)而反演定量的期望損失值，這種模型需要大力發(fā)展[19，20]。

本研究依據(jù)災(zāi)害風(fēng)險(xiǎn)系統(tǒng)理論，對(duì)河北省日光溫室大風(fēng)災(zāi)害的致災(zāi)因子危險(xiǎn)性、承災(zāi)體暴露度和脆弱性進(jìn)行分析，采用脆弱性曲線(xiàn)法構(gòu)建評(píng)估模型，開(kāi)展風(fēng)險(xiǎn)定量評(píng)估，為風(fēng)災(zāi)風(fēng)險(xiǎn)預(yù)警服務(wù)及防御措施的制定提供技術(shù)支撐。

1? 資料與方法

1.1? 資料

1）氣象數(shù)據(jù)。河北省142個(gè)國(guó)家地面氣象站（1984—2016年）和602個(gè)區(qū)域地面氣象站（2014—2016年）的逐日觀(guān)測(cè)資料，觀(guān)測(cè)指標(biāo)為日最大風(fēng)速和日極大風(fēng)速，取自河北省氣象局。地理信息數(shù)據(jù)取自河北省地理信息局。

2）承災(zāi)體數(shù)據(jù)。河北省2015年各縣域日光溫室種植面積，取自河北省農(nóng)業(yè)廳，空間分布情況見(jiàn)圖1。根據(jù)河北省地方標(biāo)準(zhǔn)[21]，日光溫室是跨度大于6 m，中高大于2 m的利用太陽(yáng)輻射增溫，有后墻、山墻和后坡，前屋面硬蓋透明覆蓋物和保溫覆蓋物的栽培設(shè)施。

3）災(zāi)情數(shù)據(jù)。實(shí)地走訪(fǎng)河北省11個(gè)地市34個(gè)縣69個(gè)村，通過(guò)問(wèn)卷調(diào)查采集大風(fēng)風(fēng)速-設(shè)施農(nóng)業(yè)損失樣本數(shù)據(jù);1984—2016年河北省各區(qū)縣風(fēng)災(zāi)統(tǒng)計(jì)數(shù)據(jù)共1 306條記錄，包括風(fēng)災(zāi)發(fā)生的地點(diǎn)、開(kāi)始和結(jié)束時(shí)間、天氣過(guò)程、受災(zāi)人口、農(nóng)作物受災(zāi)情況等信息，取自河北省氣象局和《中國(guó)氣象災(zāi)害大典（河北卷）（2008）》。對(duì)數(shù)據(jù)進(jìn)行逐條核對(duì)、除重，共梳理得到有效樣本98個(gè)。

1.2? 研究方法

1.2.1? 致災(zāi)因子危險(xiǎn)性評(píng)估? 致災(zāi)因子危險(xiǎn)性是指致災(zāi)要素的強(qiáng)度、時(shí)間和變化規(guī)模。本研究選取日極大風(fēng)速作為大風(fēng)災(zāi)害強(qiáng)度指標(biāo)，進(jìn)行大風(fēng)災(zāi)害危險(xiǎn)性評(píng)估。日極大風(fēng)速觀(guān)測(cè)時(shí)間較短，參考姬鴻麗等[22]研究方法，建立區(qū)域站與國(guó)家站、最大風(fēng)速與極大風(fēng)速的回歸方程，推算日極大風(fēng)速長(zhǎng)時(shí)間序列數(shù)據(jù)。通過(guò)高低異常值檢驗(yàn)、空間異常值檢驗(yàn)、時(shí)間一致性檢查等數(shù)據(jù)質(zhì)量控制[23，24]，最終選擇反演效果較好的區(qū)域站385個(gè)。

根據(jù)極值理論來(lái)模擬極端天氣事件分布，是基于歷史實(shí)測(cè)資料的樣本分布情況來(lái)推算極端事件發(fā)生的可能性[25]，對(duì)于客觀(guān)反映災(zāi)害強(qiáng)度具有重要意義。重現(xiàn)期是超越概率的具體形式，通過(guò)多少年一遇來(lái)體現(xiàn)，對(duì)比不同年遇型風(fēng)險(xiǎn)的大小，可直觀(guān)識(shí)別由致災(zāi)因子強(qiáng)度變化而引起的空間變化趨勢(shì)。對(duì)于連續(xù)型隨機(jī)變量X，小于任意實(shí)數(shù)x的累積概率為F（X），則重現(xiàn)期（Return Period，RP）計(jì)算方法如下[26]：

式中，LR為不同風(fēng)速下的溫室損失率（Loss ratio，%）;L為風(fēng)災(zāi)造成的損失（萬(wàn)元）;W為承災(zāi)體（日光溫室）的重置費(fèi)用（萬(wàn)元）。根據(jù)災(zāi)情調(diào)查和風(fēng)災(zāi)統(tǒng)計(jì)數(shù)據(jù)，建立致災(zāi)強(qiáng)度與因?yàn)?zāi)損失率的關(guān)系進(jìn)行大風(fēng)災(zāi)害日光溫室脆弱性評(píng)估。建立線(xiàn)性、冪、指數(shù)、對(duì)數(shù)等回歸方程，擬合“風(fēng)速-日光溫室損失率”曲線(xiàn)，根據(jù)擬合度效果，確定脆弱性曲線(xiàn)。

1.2.3? 承災(zāi)體損失風(fēng)險(xiǎn)評(píng)估? 從風(fēng)險(xiǎn)產(chǎn)生的條件來(lái)看，風(fēng)險(xiǎn)是致災(zāi)因子、承災(zāi)體暴露度和脆弱性三者共同作用的結(jié)果[28]。承災(zāi)體損失風(fēng)險(xiǎn)為由于致災(zāi)因子導(dǎo)致的損失期望值，其計(jì)算公式如下[26]：

R（t）=V（Ht）×E? （6）

式中，t為重現(xiàn)期;R（t）為風(fēng)險(xiǎn);Ht為致災(zāi)危險(xiǎn)性;V為承災(zāi)體脆弱性曲線(xiàn);E為承災(zāi)體的價(jià)值。

2? 結(jié)果與分析

2.1? 致災(zāi)因子危險(xiǎn)性評(píng)估

河北省大風(fēng)災(zāi)害致災(zāi)因子危險(xiǎn)性評(píng)估結(jié)果見(jiàn)圖2。不同重現(xiàn)期下的風(fēng)災(zāi)危險(xiǎn)性分布大體一致，并且隨著重現(xiàn)期的增大，極大風(fēng)速的致災(zāi)強(qiáng)度逐漸增加，影響范圍也在擴(kuò)大。風(fēng)災(zāi)高危險(xiǎn)區(qū)主要分布在河北省東南部、承德東北部以及西北邊緣地區(qū)。其中，東南部地區(qū)的風(fēng)災(zāi)危險(xiǎn)性最高，主要位于滄州東部、衡水北部和邯鄲東部;東北部地區(qū)的危險(xiǎn)性次之，以承德中北部為主;西北部地區(qū)以張家口西部和保定中西部為主。5年一遇極大風(fēng)速致災(zāi)強(qiáng)度的高值區(qū)位于太行山以東的衡水北部和邯鄲東部地區(qū)，零星分布有28.5 m/s（11級(jí)）以上的大風(fēng)。10年一遇全省大部分地區(qū)的極大風(fēng)速?gòu)?qiáng)度已超過(guò)20.8 m/s，風(fēng)力在9級(jí)及以上。從30年、50年到百年一遇，風(fēng)力和高值區(qū)域進(jìn)一步增大，承德中北部、張家口西部、保定以南的風(fēng)速在24.5 m/s（10級(jí)）以上，其中東南部的平原城市風(fēng)速已超過(guò)32.6 m/s，達(dá)到12級(jí)大風(fēng)水平。

根據(jù)國(guó)家標(biāo)準(zhǔn)《建筑結(jié)構(gòu)荷載規(guī)范》[27]的規(guī)定，基本風(fēng)壓采用重現(xiàn)期為50年一遇日最大風(fēng)速計(jì)算得到，見(jiàn)圖3。其空間分布情況與河北省的地形地貌趨于一致，在太行山脈以東-燕山山脈以南的平原地區(qū)和冀西北的張家口高原地區(qū)，是基本風(fēng)壓的高值區(qū)，在0.22 kN/m2以上。張家口的康保、尚義和蔚縣，廊坊的霸州、永清和安次，石家莊無(wú)極、辛集，衡水深州和安平，保定市區(qū)以及滄州海興縣，建筑結(jié)構(gòu)設(shè)計(jì)應(yīng)滿(mǎn)足基本風(fēng)壓在0.3 kN/m2以上，最高值出現(xiàn)在張家口蔚縣，達(dá)1.1 kN/m2。

2.2? 承災(zāi)體脆弱性評(píng)估

從災(zāi)情數(shù)據(jù)中選取有日光溫室損失的大風(fēng)災(zāi)害案例，計(jì)算該次風(fēng)災(zāi)事件的溫室損失率，以致災(zāi)因子強(qiáng)度為x軸，以損失率為y軸，構(gòu)建大風(fēng)風(fēng)速-日光溫室損失率對(duì)應(yīng)關(guān)系，選擇線(xiàn)性、指數(shù)、冪和對(duì)數(shù)函數(shù)分別對(duì)大風(fēng)災(zāi)害日光溫室脆弱性曲線(xiàn)進(jìn)行擬合，對(duì)結(jié)果進(jìn)行擬合優(yōu)度檢驗(yàn)，結(jié)果見(jiàn)圖4。確定系數(shù)（R2）越接近于1，擬合度越好，因此選擇R2較高且形式簡(jiǎn)單的線(xiàn)性函數(shù)表征日光溫室風(fēng)災(zāi)脆弱性。由于大風(fēng)災(zāi)害造成的損失具有較大的不確定性，為避免這種影響，采用分段函數(shù)形式表示如下。根據(jù)脆弱性曲線(xiàn)可以看出，在一定風(fēng)速范圍內(nèi)，日光溫室風(fēng)災(zāi)損失率隨著風(fēng)速的增大而增加，開(kāi)始出現(xiàn)損失的臨界風(fēng)速值為20.0 m/s。

y=0 （x≤20.0）

y=4.856 2x-97.149 （20.0

y=100 （x≥40.6）

2.3? 承災(zāi)體損失風(fēng)險(xiǎn)評(píng)估

河北省日光溫室在冀北的張家口種植較少，種植密集區(qū)主要分布在東北部的承德、唐山、秦皇島，面積在2 000 hm2以上。綜合考慮致災(zāi)因子和承災(zāi)體因素，得到河北省日光溫室大風(fēng)災(zāi)害的受損風(fēng)險(xiǎn)如圖5所示。日光溫室受損風(fēng)險(xiǎn)總體呈現(xiàn)西北低東南高的趨勢(shì)。隨著重現(xiàn)期的增加，損失范圍逐漸擴(kuò)大，損失值有所增加。重現(xiàn)期為5年和10年一遇，唐山西部、廊坊北部、保定中部、衡水北部和邯鄲東部的溫室損失風(fēng)險(xiǎn)相對(duì)較高，其他地區(qū)的地均損失在10萬(wàn)元/km2以下。承德中部和唐山東部以及邢臺(tái)東部遭遇30年一遇的風(fēng)險(xiǎn)有所增加，損失值達(dá)30萬(wàn)元/km2以上。到50年和100年一遇時(shí)，高風(fēng)險(xiǎn)區(qū)域分布在唐山-秦皇島交界、唐山西部、廊坊以及太行山以東、沿海以西的中南部地區(qū)，損失超過(guò)50萬(wàn)元/km2。

3? 小結(jié)與討論

本研究基于極值理論和概率分布函數(shù)，采用重現(xiàn)期的形式來(lái)定量評(píng)價(jià)河北省大風(fēng)致災(zāi)強(qiáng)度，充分考慮了全球氣候變化背景下極端事件的影響，而傳統(tǒng)的數(shù)據(jù)統(tǒng)計(jì)主要采用平均數(shù)，以中間數(shù)據(jù)為基礎(chǔ)進(jìn)行分析，但極端事件的影響大多分布在函數(shù)的尾端[29]。研究得到的大風(fēng)致災(zāi)高危險(xiǎn)區(qū)分布在東南部、東北部和西北部，這與孫霞等[30]基于灰色理論得到的風(fēng)災(zāi)危險(xiǎn)度分布結(jié)果一致，由于東南部的平原地區(qū)地勢(shì)低平，地表的阻擋物（如山地、丘陵等）較少，為風(fēng)的流動(dòng)和發(fā)展提供了空間，風(fēng)力作用更強(qiáng)，其中衡水北部和邯鄲東部的風(fēng)速最大，超過(guò)32.6 m/s，達(dá)到12級(jí)水平?；撅L(fēng)壓的分布與地形呈現(xiàn)一致，西北高原和東南平原區(qū)較高，對(duì)于壩上高原，其下墊面廣闊而單一，有利于氣壓中心的形成，易與周邊產(chǎn)生較大的氣壓差[31]，冷熱不均導(dǎo)致風(fēng)力較大。而太行-燕山山區(qū)，由于山的阻擋，氣流阻力大，其風(fēng)速和風(fēng)壓值均較小。

前人的風(fēng)險(xiǎn)評(píng)估方法大多建立在指標(biāo)體系的基礎(chǔ)上[32]，脆弱性曲線(xiàn)法這種定量的評(píng)估方法用于風(fēng)災(zāi)研究的不多，并且關(guān)于日光溫室風(fēng)災(zāi)研究集中于等級(jí)劃分，側(cè)重于考慮致災(zāi)因子。本研究綜合考慮致災(zāi)因子和承災(zāi)體的設(shè)施農(nóng)業(yè)風(fēng)災(zāi)定量風(fēng)險(xiǎn)評(píng)估，通過(guò)構(gòu)建風(fēng)速-損失率脆弱性曲線(xiàn)，能夠判斷風(fēng)災(zāi)對(duì)日光溫室的致災(zāi)閾值，即20.0 m/s，進(jìn)而計(jì)算出地均損失值，其空間分布與致災(zāi)因子、承災(zāi)體呈現(xiàn)不一致，總體表現(xiàn)為西北低東南高。雖然日光溫室在中南部的種植面積小于東北部，但其致災(zāi)因子危險(xiǎn)性較高，因此受風(fēng)災(zāi)影響的損失風(fēng)險(xiǎn)較高，而張家口由于承災(zāi)體暴露度較低，其風(fēng)災(zāi)風(fēng)險(xiǎn)低，充分說(shuō)明承災(zāi)體的損失風(fēng)險(xiǎn)既受致災(zāi)因子強(qiáng)度的影響，也與承災(zāi)體的暴露情況和內(nèi)在易損性密切相關(guān)[32]，這是風(fēng)險(xiǎn)存在的必要條件。

本研究為河北省設(shè)施農(nóng)業(yè)大風(fēng)災(zāi)害風(fēng)險(xiǎn)預(yù)警業(yè)務(wù)的開(kāi)展奠定了一定基礎(chǔ)，但仍存在一些不足：選取致災(zāi)因子危險(xiǎn)性評(píng)估指標(biāo)為極大風(fēng)速，即瞬時(shí)風(fēng)速，而風(fēng)災(zāi)發(fā)生還與大風(fēng)持續(xù)時(shí)長(zhǎng)有關(guān)，因此有必要進(jìn)一步參考熵權(quán)系數(shù)[23]等方法，嘗試構(gòu)建綜合性的致災(zāi)強(qiáng)度指標(biāo);此外，本研究考慮了絕大多數(shù)溫室的結(jié)構(gòu)、材料和跨度情況，并非全部，因此實(shí)地調(diào)查結(jié)果及脆弱性曲線(xiàn)的適用性存在一定局限，具體類(lèi)型還需具體分析;相同風(fēng)速不同風(fēng)向?qū)厥业钠茐淖饔靡膊煌柽M(jìn)一步研究不同風(fēng)向?qū)厥业膿p失情況，以更好地指導(dǎo)日光溫室風(fēng)險(xiǎn)預(yù)警和防御工作。

參考文獻(xiàn)：

[1] 高? 浩，黎貞發(fā)，潘學(xué)標(biāo)，等.中國(guó)設(shè)施農(nóng)業(yè)氣象業(yè)務(wù)服務(wù)現(xiàn)狀與對(duì)策[J].中國(guó)農(nóng)業(yè)氣象，2010，31（3）：402-406.

[2] 李? 楠，薛曉萍，李鴻怡.基于信息擴(kuò)散的山東省日光溫室風(fēng)災(zāi)風(fēng)險(xiǎn)評(píng)估[J].氣象與環(huán)境學(xué)報(bào)，2018，34（5）：149-155.

[3] 楊再?gòu)?qiáng)，張婷華，黃海靜，等.北方地區(qū)日光溫室氣象災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)氣象，2013，34（3）：342-349.

[4] 魏瑞江，李春強(qiáng)，康西言.河北省日光溫室低溫寡照災(zāi)害風(fēng)險(xiǎn)分析[J].自然災(zāi)害學(xué)報(bào)，2008，17（3）：56-62.

[5] 陳貴林，張進(jìn)文.河北省設(shè)施蔬菜生產(chǎn)可持續(xù)發(fā)展問(wèn)題的探討[J].沈陽(yáng)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2000，31（1）：29-31.

[6] 王? 瓊，魏瑞江，王榮英，等.河北日光溫室氣象災(zāi)害影響和氣象服務(wù)評(píng)估[J].中國(guó)農(nóng)業(yè)氣象，2014，35（6）：682-689.

[7] 劉金玉，孫愛(ài)良，董航宇.黃驊市設(shè)施蔬菜生產(chǎn)主要?dú)庀鬄?zāi)害特征及其防御措施研究[J].現(xiàn)代農(nóng)業(yè)科技，2013（20）：239-240，247.

[8] 閻? 訪(fǎng)，岳艷霞，張秉祥，等.石家莊地區(qū)大風(fēng)日數(shù)氣候變化特征[J].氣象與環(huán)境學(xué)報(bào)，2013，29（2）：81-88.

[9] 郭曉冬.災(zāi)害性天氣對(duì)日光溫室蔬菜的危害及預(yù)防措施[J].甘肅農(nóng)業(yè)科技，2008（9）：50-52.

[10] 孫? 霞，俞海洋，孫? 斌，等.河北省主要?dú)庀鬄?zāi)害時(shí)空變化的統(tǒng)計(jì)分析[J].干旱氣象，2014，32（3）：388-392.