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

    利用毫米波云雷達(dá)數(shù)據(jù)反演層云微物理參數(shù)和云內(nèi)湍流耗散率

    2020-12-11 09:09:15黃興友陸琳洪滔梅垚楊敏
    大氣科學(xué)學(xué)報(bào) 2020年5期
    關(guān)鍵詞:反演

    黃興友 陸琳 洪滔 梅垚 楊敏

    摘要 利用地基毫米波雷達(dá)進(jìn)行云參數(shù)及云內(nèi)湍流特性的探測(cè)和反演。根據(jù)云雷達(dá)回波的功率譜數(shù)據(jù),反演出大氣垂直運(yùn)動(dòng)速度和云微物理參數(shù),得到云內(nèi)湍流耗散率ε的大小和分布情況,并進(jìn)一步研究和分析了云內(nèi)空氣垂直運(yùn)動(dòng)與云微物理參數(shù)、反射率因子、多普勒速度、速度譜寬變化的關(guān)系,更好地了解云的演變情況。對(duì)2016年8月8日四川稻城的一次層狀云過(guò)程的探測(cè)和反演表明:1)粒子有效半徑隨著上升氣流的增強(qiáng)而增大,由于碰并聚合的作用,粒子數(shù)濃度也呈現(xiàn)相應(yīng)減小的趨勢(shì)。2)云內(nèi)湍流耗散率ε在云底、云頂較大,云內(nèi)較小,量級(jí)在10-8~10-2 m2·s-3,多普勒速度能譜驗(yàn)證了假設(shè)雷達(dá)探測(cè)湍渦的尺度在慣性副區(qū)的合理性。3)隨著上升氣流的增強(qiáng),云粒子的下沉運(yùn)動(dòng)相應(yīng)減小,速度譜寬相應(yīng)增大。

    關(guān)鍵詞 毫米波雷達(dá);湍流耗散率;空氣垂直運(yùn)動(dòng);云微物理參數(shù);反演

    在地球-大氣系統(tǒng)中,云是調(diào)節(jié)氣候的重要因素。云的形成、結(jié)構(gòu)和生命周期與湍流緊密相關(guān)(Feingold et al.,1999)。湍流不僅加強(qiáng)了云內(nèi)的混合,還因卷夾作用,引起地面熱量、動(dòng)量、水汽和物質(zhì)的向上輸送,可以為成云區(qū)提供水汽和云凝結(jié)核,對(duì)云的生長(zhǎng)有促進(jìn)作用,進(jìn)而改變?cè)铺匦院徒邓植记闆r。所以云內(nèi)湍流能改變?cè)苾?nèi)的溫-濕環(huán)境,影響云的生消過(guò)程或云的生命周期(Bouniol et al.,2004;Wang et al.,2006)。

    由于云在高空,直接探測(cè)云微物理參數(shù)(彭沖等,2016;楊文霞等,2018)和湍流比較困難,利用地基毫米波云雷達(dá)(邱玉琚等,2012),進(jìn)行云參數(shù)和云內(nèi)湍流耗散率的遙感探測(cè)是比較有效的途徑。除了云高云厚等宏觀信息外,根據(jù)毫米波雷達(dá)回波數(shù)據(jù),能夠反演云內(nèi)粒子大小、數(shù)濃度、冰水含量等云微物理參數(shù)(樊雅文等,2013;黃興友等,2019),因而在觀測(cè)非降水云和弱降水云方面有獨(dú)特的優(yōu)勢(shì)(Hobbs and Funk,1984;Kollias et al.,2007;宗蓉等,2014)。毫米波雷達(dá)回波信號(hào)的功率譜中包含云粒子后向散射能量和速度信息,可以用來(lái)提取云粒子下落速度、空氣垂直運(yùn)動(dòng)速度以及云內(nèi)湍流耗散率(Frisch et al.,1995;Shupeet al.,2008a,2008b,2012),可用于云物理過(guò)程和天氣演變的研究。

    云滴的自由下落速度是反演云微物理參數(shù)的關(guān)鍵參量之一,Battan(1964)提出速度譜低端法反演云粒子下落速度的方案,其核心思想是云粒子的下落速度等于垂直指向雷達(dá)獲取的多普勒速度與空氣垂直運(yùn)動(dòng)速度的差。Shupe et al.(2008a)在速度譜低端法的基礎(chǔ)上,利用云雷達(dá)功率譜數(shù)據(jù)并結(jié)合小粒子示蹤法,反演了空氣垂直運(yùn)動(dòng)速度,并認(rèn)為云粒子下落速度本身對(duì)譜寬起了決定性作用。彭亮等(2012)結(jié)合小粒子示蹤法以及云粒子下降速度W0與反射率因子Z的關(guān)系、反演了云內(nèi)空氣垂直運(yùn)動(dòng)速度,進(jìn)而可計(jì)算云粒子的下落速度。段藝萍等(2013)通過(guò)雷達(dá)功率譜數(shù)據(jù)估算層狀云內(nèi)空氣垂直運(yùn)動(dòng)速度,在此基礎(chǔ)上,反演得到了層狀云的微物理參數(shù)。

    關(guān)于云內(nèi)湍流耗散率,Bouniol et al.(2004)利用94 GHz地基雷達(dá)較高時(shí)間分辨率的多普勒速度資料,反演了云內(nèi)湍流耗散率,發(fā)現(xiàn)不同云區(qū)的耗散率有差異:在云底較大,約為10-4 m2·s-3;在云內(nèi)較小,約為10-8 m2·s-3,云頂略有增大。O'Connor et al.(2005)利用94 GHz毫米波雷達(dá)和激光雷達(dá)對(duì)降水層狀云進(jìn)行觀測(cè),在Bouniol et al.(2004)湍流耗散率方法的基礎(chǔ)上,獲得了云層的液水通量和雷達(dá)反射率因子之間的冪指數(shù)關(guān)系。Shupe et al.(2008b)利用35 GHz毫米波雷達(dá)對(duì)北極的混合相層狀云進(jìn)行觀測(cè),研究云內(nèi)空氣垂直運(yùn)動(dòng)對(duì)云微物理特性的影響,并反演了云內(nèi)湍流。

    本文參考小粒子示蹤法(Shupe et al.,2008a),利用Ka波段云雷達(dá)探測(cè)的云回波多普勒譜(功率譜)數(shù)據(jù),反演云微物理參數(shù),并借鑒Bouniol et al.(2004)的方法計(jì)算湍流耗散率,對(duì)云內(nèi)湍流及云特性參數(shù)進(jìn)行初步分析研究。

    1 毫米波雷達(dá)與觀測(cè)資料

    2016年7月至8月,中國(guó)氣象科學(xué)研究院的Ka波段云雷達(dá)在四川稻城進(jìn)行云特性的觀測(cè),回波數(shù)據(jù)包括功率譜數(shù)據(jù)和基數(shù)據(jù),本文選取8月8日的層狀云過(guò)程資料進(jìn)行研究。

    該雷達(dá)使用小功率的固態(tài)發(fā)射機(jī)體制,為增強(qiáng)其探測(cè)能力,采用脈沖壓縮技術(shù),能探測(cè)到對(duì)流層以內(nèi)的云目標(biāo),探測(cè)高度可達(dá)20 km。雷達(dá)采用垂直指向探測(cè)的方案,所以,回波信號(hào)中包含云粒子和大氣的垂直運(yùn)動(dòng)信息。雷達(dá)的主要工作參數(shù)如表1所示。

    2 反演方法

    2.1 毫米波雷達(dá)功率譜數(shù)據(jù)處理

    回波的功率譜數(shù)據(jù),是通過(guò)對(duì)云雷達(dá)的I、Q數(shù)據(jù)進(jìn)行相干積分、去直流加窗處理、快速傅里葉變換和譜平均等基本處理后得到的,可用于提取目標(biāo)粒子特性參數(shù)、粒子垂直運(yùn)動(dòng)速度和大氣垂直運(yùn)動(dòng)速度。

    回波信號(hào)中不可避免地含有噪聲,假設(shè)噪聲服從自由度為2N-FFT/K的中心χ^2分布,將頻譜平均分為K段,每一段數(shù)據(jù)個(gè)數(shù)為N-FFT/K,N-FFT為FFT點(diǎn)數(shù)。對(duì)每一段數(shù)據(jù)取平均,其中的最小值可視為整個(gè)譜的平均噪聲功率。研究表明,劃分為8段或16段最合適(Petitdidier et al.,1997;周旭輝,2011),本文中K取8,處理得到噪聲電平(Noise level)。如果直接對(duì)噪聲電平以上信號(hào)進(jìn)行積分等處理,可以得到云的回波功率、多普勒速度、譜寬(Woodman,1985;宗蓉等,2016)基數(shù)據(jù),但有一定的誤差,尤其在信噪比較低的情況下,鄭佳鋒(2016)認(rèn)為,應(yīng)該先進(jìn)行云回波信號(hào)的識(shí)別,然后再對(duì)云回波信號(hào)進(jìn)行積分等處理,獲得的基數(shù)據(jù)才更準(zhǔn)確。識(shí)別云回波信號(hào)的基本思路是:先找到最大譜峰,再與噪聲電平比較,找到噪聲電平與最大譜峰交點(diǎn)外端區(qū)域的最大值,將其作為云回波信號(hào)和噪聲的分界線(Boundary of Signal and Noise),參見圖1,并參考Shupe et al.(2008a),云回波信號(hào)應(yīng)該滿足至少有7個(gè)連續(xù)譜點(diǎn)的功率在噪聲電平以上。根據(jù)這個(gè)標(biāo)準(zhǔn),可以識(shí)別出云回波信號(hào),再對(duì)云回波譜數(shù)據(jù)進(jìn)行積分,計(jì)算出云信號(hào)的譜矩和基數(shù)據(jù)。其中,回波功率PR和信噪比RSN為:

    2.2 云微物理參數(shù)的反演

    云粒子的大小與云粒子的自然降落速度相關(guān),垂直指向的云雷達(dá)測(cè)量的徑向速度中包括云粒子的自然降落速度和大氣垂直運(yùn)動(dòng)速度。微米量級(jí)尺度的云粒子,自然降落速度很?。? cm/s以內(nèi)),而空氣垂直運(yùn)動(dòng)速度可以比云粒子降落大1~2個(gè)量級(jí)(Kollias et al.,2001)。如果云粒子足夠小,它應(yīng)該是隨大氣運(yùn)動(dòng),因此,云回波功率譜左側(cè)第一個(gè)譜點(diǎn)對(duì)應(yīng)的速度就是小云滴的下落速度,也就是大氣的垂直速度(Shupe et al.,2008a;鄭佳鋒,2016;Siebert et al.,2006),如圖1中的Wm所示。得到了大氣垂直速度Wm后,利用徑向速度VD和公式(5),可以計(jì)算出云粒子的自然降落速度Vg。

    3.1 估算云內(nèi)垂直速度和湍流耗散率

    圖3a是通過(guò)功率譜數(shù)據(jù)估算的空氣垂直運(yùn)動(dòng)速度,云內(nèi)既有上升也有下沉運(yùn)動(dòng),主要集中在-0.6~1.5 m/s,在較強(qiáng)上升氣流(>0.8 m/s)集中的地方,云頂高度有所增加,云體變厚,較弱的上升氣流對(duì)云頂高度變化的影響并不大,這只是比較微弱的關(guān)系,與Shupe et al.(2008b)的研究結(jié)果一致。圖3b是通過(guò)方程(20)計(jì)算獲得的湍流耗散率,通過(guò)對(duì)數(shù)形式log10(簡(jiǎn)寫為lg)表示,量級(jí)在10-8~10-2 m2·s-3。云底和云頂?shù)暮纳⒙瘦^大,云內(nèi)較小,這與Bouniol et al.(2004)的分析一致。對(duì)比圖3a與圖3b,可以看到,云內(nèi)上升氣流較強(qiáng)的地方,湍流耗散率相應(yīng)增大、湍流強(qiáng)度也較大,對(duì)應(yīng)的譜寬也增大(圖2c)。圖3c為多普勒速度在7.5 km高度的時(shí)間序列,該高度均為下沉氣流,速度大小在-1.5~0 m/s,NCEP資料顯示該時(shí)段的水平風(fēng)速約為7 m/s。

    圖4是計(jì)算出的能譜圖Sv(f),能譜的高頻區(qū)存在近似-5/3的斜率,不僅反映了雷達(dá)已經(jīng)探測(cè)到了湍流的慣性副區(qū),也間接驗(yàn)證了湍流能譜的-5/3指數(shù)特性,符合科爾莫戈羅夫的湍流理論(盛裴軒等,2003;Shupe et al.,2012)。

    3.2 層云微物理參數(shù)的反演

    圖5為云層內(nèi)云粒子有效半徑、液態(tài)水含量、粒子數(shù)濃度的反演結(jié)果。粒子有效半徑主要分布在20~90 μm。大粒子區(qū)對(duì)應(yīng)較強(qiáng)的回波,說(shuō)明粒子大小是影響回波強(qiáng)度的主要因素。在云內(nèi)垂直上升氣流較大的地方,往往對(duì)應(yīng)粒子尺度也大,如果較強(qiáng)上升氣流區(qū)的上部存在下沉氣流,則云滴粒子在反復(fù)抬升和降落過(guò)程中,可能逐漸增長(zhǎng),

    如在22:13—22:49,8 km高度處的垂直上升氣流較強(qiáng),其上部存在下沉氣流,粒子有效半徑較大,但粒子數(shù)濃度很小,可能因碰并作用所致。液態(tài)水含量在0.001~1.9 g/m3,集中在回波強(qiáng)度較強(qiáng)的地方;粒子數(shù)濃度通過(guò)對(duì)數(shù)形式log10(簡(jiǎn)寫為lg)表示,主要集中在105~109 m-3,區(qū)域分布較為均勻。

    3.3 云微物理參數(shù)與云內(nèi)湍流關(guān)系

    圖6反映出云內(nèi)上升氣流速度大、則云粒子大、數(shù)濃度小,符合云內(nèi)空氣垂直運(yùn)動(dòng)與云微物理參數(shù)存在一些較為直觀的聯(lián)系。圖6a中粒子有效半徑隨著上升氣流的增強(qiáng)而增大,在一定程度上說(shuō)明上升氣流對(duì)粒子的生長(zhǎng)及碰并聚合有促進(jìn)作用,與Shupe et al.(2008b)得出的研究結(jié)論一致。圖6c中粒子數(shù)濃度則呈現(xiàn)減小趨勢(shì),這可能是粒子的碰并聚合作用,小粒子結(jié)合成較大粒子,使得粒子數(shù)濃度降低的原因。本文個(gè)例圖6b中液態(tài)水含量隨著上升氣流的增加呈減小趨勢(shì),這與Shupe et al.(2008b)相似研究的結(jié)果存在較大差異,可能是由本文個(gè)例中云體不穩(wěn)定等原因引起。

    圖6d、e、f為空氣垂直運(yùn)動(dòng)與雷達(dá)基本產(chǎn)品的關(guān)系分析,圖6d隨著空氣垂直運(yùn)動(dòng)由下沉運(yùn)動(dòng)轉(zhuǎn)變成上升運(yùn)動(dòng),云粒子的下沉運(yùn)動(dòng)也在相應(yīng)減小,這一部分也是因?yàn)槎嗥绽账俣劝丝諝獯怪边\(yùn)動(dòng)速度。在圖6e中,隨著上升氣流增強(qiáng),速度譜寬相應(yīng)存在增大的趨勢(shì),這很可能是粒子在上升氣流中形成、生長(zhǎng),從而導(dǎo)致了粒子譜分布的展寬。Shupe et al.(2008b)還提出,隨著上升氣流從0 m/s開始增加,強(qiáng)度在-40 dBZ反射率因子也以0.002 m·s-1·(dBZ)-1的速率增大,但本文由于反射率因子強(qiáng)度較強(qiáng),空氣垂直運(yùn)動(dòng)的時(shí)間精度較低,與Shupe et al.(2008b)的研究條件有所差別,空氣垂直運(yùn)動(dòng)和反射率因子變化的關(guān)系并不明顯。

    4 結(jié)論

    通過(guò)毫米波雷達(dá)對(duì)層狀云個(gè)例內(nèi)部湍流的分析和云微物理特征的反演,可以初步得到一些結(jié)論:

    1)利用毫米波雷達(dá)采集的功率譜數(shù)據(jù),通過(guò)分段法計(jì)算噪聲電平,合理選取云信號(hào),能夠較為準(zhǔn)確地計(jì)算出雷達(dá)反射率因子、多普勒速度、速度譜寬等基礎(chǔ)產(chǎn)品,并反演了空氣垂直運(yùn)動(dòng)速度和云微物理參數(shù)。

    2)通過(guò)平均多普勒速度方差計(jì)算得到的湍流耗散率,量級(jí)在10-8~10-2 m2·s-3,耗散率在云底、云頂較大,云內(nèi)較小,耗散率的大值區(qū),往往對(duì)應(yīng)云內(nèi)上升氣流較強(qiáng)的區(qū)域,意味著強(qiáng)上升氣流區(qū)對(duì)應(yīng)的湍流強(qiáng)度大,同時(shí)多普勒速度能譜驗(yàn)證了雷達(dá)探測(cè)湍渦的尺度在慣性副區(qū)假設(shè)的合理性。

    3)粒子有效半徑隨著上升氣流的增強(qiáng)而增大,云滴粒子通過(guò)碰并作用增長(zhǎng),若在較強(qiáng)上升氣流區(qū)上部存在下沉區(qū),由于湍流作用的增強(qiáng),云滴粒子在反復(fù)抬升和降落中能夠更加快速的生長(zhǎng)。由于云滴粒子的碰并作用,隨著上升氣流的增加,粒子數(shù)濃度則呈現(xiàn)減小趨勢(shì);液態(tài)水含量與空氣垂直運(yùn)動(dòng)的關(guān)系還有待進(jìn)一步研究發(fā)現(xiàn)。

    4)隨著上升氣流的增強(qiáng),云粒子的下沉運(yùn)動(dòng)相應(yīng)減小,速度譜寬也相應(yīng)有所增大,但本次個(gè)例尚未發(fā)現(xiàn)上升運(yùn)動(dòng)與反射率因子較為直觀的關(guān)系。

    云微物理參數(shù)、云內(nèi)湍流耗散率以及空氣垂直運(yùn)動(dòng),都受到云頂輻射冷卻、風(fēng)切變、地表湍流熱通量、潛熱冷卻、云頂干空氣的拖拽作用等影響。而云內(nèi)湍流對(duì)云微物理特性的影響,也并非是單一的線性關(guān)系,這與云的類型、發(fā)展階段、溫濕條件、地區(qū)等都有很大的關(guān)系。因此,本文的研究方案和結(jié)果,還需要更多的觀測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證,以便深入研究云內(nèi)湍流與云微物理特性的關(guān)系。

    參考文獻(xiàn)(References)

    Battan L J,1964.Some observations of vertical velocities and precipitation sizes in a thunderstorm[J].J Appl Meteor,3(4):415-420.

    Bouniol D,Illingworth A,Hogan R,2004.Deriving turbulent kinetic energy dissipation rate within clouds using ground based radar[C]//European Conference on Radar Meteorology.Third European Conference on Radar Meteorology (ERAD).

    段藝萍,劉壽東,劉黎平,等,2013.利用云雷達(dá)反演層狀云空氣垂直速度及微物理參數(shù)的個(gè)例研究[J].科學(xué)技術(shù)與工程,13(27):7933-7940. Duan Y P,Liu S D,Liu L P,et al.,2013.A case test of retrieving stratus vertical velocity & microphysical parameters from millimeter-wave cloud radar[J].Science Technology and Engineering,13(27):7933-7940.(in Chinese).

    樊雅文,黃興友,李鋒,2013.毫米波雷達(dá)測(cè)云個(gè)例研究[J].大氣科學(xué)學(xué)報(bào),36(5):554-559. Fan Y W,Huang X Y,Li F,2013.A case study on cloud measurement with a 35GHz millimeter-wave cloud radar[J].Trans Atmos Sci,36(5):554-559.(in Chinese).

    Feingold G,F(xiàn)risch A S,Stevens B,et al.,1999.On the relationship among cloud turbulence,droplet formation and drizzle as viewed by Doppler radar,microwave radiometer and lidar[J].J Geophys Res,104(D18):22195-22203.

    Frisch A S,F(xiàn)airall C W,Snider J B,1995.Measurement of stratus cloud and drizzle parameters in ASTEX with a kα-band Doppler radar and a microwave radiometer[J].J Atmos Sci,52(16):2788-2799.

    Gossard E E,Strauch R O,Rogers R R,1990.Evolution of dropsize distributions in liquid precipitation observed by ground-based Doppler radar[J].J Atmos Oceanic Technol,7(6):815-828.

    Gossard E E,1994.Measurement of cloud droplet size spectra by Doppler radar[J].J Atmos Oceanic Technol,11(3):712-726.

    Hobbs P V,F(xiàn)unk N T,1984.Cloud and precipitation studies with a millimetre-wave radar:a pictorial overview[J].Weather,39(11):334-339.

    黃興友,蘆荀,黃勇,等,2019.層狀云微物理參數(shù)反演及其輻射效應(yīng)的個(gè)例研究[J].大氣科學(xué)學(xué)報(bào),42(5):769-777. Huang X Y,Lu X,Huang Y,et al.,2019.A case study on the microphysical parameter retrieval and radiative effects of stratus clouds[J].Trans Atmos Sci,42(5):769-777.(in Chinese).

    Kollias P,Albrecht B A,Lhermitte R,et al.,2001.Radar observations of updrafts,downdrafts,and turbulence in fair-weather cumuli[J].J Atmos Sci,58(13):1750-1766.

    Kollias P,Clothiaux E E,Miller M A,et al.,2007.Millimeter-wavelength radars:new frontier in atmospheric cloud and precipitation research[J].Bull Amer Meteor Soc,88(10):1608-1624.

    OConnor E J,Hogan R J,Illingworth A J,2005.Retrieving stratocumulus drizzle parameters using Doppler radar and lidar[J].J Appl Meteorol,44(1):14-27.

    宗蓉,劉黎平,銀燕,2014.基于飛機(jī)觀測(cè)資料的降水粒子反射率因子閾值分析[J].大氣科學(xué)學(xué)報(bào),37(4):469-475. Zong R,Liu L P,Yin Y,2014.Analysis of drizzle reflectivity threshold using in-situ measurements[J].Trans Atmos Sci,37(4):469-475.(in Chinese).

    宗蓉,郭偉,仲凌志,2016.毫米波雷達(dá)和激光雷達(dá)探測(cè)云邊界的對(duì)比分析[J].激光技術(shù),40(2):171-177. Zong R,Guo W,Zhong L Z,2016.Comparison analysis of cloud boundary detection by means of millimeter wavelength cloud radar and lidar[J].Laser Technol,40(2):171-177.(in Chinese).

    Clouds exert an important influence on the weather and precipitation.The detection of cloud microphysical parameters(particle effective radius,drop number concentration,etc.),as well as of macro-physical parameters(cloud height,cloud thickness),is essential to studying cloud physical process and cloud radiation effects.The present paper aims to perform the detection and retrieval of cloud parameters and in-cloud turbulence characteristics with a ground-based millimeter cloud radar.By applying the Doppler spectrum data,the ambient vertical air velocity and cloud microphysical parameters are solved,whereby the turbulent dissipation rate ε and its distribution are obtained.Next,the relationships of vertical air motion to cloud microphysical parameters,radar reflectivity factor,Doppler velocity and Doppler spectrum width are thoroughly studied and analyzed so as to aid in understanding the cloud evolution.Cloud radar measurements and retrievals performed on a Stratiform cloud process case at Daocheng,Sichuan on August 8,2016 reveal the following:1)The cloud particle effective radius increases as the updraft increase,and the cloud drop number concentration decreases accordingly due to the collision-coalescence process;2)The turbulent dissipation rates ε are strong the near cloud base and cloud top,and weak in the middle layer,ranging from 10-8 to 10-2 m2·s-3.The Doppler spectrum validates the hypotheses that the turbulent scale captured by cloud radar is in the inertial subrange;3)As the air updraft grows stronger,the cloud particles deposition speed is weakened,while the Doppler spectrum width is strengthened.

    millimeter-wave cloud radar;turbulent dissipation rate;vertical air motions;cloud microphysical parameters;retrieval

    doi:10.13878/j.cnki.dqkxxb.20180203003

    (責(zé)任編輯:劉菲)

    猜你喜歡
    反演
    反演對(duì)稱變換在解決平面幾何問(wèn)題中的應(yīng)用
    基于Sentinel-2遙感影像的黑土區(qū)土壤有效磷反演
    基于ADS-B的風(fēng)場(chǎng)反演與異常值影響研究
    利用錐模型反演CME三維參數(shù)
    兩種測(cè)高重力異常反演海底地形方法比較
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種穩(wěn)健的彈性阻抗反演方法
    拉普拉斯變換反演方法探討
    封閉母線接頭溫度場(chǎng)反演
    欧美变态另类bdsm刘玥| 久热这里只有精品99| 免费看不卡的av| 国产成人freesex在线| 亚洲精品久久久久久婷婷小说| 一本—道久久a久久精品蜜桃钙片| 国产成人a∨麻豆精品| 999精品在线视频| h视频一区二区三区| 日韩熟女老妇一区二区性免费视频| 男女啪啪激烈高潮av片| av一本久久久久| a级毛片黄视频| 国产一区有黄有色的免费视频| 国产一区二区三区综合在线观看 | 狂野欧美激情性xxxx在线观看| 在线 av 中文字幕| 成年女人在线观看亚洲视频| 国产成人精品无人区| 国产成人免费观看mmmm| 午夜福利影视在线免费观看| 亚洲av成人精品一区久久| 欧美日韩综合久久久久久| 亚洲天堂av无毛| 丰满迷人的少妇在线观看| 久久精品夜色国产| 久久婷婷青草| 亚洲精品乱码久久久久久按摩| 久久久国产一区二区| 精品久久久精品久久久| 乱人伦中国视频| 国产在线视频一区二区| 亚洲欧洲国产日韩| 午夜福利视频精品| videos熟女内射| 波野结衣二区三区在线| av免费在线看不卡| 亚洲精品国产av成人精品| 亚洲精品国产av蜜桃| 久久亚洲国产成人精品v| 亚洲激情五月婷婷啪啪| 黑人欧美特级aaaaaa片| 韩国av在线不卡| 午夜免费男女啪啪视频观看| 这个男人来自地球电影免费观看 | 岛国毛片在线播放| 亚洲国产av影院在线观看| 91精品国产国语对白视频| 亚洲成人一二三区av| 亚洲性久久影院| 免费黄色在线免费观看| 老女人水多毛片| 国产成人精品在线电影| 色婷婷av一区二区三区视频| 一本一本综合久久| 国产av一区二区精品久久| 精品人妻偷拍中文字幕| 99久久中文字幕三级久久日本| 能在线免费看毛片的网站| av线在线观看网站| 熟女电影av网| 久久久国产一区二区| 国产精品久久久久久精品电影小说| 国产色婷婷99| 18在线观看网站| 搡老乐熟女国产| 少妇熟女欧美另类| 久久久久国产网址| 七月丁香在线播放| 插阴视频在线观看视频| 免费黄频网站在线观看国产| 欧美精品人与动牲交sv欧美| 黑人欧美特级aaaaaa片| 国产深夜福利视频在线观看| 久久国产精品大桥未久av| 亚洲国产精品专区欧美| 欧美 亚洲 国产 日韩一| 免费黄频网站在线观看国产| 亚洲国产精品一区二区三区在线| 国产片内射在线| 99热全是精品| 国产精品一区二区三区四区免费观看| 男人操女人黄网站| 亚洲国产最新在线播放| 久久久国产欧美日韩av| 国产女主播在线喷水免费视频网站| 又粗又硬又长又爽又黄的视频| 人妻少妇偷人精品九色| 美女主播在线视频| 插阴视频在线观看视频| 一区二区日韩欧美中文字幕 | 欧美日韩视频精品一区| 欧美日韩国产mv在线观看视频| 国产日韩欧美亚洲二区| 男的添女的下面高潮视频| 精品久久国产蜜桃| 嘟嘟电影网在线观看| 国产白丝娇喘喷水9色精品| 精品久久久久久电影网| 国产精品一区www在线观看| 国产成人午夜福利电影在线观看| 精品久久久久久久久亚洲| 亚洲,一卡二卡三卡| 18禁动态无遮挡网站| 免费大片黄手机在线观看| 视频区图区小说| 国模一区二区三区四区视频| 午夜视频国产福利| 看非洲黑人一级黄片| 嫩草影院入口| 人妻 亚洲 视频| 欧美精品人与动牲交sv欧美| 天美传媒精品一区二区| 午夜福利视频在线观看免费| 18在线观看网站| 精品少妇久久久久久888优播| 久久久久久久大尺度免费视频| 免费看光身美女| 亚洲av成人精品一二三区| 一级毛片我不卡| 黄色欧美视频在线观看| 哪个播放器可以免费观看大片| 国产成人免费无遮挡视频| 国产亚洲精品第一综合不卡 | 免费久久久久久久精品成人欧美视频 | 久久精品人人爽人人爽视色| 少妇的逼好多水| 亚洲精品成人av观看孕妇| 99re6热这里在线精品视频| 3wmmmm亚洲av在线观看| 免费日韩欧美在线观看| 伦精品一区二区三区| av福利片在线| 女性被躁到高潮视频| 亚洲精华国产精华液的使用体验| 精品久久久久久电影网| 在线免费观看不下载黄p国产| 尾随美女入室| 97超视频在线观看视频| 韩国av在线不卡| 国产一区二区三区综合在线观看 | 午夜老司机福利剧场| 国产日韩欧美在线精品| 国产精品国产三级专区第一集| 亚洲,一卡二卡三卡| 国产日韩欧美亚洲二区| 国产精品偷伦视频观看了| 国精品久久久久久国模美| 午夜福利视频在线观看免费| 亚洲美女黄色视频免费看| 热re99久久国产66热| av在线老鸭窝| 观看av在线不卡| 插阴视频在线观看视频| 国产毛片在线视频| 国产精品人妻久久久影院| 久久精品夜色国产| 久久久久久久久久人人人人人人| 国产精品欧美亚洲77777| 嘟嘟电影网在线观看| 精品国产露脸久久av麻豆| 少妇的逼好多水| 成人二区视频| 国产女主播在线喷水免费视频网站| 亚洲av电影在线观看一区二区三区| 熟妇人妻不卡中文字幕| 视频在线观看一区二区三区| 最后的刺客免费高清国语| 日韩精品免费视频一区二区三区 | 99久久中文字幕三级久久日本| 国产成人aa在线观看| 日韩成人伦理影院| 视频区图区小说| 精品一区二区三区视频在线| 久久久国产精品麻豆| 国产成人精品婷婷| 美女视频免费永久观看网站| 22中文网久久字幕| 王馨瑶露胸无遮挡在线观看| 中国美白少妇内射xxxbb| 成人国产麻豆网| a级毛片黄视频| 免费看光身美女| 久久久久视频综合| 女的被弄到高潮叫床怎么办| 永久网站在线| 校园人妻丝袜中文字幕| 日韩精品免费视频一区二区三区 | 日日撸夜夜添| 蜜桃国产av成人99| a级毛片在线看网站| av在线播放精品| 久久久久久久亚洲中文字幕| 又大又黄又爽视频免费| 久久精品国产亚洲网站| 免费观看av网站的网址| 日日爽夜夜爽网站| 亚洲精华国产精华液的使用体验| 亚洲av福利一区| 妹子高潮喷水视频| 寂寞人妻少妇视频99o| 久久人人爽av亚洲精品天堂| 久久狼人影院| 色婷婷av一区二区三区视频| 久热久热在线精品观看| 亚洲精品日韩在线中文字幕| 欧美 亚洲 国产 日韩一| 高清视频免费观看一区二区| 一边亲一边摸免费视频| 精品亚洲乱码少妇综合久久| 国产 一区精品| 夜夜骑夜夜射夜夜干| 人人妻人人爽人人添夜夜欢视频| 亚洲内射少妇av| 三上悠亚av全集在线观看| 三级国产精品欧美在线观看| 王馨瑶露胸无遮挡在线观看| 国产一区二区三区综合在线观看 | 水蜜桃什么品种好| 人人妻人人澡人人爽人人夜夜| 精品久久久久久久久av| av在线app专区| 亚洲精品日韩在线中文字幕| 婷婷色综合www| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 少妇高潮的动态图| √禁漫天堂资源中文www| 内地一区二区视频在线| 久久免费观看电影| 热99国产精品久久久久久7| 免费看光身美女| 一边亲一边摸免费视频| 亚洲美女搞黄在线观看| 在线观看免费日韩欧美大片 | 最新的欧美精品一区二区| 人妻人人澡人人爽人人| 亚洲中文av在线| 欧美少妇被猛烈插入视频| 日日啪夜夜爽| 精品国产乱码久久久久久小说| a 毛片基地| 国产白丝娇喘喷水9色精品| a级毛片免费高清观看在线播放| 97在线人人人人妻| 国产欧美日韩综合在线一区二区| 性色av一级| 男女免费视频国产| 亚洲国产最新在线播放| 日韩熟女老妇一区二区性免费视频| 伦理电影免费视频| 亚洲精品久久成人aⅴ小说 | 成年av动漫网址| 久久久久久久久久久久大奶| 亚洲精品国产av成人精品| 国产成人freesex在线| 多毛熟女@视频| videos熟女内射| 日本色播在线视频| a级毛色黄片| 高清黄色对白视频在线免费看| 91精品国产国语对白视频| 日韩在线高清观看一区二区三区| 免费久久久久久久精品成人欧美视频 | 人妻系列 视频| 国产成人精品久久久久久| 免费观看无遮挡的男女| 精品人妻熟女毛片av久久网站| 久久久国产一区二区| 久久精品国产鲁丝片午夜精品| 街头女战士在线观看网站| 久久国内精品自在自线图片| 亚洲经典国产精华液单| 18禁在线播放成人免费| 97超碰精品成人国产| xxx大片免费视频| 久久久久久人妻| 国产极品天堂在线| 欧美日韩av久久| .国产精品久久| 黄色怎么调成土黄色| 精品人妻在线不人妻| 亚洲国产av新网站| 99热这里只有精品一区| 久久精品国产自在天天线| 日韩一区二区三区影片| 日韩av不卡免费在线播放| av电影中文网址| av有码第一页| 国产av国产精品国产| 少妇的逼水好多| 一本色道久久久久久精品综合| 在线亚洲精品国产二区图片欧美 | 日日啪夜夜爽| 精品一区二区免费观看| 国产在线一区二区三区精| 69精品国产乱码久久久| 亚洲成人一二三区av| 亚洲一区二区三区欧美精品| 午夜激情av网站| 成年美女黄网站色视频大全免费 | 纵有疾风起免费观看全集完整版| 中文字幕精品免费在线观看视频 | 久久久久精品久久久久真实原创| 欧美日韩亚洲高清精品| 如日韩欧美国产精品一区二区三区 | 久久婷婷青草| 亚洲内射少妇av| 又黄又爽又刺激的免费视频.| 国产精品秋霞免费鲁丝片| 国产精品一二三区在线看| 欧美精品人与动牲交sv欧美| 欧美日韩av久久| 制服丝袜香蕉在线| 99久国产av精品国产电影| 黑丝袜美女国产一区| xxxhd国产人妻xxx| 中国美白少妇内射xxxbb| 黑人猛操日本美女一级片| 在线观看一区二区三区激情| 国产欧美另类精品又又久久亚洲欧美| 老女人水多毛片| 国产伦精品一区二区三区视频9| 性色av一级| 国产精品国产三级专区第一集| 国产 精品1| 国产成人精品久久久久久| 国语对白做爰xxxⅹ性视频网站| 黄色一级大片看看| 一二三四中文在线观看免费高清| 欧美人与善性xxx| 青春草亚洲视频在线观看| 久久精品熟女亚洲av麻豆精品| 日韩大片免费观看网站| 欧美精品高潮呻吟av久久| 亚洲av男天堂| 在线观看人妻少妇| 各种免费的搞黄视频| 十八禁高潮呻吟视频| 日韩不卡一区二区三区视频在线| 亚洲人成网站在线播| 日本午夜av视频| 婷婷色av中文字幕| 69精品国产乱码久久久| 18禁观看日本| 亚洲欧美成人综合另类久久久| 中文字幕制服av| 男人操女人黄网站| 国产欧美日韩一区二区三区在线 | 国产精品麻豆人妻色哟哟久久| 亚洲伊人久久精品综合| 高清在线视频一区二区三区| 桃花免费在线播放| 少妇精品久久久久久久| av又黄又爽大尺度在线免费看| 国产亚洲av片在线观看秒播厂| 黑人高潮一二区| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美成人综合另类久久久| 中文字幕免费在线视频6| 国产一区二区在线观看日韩| 中文字幕亚洲精品专区| 人人妻人人澡人人看| 飞空精品影院首页| 人人妻人人澡人人看| 蜜桃久久精品国产亚洲av| 欧美bdsm另类| 老熟女久久久| 日日撸夜夜添| 夜夜爽夜夜爽视频| 欧美少妇被猛烈插入视频| 成人午夜精彩视频在线观看| 国产亚洲午夜精品一区二区久久| 亚洲人成网站在线观看播放| 99九九线精品视频在线观看视频| 99国产精品免费福利视频| 久久综合国产亚洲精品| 亚洲欧美成人精品一区二区| 天美传媒精品一区二区| 看十八女毛片水多多多| 毛片一级片免费看久久久久| 国产男女超爽视频在线观看| 亚洲婷婷狠狠爱综合网| 狠狠婷婷综合久久久久久88av| 亚洲精品国产av成人精品| 肉色欧美久久久久久久蜜桃| a级毛片黄视频| 三级国产精品欧美在线观看| 久久久久国产精品人妻一区二区| 国产极品天堂在线| 丝袜在线中文字幕| 人人妻人人澡人人看| 国产在线一区二区三区精| 日日撸夜夜添| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久av| 日韩成人伦理影院| 久久97久久精品| 欧美激情极品国产一区二区三区 | 老熟女久久久| 日韩三级伦理在线观看| 一区二区日韩欧美中文字幕 | 亚洲高清免费不卡视频| 一级毛片黄色毛片免费观看视频| 久久久久网色| 这个男人来自地球电影免费观看 | 99热网站在线观看| 国产日韩欧美在线精品| a级毛片免费高清观看在线播放| 亚洲精华国产精华液的使用体验| 最近中文字幕高清免费大全6| 国产探花极品一区二区| 黄色怎么调成土黄色| 少妇高潮的动态图| 夜夜骑夜夜射夜夜干| 久久99一区二区三区| 日韩 亚洲 欧美在线| 免费高清在线观看视频在线观看| 欧美日韩亚洲高清精品| 久久久久久久久大av| 国产一区有黄有色的免费视频| 大片电影免费在线观看免费| 国产成人午夜福利电影在线观看| xxxhd国产人妻xxx| 国产av精品麻豆| 99久国产av精品国产电影| 国产永久视频网站| 丰满饥渴人妻一区二区三| 五月伊人婷婷丁香| videos熟女内射| 自拍欧美九色日韩亚洲蝌蚪91| 国产色婷婷99| 青春草亚洲视频在线观看| 日本黄色日本黄色录像| 午夜老司机福利剧场| 亚洲国产精品一区二区三区在线| 秋霞在线观看毛片| 日韩中字成人| 国产日韩欧美亚洲二区| 夫妻午夜视频| 久久精品国产亚洲av涩爱| 男男h啪啪无遮挡| 精品国产露脸久久av麻豆| 中文乱码字字幕精品一区二区三区| 视频在线观看一区二区三区| 欧美日韩亚洲高清精品| 丰满少妇做爰视频| 看十八女毛片水多多多| 国产精品嫩草影院av在线观看| 交换朋友夫妻互换小说| 免费看av在线观看网站| 亚洲怡红院男人天堂| 亚洲欧美色中文字幕在线| 多毛熟女@视频| 欧美+日韩+精品| 男男h啪啪无遮挡| 日本-黄色视频高清免费观看| a 毛片基地| 亚洲美女黄色视频免费看| 自线自在国产av| 国产一区亚洲一区在线观看| 国产精品人妻久久久久久| 91在线精品国自产拍蜜月| 99精国产麻豆久久婷婷| 精品人妻熟女av久视频| 成人黄色视频免费在线看| 十八禁高潮呻吟视频| 日韩av免费高清视频| 97超碰精品成人国产| 啦啦啦视频在线资源免费观看| 亚洲成人av在线免费| 国内精品宾馆在线| 国产精品偷伦视频观看了| 高清毛片免费看| av又黄又爽大尺度在线免费看| 18+在线观看网站| 免费观看无遮挡的男女| 亚洲国产精品一区二区三区在线| 建设人人有责人人尽责人人享有的| 日韩伦理黄色片| 亚洲国产毛片av蜜桃av| 狂野欧美激情性xxxx在线观看| 九草在线视频观看| 亚洲四区av| 欧美亚洲 丝袜 人妻 在线| 成年av动漫网址| 中文精品一卡2卡3卡4更新| 久久亚洲国产成人精品v| av卡一久久| 在线看a的网站| 日本黄色片子视频| 七月丁香在线播放| 在线播放无遮挡| 国产成人a∨麻豆精品| 三级国产精品欧美在线观看| 免费看av在线观看网站| 一本—道久久a久久精品蜜桃钙片| 午夜免费男女啪啪视频观看| 毛片一级片免费看久久久久| 欧美精品人与动牲交sv欧美| 久久久久久久大尺度免费视频| 黑人高潮一二区| 尾随美女入室| 亚洲综合精品二区| 国产极品天堂在线| 亚洲少妇的诱惑av| 欧美少妇被猛烈插入视频| 少妇高潮的动态图| 丝袜脚勾引网站| 如日韩欧美国产精品一区二区三区 | 国产精品一区二区在线观看99| 日韩亚洲欧美综合| 欧美日韩视频精品一区| 一级爰片在线观看| 午夜免费鲁丝| 亚洲欧美精品自产自拍| 免费看av在线观看网站| 日韩三级伦理在线观看| 欧美+日韩+精品| 亚洲欧美清纯卡通| 超色免费av| 高清在线视频一区二区三区| 狂野欧美激情性xxxx在线观看| 国产白丝娇喘喷水9色精品| 欧美日韩av久久| av福利片在线| h视频一区二区三区| 亚洲精品av麻豆狂野| 两个人的视频大全免费| 国产精品久久久久久久久免| 各种免费的搞黄视频| 人人妻人人爽人人添夜夜欢视频| 中文字幕精品免费在线观看视频 | 亚洲精品久久午夜乱码| 永久网站在线| 久热久热在线精品观看| 日韩av免费高清视频| www.av在线官网国产| 老司机影院成人| 亚洲天堂av无毛| 黄色欧美视频在线观看| 亚洲精品av麻豆狂野| 亚洲av欧美aⅴ国产| 国产爽快片一区二区三区| 亚洲内射少妇av| 国产精品.久久久| 水蜜桃什么品种好| 亚洲国产av影院在线观看| 免费看不卡的av| 最后的刺客免费高清国语| av在线观看视频网站免费| 精品人妻一区二区三区麻豆| 搡老乐熟女国产| 国产有黄有色有爽视频| 少妇 在线观看| 精品国产一区二区久久| h视频一区二区三区| 毛片一级片免费看久久久久| 亚洲av国产av综合av卡| 国产亚洲av片在线观看秒播厂| 免费大片18禁| 人妻夜夜爽99麻豆av| 中文字幕免费在线视频6| 久久影院123| 欧美日韩一区二区视频在线观看视频在线| 99热这里只有是精品在线观看| 97超视频在线观看视频| 日韩,欧美,国产一区二区三区| 色婷婷久久久亚洲欧美| 最新中文字幕久久久久| 国产成人午夜福利电影在线观看| av电影中文网址| 午夜av观看不卡| 国产免费又黄又爽又色| 中文字幕av电影在线播放| 美女大奶头黄色视频| 亚洲av在线观看美女高潮| av天堂久久9| 免费观看在线日韩| 中文字幕亚洲精品专区| 国产片内射在线| 你懂的网址亚洲精品在线观看| 蜜桃国产av成人99| 精品久久国产蜜桃| 最近中文字幕2019免费版| 另类精品久久| 国产精品免费大片| 黑丝袜美女国产一区| 成人二区视频| 男女边吃奶边做爰视频| 色网站视频免费| 亚洲成人一二三区av| 亚洲内射少妇av| 久久国产精品大桥未久av| 免费看av在线观看网站| 成人国语在线视频| 欧美激情极品国产一区二区三区 | 久久av网站| 老司机影院成人| 九九久久精品国产亚洲av麻豆| 亚洲国产精品999| 男人爽女人下面视频在线观看| 国产一区二区在线观看av| 麻豆成人av视频| 国产 精品1| 色吧在线观看| 国产成人精品福利久久| 国产精品国产三级国产专区5o| 大码成人一级视频| 少妇高潮的动态图| 国产精品一二三区在线看| 一级片'在线观看视频| 国产高清不卡午夜福利| 免费看av在线观看网站| 国产精品欧美亚洲77777| 狠狠婷婷综合久久久久久88av|