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

    基于共形網(wǎng)格技術(shù)的ATEM三維數(shù)值模擬

    2020-03-02 10:05:20嵇艷鞠任桂瑩關(guān)珊珊黎東升
    關(guān)鍵詞:共形電導(dǎo)率電場

    嵇艷鞠,任桂瑩,關(guān)珊珊?,黎東升

    (1.吉林大學(xué) 地球信息探測儀器教育部重點實驗室,吉林長春130026;2.吉林大學(xué) 儀器科學(xué)與電氣工程學(xué)院,吉林 長春130026;3.吉林大學(xué)國家地球物理探測儀器工程技術(shù)研究中心,吉林長春130026)

    航空瞬變電磁方法(Airborne Transient Electromagnetic Method,ATEM)采用飛行器搭載電磁系統(tǒng)對大地的二次場進行快速測量,通過分析二次場獲得地下電阻率分布,從而了解地質(zhì)體的構(gòu)造.相比于地面瞬變電磁法,航空瞬變電磁具有探測速度快,范圍廣,可在短時間內(nèi)獲得較多信息等特點,更適用于森林、沼澤、割裂地形、山區(qū)等地形復(fù)雜區(qū)域的勘探.航空瞬變電磁法在國內(nèi)外發(fā)展較為迅速,在淺覆蓋區(qū)尋找金屬礦、隱蔽洞體等典型目標(biāo)體探測方面起到了重要的作用[1-4].常用的瞬變電磁三維正演模擬方法包括有限差分法、積分方程法、有限元方法等,其中有限差分法具有數(shù)學(xué)表達(dá)簡單、直觀高效等特點,因此發(fā)展得較為成熟.Wang等人[5]采用有限差分法進行瞬變電磁三維數(shù)值模擬.Commer等[6]2004年利用并行有限差分方法實現(xiàn)了電性源瞬變電磁三維數(shù)值模擬,2015年對時間步長進行了改進以減少計算時間[7].Li等人[8]基于有限差分法研究了含水結(jié)構(gòu)的瞬變電磁三維數(shù)值模擬.隨著瞬變電磁探測技術(shù)的發(fā)展,地下復(fù)雜三維異常體結(jié)構(gòu)的精細(xì)化勘探得到關(guān)注,采用傳統(tǒng)的時域有限差分(Finite Difference Time Domain,F(xiàn)DTD)方法進行地下不規(guī)則異常體的數(shù)值模擬時,常常會由于邊界處網(wǎng)格的階梯狀網(wǎng)格近似導(dǎo)致較大誤差,故如何處理復(fù)雜異常體的邊界問題以及提高三維不規(guī)則體的計算精度,成為地球物理勘探中的熱點問題.Cao等人[9]提出了自適應(yīng)有限差分方法,采用局部細(xì)化的方式提高了三維不規(guī)則異常體的模擬精度,但增加了計算的復(fù)雜性.Chang等人[10]采用非均勻網(wǎng)格方式對煤礦充水區(qū)域的全波形瞬變電磁響應(yīng)進行了數(shù)值模擬,但網(wǎng)格剖分較難.徐巖等人[11]在瞬變電磁三維正演模擬時,將起伏地形或不規(guī)則體表面網(wǎng)格進行了精細(xì)剖分,雖然提高了計算精度,但會增加計算量.

    共形網(wǎng)格技術(shù)是近年來新興的一種數(shù)值模擬技術(shù),它基于有限差分方法,對目標(biāo)體邊界進行特殊處理,以減小傳統(tǒng)FDTD方法所產(chǎn)生的階梯狀網(wǎng)格近似誤差,提高了有限差分方法的計算精度,在不增加網(wǎng)格數(shù)量的基礎(chǔ)上,提高了計算效率.Ilinca等人[12]將共形網(wǎng)格技術(shù)應(yīng)用于材料物理學(xué)中,解決了材料復(fù)雜界面的熱傳遞問題.Wang等人[13]提出了一種改進的時域有限差分共形方案,可快速預(yù)測雷達(dá)截面以及復(fù)雜三維結(jié)構(gòu)的感應(yīng)電流分布.Guo等人[14]利用改進的共形FDTD網(wǎng)格劃分算法模擬復(fù)雜三維模型.武超[15]將共形網(wǎng)格技術(shù)應(yīng)用于無線電物理方面,研究了復(fù)雜目標(biāo)的電磁特性計算.范宜仁等人[16]將共形網(wǎng)格技術(shù)應(yīng)用于隨鉆電磁波測井響應(yīng)數(shù)值模擬中,為隨鉆快速反演提供支持.Nicolas等人[17]提出了一種基于六面體網(wǎng)格的共形方案,并在材料邊緣處進行特定分裂,實現(xiàn)了細(xì)化和共形.Wei等人[18]擴展了共形時域有限差分方法,將該方法應(yīng)用于地震中的彈性波計算,并取得了很好的結(jié)果.Cabello等人[19]將共形網(wǎng)格技術(shù)應(yīng)用于雷達(dá)散射截面的計算中,取得了較好的效果.但在三維地質(zhì)體探測中,如何采用共形網(wǎng)格技術(shù)更好地擬合曲面異常體的邊界問題還有待研究.

    本文在傳統(tǒng)FDTD方法的基礎(chǔ)上,結(jié)合共形網(wǎng)格技術(shù)開展了航空瞬變電磁的數(shù)值模擬研究.基于介質(zhì)參數(shù)等效的思想,更新了傳統(tǒng)FDTD方法的電場表達(dá)式;對于具有解析式的典型目標(biāo)體,討論了基于解析式建模的FDTD共形網(wǎng)格生成過程,采用共形網(wǎng)格技術(shù)實現(xiàn)了地下曲面異常體的三維數(shù)值模擬,提高了傳統(tǒng)FDTD方法對曲面異常體進行數(shù)值模擬時的精度.

    1 共形網(wǎng)格技術(shù)

    利用常規(guī)FDTD方法進行數(shù)值模擬時,Yee元胞僅可以賦一個電導(dǎo)率值,在模擬曲面異常體邊界時勢必會帶來由于階梯狀網(wǎng)格近似而產(chǎn)生的誤差.圖1為球形異常體二維截面階梯狀網(wǎng)格近似圖,介質(zhì)1表示目標(biāo)體內(nèi)部,介質(zhì)2表示目標(biāo)體外部,其中陰影區(qū)域表示位于目標(biāo)體表面附近的網(wǎng)格,為了減小誤差,在附近網(wǎng)格處,采用介質(zhì)參數(shù)等效技術(shù),更新傳統(tǒng)FDTD迭代方程,這種技術(shù)叫做共形網(wǎng)格技術(shù),而目標(biāo)體附近的網(wǎng)格被稱作共形網(wǎng)格.

    圖1 球形異常體二維截面階梯狀網(wǎng)格近似圖Fig.1 Spherical two-dimensional section step approximation

    在計算時,從曲面異常體的邊界出發(fā),確定共形網(wǎng)格,基于線性加權(quán)平均思想,重新計算邊界共形網(wǎng)格處電場的等效電導(dǎo)率值,將等效電導(dǎo)率代入相應(yīng)位置的電場表達(dá)式,進行更精確的數(shù)值模擬[20].其中,對于具有解析式的異常體模型,可直接根據(jù)其解析式以及整個計算區(qū)域網(wǎng)格的坐標(biāo),得到具體共形網(wǎng)格的坐標(biāo)位置以及需要共形處理的電場所在棱邊位置,為異常體表面的共形計算提供參數(shù).

    2 基于共形網(wǎng)格技術(shù)的電場迭代方程推導(dǎo)

    圖2為航空瞬變電磁法的工作原理示意圖.當(dāng)發(fā)射線圈中通入發(fā)射電流時,在整個空間產(chǎn)生一次場,發(fā)射電流的瞬間變化在大地中形成渦流并產(chǎn)生二次場,接收線圈接收到二次場后以感應(yīng)電壓的形式記錄下來,文中的數(shù)值模擬采用發(fā)射線圈與接收線圈共中心的方式.

    圖2 航空瞬變電磁工作原理圖Fig.2 Schematic diagram of airborne transient electromagnetic working principle

    在用有限差分方法對三維目標(biāo)體進行數(shù)值模擬時,基于Maxwell旋度方程進行差分離散,式(1)和(2)為時域Maxwell旋度方程.常規(guī)FDTD方法在數(shù)值模擬時,需要將整個計算區(qū)域離散成規(guī)則的六面體,即近似為圖3所示的Yee元胞.電場和磁場在時間上相差半個步長交替取樣,Maxwell旋度方程離散后推導(dǎo)得到顯式差分方程.擬合復(fù)雜模型的曲面邊界時,在處理介質(zhì)共形中,只涉及到電導(dǎo)率的變化,因此只需更新Maxwell方程組的電場迭代方程[21].

    式中:e(r,t)為電場強度;b(r,t)為磁感應(yīng)強度;h(r,t)為磁場強度;j(r,t)為傳導(dǎo)電流密度;ε(r)為介電常數(shù).

    圖3 FDTD中的Yee元胞Fig.3 The Yee cell of FDTD

    基于共形網(wǎng)格技術(shù),將傳統(tǒng)電場迭代方程中的電導(dǎo)率進行等效替代.圖4為共形網(wǎng)格中電導(dǎo)率的等效.圖4(a)為三維FDTD共形網(wǎng)格示意圖,網(wǎng)格被目標(biāo)體分為介質(zhì)1和介質(zhì)2兩部分,由于介質(zhì)不同,電導(dǎo)率存在差異,因此,需求出共形網(wǎng)格中各電場位置的等效電導(dǎo)率.圖4(b)為三維共形網(wǎng)格在xoy面上的投影,介質(zhì)1的電導(dǎo)率為σ1,介質(zhì)2的電導(dǎo)率為σ2,介質(zhì)1與介質(zhì)2的分界面將一條棱邊分為兩部分,因電場分布在棱邊上,故相應(yīng)棱邊上的等效電導(dǎo)率值由分界面所截長度的線性加權(quán)平均得到.式(3)為A邊上x方向電場的等效電導(dǎo)率.

    以x方向為例,將式(3)代入傳統(tǒng)電場迭代方程中,得式(4)所示更新后的電場方程.y、z方向同理可得.

    圖4 共形網(wǎng)格中電導(dǎo)率的等效Fig.4 Equivalent of conductivity in a conformal grid

    3 共形參數(shù)確定方法

    通過分析三維共形網(wǎng)格可知,在利用共形網(wǎng)格技術(shù)進行數(shù)值模擬時,其關(guān)鍵是得到需要共形的網(wǎng)格位置以及網(wǎng)格共形參數(shù)lz1和lz2.本文以球形異常體表面的共形網(wǎng)格為例,分析說明具有解析式目標(biāo)體的共形參數(shù)確定方法.

    若交點z滿足式(5),則可判斷z所在的網(wǎng)格編號為k:

    設(shè)zmin所在的網(wǎng)格編號為kmin,zmax所在的網(wǎng)格編號為kmax,lz1k為z方向編號為k的網(wǎng)格線被目標(biāo)體所截的內(nèi)部長度,lz2k為z方向編號為k的網(wǎng)格線被目標(biāo)體所截的外部長度,由此可得z方向的共形網(wǎng)格坐標(biāo)及共形參數(shù),如式(6)和式(7)所示.

    圖5 球形異常體共形網(wǎng)格參數(shù)示意圖Fig.5 Schematic diagram of spherical conformal mesh parameters

    同理,可分別得到x方向、y方向的共形網(wǎng)格坐標(biāo)及共形參數(shù).

    4 模型計算及分析

    4.1 精度驗證

    為了驗證共形網(wǎng)格技術(shù)的有效性,建立5 980 m×5 980 m×3 910 m的包含吸收邊界條件的計算區(qū)域,大地被剖分為127×127×63個非均勻網(wǎng)格.采用發(fā)射與接收共中心的測量裝置,發(fā)射線圈位于整個計算區(qū)域的中心,距地面高度為30 m,發(fā)射波形為負(fù)階躍波,電流大小為300 A,其中三維異常體為三棱柱模型,位于發(fā)射線圈的正下方,截面為等腰三角形,腰長為,底邊長為120 m,三棱柱高度為300 m,三維異常體埋深為250 m,異常體電導(dǎo)率為1 S/m,圍巖電導(dǎo)率為0.01 S/m.圖6為三棱柱模型設(shè)置示意圖.

    圖6 三棱柱模型設(shè)置示意圖Fig.6 Schematic diagram of the setting of the triangular prism model

    在數(shù)值計算時,采用加載初始場的方式代替實際探測中的發(fā)射源,采用數(shù)字濾波算法解決計算式中含有的一階貝塞爾函數(shù)積分問題,之后利用線性濾波方法進行頻時轉(zhuǎn)換,得到時間域的初始場值[23-24].

    為了驗證采用共形網(wǎng)格技術(shù)進行數(shù)值模擬的有效性,將階梯狀網(wǎng)格近似結(jié)果及共形網(wǎng)格處理后結(jié)果與有限元方法進行了對比.三棱柱模型截面電導(dǎo)率階梯狀網(wǎng)格近似賦值示意圖如圖7(a)所示,其中虛線內(nèi)部近似為異常體電導(dǎo)率值.在利用共形網(wǎng)格技術(shù)處理時,需要將共形網(wǎng)格內(nèi)部電場處的電導(dǎo)率作等效,截面上需共形處理的電場位置示意圖如圖7(b)所示,以三棱柱截面左側(cè)為例,箭頭所在位置表示x、y方向需共形的電場位置坐標(biāo).

    圖7 三棱柱模型截面共形處理示意圖Fig.7 Schematic diagram of a conformal mesh of a triangular prism model cross section

    圖8(a)為階梯狀網(wǎng)格近似和共形網(wǎng)格技術(shù)處理后的電磁響應(yīng)與有限元方法的電磁響應(yīng)對比圖.由于有限元方法采用四面體模擬地下異常體,故該方法可精確剖分三棱柱模型,因此將有限元方法對三棱柱模型的數(shù)值模擬作為基準(zhǔn)解.由圖8(a)可知,3種方法得出的曲線一致性較好.從圖8(b)中可以看出,在1 ms之前,兩種處理方法的相對誤差相差不大,而在1 ms后,可看出利用共形網(wǎng)格技術(shù)處理后的相對誤差小于階梯狀網(wǎng)格近似方法處理的相對誤差,減小了傳統(tǒng)FDTD方法在目標(biāo)剖分過程中引入的階梯狀網(wǎng)格近似誤差,驗證了共形網(wǎng)格技術(shù)的有效性.由于在對三棱柱模型進行階梯狀網(wǎng)格近似時,采用的最小網(wǎng)格為10 m,因此兩者之間的誤差相差不大,當(dāng)增大最小剖分網(wǎng)格尺寸后,經(jīng)共形網(wǎng)格技術(shù)處理后的數(shù)值模擬優(yōu)勢會更加明顯.圖8(c)為空中垂直磁場的等值線分布.

    圖8 三棱柱模型計算結(jié)果圖Fig.8 Calculation results of the triangular prism model

    4.2 球形異常體模型的數(shù)值模擬

    采用與4.1節(jié)中同樣的系統(tǒng)模型參數(shù),建立4 060 m×4 060 m×2 950 m的包含吸收邊界層的計算區(qū)域,大地被剖分為111×111×77個網(wǎng)格.初始場及時間步長的選取方式與4.1節(jié)相同.在發(fā)射源的正下方設(shè)置球形異常體模型,半徑為100 m,埋深為100 m,剖分異常體的網(wǎng)格步長為10 m,計算區(qū)域示意圖如圖9(a)所示.以球形異常體的最大橫截面為例,作如圖9(b)所示的球形異常體截面階梯狀網(wǎng)格近似示意圖,并將與球形異常體邊界有交點的棱邊處的電場進行共形處理.其中,異常體電導(dǎo)率為1 S/m,圍巖電導(dǎo)率為0.01 S/m.球形異常體模型基準(zhǔn)解采用5 m×5 m×5 m的小網(wǎng)格進行剖分.計算區(qū)域不變,將球形異常體附近網(wǎng)格剖分成5 m×5 m×5 m的小網(wǎng)格,此時大地被剖分為163×163×160個網(wǎng)格,將此時的電磁響應(yīng)作為基準(zhǔn)解進行分析.圖10(a)為階梯狀網(wǎng)格近似和共形網(wǎng)格技術(shù)處理后的電磁響應(yīng)與基準(zhǔn)解的電磁響應(yīng)對比圖.圖10(b)為兩種處理方法的相對誤差曲線,從圖中可以看出在0.4 ms之前,兩種處理方法的誤差相差不大,而在0.4 ms后,共形處理的相對誤差明顯小于階梯狀網(wǎng)格近似處理的相對誤差,共形網(wǎng)格技術(shù)在處理復(fù)雜異常體邊界問題時具有明顯優(yōu)勢,相對誤差減小了1.65%.圖10(c)為空中垂直磁場的等值線分布圖.

    圖9 球形異常體模型設(shè)置示意圖Fig.9 Schematic diagram of the sphere model setting

    圖10 球形異常體模型計算結(jié)果圖Fig.10 Calculation results of sphere model

    圖11為球形異常體電磁響應(yīng)切片圖,異常體由于埋深較淺,對源的影響在早期較為明顯,隨著時間的推移,電磁響應(yīng)逐漸向下、向外擴散,且響應(yīng)逐漸減小.

    圖11 球形異常體電磁響應(yīng)切片圖Fig.11 Electromagnetic response slice diagram of spherical anomalous body

    4.3 橢球形異常體模型的數(shù)值模擬

    系統(tǒng)模型參數(shù)的選取方式均與4.1節(jié)相同.計算區(qū)域及網(wǎng)格剖分?jǐn)?shù)與4.2節(jié)相同.在發(fā)射源的正下方放置一個半長軸a=150 m,半短軸,埋深為100 m的橢球形異常體,采用基于解析式建模的共形網(wǎng)格技術(shù)求解橢球形異常體的時間域電磁響應(yīng).圖12(a)為計算區(qū)域示意圖,圖12(b)為橢球截面階梯狀網(wǎng)格近似示意圖.其中,異常體電導(dǎo)率為1 S/m,圍巖電導(dǎo)率為0.01 S/m.將橢球形異常體附近網(wǎng)格剖分成5 m×5 m×5 m的小網(wǎng)格,網(wǎng)格剖分?jǐn)?shù)同4.2節(jié)中的描述,以此時的電磁響應(yīng)作為基準(zhǔn)解進行分析.圖13(a)為階梯狀網(wǎng)格近似和共形網(wǎng)格技術(shù)處理后的電磁響應(yīng)與基準(zhǔn)解的電磁響應(yīng)對比圖.圖13(b)為兩種處理方法的相對誤差曲線圖,從圖中可以看出在0.5 ms之前,兩種處理方法的誤差相差不大,而在0.5 ms后共形處理的相對誤差明顯小于階梯狀網(wǎng)格近似處理的相對誤差,其中相對誤差減小了4.64%.與圖10(b)中曲線進行比較,在處理橢球形異常體時,共形網(wǎng)格處理的相對誤差相對于階梯狀網(wǎng)格近似處理的相對誤差減小較為明顯.通過分析圖9(b)、圖12(b)可知,需要共形處理的棱邊電場越多,共形處理的效果越明顯,精度就越高.

    圖12 橢球形異常體模型設(shè)置示意圖Fig.12 Schematic diagram of the ellipsoid model setting

    圖13 橢球形異常體模型計算結(jié)果圖Fig.13 Calculation results of ellipsoid model

    圖14為橢球形異常體電磁響應(yīng)切片圖.通過比較球形與橢球形異常體在t=0.5 ms,y=1 980 m的電磁響應(yīng)切片圖,在電磁場擴散的過程中,能夠反映出三維異常體的近似形狀,并且橢球形異常體的擴散速度要高于球形異常體的電磁場擴散速度.圖15為t=7 ms時,z分別為150、200、250 m時的橢球形異常體感應(yīng)電動勢等值線分布圖.

    本文分別計算了三棱柱、球形異常體、橢球形異常體作為三維異常體時的時間域電磁響應(yīng),在相同運行環(huán)境下,計算效率如表1所示.其中三棱柱模型利用共形網(wǎng)格技術(shù)以及階梯狀網(wǎng)格近似方法的運行時間分別為2 460 s和2 214 s;球形異常體模型利用共形網(wǎng)格技術(shù)以及階梯狀近似方法的運行時間分別為2 743 s和1 447 s;橢球形異常體模型利用共形網(wǎng)格技術(shù)以及階梯狀網(wǎng)格近似方法的運行時間分別為2 909 s和1 564 s.由于三棱柱模型的共形參數(shù)確定只涉及到二維平面,故相對容易獲得,因此與階梯狀網(wǎng)格近似方法的運行時間相比較相差不多,耗時為階梯狀網(wǎng)格近似處理的1.11倍.球形異常體模型和橢球形異常體模型為具有解析式的目標(biāo)體,需根據(jù)解析式確定共形參數(shù),因此在電場和磁場的每一次迭代計算中都需要判斷x、y、z三個方向共形網(wǎng)格的位置,即異常體表面與網(wǎng)格相交位置,以及交點所分相應(yīng)棱邊的內(nèi)外邊長.并將共形參數(shù)代入迭代式中進行計算,因此與三棱柱模型相比要更耗時,與階梯狀網(wǎng)格近似相比,其時間比值分別為1.90倍和1.86倍.由此可知,共形網(wǎng)格技術(shù)可在保證計算效率的同時提高計算精度.

    圖14 橢球形異常體電磁響應(yīng)切片圖Fig.14 Electromagnetic response slice diagram of ellipsoidal anomalous body

    圖15 t=7 ms時橢球形異常體感應(yīng)電動勢等值線分布圖Fig.15 Contour distribution of ellipsoidal anomalous electromagnetic response at t=7 ms

    表1 共形網(wǎng)格技術(shù)及階梯狀網(wǎng)格近似計算效率對比表Tab.1 The calculation efficiency comparison of conformal technology and stepped grid approximation

    5 結(jié)論

    本文基于傳統(tǒng)FDTD方法,結(jié)合了共形網(wǎng)格技術(shù),采用介質(zhì)參數(shù)等效思想,求得了線性加權(quán)平均電導(dǎo)率,更新了電場迭代方程;研究了有解析式的三維目標(biāo)體共形網(wǎng)格技術(shù);采用三棱柱模型驗證了共形網(wǎng)格技術(shù)的有效性,提高了地下曲面三維異常體的數(shù)值模擬精度.主要結(jié)論如下:

    1)在不改變網(wǎng)格大小及數(shù)量的基礎(chǔ)上,共形網(wǎng)格技術(shù)僅對電導(dǎo)率進行線性加權(quán)平均處理,并將計算求得的等效電導(dǎo)率代入傳統(tǒng)FDTD電場表達(dá)式,更新了電場迭代方程.共形網(wǎng)格技術(shù)在保證計算效率的條件下,與階梯狀網(wǎng)格近似方法相比提高了曲面三維異常體的計算精度.

    2)可通過三維坐標(biāo)投影法確定出共形網(wǎng)格位置,進而判斷出所截網(wǎng)格棱邊位置,以確定出曲面共形參數(shù).

    3)在共形網(wǎng)格處理中,網(wǎng)格棱邊與異常體邊界的交點越多,即需要共形處理的電場越多,共形網(wǎng)格技術(shù)的效果越明顯,精度提高越多.

    本文目前開展的工作主要是針對有解析式的三維目標(biāo)模型進行共形網(wǎng)格技術(shù)的研究,而對于無法寫出解析式的任意復(fù)雜模型,則應(yīng)分以下幾點進行共形網(wǎng)格技術(shù)處理:1)通過3ds Max等建模軟件建立三維目標(biāo)模型.2)根據(jù)目標(biāo)模型,生成.stl文件,導(dǎo)出三角面元數(shù)據(jù).3)由三角面元數(shù)據(jù),通過編程生成Yee網(wǎng)格坐標(biāo)以及共形參數(shù),并導(dǎo)入3ds Max中進行可視化.4)結(jié)合有限差分方法與共形網(wǎng)格技術(shù)進行計算,最終實現(xiàn)任意三維復(fù)雜模型的高精度數(shù)值模擬.

    猜你喜歡
    共形電導(dǎo)率電場
    具有共形能力的阻抗可調(diào)天線
    巧用對稱法 妙解電場題
    基于共形超表面的波束聚焦研究
    共形雙曲度量的孤立奇點
    基于比較測量法的冷卻循環(huán)水系統(tǒng)電導(dǎo)率檢測儀研究
    低溫脅迫葡萄新梢電導(dǎo)率和LT50值的研究
    電場強度單個表達(dá)的比較
    電場中六個常見物理量的大小比較
    高電導(dǎo)率改性聚苯胺的合成新工藝
    感生電場與動生電場的等效性探究
    物理與工程(2014年5期)2014-02-27 11:23:20
    久久人人爽av亚洲精品天堂 | 久久精品久久久久久久性| 成人毛片a级毛片在线播放| 美女高潮的动态| 99久久精品国产国产毛片| 色吧在线观看| 只有这里有精品99| 欧美高清成人免费视频www| 人妻少妇偷人精品九色| av天堂中文字幕网| 午夜免费观看性视频| 国产精品国产三级国产专区5o| 91精品伊人久久大香线蕉| 久久久久久久午夜电影| 欧美日韩在线观看h| 又爽又黄a免费视频| 日本一二三区视频观看| 一级毛片 在线播放| 2022亚洲国产成人精品| 九九在线视频观看精品| 久久久欧美国产精品| 高清毛片免费看| 欧美bdsm另类| 亚洲最大成人中文| 亚洲国产精品国产精品| 久久午夜福利片| 噜噜噜噜噜久久久久久91| 亚洲国产欧美在线一区| 看非洲黑人一级黄片| av在线蜜桃| 日韩亚洲欧美综合| 丰满人妻一区二区三区视频av| 亚洲av中文av极速乱| 亚洲av二区三区四区| 国语对白做爰xxxⅹ性视频网站| 国产高清国产精品国产三级 | 亚洲人成网站在线观看播放| 美女国产视频在线观看| 久久久亚洲精品成人影院| 国产精品一区二区在线观看99| 不卡视频在线观看欧美| 日日摸夜夜添夜夜添av毛片| 高清视频免费观看一区二区| 免费人成在线观看视频色| 嘟嘟电影网在线观看| 天天躁日日操中文字幕| 一边亲一边摸免费视频| 肉色欧美久久久久久久蜜桃 | 国产色婷婷99| 国产精品久久久久久久久免| 国产成人a∨麻豆精品| 日本wwww免费看| 免费av观看视频| 久久久久久久大尺度免费视频| 欧美成人精品欧美一级黄| 涩涩av久久男人的天堂| 人人妻人人看人人澡| 在线精品无人区一区二区三 | 高清在线视频一区二区三区| 亚洲av成人精品一区久久| 国产爽快片一区二区三区| 国产精品麻豆人妻色哟哟久久| 韩国高清视频一区二区三区| 亚洲欧美日韩另类电影网站 | 国内少妇人妻偷人精品xxx网站| 日本黄大片高清| 2022亚洲国产成人精品| 蜜臀久久99精品久久宅男| 免费观看a级毛片全部| 久久精品国产自在天天线| 日韩大片免费观看网站| 男女边吃奶边做爰视频| 校园人妻丝袜中文字幕| 国精品久久久久久国模美| 成人无遮挡网站| 久久午夜福利片| 国产高清有码在线观看视频| 麻豆国产97在线/欧美| 3wmmmm亚洲av在线观看| 精品人妻视频免费看| 色综合色国产| 精品一区二区三卡| 丝瓜视频免费看黄片| 简卡轻食公司| 午夜福利视频1000在线观看| 久久6这里有精品| 噜噜噜噜噜久久久久久91| 久久人人爽人人爽人人片va| 免费观看a级毛片全部| 欧美精品人与动牲交sv欧美| .国产精品久久| 日韩视频在线欧美| 在线亚洲精品国产二区图片欧美 | 日韩,欧美,国产一区二区三区| 中文字幕免费在线视频6| 久久精品久久久久久久性| 国产有黄有色有爽视频| 黑人高潮一二区| 国产精品不卡视频一区二区| 嫩草影院入口| 免费不卡的大黄色大毛片视频在线观看| 欧美精品一区二区大全| 日韩,欧美,国产一区二区三区| 内地一区二区视频在线| 亚洲精品一二三| 欧美日本视频| 午夜激情久久久久久久| 国产午夜精品久久久久久一区二区三区| 国产精品不卡视频一区二区| 成人鲁丝片一二三区免费| 亚洲精品国产成人久久av| 91aial.com中文字幕在线观看| 伦精品一区二区三区| 香蕉精品网在线| 一二三四中文在线观看免费高清| 黄片无遮挡物在线观看| 免费大片黄手机在线观看| av在线app专区| 成人特级av手机在线观看| 国产一区二区三区av在线| 亚洲经典国产精华液单| 日本一本二区三区精品| 啦啦啦啦在线视频资源| 亚洲精品影视一区二区三区av| 国产成人午夜福利电影在线观看| 九九久久精品国产亚洲av麻豆| 建设人人有责人人尽责人人享有的 | 亚洲精品乱码久久久久久按摩| 精品视频人人做人人爽| 免费黄色在线免费观看| 日日摸夜夜添夜夜爱| 少妇丰满av| 夜夜看夜夜爽夜夜摸| av专区在线播放| 性色av一级| 男男h啪啪无遮挡| 日本猛色少妇xxxxx猛交久久| 国产乱人视频| 极品少妇高潮喷水抽搐| 少妇 在线观看| av国产久精品久网站免费入址| 国产精品国产三级国产av玫瑰| 国产成人精品福利久久| 久久久a久久爽久久v久久| 六月丁香七月| 久久久精品94久久精品| 亚洲av成人精品一二三区| 99视频精品全部免费 在线| 看黄色毛片网站| 久久精品国产亚洲av天美| 午夜老司机福利剧场| 麻豆国产97在线/欧美| 日本午夜av视频| 亚洲第一区二区三区不卡| 国产探花极品一区二区| 别揉我奶头 嗯啊视频| 水蜜桃什么品种好| 欧美激情在线99| av又黄又爽大尺度在线免费看| 精品熟女少妇av免费看| 国产色婷婷99| 亚洲婷婷狠狠爱综合网| 18禁裸乳无遮挡动漫免费视频 | 欧美激情久久久久久爽电影| 国产成人精品一,二区| 少妇人妻精品综合一区二区| 99热国产这里只有精品6| 身体一侧抽搐| 欧美一级a爱片免费观看看| 亚洲精品久久午夜乱码| 美女主播在线视频| 两个人的视频大全免费| 丝袜喷水一区| 禁无遮挡网站| 国产高清国产精品国产三级 | 一区二区三区精品91| 久久久亚洲精品成人影院| 亚洲精品日韩在线中文字幕| 日韩人妻高清精品专区| 黄片无遮挡物在线观看| 中文欧美无线码| 搡女人真爽免费视频火全软件| 少妇人妻久久综合中文| 国产精品av视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久久久电影网| 天堂中文最新版在线下载 | 内地一区二区视频在线| 香蕉精品网在线| av女优亚洲男人天堂| 少妇 在线观看| 观看免费一级毛片| 国产探花极品一区二区| 国内少妇人妻偷人精品xxx网站| 亚洲美女视频黄频| 色视频在线一区二区三区| 欧美xxxx性猛交bbbb| 国产精品麻豆人妻色哟哟久久| 伦精品一区二区三区| 久久精品人妻少妇| xxx大片免费视频| 性插视频无遮挡在线免费观看| 少妇 在线观看| 国产欧美日韩一区二区三区在线 | 尾随美女入室| www.av在线官网国产| 女人久久www免费人成看片| 黄色一级大片看看| 黄色欧美视频在线观看| 亚洲久久久久久中文字幕| 乱系列少妇在线播放| 五月天丁香电影| 综合色av麻豆| 久久99热这里只频精品6学生| 大片免费播放器 马上看| 男人和女人高潮做爰伦理| 在现免费观看毛片| 亚洲成人一二三区av| 涩涩av久久男人的天堂| 免费人成在线观看视频色| 另类亚洲欧美激情| 成年av动漫网址| 啦啦啦中文免费视频观看日本| 青青草视频在线视频观看| 天天躁日日操中文字幕| 亚洲va在线va天堂va国产| 老司机影院成人| 插阴视频在线观看视频| 日本猛色少妇xxxxx猛交久久| 久久97久久精品| 国产成人a∨麻豆精品| 91久久精品国产一区二区成人| 日日摸夜夜添夜夜爱| 亚洲av免费高清在线观看| 少妇 在线观看| 久久精品夜色国产| 狂野欧美白嫩少妇大欣赏| 久久久色成人| 又爽又黄无遮挡网站| 精品亚洲乱码少妇综合久久| 下体分泌物呈黄色| 少妇人妻精品综合一区二区| 亚洲精品国产av蜜桃| 亚洲人成网站在线播| 午夜精品国产一区二区电影 | 国产成人免费观看mmmm| 国产精品一及| 一级毛片 在线播放| 一级毛片电影观看| 亚州av有码| 欧美日韩一区二区视频在线观看视频在线 | 成人毛片a级毛片在线播放| 街头女战士在线观看网站| 精品久久久久久久久亚洲| 亚洲av.av天堂| 久久久成人免费电影| 男女国产视频网站| 欧美成人a在线观看| 99久久精品国产国产毛片| 色婷婷久久久亚洲欧美| 国产亚洲av片在线观看秒播厂| 白带黄色成豆腐渣| 亚洲国产高清在线一区二区三| 免费av毛片视频| 久久人人爽av亚洲精品天堂 | 自拍偷自拍亚洲精品老妇| 日本wwww免费看| 国产精品国产三级国产av玫瑰| 亚洲美女视频黄频| 日韩人妻高清精品专区| 少妇被粗大猛烈的视频| 国产国拍精品亚洲av在线观看| 搡老乐熟女国产| 一级片'在线观看视频| 日本色播在线视频| 日日撸夜夜添| 夜夜看夜夜爽夜夜摸| 欧美成人精品欧美一级黄| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜爱| 午夜福利视频1000在线观看| 一本色道久久久久久精品综合| 黑人高潮一二区| 欧美少妇被猛烈插入视频| 69人妻影院| 18禁裸乳无遮挡免费网站照片| 国产精品国产三级国产专区5o| 中文字幕制服av| 亚洲欧美日韩另类电影网站 | 国产伦精品一区二区三区视频9| 国产亚洲91精品色在线| 成人漫画全彩无遮挡| 国产精品女同一区二区软件| 国产亚洲av嫩草精品影院| 狂野欧美激情性bbbbbb| 日本三级黄在线观看| 亚洲精品日本国产第一区| 嫩草影院入口| 三级国产精品片| 亚洲精品乱码久久久v下载方式| 国产在线男女| 看十八女毛片水多多多| 亚洲激情五月婷婷啪啪| 亚洲怡红院男人天堂| 国产精品女同一区二区软件| 大片免费播放器 马上看| 在线 av 中文字幕| 六月丁香七月| 深爱激情五月婷婷| 久久6这里有精品| av卡一久久| 精品久久久久久久久av| 少妇的逼好多水| 国产一区二区亚洲精品在线观看| 在线观看免费高清a一片| 亚洲高清免费不卡视频| 亚洲国产最新在线播放| 国产中年淑女户外野战色| 观看免费一级毛片| 不卡视频在线观看欧美| 国产男女超爽视频在线观看| 在线观看一区二区三区| 久久6这里有精品| 中文字幕制服av| 在线精品无人区一区二区三 | 18+在线观看网站| 看免费成人av毛片| 国产免费一级a男人的天堂| 丰满乱子伦码专区| 国产 精品1| 99热这里只有精品一区| 亚洲经典国产精华液单| 中国三级夫妇交换| 日本一二三区视频观看| 中文字幕制服av| 不卡视频在线观看欧美| 精品人妻视频免费看| a级毛片免费高清观看在线播放| 亚洲av免费高清在线观看| 中文字幕亚洲精品专区| 美女脱内裤让男人舔精品视频| 草草在线视频免费看| 最近中文字幕高清免费大全6| 99久国产av精品国产电影| 97超碰精品成人国产| 国产毛片a区久久久久| 亚洲av国产av综合av卡| xxx大片免费视频| 久久ye,这里只有精品| 亚洲欧美日韩另类电影网站 | 久久久精品94久久精品| 免费看av在线观看网站| 亚洲国产av新网站| 亚洲在久久综合| 国产成人免费观看mmmm| 欧美成人一区二区免费高清观看| 国产精品麻豆人妻色哟哟久久| 亚洲电影在线观看av| 亚洲欧洲国产日韩| 好男人在线观看高清免费视频| 免费少妇av软件| 亚洲真实伦在线观看| 精品久久久久久久久av| 国产精品蜜桃在线观看| 九九在线视频观看精品| 亚洲av免费在线观看| 国产精品一区www在线观看| 久久久欧美国产精品| 国内少妇人妻偷人精品xxx网站| 久久久久久久精品精品| 国产精品蜜桃在线观看| www.av在线官网国产| 国产男人的电影天堂91| 亚洲欧美日韩另类电影网站 | av在线天堂中文字幕| 99热这里只有是精品在线观看| 一级毛片aaaaaa免费看小| 日产精品乱码卡一卡2卡三| 嫩草影院新地址| 亚洲美女搞黄在线观看| 午夜福利高清视频| 久久久精品欧美日韩精品| 天堂俺去俺来也www色官网| 色网站视频免费| 日韩三级伦理在线观看| 久久久久性生活片| 国产毛片在线视频| 亚洲av.av天堂| 日韩亚洲欧美综合| 色5月婷婷丁香| 亚洲成人中文字幕在线播放| 男插女下体视频免费在线播放| 国产欧美日韩精品一区二区| av播播在线观看一区| 69人妻影院| 国产乱人偷精品视频| 男女边摸边吃奶| 国产精品精品国产色婷婷| av在线app专区| 制服丝袜香蕉在线| 精品久久国产蜜桃| 中文字幕免费在线视频6| 国产精品久久久久久久电影| 久久精品久久精品一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国产午夜福利久久久久久| 视频中文字幕在线观看| av在线老鸭窝| 日本与韩国留学比较| 国产片特级美女逼逼视频| kizo精华| 哪个播放器可以免费观看大片| 最后的刺客免费高清国语| 国产成人精品一,二区| 国产女主播在线喷水免费视频网站| 久久久精品94久久精品| 建设人人有责人人尽责人人享有的 | 99热这里只有是精品50| 日本免费在线观看一区| 夫妻午夜视频| 亚洲国产成人一精品久久久| 亚洲精品aⅴ在线观看| av在线天堂中文字幕| 在线亚洲精品国产二区图片欧美 | 韩国高清视频一区二区三区| 国产精品一区二区三区四区免费观看| 成人国产麻豆网| 亚洲欧美一区二区三区国产| 久久久久久久久久成人| 少妇人妻 视频| 少妇熟女欧美另类| 一个人看视频在线观看www免费| 婷婷色综合大香蕉| 久久午夜福利片| 中文在线观看免费www的网站| 久久鲁丝午夜福利片| 少妇丰满av| 国产精品99久久99久久久不卡 | 欧美最新免费一区二区三区| 99热网站在线观看| 秋霞在线观看毛片| 三级经典国产精品| 国产男女内射视频| 国产久久久一区二区三区| 国产成人午夜福利电影在线观看| 亚洲精品久久久久久婷婷小说| 夫妻午夜视频| 少妇 在线观看| 欧美xxⅹ黑人| 亚洲av一区综合| 涩涩av久久男人的天堂| 精华霜和精华液先用哪个| 日韩免费高清中文字幕av| 成人一区二区视频在线观看| 久久久成人免费电影| 少妇 在线观看| 女人久久www免费人成看片| 国产亚洲最大av| 久久久久久久国产电影| 99精国产麻豆久久婷婷| 综合色av麻豆| 蜜桃亚洲精品一区二区三区| 高清午夜精品一区二区三区| 久久影院123| 午夜精品国产一区二区电影 | 午夜免费男女啪啪视频观看| 久久久午夜欧美精品| 国产成人精品一,二区| 80岁老熟妇乱子伦牲交| 国产成年人精品一区二区| 蜜桃亚洲精品一区二区三区| 久久久久久久久久人人人人人人| 在线观看国产h片| 国产伦精品一区二区三区四那| 亚洲精品乱码久久久v下载方式| 免费高清在线观看视频在线观看| 国产精品久久久久久精品电影小说 | 搡老乐熟女国产| 欧美日韩亚洲高清精品| 大话2 男鬼变身卡| 成人无遮挡网站| 日韩在线高清观看一区二区三区| 久久久久久伊人网av| 男人狂女人下面高潮的视频| 国产免费福利视频在线观看| 日韩成人av中文字幕在线观看| 亚洲人成网站在线播| 麻豆久久精品国产亚洲av| 国产精品精品国产色婷婷| 成年av动漫网址| 日韩欧美精品v在线| 禁无遮挡网站| 欧美+日韩+精品| 国产亚洲最大av| 熟妇人妻不卡中文字幕| 在线看a的网站| 日本一二三区视频观看| 80岁老熟妇乱子伦牲交| 精品久久久久久久久亚洲| 少妇 在线观看| 日韩强制内射视频| 三级经典国产精品| 国产午夜福利久久久久久| 国产伦精品一区二区三区视频9| 国产亚洲av嫩草精品影院| 国语对白做爰xxxⅹ性视频网站| 51国产日韩欧美| av专区在线播放| 伦精品一区二区三区| av网站免费在线观看视频| 美女国产视频在线观看| 99热这里只有是精品在线观看| 欧美日韩综合久久久久久| 色哟哟·www| 91精品国产九色| 国产高清不卡午夜福利| 亚洲美女视频黄频| 亚洲欧美清纯卡通| 少妇的逼水好多| 视频区图区小说| 欧美日韩国产mv在线观看视频 | 性色av一级| 国产毛片在线视频| 日韩欧美精品v在线| 精品少妇久久久久久888优播| 亚洲图色成人| 成人毛片60女人毛片免费| 一级毛片电影观看| 亚洲精品乱码久久久v下载方式| 肉色欧美久久久久久久蜜桃 | 亚洲在线观看片| 日本午夜av视频| 建设人人有责人人尽责人人享有的 | 亚洲av中文av极速乱| 国产亚洲av嫩草精品影院| 成人免费观看视频高清| 国产国拍精品亚洲av在线观看| 国产综合懂色| av播播在线观看一区| 99久久精品国产国产毛片| 国产精品一区二区在线观看99| 久久精品综合一区二区三区| 下体分泌物呈黄色| 色综合色国产| 啦啦啦在线观看免费高清www| 亚洲欧美中文字幕日韩二区| 成人亚洲精品一区在线观看 | 久久久久久久国产电影| 国产永久视频网站| 少妇猛男粗大的猛烈进出视频 | 国产女主播在线喷水免费视频网站| 能在线免费看毛片的网站| av在线亚洲专区| 亚洲怡红院男人天堂| 男人和女人高潮做爰伦理| 欧美最新免费一区二区三区| 亚洲最大成人av| 亚洲av中文字字幕乱码综合| 九九久久精品国产亚洲av麻豆| 国产亚洲一区二区精品| 欧美日韩综合久久久久久| 一级二级三级毛片免费看| 极品少妇高潮喷水抽搐| 国内揄拍国产精品人妻在线| 亚洲人与动物交配视频| 网址你懂的国产日韩在线| 女人十人毛片免费观看3o分钟| av专区在线播放| 99热国产这里只有精品6| 国产v大片淫在线免费观看| 纵有疾风起免费观看全集完整版| 69av精品久久久久久| 久久鲁丝午夜福利片| 亚洲内射少妇av| av在线app专区| 日日撸夜夜添| 免费看av在线观看网站| 三级国产精品欧美在线观看| 午夜激情久久久久久久| 国产淫语在线视频| 一级毛片 在线播放| 国产综合懂色| 免费不卡的大黄色大毛片视频在线观看| 国产日韩欧美亚洲二区| 日韩一本色道免费dvd| 性插视频无遮挡在线免费观看| 亚洲国产精品国产精品| 神马国产精品三级电影在线观看| 免费观看无遮挡的男女| 欧美高清成人免费视频www| 日产精品乱码卡一卡2卡三| 麻豆成人午夜福利视频| 白带黄色成豆腐渣| 欧美性猛交╳xxx乱大交人| 亚洲丝袜综合中文字幕| 男女边吃奶边做爰视频| 亚洲av成人精品一区久久| 中国三级夫妇交换| 啦啦啦啦在线视频资源| 国产视频内射| 少妇人妻 视频| 成人一区二区视频在线观看| 综合色丁香网| 美女脱内裤让男人舔精品视频| 久久久久网色| 亚洲欧美精品自产自拍| 男人爽女人下面视频在线观看| 夫妻性生交免费视频一级片| 欧美激情在线99| 国产午夜精品一二区理论片| 亚洲av免费在线观看| 国产免费又黄又爽又色| 日本wwww免费看| 热re99久久精品国产66热6| 亚洲欧美精品专区久久|