• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于三維多上行先導(dǎo)模型的高建筑物年雷擊次數(shù)分析

    2023-07-31 05:07:02樊佳樂(lè)譚涌波余駿皓林雨荷
    科學(xué)技術(shù)與工程 2023年20期
    關(guān)鍵詞:電位次數(shù)建筑物

    樊佳樂(lè),譚涌波*,余駿皓,林雨荷

    (1.南京信息工程大學(xué)大氣物理學(xué)院,南京 210044;2.南京信息工程大學(xué)中國(guó)氣象局氣溶膠與云降水重點(diǎn)開(kāi)放實(shí)驗(yàn)室,南京 210044)

    隨著社會(huì)和科技的發(fā)展,高建筑物越來(lái)越多,高建筑物的雷擊災(zāi)害逐年頻發(fā),對(duì)建筑物的雷電防護(hù)和雷擊災(zāi)害風(fēng)險(xiǎn)評(píng)估也越發(fā)被重視。年預(yù)計(jì)雷擊次數(shù)是確定建筑物防雷類別和風(fēng)險(xiǎn)評(píng)估中最主要的參數(shù)之一,對(duì)雷擊次數(shù)的準(zhǔn)確計(jì)算是提高建筑物安全性,準(zhǔn)確評(píng)估雷擊風(fēng)險(xiǎn)的基礎(chǔ)[1-2]。目前關(guān)于建筑物雷擊次數(shù)的工程計(jì)算主要有ProtectionofStructuresAgainstLightning:Part1:GeneralPrinciples(IEC61024-1)[3]、ProtectionAgainstLightning:Part2:RiskManagementforStructuresandServices(IEC 62305-2[4],下文簡(jiǎn)稱IEC標(biāo)準(zhǔn))給出的具體計(jì)算公式;除此之外,一些學(xué)者通過(guò)閃電觀測(cè)和數(shù)值模擬分析不同情況下雷擊特征,從理論的角度研究影響建筑物雷擊概率的決定性因素[5-6]。

    國(guó)際電工委員會(huì)(International Electrotechnical Commission,IEC)標(biāo)準(zhǔn)中給出了建筑物雷擊次數(shù)的經(jīng)驗(yàn)公式,利用建筑物等效截收面積與所處地區(qū)的年平均地閃密度兩者相乘得到建筑物年預(yù)計(jì)雷擊次數(shù)。根據(jù)定義等效截收面積為與建筑物遭受相同雷擊次數(shù)的等效面積,利用建筑物自身的長(zhǎng)寬高計(jì)算得出。ProtectionofStructuresAgainstLightning:Part1:GeneralPrinciples(IEC61024-1)[3]首次給出建筑物年預(yù)計(jì)雷擊次數(shù)的計(jì)算公式,在之后的ProtectionAgainstLightning:Part2:RiskManagementforStructuresandServices(IEC 62305-2)[4]中對(duì)年平均地閃密度的估算進(jìn)行了修改,等效截收面積的計(jì)算沒(méi)有修改。中國(guó)《建筑物防雷設(shè)計(jì)規(guī)范》(GB 50057—2010[7],下文簡(jiǎn)稱國(guó)標(biāo))中所采用的雷擊次數(shù)計(jì)算方法參照了IEC標(biāo)準(zhǔn)中的經(jīng)驗(yàn)公式。相較于IEC標(biāo)準(zhǔn),國(guó)標(biāo)中區(qū)分了高于100 m和低于100 m建筑物等效截收面積的計(jì)算,使得對(duì)建筑物雷擊次數(shù)的計(jì)算更加細(xì)致。不管是IEC標(biāo)準(zhǔn)還是中國(guó)國(guó)標(biāo)關(guān)于建筑物等效截收面積的計(jì)算公式中均考慮了建筑物自身幾何特征對(duì)雷擊次數(shù)的影響,但并沒(méi)有考慮自然閃電的極性,強(qiáng)度以及空間形態(tài)等特征在建筑物雷擊事件中所起到的關(guān)鍵性作用[8-9]。標(biāo)準(zhǔn)中的經(jīng)驗(yàn)公式在一定程度上方便了工程應(yīng)用,但與閃電臨近地面引發(fā)多上行先導(dǎo)(upward leader,UL)并選擇和其中一個(gè)或多個(gè)先導(dǎo)連接從而形成回?fù)暨^(guò)程的物理過(guò)程還存在一定的差異。

    大量觀測(cè)表明,閃電更容易擊中高建筑物,建筑物的高度是影響雷擊次數(shù)的重要因素。目前,對(duì)于高建筑物雷擊特征的觀測(cè)主要有光學(xué)觀測(cè)和閃電定位技術(shù)兩種手段。Eriksson[10]根據(jù)不同高度建筑物觀測(cè)到的閃電數(shù)目推導(dǎo)出雷電擊中建筑物的概率公式,并通過(guò)年平均地閃密度給出建筑物年雷擊次數(shù)的經(jīng)驗(yàn)公式。吳珊珊等[11]利用廣州高建筑物雷電觀測(cè)站6年的光學(xué)觀測(cè)資料進(jìn)行分析,發(fā)現(xiàn)在119次下行地閃中,廣州塔附近1 km區(qū)域內(nèi)未觀測(cè)到地閃,1~3 km范圍內(nèi)隨著距離增加下行地閃密度逐漸增加,表明廣州塔吸引了周?chē)南滦械亻W,且隨著距離的增加吸引作用逐漸減弱。光學(xué)觀測(cè)能明確閃電擊中建筑物的情況,但由于樣本量的有限,導(dǎo)致對(duì)建筑物年雷擊次數(shù)的研究相對(duì)較少。除此之外,一些學(xué)者通過(guò)對(duì)閃電定位數(shù)據(jù)的分析來(lái)評(píng)估高建筑物對(duì)地閃活動(dòng)的影響。Ngqungqa[12]通過(guò)地閃定位數(shù)據(jù)分析南非兩座鐵塔附近的閃電活動(dòng),發(fā)現(xiàn)鐵塔周?chē)霃?.5 km內(nèi)的地閃密度大于2.5~10 km圓環(huán)區(qū)域的地閃密度,表明鐵塔吸引了閃電。Zhang等[13]分析了廣州塔附近的地閃定位數(shù)據(jù),發(fā)現(xiàn)廣州塔近1 km半徑范圍內(nèi)閃電密度顯著增加,而半徑1~4 km區(qū)域明顯減小,推測(cè)廣州塔吸引了一些下行負(fù)地閃。閃電定位的數(shù)據(jù)充足,卻由于精度較差的原因,導(dǎo)致難以明確建筑物是否遭雷擊,也較難給出建筑物年雷擊次數(shù)的經(jīng)驗(yàn)公式。

    模擬仿真是觀測(cè)研究很好的補(bǔ)充手段之一。許多學(xué)者利用放電模型模擬地閃的發(fā)展與連接過(guò)程,通過(guò)敏感性實(shí)驗(yàn),研究各種特征參量在建筑物雷擊事件中的影響。Zon等[14]利用二維隨機(jī)模型模擬了閃電的起始,下行先導(dǎo)(downward leader,DL)的梯階傳播以及始發(fā)于地面建筑物的UL,并研究了建筑物的雷擊情況,發(fā)現(xiàn)建筑物越高雷擊概率越大,閃電始發(fā)距離越遠(yuǎn)建筑物雷擊概率越小。Petrov等[15]利用三維隨機(jī)模型模擬閃電放電,探討了閃電極性以及建筑物幾何特性對(duì)建筑物雷擊概率的影響。譚涌波等[8]利用二維隨機(jī)模型,通過(guò)改變閃電空間形態(tài),研究閃電形態(tài)的差異對(duì)雷擊位置,UL長(zhǎng)度以及連接過(guò)程形態(tài)等參數(shù)的影響。這些研究在閃電始發(fā)位置相對(duì)固定的前提下,給出了影響雷擊概率的關(guān)鍵性因子,但沒(méi)有進(jìn)一步給出建筑物年雷擊次數(shù)的計(jì)算公式。還有一些研究利用拓展的電氣幾何模型計(jì)算多建筑物不同位置的雷擊概率以及附近區(qū)域的雷擊點(diǎn)分布。結(jié)果表明,建筑物拐角和避雷針更容易被閃電擊中[16-17]。這些研究對(duì)更復(fù)雜的建筑物群進(jìn)行研究,但未能考慮自然閃電的不同特征對(duì)雷擊概率的影響。Jiang等[18]利用二維隨機(jī)模型模擬負(fù)地閃的發(fā)展與連接過(guò)程,探討了高建筑物對(duì)雷擊點(diǎn)的影響,分析了高建筑物的保護(hù)距離,給出了建筑物下行負(fù)地閃頻次的公式。由于模式為二維隨機(jī)模型,與真實(shí)三維空間存在一定的差異。

    綜上所述,建筑物年雷擊次數(shù)是建筑物安全性能評(píng)估的重要指標(biāo)之一,IEC標(biāo)準(zhǔn)和中國(guó)國(guó)標(biāo)都給出了各自的經(jīng)驗(yàn)計(jì)算公式,但均未考慮自然閃電特性對(duì)這個(gè)關(guān)鍵因子的可能影響,而已有的一些科研工作還是以定性的理論探討為主。如何將建筑物的幾何特性與自然閃電的傳播特征相結(jié)合,在盡可能還原真實(shí)雷擊物理過(guò)程的前提下,給出建筑物年雷擊次數(shù)的計(jì)算公式,仍需要深入研究。鑒于此,在最新搭建的三維高分辨率多先導(dǎo)隨機(jī)模型基礎(chǔ)上進(jìn)行修改,通過(guò)改變建筑物高度,閃電始發(fā)的位置以及閃電強(qiáng)度進(jìn)行大量模擬,分析不同情況下建筑物的雷擊概率和年雷擊次數(shù)。在一定程度上更精細(xì)的計(jì)算了建筑物的年雷擊次數(shù),為高建筑物風(fēng)險(xiǎn)評(píng)估提供可靠依據(jù)。

    1 模式介紹

    閃電連接過(guò)程的模型涉及到建筑物的高度,寬度等特征,因此需要較高的分辨率才能更精準(zhǔn)的還原建筑物的特征以及閃電的空間形態(tài)。而DL每往前延伸一步,都需通過(guò)迭代技術(shù)對(duì)空間電位進(jìn)行重置計(jì)算[19],需要耗費(fèi)大量的計(jì)算資源,這也是目前難以實(shí)現(xiàn)高分辨率三維閃電連接過(guò)程模擬的重要原因[20]。本課題組基于GPU并行計(jì)算解決了三維高分辨率泊松方程的求解問(wèn)題[20],并在此基礎(chǔ)上搭建了三維多先導(dǎo)隨機(jī)閃電連接過(guò)程的模式框架[21],考慮了DL在傳播過(guò)程中毛細(xì)分支的退電離過(guò)程,讓模擬得到的先導(dǎo)能更好地再現(xiàn)真實(shí)閃電的空間多分叉且蜿蜒延伸的空間形態(tài)。在已有的三維高分辨率隨機(jī)多先導(dǎo)連接模型的基礎(chǔ)上,通過(guò)敏感性實(shí)驗(yàn),探討孤立建筑物年雷擊次數(shù)的影響因子,最終給出計(jì)算公式(模式的基本框架可參考文獻(xiàn)[22])。

    本文模型不考慮云中放電的部分,只截取閃電到達(dá)近地面附近的發(fā)展情況,如圖1所示。考慮到自然界中負(fù)地閃占所有地閃過(guò)程的90%以上[23],只考慮負(fù)地閃所引發(fā)的上行正先導(dǎo)以及它們之間的回?fù)暨^(guò)程。在模擬域頂部設(shè)置一段初始的DL,并賦予這段DL一個(gè)初始的電位,作為閃電在近地面發(fā)展的起始條件??紤]到實(shí)際的閃電會(huì)從空間中任意位置始發(fā),對(duì)初始DL在模擬域頂部的位置做出一些調(diào)整(圖2)。模式中DL的傳播是由通道點(diǎn)與周?chē)h(huán)境點(diǎn)的電位差決定,當(dāng)電位差超過(guò)150 kV/m時(shí)[18],閃電先導(dǎo)發(fā)展。根據(jù)觀測(cè)經(jīng)驗(yàn),DL可以產(chǎn)生很多分支[24-25],在模型中允許DL每一個(gè)發(fā)展點(diǎn)都可以向周?chē)狞c(diǎn)發(fā)展,再根據(jù)周?chē)c(diǎn)的電場(chǎng)值的權(quán)重隨機(jī)選取待發(fā)展點(diǎn)[20]。當(dāng)DL到達(dá)近地面時(shí),會(huì)從建筑物尖端或者地面引發(fā)多個(gè)UL,且UL不會(huì)有明顯的分支,本文模型允許滿足UL起始閾值180 kV/m時(shí)的建筑物或者地面上的多個(gè)格點(diǎn)始發(fā),判斷頭部點(diǎn)與周?chē)h(huán)境點(diǎn)的電位梯度,當(dāng)滿足傳播閾值150 kV/m時(shí)[18],隨機(jī)選取一個(gè)待發(fā)展點(diǎn)發(fā)展且只允許UL頭部發(fā)展。根據(jù)觀測(cè)經(jīng)驗(yàn),設(shè)置正負(fù)先導(dǎo)的速度比為1∶6[26-27]。先導(dǎo)每發(fā)展一步都會(huì)對(duì)全域電位重解,當(dāng)DL和UL之間兩點(diǎn)的電位差大于連接閾值500 kV/m時(shí)發(fā)生連接[18,28]。實(shí)際雷暴過(guò)程中,由于雷暴云中存在大量電荷使得環(huán)境電場(chǎng)增強(qiáng),故設(shè)置背景電場(chǎng)為-10 kV/m。在模型中,建筑物、地面、先導(dǎo)通道以及模擬域上邊界均滿足Dirichlet邊界條件,模擬域的側(cè)邊界滿足Neumann邊界條件。

    藍(lán)線表示DL通道;紅線表示UL通道結(jié)構(gòu)圖1 模型示意圖Fig.1 Schematic diagram of the model

    圖2 模擬區(qū)域示意圖Fig.2 Schematic diagram of simulation area

    選取模擬域大小為2 km×2 km×1.5 km,空間分辨率為10 m×10 m×10 m。初始DL長(zhǎng)度設(shè)為50 m,為了研究不同強(qiáng)度的閃電在建筑物雷擊事件中的影響,初始下行負(fù)先導(dǎo)內(nèi)部初始電位設(shè)為-20、-30、-40 MV來(lái)模擬不同強(qiáng)度的負(fù)地閃到達(dá)近地面的情況[9,29]。圖2展示了模擬區(qū)域示意圖。大量研究表明,建筑物高度是影響雷擊事件的重要因素[11],將建筑物的長(zhǎng)寬設(shè)為40 m,改變建筑物高度(100~600 m,間隔為100 m),并將其固定在模擬域中心位置,圖2中不同高度用不同顏色線條進(jìn)行標(biāo)注??紤]到閃電會(huì)從建筑物上空任意位置始發(fā),設(shè)置建筑物與初始DL位置的水平距離范圍為0~600 m,間隔為100 m,形成了6個(gè)半徑相差100 m同心圓(圖2)。在選取DL起始位置時(shí)只能取整數(shù)點(diǎn),不能取到一個(gè)完整的圓,因此設(shè)置每個(gè)圓的半徑誤差在10 m內(nèi)來(lái)取到更多的點(diǎn),圖2中展示了在一組模擬中取到的DL起始點(diǎn)。主要探討閃電強(qiáng)度(3種),閃電初始位置相對(duì)于建筑物的水平距離(7種),以及建筑物高度(6種)對(duì)雷擊概率的影響,總共有126種方案,每種方案模擬300次。來(lái)獲得統(tǒng)計(jì)上有意義的結(jié)果。

    2 結(jié)果與分析

    2.1 模擬結(jié)果

    隨機(jī)放電模型能夠較好地再現(xiàn)先導(dǎo)通道的蜿蜒延伸、多分叉的空間形態(tài)[20],在建筑物雷擊過(guò)程的模擬中,能更好地仿真DL幾何特征對(duì)連接過(guò)程的影響,有利于探討閃電的多樣性對(duì)建筑物年雷擊次數(shù)的影響。從37 800個(gè)模擬結(jié)果中選取了6次典型的雷擊個(gè)例(圖3),圖3(a)為建筑物只始發(fā)一個(gè)UL,閃電擊地位置唯一確定;圖3(b)、圖3(c)為建筑物始發(fā)了多個(gè)UL集中建筑物頂端;圖3(d)為DL繞過(guò)建筑物頂端始發(fā)的UL,擊中建筑物側(cè)面從而形成側(cè)擊現(xiàn)象;圖3(e)為閃電直接擊中地面的情形;圖3(f)為建筑物雖然始發(fā)UL,但并未被擊中的情況。而這些模擬的典型雷擊個(gè)例在以往的觀測(cè)研究中均有報(bào)道[25,30-31],從而驗(yàn)證了LAMM能較好地仿真各種雷擊事件。與DL連接的UL的起始點(diǎn)為閃電的接地點(diǎn),根據(jù)接地點(diǎn)位置的不同,將模擬實(shí)驗(yàn)結(jié)果分為擊中建筑物和未擊中建筑物兩種情形。其中,定義擊中建筑物的情形為建筑物觸發(fā)單個(gè)UL并發(fā)生連接[圖3(a)]、觸發(fā)多個(gè)UL且與其中一個(gè)發(fā)生連接[圖3(b)、圖3(c)]以及擊中建筑物側(cè)面[圖3(d)]。定義未擊中建筑物的情形分別為建筑物未觸發(fā)UL直接擊中地面[圖3(e)]和建筑物觸發(fā)UL但未產(chǎn)生連接擊中地面[圖3(f)]。在未擊中建筑物的個(gè)例中,一般建筑物較矮時(shí)不易觸發(fā)上行先導(dǎo),而隨著建筑物的增高,上行先導(dǎo)觸發(fā)概率增大,但其并不與DL連接,從而形成上行未連接先導(dǎo)。

    黑色長(zhǎng)方為建筑物;藍(lán)線表示DL通道結(jié)構(gòu);紅線表示UL通道結(jié)構(gòu)圖3 模擬結(jié)果示意圖。Fig.3 Schematic diagram of simulation results

    在模擬過(guò)程中,建筑物高度和閃電始發(fā)位置的改變,都會(huì)出現(xiàn)不同的雷擊特征。如圖3(a)~圖3(c)所示,DL初始電位均為-30 MV,建筑物高度分別為100、300、600 m,閃電始發(fā)位置距建筑物水平距離均為100 m。隨著建筑物高度的增加,觸發(fā)多個(gè)UL的概率明顯增大。圖3(c)~圖3(f)中,建筑物高度均為600 m,DL初始電位均為-30 MV,閃電始發(fā)位置距建筑物水平距離分別為100、400、600 m。建筑物高度相同時(shí),隨著閃電始發(fā)位置的遠(yuǎn)離,偶爾會(huì)發(fā)生側(cè)擊現(xiàn)象,未擊中建筑物情況的次數(shù)也呈現(xiàn)出增多的趨勢(shì)。從圖3給出的幾次典型個(gè)例中,不難發(fā)現(xiàn)建筑物高度以及閃電始發(fā)距離都會(huì)對(duì)雷擊事件產(chǎn)生影響。

    為了分析不同強(qiáng)度的閃電對(duì)建筑物雷擊事件的影響,通過(guò)改變DL初始電位進(jìn)行模擬。以600 m的建筑物為例,圖4給出了不同強(qiáng)度的雷擊過(guò)程模擬結(jié)果。圖4(a)、圖4(c)和圖4(b)、圖4(d)的初始電位分別為-20 MV和-40 MV,閃電始發(fā)位置距建筑物水平距離分別為200 m和400 m。從圖4中不難發(fā)現(xiàn),不同強(qiáng)度的閃電在模擬地閃過(guò)程中表現(xiàn)出了明顯的區(qū)別:①閃電初始電位越大,DL的分支越多,閃電發(fā)展越旺盛;②閃電強(qiáng)度越大,越容易觸發(fā)UL,UL發(fā)展的長(zhǎng)度越長(zhǎng);③隨著閃電始發(fā)位置的遠(yuǎn)離,初始電位越小,UL觸發(fā)時(shí)間較晚,使得閃電較難擊中建筑物。

    黑色長(zhǎng)方為建筑物;藍(lán)線表示DL通道結(jié)構(gòu);紅線表示UL通道結(jié)構(gòu)圖4 不同強(qiáng)度閃電的雷擊過(guò)程模擬結(jié)果Fig.4 Simulation results of lightning strike process with different intensities

    2.2 建筑物雷擊概率與年雷擊次數(shù)分析

    固定DL初始電位為-40 MV,探討建筑物高度和閃電始發(fā)位置對(duì)雷擊概率以及年雷擊次數(shù)的影響。根據(jù)式(1)計(jì)算建筑物雷擊概率D[h,r],定義雷擊概率是一個(gè)相對(duì)量,其數(shù)值受到閃電始發(fā)位置r的影響。

    (1)

    式(1)中:h為建筑物高度;r為閃電起始位置與建筑物的水平距離;n為當(dāng)建筑物高度為h,閃電始發(fā)位置為r時(shí),閃電擊中建筑物的次數(shù);N為該組模擬的總次數(shù),N=300。

    圖5展現(xiàn)了當(dāng)DL起始位置與建筑物水平距離為0,100,200,…,600 m時(shí),各高度建筑物的雷擊概率,其具體數(shù)值用十字符號(hào)標(biāo)出。由圖5可知,隨著建筑物高度的增加,雷擊概率呈現(xiàn)出增大的趨勢(shì),即建筑物越高,雷擊概率越大。當(dāng)閃電從建筑物正上方始發(fā)時(shí),100 m和200 m高度的建筑物,雷擊概率分別為0.86和0.99,其他高度建筑物的雷擊概率為1。隨著閃電始發(fā)位置逐漸遠(yuǎn)離建筑物,不同高度建筑物的雷擊概率呈現(xiàn)出明顯差異:其中,100 m和200 m建筑物對(duì)應(yīng)的雷擊概率,在閃電始發(fā)位置從建筑物正上方遠(yuǎn)離時(shí),立刻出現(xiàn)遞減趨勢(shì),而300 m和400 m的建筑物在閃電始發(fā)距離大于200 m時(shí),其雷擊概率才出現(xiàn)較明顯的遞減趨勢(shì),500 m和600 m的建筑物對(duì)應(yīng)的雷擊概率則分別在閃電始發(fā)距離大于300 m和400 m時(shí)呈現(xiàn)出遞減趨勢(shì)。

    由上述分析可知,除了建筑物高度以外,閃電的起始位置也是影響建筑物雷擊概率的重要影響因素,在各個(gè)高度組下的建筑物雷擊概率均隨閃電始發(fā)距離r的增大而減小。如圖5所示,隨著閃電起始位置逐漸偏離建筑物,各雷擊概率曲線間的間距先增大后減小,這是因?yàn)楫?dāng)閃電在離建筑物水平距離很近或很遠(yuǎn)處始發(fā)時(shí),100~600 m建筑物的雷擊概率都將趨近于1或0,此時(shí)建筑物高度對(duì)雷擊概率的影響不明顯。值得注意的是,閃電起始位置只對(duì)雷擊概率產(chǎn)生影響,而建筑物的年雷擊次數(shù)則不受其影響。事實(shí)上,建筑物年雷擊次數(shù)受區(qū)域地閃密度的影響,這需要進(jìn)行全區(qū)域的地閃活動(dòng)分析,后文中將會(huì)具體說(shuō)明。

    隨著閃電起始位置的變化,雖然不同高度建筑物對(duì)應(yīng)的雷擊概率在數(shù)值上表現(xiàn)出明顯差異,但其變化趨勢(shì)卻呈現(xiàn)出相似的特征。即當(dāng)閃電在建筑物一定范圍內(nèi)起始,不同高度建筑物的雷擊概率都趨近于1,一旦閃電起始位置超過(guò)這一范圍,雷擊概率則會(huì)快速下降,而當(dāng)閃電起始位置距建筑物很遠(yuǎn)時(shí),雷擊概率都趨近于0。根據(jù)上述共有特征,選取式(2)對(duì)每個(gè)高度的建筑物的雷擊概率進(jìn)行線性擬合,擬合結(jié)果如圖5中曲線所示,得出特定建筑物高度對(duì)應(yīng)的雷擊概率曲線fh(r)。表1給出式(2)中的參數(shù)a、b、c的數(shù)值以及擬合優(yōu)度R2。

    表1 雷擊概率擬合曲線公式各參數(shù)值及擬合優(yōu)度計(jì)算值Table 1 Parameter values and goodness of fit calculation values of lightning strike probability fitting curve formula

    (2)

    +表示當(dāng)DL起始位置與建筑物水平距離為 0,100,200,…,600 m 時(shí)不同高度建筑物雷擊概率的具體數(shù)值圖5 DL電位為-40 MV時(shí),不同高度建筑物的雷擊概率隨DL起始點(diǎn)水平距離的變化曲線Fig.5 Lightning strike probability of buildings with different heights varies with the horizontal distance of DL starting point when DL potential is -40 MV

    為分析閃電在建筑物上空一定范圍內(nèi)隨機(jī)發(fā)生時(shí)雷擊建筑物的情況,基于三維多上行先導(dǎo)模型的模擬結(jié)果,采用式(3)對(duì)建筑物雷擊概率進(jìn)行積分。考慮到閃電始發(fā)位置很遠(yuǎn)時(shí)雷擊概率都接近于0,并不會(huì)對(duì)積分后的數(shù)值造成影響,選取建筑物周?chē)?.5 km區(qū)域的雷擊概率進(jìn)行積分,得到閃電年平均密度為1次/km2/a時(shí),不同高度建筑物的年雷擊次數(shù)Nh。

    (3)

    式(3)中:Nh為當(dāng)區(qū)域閃電年平均密度為1次/km2/a時(shí),高度為h的建筑物的年雷擊次數(shù);r為DL起始位置距建筑物的水平距離,km;fh(r)為距建筑物不同距離始發(fā)閃電時(shí)建筑物的雷擊概率,其具體計(jì)算見(jiàn)式(2)。

    根據(jù)式(3)計(jì)算100~600 m高度建筑物對(duì)應(yīng)的閃電年雷擊次數(shù)分別為0.39、0.66、0.86、1.33、1.51、1.94,此時(shí)DL初始電位為-40 MV。顯然,年雷擊次數(shù)與建筑物高度呈正相關(guān),其中600 m建筑物被負(fù)地閃擊中的次數(shù)是100 m建筑物的5倍。根據(jù)計(jì)算值,如圖5所示進(jìn)行線性擬合,得出當(dāng)閃電年平均密度為1次/km2/a時(shí),利用建筑物高度h計(jì)算的年雷擊次數(shù)N的經(jīng)驗(yàn)公式(適用100~600 m高的建筑物)。

    N=1.29×10-6h2+2.18×10-3h+0.16

    (4)

    圖6對(duì)比了由式(4)以及國(guó)標(biāo)計(jì)算的建筑物年雷擊次數(shù)??梢钥闯?所得到的建筑物年雷擊次數(shù)曲線與國(guó)標(biāo)中給出的具有相同的變化趨勢(shì),建筑物高度與建筑物年雷擊次數(shù)呈明顯的正相關(guān)。并且隨著高度的增加曲線的遞增速率均增大,這是由于閃電從較遠(yuǎn)處始發(fā)時(shí),矮建筑物的雷擊概率為0,而高建筑物仍有被擊中的情況,通過(guò)積分得到的年雷擊次數(shù)會(huì)隨著建筑物高度的增加而增大。值得注意的是,雖然各高度組建筑物的年雷擊次數(shù)計(jì)算結(jié)果均略大于國(guó)標(biāo)的計(jì)算結(jié)果,但其在一定程度上更符合觀測(cè)經(jīng)驗(yàn)。Ngqungqa[12]通過(guò)觀測(cè)得出100 m建筑物的年雷擊次數(shù)為0.29(本文取0.39),而根據(jù)廣州塔(600 m)附近1 km范圍內(nèi)未觀測(cè)到閃電的觀測(cè)事實(shí)[11],推測(cè)當(dāng)閃電年平均密度為1次/km2/a時(shí),廣州塔的年雷擊次數(shù)約為3(本文取1.94)。廣州塔是一個(gè)尖頂?shù)慕ㄖ?當(dāng)雷暴云存在時(shí),建筑物頂端的電場(chǎng)畸變更大,使其更容易遭受雷擊。所探討的建筑物都為平頂,計(jì)算值略小具有一定的合理性。

    圖6 DL在高度1.5 km處始發(fā),且閃電年平均密度為1次/km2/a時(shí),不同高度建筑物的年雷擊次數(shù)Fig.6 Annual lightning strikes of buildings at different heights when DL starts at a height of 1.5 km and the average annual lightning density is 1 strike /km2/a

    2.3 改變閃電強(qiáng)弱分析

    觀測(cè)和模擬研究表明,閃電的強(qiáng)弱會(huì)影響閃擊距離等參數(shù)。不同強(qiáng)度的閃電對(duì)建筑物雷擊概率的影響不可忽略。由于模型中未模擬回?fù)暨^(guò)程,因此無(wú)法如觀測(cè)中利用回?fù)舴逯惦娏髋袛嚅W電的強(qiáng)弱,而是通過(guò)設(shè)置不同的DL初始電位,使得閃電通道與地面形成不同的電勢(shì)差,以此間接體現(xiàn)閃電的強(qiáng)弱。

    設(shè)置DL初始電位為-20 MV和-30 MV,統(tǒng)計(jì)閃電從距離建筑物不同水平位置始發(fā)時(shí),各個(gè)高度建筑物對(duì)應(yīng)的雷擊概率,并利用式(2)對(duì)所得數(shù)據(jù)進(jìn)行線性擬合。表2、表3分別給出了這兩個(gè)初始電位對(duì)應(yīng)的擬合公式中的各參數(shù)數(shù)值以及擬合優(yōu)度,并且根據(jù)擬合公式繪制不同高度建筑物的雷擊概率曲線,如圖7所示。

    表2 DL初始電位-20 MV,雷擊概率擬合曲線公式各參數(shù)值及擬合優(yōu)度計(jì)算值Table 2 DL initial potential -20 MV,lightning strike probability fitting curve formula each parameter value and goodness of fit calculation value

    表3 DL初始電位-30 MV,雷擊概率擬合曲線公式各參數(shù)值及擬合優(yōu)度計(jì)算值Table 3 DL initial potential -30MV,lightning strike probability fitting curve formula parameter values and goodness of fit calculation values

    +表示當(dāng)DL起始位置與建筑物水平距離為 0,100,200,…,600 m 時(shí)不同高度建筑物雷擊概率的具體數(shù)值圖7 DL電位為-20 MV和-30 MV時(shí),不同高度建筑物的雷擊概率隨DL起始點(diǎn)水平距離的變化曲線Fig.7 Lightning strike probability of buildings with different heights varies with horizontal distance of DL starting point when DL potential is -20 MV and -30 MV

    圖7展示了不同閃電強(qiáng)度下的雷擊概率曲線,各圖的變化趨勢(shì)與圖5的變化趨勢(shì)相同。對(duì)比圖7(a)和圖7(b)可知,當(dāng)建筑物在300 m以下,圖7(a)中的每個(gè)初始位置對(duì)應(yīng)的雷擊概率都小于圖7(b)中的雷擊概率,在建筑物高度小于300 m以下時(shí)尤其明顯。而400~600 m建筑物的雷擊概率曲線明顯比圖7(b)中的曲線更早衰減。

    根據(jù)式(1)分別計(jì)算DL初始電位為-20 MV及-30 MV條件下的建筑物年雷擊次數(shù),并如圖8進(jìn)行線性擬合,得出利用建筑物高度計(jì)算的年雷擊次數(shù)N的經(jīng)驗(yàn)公式。式(5)、式(6)分別是當(dāng)DL初始電位為-20 MV和-30 MV時(shí),特定建筑物高度的年雷擊次數(shù)N′的經(jīng)驗(yàn)公式(適用100~600 m高的建筑物)。

    +表示各個(gè)高度建筑物的年雷擊次數(shù)圖8 DL初始電位為-20、-30、-40 MV,在高度1.5 km處始發(fā)時(shí),且閃電年平均密度為1次/km2/a時(shí),不同高度建筑物的年雷擊次數(shù)Fig.8 Annual lightning strikes of buildings at different heights when DL initial potentials are -20,-30,-40 MV,starting at a height of 1.5 km,and the annual average density of lightning is 1 strike/km2/a

    N′=3.57×10-7h2+1×10-3h+0.02

    (5)

    N′=4.82×10-7h2+1.57×10-3h+0.08

    (6)

    如圖8所示,當(dāng)DL初始電位為-20 MV時(shí),100~600 m的建筑物的年雷擊次數(shù)分別為0.12、0.23、0.36、0.47、0.6、0.75,其中600 m建筑物年雷擊次數(shù)是100 m建筑物的6.25倍。當(dāng)DL初始電位為-30 MV時(shí),各高度建筑物對(duì)應(yīng)的年雷擊次數(shù)為0.22、0.43、0.65、0.77、0.92、1.24,其中600 m建筑物年雷擊次數(shù)是100 m的5.6倍。隨著閃電強(qiáng)度的減弱,高建筑物的年雷擊次數(shù)與矮建筑物的年雷擊次數(shù)比值不斷增大,這表示低矮建筑物的年雷擊次數(shù)衰減的程度更大。三組電位對(duì)應(yīng)的年雷擊次數(shù)均與建筑物高度呈正相關(guān),且隨著建筑物高度增加,年雷擊次數(shù)遞增速率不斷增大。隨著閃電強(qiáng)度的減弱,年雷擊次數(shù)隨建筑物高度增加,其增速明顯減小。

    通過(guò)上述分析,不難發(fā)現(xiàn)建筑物高度和閃電強(qiáng)度都是影響年雷擊次數(shù)的重要因子,二者與年雷擊次數(shù)之間呈現(xiàn)明顯的正相關(guān)。為了體現(xiàn)上述因子對(duì)年雷擊次數(shù)的影響,從而得出一個(gè)改進(jìn)的建筑物年雷擊次數(shù)計(jì)算公式,對(duì)式(4)~式(6)進(jìn)行整合。即在不改變建筑物年雷擊次數(shù)變化趨勢(shì)的前提下,將閃電的強(qiáng)度引入公式。式(7)為利用DL初始電位和建筑物高度h共同計(jì)算的建筑物年雷擊次數(shù)經(jīng)驗(yàn)公式。

    N′=(0.17x+8)×10-7h2-(0.06x+0.18)×10-3h-(0.007x+0.12)

    (7)

    式(7)中:x為DL初始電位(適用閃電初始電位為-20~-40 MV);h為建筑物高度(適用于100~600 m建筑物)。

    3 結(jié)論

    (1)基于三維多上行先導(dǎo)模型的模擬,可知建筑物高度,閃電的始發(fā)位置以及閃電的強(qiáng)度都是影響建筑物雷擊概率的重要因素,隨著閃電始發(fā)位置遠(yuǎn)離建筑物,不同高度的建筑物呈現(xiàn)出不同的衰減趨勢(shì)。

    (2)通過(guò)對(duì)不同高度建筑物的雷擊概率進(jìn)行積分,得出了利用建筑物高度計(jì)算的建筑物年雷擊次數(shù)經(jīng)驗(yàn)公式。本文計(jì)算結(jié)果與觀測(cè)結(jié)果較為接近。相較于國(guó)標(biāo)中的公式,更多考慮了閃電本身對(duì)建筑物年雷擊次數(shù)的影響。改變閃電強(qiáng)度進(jìn)行模擬計(jì)算后,得出建筑物的高度和閃電強(qiáng)度都是影響建筑物年雷擊次數(shù)的重要因素。

    (3)本文模型利用閃電初始電位判斷閃電強(qiáng)度,而觀測(cè)中利用回?fù)舴逯惦娏髋袛?如何結(jié)合回?fù)舴逯惦娏髋袛嗄P椭械拈W電強(qiáng)度,需要通過(guò)更多的觀測(cè)資料來(lái)完善模型中的參數(shù),使模擬更接近真實(shí)的雷擊情況。

    猜你喜歡
    電位次數(shù)建筑物
    機(jī)場(chǎng)航站樓年雷擊次數(shù)計(jì)算
    電位滴定法在食品安全檢測(cè)中的應(yīng)用
    2020年,我國(guó)汽車(chē)召回次數(shù)同比減少10.8%,召回?cái)?shù)量同比增長(zhǎng)3.9%
    鄰近既有建筑物全套管回轉(zhuǎn)鉆機(jī)拔樁技術(shù)
    一類無(wú)界算子的二次數(shù)值域和譜
    描寫(xiě)建筑物的詞語(yǔ)
    依據(jù)“次數(shù)”求概率
    電鍍廢水處理中的氧化還原電位控制
    淺談等電位聯(lián)結(jié)
    火柴游戲
    少妇 在线观看| 女同久久另类99精品国产91| 欧美亚洲日本最大视频资源| 国产一卡二卡三卡精品| 亚洲国产精品成人综合色| 中文字幕人妻熟女乱码| 亚洲人成网站高清观看| 少妇 在线观看| 美女免费视频网站| 亚洲精品中文字幕一二三四区| 在线观看日韩欧美| 成人一区二区视频在线观看| 黄色 视频免费看| 国产不卡一卡二| 视频区欧美日本亚洲| 一区二区三区精品91| 亚洲男人的天堂狠狠| 美国免费a级毛片| 国产高清视频在线播放一区| 亚洲精品美女久久久久99蜜臀| cao死你这个sao货| 亚洲avbb在线观看| 天天一区二区日本电影三级| 国产成人系列免费观看| 91成年电影在线观看| 久久性视频一级片| 50天的宝宝边吃奶边哭怎么回事| 欧美三级亚洲精品| 亚洲第一青青草原| 成年人黄色毛片网站| 色尼玛亚洲综合影院| 99国产精品一区二区三区| 12—13女人毛片做爰片一| 免费av毛片视频| 满18在线观看网站| 亚洲一区高清亚洲精品| 欧美性长视频在线观看| 日本 av在线| 黄色视频,在线免费观看| 久久久国产成人免费| 99热这里只有精品一区 | 淫秽高清视频在线观看| bbb黄色大片| 最近最新中文字幕大全电影3 | 日韩视频一区二区在线观看| 十分钟在线观看高清视频www| 99国产精品99久久久久| 韩国精品一区二区三区| 日韩视频一区二区在线观看| e午夜精品久久久久久久| 99精品在免费线老司机午夜| 久久婷婷成人综合色麻豆| 黄片小视频在线播放| 欧美在线一区亚洲| 高清在线国产一区| 动漫黄色视频在线观看| 国产午夜精品久久久久久| 欧美乱妇无乱码| 国产高清videossex| 成人国产综合亚洲| 男女做爰动态图高潮gif福利片| 久久精品国产亚洲av高清一级| 精品熟女少妇八av免费久了| 国产精品一区二区免费欧美| 999精品在线视频| 欧美激情 高清一区二区三区| 18美女黄网站色大片免费观看| 91麻豆av在线| 婷婷亚洲欧美| 久久99热这里只有精品18| 久久久久久久精品吃奶| 色尼玛亚洲综合影院| 美女国产高潮福利片在线看| 国产欧美日韩一区二区三| 男女之事视频高清在线观看| 老司机午夜福利在线观看视频| 亚洲国产日韩欧美精品在线观看 | 99精品久久久久人妻精品| 我的亚洲天堂| 少妇粗大呻吟视频| 久久香蕉国产精品| 美女 人体艺术 gogo| 看片在线看免费视频| 亚洲精品中文字幕在线视频| 19禁男女啪啪无遮挡网站| 免费看美女性在线毛片视频| 精品久久久久久成人av| 黑人巨大精品欧美一区二区mp4| 男女做爰动态图高潮gif福利片| 69av精品久久久久久| 午夜福利成人在线免费观看| 女性被躁到高潮视频| 午夜a级毛片| 精品免费久久久久久久清纯| 校园春色视频在线观看| 黄网站色视频无遮挡免费观看| 一进一出好大好爽视频| 久久久久九九精品影院| 久久青草综合色| 欧美一区二区精品小视频在线| 久久久水蜜桃国产精品网| 在线观看免费午夜福利视频| 老司机深夜福利视频在线观看| 久久久国产欧美日韩av| 精品久久久久久久末码| 色综合欧美亚洲国产小说| 一本精品99久久精品77| 久久香蕉精品热| 好男人电影高清在线观看| 亚洲成人国产一区在线观看| 女警被强在线播放| 欧美精品啪啪一区二区三区| 99热只有精品国产| 中文字幕精品亚洲无线码一区 | 欧美激情 高清一区二区三区| 精品一区二区三区视频在线观看免费| 亚洲av成人av| 热99re8久久精品国产| 波多野结衣巨乳人妻| 91字幕亚洲| 欧美性猛交黑人性爽| 最新在线观看一区二区三区| 无限看片的www在线观看| 日本免费一区二区三区高清不卡| 丝袜美腿诱惑在线| 天堂影院成人在线观看| 伊人久久大香线蕉亚洲五| a级毛片a级免费在线| 亚洲免费av在线视频| 超碰成人久久| 国内揄拍国产精品人妻在线 | 午夜免费鲁丝| 午夜福利在线在线| 99国产精品一区二区蜜桃av| 九色国产91popny在线| 亚洲电影在线观看av| 老汉色av国产亚洲站长工具| 性欧美人与动物交配| 老熟妇乱子伦视频在线观看| 在线天堂中文资源库| 亚洲欧美激情综合另类| 精品少妇一区二区三区视频日本电影| 久久久久久久午夜电影| 久久久久国产一级毛片高清牌| 黄色成人免费大全| 国产精品爽爽va在线观看网站 | 黄色 视频免费看| 欧美性猛交╳xxx乱大交人| 亚洲第一欧美日韩一区二区三区| 黄网站色视频无遮挡免费观看| 女生性感内裤真人,穿戴方法视频| 国产视频内射| 日本黄色视频三级网站网址| 日韩欧美免费精品| 亚洲人成网站高清观看| 成人三级黄色视频| 成人手机av| 啦啦啦韩国在线观看视频| 波多野结衣高清无吗| 极品教师在线免费播放| 欧美久久黑人一区二区| 男女午夜视频在线观看| 中国美女看黄片| 丁香六月欧美| 男女做爰动态图高潮gif福利片| 欧美乱妇无乱码| 热re99久久国产66热| 国产又爽黄色视频| 午夜免费激情av| 一级毛片高清免费大全| 国产主播在线观看一区二区| 久久久精品国产亚洲av高清涩受| 国产精品99久久99久久久不卡| 岛国视频午夜一区免费看| 黄色片一级片一级黄色片| 99国产精品一区二区三区| 国产精品,欧美在线| 亚洲人成电影免费在线| 手机成人av网站| 日本精品一区二区三区蜜桃| 热99re8久久精品国产| 亚洲va日本ⅴa欧美va伊人久久| 97人妻精品一区二区三区麻豆 | 亚洲一区二区三区不卡视频| 一进一出好大好爽视频| 国产视频内射| 国产99白浆流出| 亚洲狠狠婷婷综合久久图片| 中文在线观看免费www的网站 | 日韩视频一区二区在线观看| 亚洲第一电影网av| 少妇粗大呻吟视频| 久久亚洲真实| 热re99久久国产66热| 国产又色又爽无遮挡免费看| 一本一本综合久久| www日本在线高清视频| 黄色女人牲交| 午夜激情福利司机影院| 亚洲精品一卡2卡三卡4卡5卡| 国产午夜精品久久久久久| 成人三级做爰电影| 久热爱精品视频在线9| 国产私拍福利视频在线观看| 首页视频小说图片口味搜索| 老司机在亚洲福利影院| 深夜精品福利| 欧美黑人欧美精品刺激| 国产精品亚洲av一区麻豆| 一区二区三区国产精品乱码| 日韩 欧美 亚洲 中文字幕| 99在线视频只有这里精品首页| 啦啦啦韩国在线观看视频| 久久久久久九九精品二区国产 | 一夜夜www| 亚洲一区高清亚洲精品| 国产精品98久久久久久宅男小说| 国产一卡二卡三卡精品| 91麻豆av在线| 一区二区三区激情视频| 久久精品影院6| 欧美日韩亚洲国产一区二区在线观看| 熟女电影av网| 欧美激情极品国产一区二区三区| 伦理电影免费视频| 亚洲成人久久性| 97碰自拍视频| 欧美性猛交╳xxx乱大交人| www国产在线视频色| 99国产精品一区二区三区| 久久人人精品亚洲av| 99riav亚洲国产免费| 制服人妻中文乱码| 亚洲精品中文字幕在线视频| 中国美女看黄片| 国产亚洲精品综合一区在线观看 | 国产精品二区激情视频| 可以在线观看的亚洲视频| 中文亚洲av片在线观看爽| 国产精品久久久人人做人人爽| 免费看十八禁软件| 国产日本99.免费观看| 欧美乱码精品一区二区三区| 亚洲三区欧美一区| 亚洲国产欧洲综合997久久, | 黄片小视频在线播放| 中文字幕人妻熟女乱码| 美女午夜性视频免费| 色播亚洲综合网| 国产久久久一区二区三区| 少妇熟女aⅴ在线视频| 国产蜜桃级精品一区二区三区| 免费看十八禁软件| 真人做人爱边吃奶动态| 亚洲三区欧美一区| 69av精品久久久久久| 中文在线观看免费www的网站 | 国产aⅴ精品一区二区三区波| 日本五十路高清| 午夜亚洲福利在线播放| 久久精品人妻少妇| 午夜福利免费观看在线| 好男人在线观看高清免费视频 | 1024香蕉在线观看| 国产色视频综合| 成人永久免费在线观看视频| 色综合亚洲欧美另类图片| 欧美黄色淫秽网站| 在线观看免费视频日本深夜| 亚洲欧美一区二区三区黑人| 国产精品电影一区二区三区| 国产精品久久久av美女十八| 日韩高清综合在线| 一区二区三区国产精品乱码| 免费在线观看完整版高清| 国产真实乱freesex| 午夜福利在线观看吧| 青草久久国产| 桃色一区二区三区在线观看| a级毛片a级免费在线| 色播亚洲综合网| 久久中文字幕人妻熟女| 好男人在线观看高清免费视频 | 亚洲自偷自拍图片 自拍| 国产亚洲av高清不卡| 国产99白浆流出| www日本黄色视频网| 免费高清在线观看日韩| 男女床上黄色一级片免费看| 亚洲熟妇中文字幕五十中出| 国产伦人伦偷精品视频| 给我免费播放毛片高清在线观看| 国产成人av教育| 欧美日韩亚洲国产一区二区在线观看| 免费观看精品视频网站| 亚洲中文日韩欧美视频| 在线观看www视频免费| 老汉色∧v一级毛片| 亚洲第一av免费看| 国产蜜桃级精品一区二区三区| 亚洲国产欧洲综合997久久, | 亚洲国产精品久久男人天堂| 国产高清videossex| 啦啦啦观看免费观看视频高清| 夜夜躁狠狠躁天天躁| 色av中文字幕| 亚洲免费av在线视频| 国产麻豆成人av免费视频| 久久久久久亚洲精品国产蜜桃av| 人妻久久中文字幕网| 国产精品 欧美亚洲| 国产又黄又爽又无遮挡在线| 99精品在免费线老司机午夜| 50天的宝宝边吃奶边哭怎么回事| svipshipincom国产片| 久久久久久久久免费视频了| 中文字幕久久专区| 男男h啪啪无遮挡| 两性午夜刺激爽爽歪歪视频在线观看 | 精品不卡国产一区二区三区| 少妇 在线观看| 国产在线观看jvid| av欧美777| 少妇被粗大的猛进出69影院| 国产成人精品久久二区二区91| 久久热在线av| 一二三四在线观看免费中文在| 午夜福利18| 啪啪无遮挡十八禁网站| 日本撒尿小便嘘嘘汇集6| 久久国产亚洲av麻豆专区| 亚洲av中文字字幕乱码综合 | 麻豆久久精品国产亚洲av| 亚洲男人天堂网一区| or卡值多少钱| 亚洲色图 男人天堂 中文字幕| 免费在线观看黄色视频的| √禁漫天堂资源中文www| 中文资源天堂在线| 亚洲午夜理论影院| 午夜福利欧美成人| 国产一区二区三区在线臀色熟女| 精品久久久久久久末码| 国产国语露脸激情在线看| 中文字幕人妻熟女乱码| 日日干狠狠操夜夜爽| 亚洲欧美激情综合另类| 亚洲国产中文字幕在线视频| 日韩中文字幕欧美一区二区| 久久久精品国产亚洲av高清涩受| 变态另类成人亚洲欧美熟女| 在线观看www视频免费| 波多野结衣av一区二区av| 国产黄a三级三级三级人| 一级片免费观看大全| 校园春色视频在线观看| 又黄又粗又硬又大视频| 亚洲自拍偷在线| 日韩大码丰满熟妇| 亚洲欧美激情综合另类| x7x7x7水蜜桃| 国产熟女xx| 久久精品人妻少妇| 99国产精品一区二区蜜桃av| 在线看三级毛片| 午夜免费鲁丝| 97碰自拍视频| 国产av不卡久久| 日韩av在线大香蕉| 最新美女视频免费是黄的| 欧美性长视频在线观看| www.自偷自拍.com| 国产三级在线视频| 久久久久久大精品| 久久久久九九精品影院| 中文字幕av电影在线播放| 国产精品一区二区精品视频观看| 午夜福利18| 夜夜爽天天搞| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看视频国产中文字幕亚洲| 婷婷六月久久综合丁香| 嫩草影院精品99| a在线观看视频网站| 热re99久久国产66热| 日日摸夜夜添夜夜添小说| www.熟女人妻精品国产| 桃红色精品国产亚洲av| 两性夫妻黄色片| 后天国语完整版免费观看| 国产精品久久视频播放| 日韩av在线大香蕉| а√天堂www在线а√下载| 久久久久久久午夜电影| 日韩欧美一区视频在线观看| 亚洲精品久久国产高清桃花| 性欧美人与动物交配| 日日夜夜操网爽| 日韩精品中文字幕看吧| 99re在线观看精品视频| 国产精品自产拍在线观看55亚洲| 精品国产乱码久久久久久男人| 国产精华一区二区三区| 人人妻人人澡人人看| 成人18禁在线播放| 香蕉久久夜色| 国产一级毛片七仙女欲春2 | 一边摸一边做爽爽视频免费| 亚洲电影在线观看av| 村上凉子中文字幕在线| 免费搜索国产男女视频| 99久久国产精品久久久| 怎么达到女性高潮| 国产精品1区2区在线观看.| 国产人伦9x9x在线观看| 亚洲成国产人片在线观看| 国产亚洲av嫩草精品影院| 免费在线观看日本一区| 一本大道久久a久久精品| 中文亚洲av片在线观看爽| 男人舔奶头视频| 岛国视频午夜一区免费看| 欧美黑人巨大hd| 麻豆av在线久日| 欧美黑人欧美精品刺激| 天天添夜夜摸| tocl精华| 欧美zozozo另类| 天堂影院成人在线观看| 看黄色毛片网站| 麻豆国产av国片精品| 色av中文字幕| 国产精品电影一区二区三区| 成人一区二区视频在线观看| 欧美性长视频在线观看| 动漫黄色视频在线观看| 人人妻人人澡欧美一区二区| 他把我摸到了高潮在线观看| 亚洲成人精品中文字幕电影| 亚洲激情在线av| 欧美激情高清一区二区三区| 麻豆av在线久日| 色精品久久人妻99蜜桃| 国产又黄又爽又无遮挡在线| 韩国av一区二区三区四区| 免费在线观看黄色视频的| 亚洲午夜理论影院| 老司机在亚洲福利影院| 天堂影院成人在线观看| 免费观看精品视频网站| 黄色丝袜av网址大全| 成人一区二区视频在线观看| 国产黄片美女视频| 久久这里只有精品19| 亚洲成人久久爱视频| 黄色a级毛片大全视频| 久久国产精品影院| 亚洲自拍偷在线| 亚洲熟妇中文字幕五十中出| tocl精华| 大型av网站在线播放| 国产成人啪精品午夜网站| 久久久久国产一级毛片高清牌| 精品无人区乱码1区二区| 国产精品一区二区免费欧美| 老司机深夜福利视频在线观看| 手机成人av网站| 一a级毛片在线观看| 国产精品日韩av在线免费观看| 亚洲激情在线av| 国产精品久久视频播放| 性欧美人与动物交配| 免费一级毛片在线播放高清视频| 999精品在线视频| 国产精品免费一区二区三区在线| 九色国产91popny在线| 国语自产精品视频在线第100页| 欧美日韩福利视频一区二区| 欧美绝顶高潮抽搐喷水| 亚洲全国av大片| 亚洲国产看品久久| 欧美激情高清一区二区三区| 黄色片一级片一级黄色片| 日本一本二区三区精品| 午夜视频精品福利| 亚洲熟妇熟女久久| 操出白浆在线播放| 亚洲专区国产一区二区| 国产精品久久久久久精品电影 | 国产麻豆成人av免费视频| 91老司机精品| 亚洲人成伊人成综合网2020| 婷婷精品国产亚洲av| avwww免费| 久久性视频一级片| 国产精品永久免费网站| 精品欧美国产一区二区三| 99热6这里只有精品| 一本精品99久久精品77| 12—13女人毛片做爰片一| 午夜福利免费观看在线| 亚洲成人久久爱视频| 99国产精品99久久久久| 少妇粗大呻吟视频| 色综合站精品国产| 国产又黄又爽又无遮挡在线| 久久精品国产综合久久久| 在线永久观看黄色视频| 丝袜人妻中文字幕| 非洲黑人性xxxx精品又粗又长| 淫妇啪啪啪对白视频| 亚洲午夜精品一区,二区,三区| 国产高清激情床上av| 51午夜福利影视在线观看| 看片在线看免费视频| 一区二区三区高清视频在线| 午夜福利欧美成人| 一级片免费观看大全| 亚洲人成网站高清观看| 99re在线观看精品视频| 男女视频在线观看网站免费 | 午夜精品在线福利| 日韩欧美国产在线观看| 欧美三级亚洲精品| 亚洲精品美女久久久久99蜜臀| 精品日产1卡2卡| 亚洲五月色婷婷综合| 国产成人精品久久二区二区91| 级片在线观看| 午夜日韩欧美国产| 午夜免费成人在线视频| 久久精品aⅴ一区二区三区四区| 在线观看午夜福利视频| 亚洲欧洲精品一区二区精品久久久| 国产av在哪里看| 亚洲欧美一区二区三区黑人| 国产精品九九99| 级片在线观看| 亚洲专区中文字幕在线| av有码第一页| 国产97色在线日韩免费| 搡老熟女国产l中国老女人| 丝袜在线中文字幕| 天天添夜夜摸| 老司机福利观看| 女性被躁到高潮视频| 精品国产超薄肉色丝袜足j| 国产一区二区激情短视频| 欧美丝袜亚洲另类 | 亚洲av第一区精品v没综合| 午夜激情福利司机影院| 色婷婷久久久亚洲欧美| 一a级毛片在线观看| 香蕉国产在线看| 亚洲男人天堂网一区| 99精品久久久久人妻精品| 亚洲精品av麻豆狂野| av在线播放免费不卡| 老鸭窝网址在线观看| 制服诱惑二区| 一本精品99久久精品77| 亚洲全国av大片| 久久久久久国产a免费观看| 亚洲成a人片在线一区二区| 国产成人精品久久二区二区免费| 精华霜和精华液先用哪个| 国产成人精品久久二区二区91| 国产精品九九99| 欧美日韩一级在线毛片| 国产免费男女视频| 成在线人永久免费视频| 欧美一级毛片孕妇| av免费在线观看网站| 亚洲熟妇熟女久久| 国产成人一区二区三区免费视频网站| 午夜日韩欧美国产| 久久精品国产亚洲av香蕉五月| 97超级碰碰碰精品色视频在线观看| 夜夜躁狠狠躁天天躁| 91麻豆精品激情在线观看国产| 亚洲一码二码三码区别大吗| 黄片大片在线免费观看| 久久青草综合色| 久久人人精品亚洲av| 后天国语完整版免费观看| 色综合婷婷激情| 欧美一区二区精品小视频在线| 亚洲性夜色夜夜综合| av有码第一页| 久久久久久久久中文| 亚洲精品国产区一区二| 制服诱惑二区| 中文字幕高清在线视频| 无限看片的www在线观看| 亚洲第一青青草原| 国产成人精品久久二区二区91| 91大片在线观看| 国产久久久一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 18美女黄网站色大片免费观看| av片东京热男人的天堂| 国产亚洲av高清不卡| 一区二区三区精品91| 亚洲自偷自拍图片 自拍| 啦啦啦免费观看视频1| 女警被强在线播放| 变态另类成人亚洲欧美熟女| 麻豆成人av在线观看| 亚洲中文av在线| 午夜成年电影在线免费观看| 日韩欧美 国产精品| 十八禁人妻一区二区| 成人18禁高潮啪啪吃奶动态图| 男女做爰动态图高潮gif福利片| 在线观看日韩欧美|