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

    時(shí)間域航空電磁擴(kuò)散特征和成像深度研究

    2016-09-29 06:48:44殷長春邱長凱劉云鶴蔡晶
    地球物理學(xué)報(bào) 2016年8期
    關(guān)鍵詞:磁偶極子煙圈感應(yīng)電流

    殷長春,邱長凱,劉云鶴,蔡晶

    吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院,長春 130026

    ?

    時(shí)間域航空電磁擴(kuò)散特征和成像深度研究

    殷長春,邱長凱,劉云鶴,蔡晶

    吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院,長春130026

    為了研究時(shí)間域航空電磁系統(tǒng)擴(kuò)散特征,基于連續(xù)性邊界條件將電磁場向下延拓得到地下介質(zhì)中各點(diǎn)處的頻率域響應(yīng),通過傅里葉變換將其變換到時(shí)間域.計(jì)算階躍波形發(fā)射時(shí)垂直磁偶極子和水平磁偶極子在均勻半空間介質(zhì)和兩層介質(zhì)中產(chǎn)生的電磁場,通過電流密度矢量圖和隨時(shí)間變化的電流密度等值線展示電磁場在地下介質(zhì)中的擴(kuò)散過程.電磁場在地下的擴(kuò)散受電阻率結(jié)構(gòu)的影響.在良導(dǎo)地區(qū)電磁場擴(kuò)散慢、衰減快,而在高阻地區(qū)電磁場擴(kuò)散快、衰減慢.垂直磁偶極子在地下產(chǎn)生的感應(yīng)電流形成一個(gè)電流環(huán)隨時(shí)間向下、向外擴(kuò)散.水平磁偶極子在地下產(chǎn)生的感應(yīng)電流形成兩個(gè)互相疊加的電流環(huán).時(shí)間域電磁場擴(kuò)散代表著真正物理意義上的電磁擴(kuò)散,因?yàn)樗故玖说叵码姶艌隹臻g分布隨時(shí)間的變化和場強(qiáng)的衰減.基于對(duì)感應(yīng)電流環(huán)的研究,我們探究了擴(kuò)散深度和時(shí)間域航空電磁數(shù)據(jù)成像深度的關(guān)系.時(shí)間域電磁場擴(kuò)散規(guī)律的研究不僅加深了對(duì)電磁擴(kuò)散特征的理解,更有助于提高航空電磁數(shù)據(jù)解釋水平.

    航空電磁;電磁擴(kuò)散;煙圈效應(yīng);成像深度

    To study the EM diffusion of a time-domain AEM system,we first calculate the frequency-domain EM field inside the earth based on the continuity boundary condition and then transform it into time-domain via a Fourier transform.We calculate the EM fields induced by a step pulse for a homogeneous half-space and a two-layered earth model and display the EM diffusion in the earth as 2D vectors or 3D animated time-varying contours (like the "smoke ring").These time-domain current rings display the true EM diffusion,because they demonstrate the spatial variation and decay of underground EM fields with time.

    Study of the EM current rings shows that EM diffusion is strongly influenced by the resistivity structure inside the earth.In a conductive earth,the current ring diffuses slowly but decays fast,while in a resistive earth it diffuses fast but decays slowly.The induced current for a vertical magnetic dipole forms a single current ring that propagates with time outward and downward,while for a horizontal magnetic dipole,the induced current forms two stacked current rings in the underground,diffusing with time into the deep earth.Based on the EM current ring,we find that the imaging depth for time-domain AEM data is about 0.55 times the EM diffusion depth.

    By showing the underground induced current as vectors and 3D time-varying animated contours,the EM diffusion and smoke ring effect are clearly observed.An animated EM smoke ring offers more information than a static contour of current density.The relation between imaging depth and diffusion depth is further confirmed by the smoke rings.Study on EM smoke ring can not only offer an insight into EM diffusion in the earth,but can also assist in AEM data interpretation.

    1 引言

    隨著航空電磁法在煤礦、石油和天然氣勘探、地下水以及環(huán)境和工程調(diào)查中得到越來越廣泛的應(yīng)用(雷棟等,2006;殷長春等,2015),航空電磁數(shù)據(jù)成像和解釋變得越來越重要.然而在航空電磁法成像技術(shù)中,一些參數(shù)比如成像深度只能通過經(jīng)驗(yàn)得到(Huang and Rudd,2008),而不能得到明確的解釋.本文通過展示電磁場在地下傳播和擴(kuò)散的過程來研究電磁場擴(kuò)散特征,并基于電磁擴(kuò)散的概念確定時(shí)間域航空電磁數(shù)據(jù)的成像深度.

    根據(jù)安培定律,通過時(shí)間域航空電磁系統(tǒng)發(fā)射線圈的電流產(chǎn)生一次磁場.當(dāng)發(fā)射電流被瞬時(shí)關(guān)斷,一次磁場消失,結(jié)合電磁感應(yīng)定律,地下產(chǎn)生了感應(yīng)電流(渦流)來阻止一次磁場的變化.感應(yīng)電流隨時(shí)間在地下傳播,分布變得越來越廣,其強(qiáng)度衰減,稱之為電磁擴(kuò)散.

    時(shí)間域航空電磁擴(kuò)散的可視化對(duì)于理解航空電磁系統(tǒng)的瞬變過程十分重要.感應(yīng)電流在地下的分布和擴(kuò)散特征是估算航空電磁系統(tǒng)的勘探深度(陳斌等,2014)、影響范圍(Beamish,2003;Yin et al.,2014)和趨膚深度(Beamish,2004)的有效手段,可進(jìn)一步為航空電磁系統(tǒng)設(shè)計(jì)和數(shù)據(jù)處理提供理論支撐.

    很多地球物理學(xué)家已經(jīng)對(duì)地面瞬變電磁系統(tǒng)的電磁擴(kuò)散現(xiàn)象做了研究.Lewis和Lee(1978)研究了位于均勻半空間地表的大回線產(chǎn)生的瞬變電場;Nabighian(1979)提出階躍波激發(fā)時(shí)在均勻半空間或者層狀大地產(chǎn)生的瞬變電磁場可以用一個(gè)和發(fā)射線圈形狀相同的、向下向外擴(kuò)散的電流環(huán)表示,組成一個(gè)“煙圈”;Hoversten和Morrison(1982)的研究表明回線線圈在四層介質(zhì)中激發(fā)的感應(yīng)電流形成一個(gè)單獨(dú)的煙圈;Wang(2002)驗(yàn)證了在層狀橫向各向異性介質(zhì)中的煙圈效應(yīng);閆述等(2002)通過二維正演模擬展示了二維靜態(tài)的感應(yīng)電流等值線.除此之外,還有一些關(guān)于頻率域航空電磁系統(tǒng)電磁場擴(kuò)散的研究.Yin和Hodges(2005)探究了頻率域電磁系統(tǒng)在均勻半空間介質(zhì)中的電磁擴(kuò)散,進(jìn)而將模型擴(kuò)展到復(fù)雜的各向異性和二、三維模型(Yin and Hodges,2007).然而,需要指出的是,截至目前關(guān)于時(shí)間域電磁場的擴(kuò)散和煙圈效應(yīng)多是基于地面電磁系統(tǒng)的靜態(tài)描述,電磁擴(kuò)散的動(dòng)態(tài)特征沒有得到真正體現(xiàn).Yin和Hodges(2007)把頻率域電磁場乘以時(shí)諧因子eiω t,然后展示了三維動(dòng)態(tài)的電磁場等值線和矢量圖,由于其假設(shè)電磁場隨相位變化(而不是隨時(shí)間變化),因此沒有真正展示電磁場的擴(kuò)散和煙圈效應(yīng).

    成像方法簡單、有效,已經(jīng)成為航空電磁數(shù)據(jù)處理的基本流程之一.然而,成像深度的確定卻困擾了地球物理學(xué)家很多年.航空電磁成像方法將利用轉(zhuǎn)換算法得到的視電阻率和相關(guān)的成像深度聯(lián)系起來以實(shí)現(xiàn)電導(dǎo)率深度成像(Huang and Rudd,2008;李永興等,2010;朱凱光等,2010;陳小紅和段奶軍,2012;毛立峰,2013).多年來,成像深度依據(jù)經(jīng)驗(yàn)近似為擴(kuò)散深度的一半,然而沒有人能解釋成像深度確定的真正原因.

    本文通過展示電磁場在地下隨時(shí)間的傳播過程研究時(shí)間域電磁場的擴(kuò)散特征和煙圈效應(yīng).我們首先根據(jù)Yin和Maurer(2001)提出的基于連續(xù)性邊界條件的算法計(jì)算磁偶極子源在地下介質(zhì)中產(chǎn)生的頻率域電磁場,然后利用傅里葉變換將其變換到時(shí)間域.通過展示地下感應(yīng)電流隨時(shí)間的動(dòng)態(tài)變化,可以清晰地看到電磁場隨時(shí)間的擴(kuò)散過程和煙圈效應(yīng).最后,本文嘗試通過煙圈效應(yīng)的概念來建立航空電磁數(shù)據(jù)成像深度和擴(kuò)散深度的關(guān)系.我們考慮垂直和水平磁偶極子兩種發(fā)射源,以均勻半空間模型和兩層介質(zhì)模型為例來探究時(shí)間域航空電磁系統(tǒng)的擴(kuò)散特征.

    2 正演計(jì)算方法

    圖1表示位于層狀各向同性介質(zhì)上方的航空電磁系統(tǒng)發(fā)射源.為簡單起見,我們在圖中僅給出垂直磁偶極子.定義直角坐標(biāo)系的原點(diǎn)位于發(fā)射源在地表投影處,x軸和y軸位于地表,z軸垂直向下.偶極子位于(0,0,-h0)處,h0代表偶極子的高度.

    圖1 層狀介質(zhì)上方的垂直磁偶極子模型Fig.1 A vertical magnetic dipole over a layered earth

    根據(jù)Weidelt(1991),對(duì)于水平磁偶極子,空中(r,φ,z)處徑向、φ方向和垂向電場分別為

    (1)

    而垂直磁偶極子的電場φ分量為

    (2)

    (3)

    其中J0和J1分別是0階和1階貝塞爾函數(shù),η(λ)是和波數(shù)λ以及地電參數(shù)有關(guān)的核函數(shù)(Weidelt,1991).

    根據(jù)電磁場在地表和層界面處的連續(xù)性,我們可以采用向下延拓算法得到地下各點(diǎn)處的電磁響應(yīng).結(jié)合Yin和Hodges(2005)的算法,水平磁偶極子在第l層(r,φ,z)處的電場為

    Ez=0,

    (4)

    而對(duì)于垂直磁偶極子

    (5)

    式中Fl(λ)是依據(jù)連續(xù)性邊界條件得到的核函數(shù),其實(shí)質(zhì)是振幅系數(shù).

    在求得柱坐標(biāo)下的電磁場之后,將其投影到直角坐標(biāo)系,并利用歐姆定律J=σE得到地下的電流密度J,式中σ代表電導(dǎo)率.

    如果已知頻率域場值F(ω),可通過傅里葉變換將其轉(zhuǎn)換到時(shí)間域,得到t時(shí)刻的階躍電磁響應(yīng)(殷長春等,2013):

    (6)

    對(duì)式(6)進(jìn)行變換得到

    (7)

    式中J1/2是0.5階貝塞爾函數(shù).(7)式可以通過漢克爾變換濾波算法進(jìn)行計(jì)算.

    3 時(shí)間域電磁場的擴(kuò)散特征

    為了對(duì)時(shí)間域航空電磁系統(tǒng)的電磁場擴(kuò)散進(jìn)行可視化,本文計(jì)算階躍波激勵(lì)時(shí)從斷電0.002 ms到5.0 ms共900個(gè)時(shí)刻地下各點(diǎn)的感應(yīng)電流,然后展示電流密度矢量和隨時(shí)間變化的三維電流密度等值線.本文首先選擇Yin和Hodges(2007)使用的電阻率為100 Ωm的均勻半空間模型來比較時(shí)間域和頻率域電磁場擴(kuò)散過程的差異.圖2表示斷電0.5 ms時(shí)刻磁偶極子源在地下產(chǎn)生的感應(yīng)電流密度分布(建議下載附件中的相關(guān)文件,以便更好地觀察電流密度的動(dòng)態(tài)分布特征).其中,圖2a和圖2b分別表示垂直磁偶極子和水平磁偶極子源在地表產(chǎn)生的感應(yīng)電流密度矢量,箭頭的中心點(diǎn)代表電流密度的計(jì)算點(diǎn)位置,箭頭長度和電流密度的強(qiáng)度成正比;圖2c和圖2d分別表示垂直和水平磁偶極子源在地下產(chǎn)生的感應(yīng)電流密度的三維分布.為了能清晰地看到感應(yīng)電流在地下的分布形態(tài),我們切掉了圖2c和圖2d的右下角部分.

    從圖2可以看出:1)對(duì)于垂直磁偶極子,地下感應(yīng)電流形成一個(gè)圓環(huán),對(duì)應(yīng)于垂直磁偶極子的磁力線垂直注入地下;2)對(duì)于水平磁偶極子,地下感應(yīng)電流形成兩個(gè)疊加的電流環(huán),對(duì)應(yīng)于水平磁偶極子的磁力線注入地下或從地下流出;3)在層狀各向同性介質(zhì)中沒有垂向感應(yīng)電流,垂直和水平磁偶極子源在地下產(chǎn)生的感應(yīng)電流都水平流動(dòng).

    本文利用Yin和Hodges(2007)的算法計(jì)算了頻率為f=900 Hz的垂直磁偶極子源在電阻率為100 Ωm的均勻半空間中產(chǎn)生的電流密度x分量的實(shí)部,并和時(shí)間域感應(yīng)電流密度x分量相比較,如圖3所示.對(duì)比圖3a和圖3b可以發(fā)現(xiàn),時(shí)間域感應(yīng)電流隨著時(shí)間增加不斷向下向外擴(kuò)散、強(qiáng)度衰減,而頻率域感應(yīng)電流周期性向外擴(kuò)散,其強(qiáng)度和極性也發(fā)生周期性變化.因此,我們得出結(jié)論,時(shí)間域感應(yīng)電流的擴(kuò)散過程表征了真正物理意義上的電磁擴(kuò)散,而頻率域這種周期性 “擴(kuò)散”現(xiàn)象是由于頻率域電磁場乘以了相位eiω t的緣故(Yin and Hodges,2005).

    為進(jìn)一步研究電磁場在層狀介質(zhì)模型中的擴(kuò)散特征,本文設(shè)計(jì)了一個(gè)兩層地電模型:第一層電阻率為50 Ωm、厚度為200 m,第二層電阻率為100 Ωm.

    圖2 斷電0.5 ms時(shí)刻磁偶極子源在均勻半空間產(chǎn)生的感應(yīng)電流分布(a)垂直磁偶極子和(b)水平磁偶極子在地表產(chǎn)生的電流密度矢量分布圖;(c)垂直磁偶極子和(d)水平磁偶極子在地下產(chǎn)生的感應(yīng)電流密度.Fig.2 Distribution of current density inside the earth for a vertical and a horizontal magnetic dipole over a homogeneous half-space at the decay time of 0.5 msCurrent density at the earth surface for (a)vertical magnetic dipole and (b)horizontal magnetic dipole;Current density inside the earth for (c)vertical magnetic dipole and (d)horizontal magnetic dipole.

    圖3 垂直磁偶極子在電阻率為100 Ωm的均勻半空間產(chǎn)生的感應(yīng)電流(a)頻率域:頻率f=900 Hz;(b)時(shí)間域:時(shí)間t=0.5 ms.Fig.3 Current density in contours for a vertical magnetic dipole over a homogeneous half-space of 100 Ωm(a)900 Hz for frequency-domain;(b)Time-domain at the decay time of 0.5 ms.

    圖4 斷電后0.5 ms時(shí)刻垂直磁偶極子在兩層介質(zhì)中產(chǎn)生的感應(yīng)電流分布Fig.4 Current density in contours for a vertical magnetic dipole over a two-layer earth model at the decay time of 0.5 ms

    垂直磁偶極子位于地表上方30 m.圖4表示垂直磁偶極子在兩層介質(zhì)中產(chǎn)生的感應(yīng)電流分布.從圖4中可以看出:1)電流密度分布和電阻率界面緊密相關(guān),在第一層和第二層的邊界上存在電流密度的突變,從電流密度分布可以清晰地分辨出地下電性分界面.產(chǎn)生這種電流密度突變的物理原因很簡單.事實(shí)上,在地下電性分界面處垂向電流連續(xù)(都是0),但是橫向電流密度不連續(xù),導(dǎo)致了總電流密度發(fā)生突變;2)觀察電磁場隨時(shí)間的變化可以進(jìn)一步發(fā)現(xiàn),電磁場在表層良導(dǎo)層傳播速度慢,這意味著在航空電磁勘探中,需要觀測較晚的電磁信號(hào)才能探測到良導(dǎo)蓋層下的目標(biāo)體.為便于直觀研究,本文中我們僅討論了半空間和兩層介質(zhì)的情況.然而,需要指出的是,無論地下結(jié)構(gòu)多么復(fù)雜,本文中所描述的電磁場在低阻和高阻介質(zhì)中的擴(kuò)散規(guī)律是不變的.

    4 煙圈效應(yīng)

    (8)

    選取電流環(huán)的半徑r′是擴(kuò)散深度的0.55倍(下文中將會(huì)證明),則

    (9)

    因此,我們?nèi)U(kuò)散深度平面最大感應(yīng)電流密度的三分之一作為水平磁偶極子電流環(huán)的閾值.

    圖5 水平磁偶極子感應(yīng)電流環(huán)中心和邊界電流分布特征Fig.5 Center current and edge current for current rings of a horizontal magnetic dipole

    圖6給出斷電0.5 ms時(shí)刻垂直磁偶極子和水平磁偶極子在電阻率為100 Ωm的均勻半空間地下產(chǎn)生的電流環(huán).由圖6可以看出:1)垂直磁偶極子的感應(yīng)電流在地下只形成一個(gè)電流環(huán)(煙圈);2)水平磁偶極子的感應(yīng)電流在地下形成兩個(gè)疊加的煙圈,其中心位于偶極子正下方;3)隨著時(shí)間增加,垂直和水平磁偶極子的電流環(huán)都向下和向外傳播,逐漸變得更大,同時(shí)振幅衰減,體現(xiàn)出感應(yīng)電流隨時(shí)間的擴(kuò)散特征.

    圖6 斷電0.5 ms時(shí)刻(a)垂直磁偶極子和(b)水平磁偶極子在均勻半空間地下產(chǎn)生的電流環(huán)Fig.6 Current rings for (a)a vertical magnetic dipole and (b)a horizontal magnetic dipole over a homogeneous half-space of 100 Ωm at the decay time of 0.5 ms

    圖7 斷電1.0 ms時(shí)刻垂直磁偶極子在兩層介質(zhì)中產(chǎn)生的電流環(huán)Fig.7 Current ring for a vertical magnetic dipole over a two-layer earth model at the decay time of 1.0 ms

    圖7表示斷電1.0 ms時(shí)垂直磁偶極子在兩層介質(zhì)中產(chǎn)生的電流環(huán).兩層介質(zhì)模型的參數(shù)同圖4.從圖7可以看出:1)即使對(duì)于層狀模型,也只存在一個(gè)電流環(huán).這和Yin和Hodges(2005,2007)所描述的頻率域感應(yīng)電流“煙圈”發(fā)生周期性變化完全不同;2)斷電后早期煙圈只在第一層傳播,隨著時(shí)間推移,向下傳播到第二層,界面處感應(yīng)電流連續(xù)性遭到破壞;3)在地表良導(dǎo)層,煙圈傳播慢,衰減較快;在基底高阻層,煙圈傳播快,衰減較慢;4)電磁場先在地表良導(dǎo)層中傳播,經(jīng)過一定時(shí)間后越過界面進(jìn)入第二層.由于上下兩層中電磁場傳播時(shí)間的差異,造成下半空間煙圈范圍較小.

    5 成像深度

    在航空電磁法領(lǐng)域,成像是一種處理航空電磁數(shù)據(jù)簡單有效的技術(shù).常用的方法包括電導(dǎo)率深度成像(CDI)和電導(dǎo)率深度轉(zhuǎn)換(CDT).為了對(duì)地下結(jié)構(gòu)進(jìn)行成像,通常將產(chǎn)生同樣響應(yīng)的均勻半空間模型的電阻率作為視電阻率,把擴(kuò)散深度的一半近似作為成像深度,構(gòu)造視電阻率-深度剖面.然而,到目前為止沒有人能解釋成像深度是如何確定的.本文從煙圈效應(yīng)出發(fā),建立最大電流位置和擴(kuò)散深度的關(guān)系,并基于這個(gè)關(guān)系確定航空電磁的成像深度.

    6 結(jié)論

    本文將磁偶極子源在地下產(chǎn)生的感應(yīng)電流可視化為矢量圖和三維動(dòng)態(tài)的等值面,從而能夠清晰地觀察時(shí)間域電磁場的擴(kuò)散過程和煙圈效應(yīng).垂直磁偶極子在均勻半空間介質(zhì)中產(chǎn)生的感應(yīng)電流形成一個(gè)電流環(huán),而水平磁偶極子在均勻半空間地下產(chǎn)生的感應(yīng)電流形成兩個(gè)疊加的電流環(huán).在層狀介質(zhì)中,感應(yīng)電流形成的電流環(huán)從表層向下傳播擴(kuò)散到底層.頻率域電磁擴(kuò)散呈現(xiàn)周期性特征,而時(shí)間域電磁場隨著時(shí)間向下、向外傳播,范圍變大,強(qiáng)度變?nèi)?因此,時(shí)間域擴(kuò)散是真正意義上的電磁擴(kuò)散.在良導(dǎo)層電磁場傳播慢、衰減快,而在高阻層電磁場傳播快、衰減慢.通過分析最大感應(yīng)電流所在深度和擴(kuò)散深度的關(guān)系,本文確定了時(shí)間域航空電磁法的成像深度是擴(kuò)散深度的0.55倍,有效地揭示了傳統(tǒng)時(shí)間域航空電磁數(shù)據(jù)成像中選擇的經(jīng)驗(yàn)參數(shù)的原因.由于三維動(dòng)態(tài)電流分布和煙圈效應(yīng)比靜態(tài)的電流密度平面等值線能提供更多的電磁擴(kuò)散信息,本文研究可為航空電磁系統(tǒng)設(shè)計(jì)(影響范圍、分辨率和勘探深度等)提供理論依據(jù),這將是我們未來的研究重點(diǎn).

    圖8 對(duì)于不同均勻半空間模型成像深度和擴(kuò)散深度的關(guān)系(a)垂直磁偶極子;(b)水平磁偶極子.Fig.8 Relationship between imaging depth and diffusion depth for different homogeneous half-space models(a)Vertical magnetic dipole;(b)Horizontal magnetic dipole.

    致謝感謝兩位評(píng)審專家提出的寶貴的修改意見.同時(shí)感謝編輯部的辛苦工作.

    References

    Beamish D.2003.Airborne EM footprints.Geophysical Prospecting,51(1):49-60.Beamish D.2004.Airborne EM skin depths.Geophysical Prospecting,52(5):439-449.Chen B,Mao L F,Liu G D.2014.The estimated prospecting depth of CHTEM-I system by the method of diffusion electric field.Chinese J.Geophys.(in Chinese),57(1):303-309,doi:10.6038/cjg20140126.

    Chen X H,Duan N J.2012.Study on fast imaging of airborne time-domain electromagnetic data.Progress in Geophys.(in Chinese),27(5):2123-2127,doi:10.6038/j.issn.1004-2903.2012.05.037.Hoversten G M,Morrison H F.1982.Transient fields of a current loop source above a layered earth.Geophysics,47(7):1068-1077,doi:10.1190/1.1441370.Huang H P,Rudd J.2008.Conductivity-depth imaging of helicopter-borne TEM data based on a pseudolayer half-space model.Geophysics,73(3):F115-F120.

    Lei D,Hu X Y,Zhang S F.2006.Development status of airborne electromagnetic.Contributions to Geology and Mineral Resources Research (in Chinese),21(1):40-44.Lewis R,Lee T.1978.The transient electric fields about a loop on a halfspace.Exploration Geophysics,9(4):173-177,doi:10.1071/EG978173.

    Li Y X,Qiang J K,Tang J T.2010.A research on 1-D forward and inverse airborne transient electromagnetic method.Chinese J.Geophys.(in Chinese),53(3):751-759,doi:10.3969/j.issn.0001-5733.2010.03.031.

    Mao L F.2013.Conductivity-depth imaging algorithm for central-loop helicopter TEM.CT Theory and Applications (in Chinese),22(3):429-437.Nabighian M N.1979.Quasi-static transient response of a conducting half-space-An approximate representation.Geophysics,44(10):1700-1705,doi:10.1190/1.1440931.Wang T.2002.The electromagnetic smoke ring in a transversely isotropic medium.Geophysics,67(6):1779-1789,doi:10.1190/1.1527078.

    Weidelt P.1991.Introduction into electromagnetic sounding.Lecture manuscript,Technical University of Braunschweig.

    Yan S,Chen M S,Fu J M.2002.Direct time-domain numerical analysis of transient electromagnetic fields.Chinese J.Geophys.(in Chinese),45(2):275-284.Yin C,Maurer H M.2001.Electromagnetic induction in a layered earth with arbitrary anisotropy.Geophysics,66(5):1405-1416.

    Yin C C,Hodges G.2005.Four dimensional visualization of EM fields for a helicopter EM system.//75th Annual International Meeting,SEG,Expanded Abstracts,595-598.Yin C C,Hodges G.2007.3D animated visualization of EM diffusion for a frequency-domain helicopter EM system.Geophysics,72(1):F1-F7,doi:10.1190/1.2374706.Yin C C,Huang W,Ben F.2013.The full-time electromagnetic modeling for time-domain airborne electromagnetic systems.Chinese J.Geophys.(in Chinese),56(9):3153-3162,doi:10.6038/cjg20130928.

    Yin C C,Huang X,Liu Y H,et al.2014.Footprint for frequency-domain airborne electromagnetic systems.Geophysics,79(6):E243-E254.

    Yin C C,Zhang B,Liu Y H,et al.2015.Review on airborne EM technology and developments.Chinese J.Geophys.(in Chinese),58(8):2637-2653,doi:10.6038/cjg20150804.Zhu K G,Lin J,Han Y H,et al.2010.Research on conductivity depth imaging of time domain helicopter-borne electromagnetic data based on neural network.Chinese J.Geophys.(in Chinese),53(3):743-750,doi:10.3969/j.issn.0001-5733.2010.03.030.

    附中文參考文獻(xiàn)

    陳斌,毛立峰,劉光鼎.2014.用擴(kuò)散電場法估算CHTEM-I系統(tǒng)的探測深度.地球物理學(xué)報(bào),57(1):303-309,doi:10.6038/cjg20140126.

    陳小紅,段奶軍.2012.時(shí)間域航空電磁快速成像研究.地球物理學(xué)進(jìn)展,27(5):2123-2127,doi:10.6038/j.issn.1004-2903.2012.05.037.

    雷棟,胡祥云,張素芳.2006.航空電磁法的發(fā)展現(xiàn)狀.地質(zhì)找礦論叢,21(1):40-44.

    李永興,強(qiáng)建科,湯井田.2010.航空瞬變電磁法一維正反演研究.地球物理學(xué)報(bào),53(3):751-759,doi:10.3969/j.issn.0001-5733.2010.03.031.

    毛立峰.2013.中心回線式直升機(jī)TEM資料的電導(dǎo)率-深度成像方法.CT理論與應(yīng)用研究,22(3):429-437.

    閆述,陳明生,傅君眉.2002.瞬變電磁場的直接時(shí)域數(shù)值分析.地球物理學(xué)報(bào),45(2):275-284.殷長春,黃威,賁放.2013.時(shí)間域航空電磁系統(tǒng)瞬變?nèi)珪r(shí)響應(yīng)正演模擬.地球物理學(xué)報(bào),56(9):3153-3162,doi:10.6038/cjg20130928.殷長春,張博,劉云鶴等.2015.航空電磁勘查技術(shù)發(fā)展現(xiàn)狀及展望.地球物理學(xué)報(bào),58(8):2637-2653,doi:10.6038/cjg20150804.朱凱光,林君,韓悅慧等.2010.基于神經(jīng)網(wǎng)絡(luò)的時(shí)間域直升機(jī)電磁數(shù)據(jù)電導(dǎo)率深度成像.地球物理學(xué)報(bào),53(3):743-750,doi:10.3969/j.issn.0001-5733.2010.03.030.

    (本文編輯胡素芳)

    Time-domain electromagnetic diffusion and imaging depth for airborne electromagnetic data

    YIN Chang-Chun,QIU Chang-Kai,LIU Yun-He,CAI Jing

    College of Geo-exploration Science and Technology,Jilin University,Changchun 130026,China

    With airborne electromagnetic (AEM)systems being widely used in the mineral,oil and gas,ground water,environmental and engineering exploration,AEM data interpretation and imaging become more and more important.However,there are imaging techniques in AEM community in which parameters like imaging depth can only be determined empirically but can′t be solved physically.In this paper,we discuss a fundamental phenomenon of the electromagnetic (EM)diffusion by showing how the EM fields propagate and diffuse in the earth and based on this concept,we try to determine parameters like the imaging depth for AEM data imaging.

    Airborne electromagnetic;Electromagnetic diffusion;Smoke ring;Imaging depth

    殷長春,邱長凱,劉云鶴等.2016.時(shí)間域航空電磁擴(kuò)散特征和成像深度研究.地球物理學(xué)報(bào),59(8):3079-3086,

    10.6038/cjg20160830.

    Yin C C,Qiu C K,Liu Y H,et al.2016.Time-domain electromagnetic diffusion and imaging depth for airborne electromagnetic data.Chinese J.Geophys.(in Chinese),59(8):3079-3086,doi:10.6038/cjg20160830.

    國家自然科學(xué)基金重點(diǎn)項(xiàng)目(41530320)和面上項(xiàng)目(41274121),國家青年基金項(xiàng)目(41404093),國家重大科研裝備研究項(xiàng)目(ZDYZ2012-1-03,20130523MTEM05)聯(lián)合資助.

    殷長春,男,1965年生,教授,主要從事電磁勘探理論,特別是航空和海洋電磁方面的研究.E-mail:yinchangchun@jlu.edu.cn

    10.6038/cjg20160830

    P631

    2015-11-15,2016-03-30收修定稿

    猜你喜歡
    磁偶極子煙圈感應(yīng)電流
    楞次定律推論的巧妙應(yīng)用
    金屬軌道的形狀影響感應(yīng)電流的大小嗎
    第23和24太陽活動(dòng)周高緯地磁感應(yīng)電流分布特性
    煙圈下落不明了
    青春(2019年5期)2019-10-20 09:25:52
    張偉大的詩
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    磁偶極子跟蹤的漸進(jìn)貝葉斯濾波方法
    老爸的煙圈
    異頻法線路參數(shù)測試中工頻感應(yīng)電流的影響分析
    老爸的煙圈
    亚洲激情在线av| 女同久久另类99精品国产91| 欧美另类亚洲清纯唯美| 亚洲一区高清亚洲精品| 一进一出抽搐动态| 久久午夜亚洲精品久久| 99riav亚洲国产免费| 国产精品 国内视频| 亚洲精品久久国产高清桃花| 丝袜人妻中文字幕| 99在线视频只有这里精品首页| 后天国语完整版免费观看| 欧美日韩一级在线毛片| 中文字幕精品亚洲无线码一区| 久久久国产成人免费| 别揉我奶头~嗯~啊~动态视频| 青草久久国产| 亚洲五月婷婷丁香| 国产精品 国内视频| 欧美黑人巨大hd| 最近最新中文字幕大全免费视频| 日本a在线网址| 国产成人精品久久二区二区91| 日本免费a在线| 免费观看人在逋| 国内毛片毛片毛片毛片毛片| 91成年电影在线观看| 精品久久久久久久毛片微露脸| 国产亚洲欧美在线一区二区| 国产精品免费一区二区三区在线| 亚洲av美国av| 18禁国产床啪视频网站| 国产精品一及| 久久99热这里只有精品18| 久久精品aⅴ一区二区三区四区| 亚洲国产看品久久| 国产不卡一卡二| 国产真人三级小视频在线观看| 国产成人系列免费观看| 一区二区三区高清视频在线| 久久久久久九九精品二区国产 | 日韩 欧美 亚洲 中文字幕| 欧美av亚洲av综合av国产av| 此物有八面人人有两片| 亚洲美女黄片视频| 国产男靠女视频免费网站| 亚洲精品在线观看二区| 国产三级中文精品| 日韩免费av在线播放| 国产精品98久久久久久宅男小说| 麻豆成人av在线观看| 999精品在线视频| 亚洲成人久久性| 大型黄色视频在线免费观看| 亚洲真实伦在线观看| 色在线成人网| 日韩欧美精品v在线| 一a级毛片在线观看| av欧美777| 一进一出好大好爽视频| 88av欧美| 中文字幕人妻丝袜一区二区| 国产欧美日韩一区二区三| 婷婷六月久久综合丁香| 国产精品久久久人人做人人爽| 在线观看一区二区三区| 免费在线观看完整版高清| 亚洲美女视频黄频| 国产蜜桃级精品一区二区三区| 久久精品影院6| 90打野战视频偷拍视频| 我要搜黄色片| 久久中文字幕一级| 亚洲人成77777在线视频| 18禁美女被吸乳视频| 99久久无色码亚洲精品果冻| 亚洲五月天丁香| 男插女下体视频免费在线播放| 操出白浆在线播放| 国产69精品久久久久777片 | 亚洲人与动物交配视频| 很黄的视频免费| 人妻久久中文字幕网| 久久伊人香网站| 高潮久久久久久久久久久不卡| 淫妇啪啪啪对白视频| 免费在线观看日本一区| 老熟妇乱子伦视频在线观看| 99国产精品一区二区三区| 在线a可以看的网站| 精品国产超薄肉色丝袜足j| 国产午夜福利久久久久久| 欧美成人午夜精品| 曰老女人黄片| 久热爱精品视频在线9| 全区人妻精品视频| 亚洲熟女毛片儿| 欧美又色又爽又黄视频| 成人三级黄色视频| 少妇被粗大的猛进出69影院| 午夜免费观看网址| 19禁男女啪啪无遮挡网站| 免费观看人在逋| 国模一区二区三区四区视频 | 国产黄色小视频在线观看| 欧美黑人巨大hd| 在线观看免费日韩欧美大片| 给我免费播放毛片高清在线观看| 国产欧美日韩一区二区精品| 亚洲成a人片在线一区二区| 特大巨黑吊av在线直播| 中文字幕最新亚洲高清| 国产av麻豆久久久久久久| 香蕉久久夜色| 99国产精品一区二区三区| 特级一级黄色大片| 久久精品亚洲精品国产色婷小说| 亚洲一区二区三区不卡视频| 精品福利观看| 国产午夜福利久久久久久| 大型黄色视频在线免费观看| 欧美最黄视频在线播放免费| av在线天堂中文字幕| 国产高清有码在线观看视频 | 久久久久久国产a免费观看| 免费在线观看影片大全网站| 妹子高潮喷水视频| 18禁观看日本| 国产又色又爽无遮挡免费看| 久久精品国产亚洲av香蕉五月| 中文字幕人妻丝袜一区二区| 听说在线观看完整版免费高清| 亚洲国产欧美一区二区综合| av福利片在线| 蜜桃久久精品国产亚洲av| 天堂av国产一区二区熟女人妻 | 中文字幕人成人乱码亚洲影| 久久精品影院6| 人成视频在线观看免费观看| 欧美色视频一区免费| 国产97色在线日韩免费| 亚洲第一电影网av| 99久久99久久久精品蜜桃| 成年版毛片免费区| 嫩草影视91久久| 国产又色又爽无遮挡免费看| 露出奶头的视频| 色播亚洲综合网| 成人av在线播放网站| 又爽又黄无遮挡网站| 亚洲欧洲精品一区二区精品久久久| 亚洲午夜理论影院| 小说图片视频综合网站| 欧美日本亚洲视频在线播放| 18禁国产床啪视频网站| 露出奶头的视频| 俺也久久电影网| 可以免费在线观看a视频的电影网站| 最近最新中文字幕大全电影3| 亚洲自偷自拍图片 自拍| 精品少妇一区二区三区视频日本电影| 国产爱豆传媒在线观看 | 国产日本99.免费观看| 国产精品 欧美亚洲| 黄色毛片三级朝国网站| 欧美精品啪啪一区二区三区| 好男人电影高清在线观看| 丰满的人妻完整版| 成在线人永久免费视频| 男插女下体视频免费在线播放| 欧美成人免费av一区二区三区| 亚洲 欧美一区二区三区| 亚洲专区国产一区二区| 午夜福利18| 亚洲av美国av| 欧美成人性av电影在线观看| 午夜精品在线福利| 性欧美人与动物交配| 中文在线观看免费www的网站 | 久久久国产精品麻豆| 1024视频免费在线观看| 精品一区二区三区四区五区乱码| 熟女少妇亚洲综合色aaa.| 久久久久免费精品人妻一区二区| 久久久久久免费高清国产稀缺| 亚洲成人中文字幕在线播放| 久久精品国产清高在天天线| 色综合站精品国产| 色精品久久人妻99蜜桃| 99久久精品国产亚洲精品| 欧美日韩一级在线毛片| xxx96com| 日本熟妇午夜| 久久欧美精品欧美久久欧美| 在线免费观看的www视频| 欧美日韩乱码在线| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇的丰满在线观看| 99国产精品99久久久久| 三级男女做爰猛烈吃奶摸视频| 欧美+亚洲+日韩+国产| 亚洲国产欧洲综合997久久,| 成人午夜高清在线视频| 一进一出抽搐动态| 九色国产91popny在线| 白带黄色成豆腐渣| 久久人人精品亚洲av| 国产成年人精品一区二区| 一个人免费在线观看电影 | 久久久久国内视频| 亚洲欧美激情综合另类| 成人一区二区视频在线观看| 国产精品久久久久久久电影 | 日韩欧美国产在线观看| 岛国在线观看网站| 日本五十路高清| 亚洲精品国产精品久久久不卡| 国产亚洲av嫩草精品影院| 午夜成年电影在线免费观看| 一本综合久久免费| 国产爱豆传媒在线观看 | 日本在线视频免费播放| 欧美乱码精品一区二区三区| 精品国产超薄肉色丝袜足j| 国产主播在线观看一区二区| 99国产综合亚洲精品| 国产熟女午夜一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人一区二区三| 不卡一级毛片| 国产一区在线观看成人免费| 嫩草影院精品99| 欧美日韩福利视频一区二区| 国产又黄又爽又无遮挡在线| 欧美精品亚洲一区二区| www日本黄色视频网| 欧美日韩中文字幕国产精品一区二区三区| 十八禁人妻一区二区| 一卡2卡三卡四卡精品乱码亚洲| netflix在线观看网站| 12—13女人毛片做爰片一| 少妇被粗大的猛进出69影院| 亚洲精华国产精华精| 老司机福利观看| 在线观看舔阴道视频| 999久久久国产精品视频| 国产私拍福利视频在线观看| 99国产精品一区二区蜜桃av| 97超级碰碰碰精品色视频在线观看| 精品高清国产在线一区| 少妇粗大呻吟视频| 精品乱码久久久久久99久播| 亚洲av熟女| 人妻久久中文字幕网| 国产亚洲精品综合一区在线观看 | 两性夫妻黄色片| 正在播放国产对白刺激| 亚洲欧美精品综合一区二区三区| 18禁黄网站禁片免费观看直播| 亚洲人成77777在线视频| 亚洲欧美日韩高清专用| svipshipincom国产片| 欧美精品啪啪一区二区三区| 亚洲自拍偷在线| 91av网站免费观看| 国产区一区二久久| 国产亚洲精品一区二区www| 精品人妻1区二区| 床上黄色一级片| 男插女下体视频免费在线播放| 亚洲avbb在线观看| www日本在线高清视频| 欧美日本亚洲视频在线播放| 午夜福利免费观看在线| 99精品久久久久人妻精品| 国产探花在线观看一区二区| 精品电影一区二区在线| 日本精品一区二区三区蜜桃| av超薄肉色丝袜交足视频| 欧美成人免费av一区二区三区| 国产成+人综合+亚洲专区| ponron亚洲| 黄色视频不卡| 日本一区二区免费在线视频| 欧美最黄视频在线播放免费| 999久久久精品免费观看国产| 美女黄网站色视频| 午夜福利视频1000在线观看| 日韩欧美国产在线观看| 一级黄色大片毛片| 一区福利在线观看| 亚洲一码二码三码区别大吗| 成年女人毛片免费观看观看9| 国内精品一区二区在线观看| 精品人妻1区二区| 久久久水蜜桃国产精品网| 别揉我奶头~嗯~啊~动态视频| 窝窝影院91人妻| www.自偷自拍.com| 国产黄a三级三级三级人| 真人一进一出gif抽搐免费| 日韩欧美国产在线观看| 亚洲国产高清在线一区二区三| 18禁国产床啪视频网站| 一边摸一边抽搐一进一小说| 精品日产1卡2卡| 欧美激情久久久久久爽电影| 亚洲精品一区av在线观看| 亚洲性夜色夜夜综合| e午夜精品久久久久久久| 人妻丰满熟妇av一区二区三区| 亚洲中文字幕日韩| 国产精品久久久久久人妻精品电影| 国产成人一区二区三区免费视频网站| 精品高清国产在线一区| 嫩草影视91久久| 成在线人永久免费视频| 99在线视频只有这里精品首页| 免费在线观看成人毛片| 亚洲国产欧美网| 国产乱人伦免费视频| 久久亚洲精品不卡| 嫁个100分男人电影在线观看| 免费在线观看日本一区| 久久欧美精品欧美久久欧美| 久久久久久大精品| 久久久国产成人精品二区| 国产黄片美女视频| 91麻豆精品激情在线观看国产| 久久精品91无色码中文字幕| 最新美女视频免费是黄的| 日本黄大片高清| 久久久国产精品麻豆| 性欧美人与动物交配| 757午夜福利合集在线观看| 久久精品国产亚洲av香蕉五月| 给我免费播放毛片高清在线观看| 免费一级毛片在线播放高清视频| 色在线成人网| 麻豆av在线久日| 国产探花在线观看一区二区| 国产三级在线视频| 久久天躁狠狠躁夜夜2o2o| 久久久久久久久免费视频了| 日韩免费av在线播放| 少妇裸体淫交视频免费看高清 | 国产午夜精品论理片| 日韩精品免费视频一区二区三区| 在线观看免费日韩欧美大片| 丁香六月欧美| 草草在线视频免费看| 一本一本综合久久| 9191精品国产免费久久| 中文字幕人成人乱码亚洲影| cao死你这个sao货| 黑人欧美特级aaaaaa片| 91字幕亚洲| 亚洲av第一区精品v没综合| 精品久久久久久成人av| 午夜两性在线视频| 色综合站精品国产| 97超级碰碰碰精品色视频在线观看| 国产精品一区二区三区四区免费观看 | 真人做人爱边吃奶动态| 久久久久亚洲av毛片大全| 亚洲人成77777在线视频| 性欧美人与动物交配| 一本综合久久免费| 亚洲国产看品久久| 桃红色精品国产亚洲av| 真人做人爱边吃奶动态| av国产免费在线观看| av超薄肉色丝袜交足视频| 可以在线观看毛片的网站| 色在线成人网| 久久天躁狠狠躁夜夜2o2o| 欧美一级a爱片免费观看看 | 欧美色视频一区免费| 日韩欧美在线乱码| 午夜福利欧美成人| 国产精华一区二区三区| 亚洲中文日韩欧美视频| 国产一区二区在线观看日韩 | 亚洲中文日韩欧美视频| 国产男靠女视频免费网站| 免费搜索国产男女视频| 淫妇啪啪啪对白视频| 搡老岳熟女国产| 日本a在线网址| 亚洲人成网站高清观看| 男人舔奶头视频| 99国产精品一区二区蜜桃av| 男女那种视频在线观看| 久久精品人妻少妇| 国产精品久久久久久人妻精品电影| 日本一本二区三区精品| 嫩草影视91久久| 一级片免费观看大全| 国产av麻豆久久久久久久| 无遮挡黄片免费观看| 18美女黄网站色大片免费观看| 精品第一国产精品| 国产av在哪里看| 久久精品91无色码中文字幕| 可以在线观看毛片的网站| 亚洲美女视频黄频| 成人午夜高清在线视频| 天堂影院成人在线观看| 身体一侧抽搐| 免费人成视频x8x8入口观看| 精品国产乱子伦一区二区三区| 亚洲成人精品中文字幕电影| 在线观看免费日韩欧美大片| 久久久久久久久中文| 在线永久观看黄色视频| 欧美日韩福利视频一区二区| 曰老女人黄片| 五月玫瑰六月丁香| 夜夜夜夜夜久久久久| 国产乱人伦免费视频| 国产一区二区在线av高清观看| e午夜精品久久久久久久| 变态另类成人亚洲欧美熟女| 国产精品爽爽va在线观看网站| svipshipincom国产片| 99国产极品粉嫩在线观看| 欧美黑人巨大hd| 天堂av国产一区二区熟女人妻 | 在线观看美女被高潮喷水网站 | 亚洲av日韩精品久久久久久密| 99久久精品国产亚洲精品| 国产视频一区二区在线看| 国产一区二区三区在线臀色熟女| 国产精品永久免费网站| 神马国产精品三级电影在线观看 | 在线a可以看的网站| 淫秽高清视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 看黄色毛片网站| 国产午夜精品论理片| 欧美黄色片欧美黄色片| 久久久国产成人精品二区| 悠悠久久av| 久久天躁狠狠躁夜夜2o2o| 欧美高清成人免费视频www| 亚洲成人国产一区在线观看| 中出人妻视频一区二区| 人妻夜夜爽99麻豆av| 人成视频在线观看免费观看| 9191精品国产免费久久| 欧美黑人欧美精品刺激| 婷婷丁香在线五月| 12—13女人毛片做爰片一| 中亚洲国语对白在线视频| netflix在线观看网站| 欧美黑人欧美精品刺激| 成人三级做爰电影| 男人舔奶头视频| 香蕉国产在线看| 曰老女人黄片| 国产精品自产拍在线观看55亚洲| 一本精品99久久精品77| 国内揄拍国产精品人妻在线| 亚洲成a人片在线一区二区| 999精品在线视频| av超薄肉色丝袜交足视频| 国产熟女xx| 99在线人妻在线中文字幕| 级片在线观看| 中文字幕人成人乱码亚洲影| 精品无人区乱码1区二区| 12—13女人毛片做爰片一| 国产精品久久久久久人妻精品电影| 久久久久久久久中文| 亚洲精品色激情综合| 久久伊人香网站| 在线观看一区二区三区| 欧美色欧美亚洲另类二区| 日韩欧美 国产精品| a在线观看视频网站| 日韩精品中文字幕看吧| 黄色视频,在线免费观看| 可以在线观看毛片的网站| 日韩欧美三级三区| 美女扒开内裤让男人捅视频| 国产精品久久久久久亚洲av鲁大| 亚洲五月婷婷丁香| 岛国在线观看网站| 最近视频中文字幕2019在线8| 国产av麻豆久久久久久久| 在线免费观看的www视频| 熟女少妇亚洲综合色aaa.| 日韩欧美精品v在线| 热99re8久久精品国产| 在线观看www视频免费| 最近最新免费中文字幕在线| 国内精品一区二区在线观看| 国内毛片毛片毛片毛片毛片| 精品久久久久久久人妻蜜臀av| 国产精品一区二区三区四区久久| 午夜成年电影在线免费观看| 成人三级黄色视频| 久久久久性生活片| 12—13女人毛片做爰片一| 日韩大码丰满熟妇| 岛国在线观看网站| 日韩欧美免费精品| 亚洲国产日韩欧美精品在线观看 | 久久久久免费精品人妻一区二区| 久久久久久国产a免费观看| 免费在线观看日本一区| 免费电影在线观看免费观看| 丰满人妻熟妇乱又伦精品不卡| 欧美在线黄色| 免费在线观看亚洲国产| 欧美日韩黄片免| 日日干狠狠操夜夜爽| 日韩国内少妇激情av| 精品国产乱码久久久久久男人| 欧美高清成人免费视频www| 一级a爱片免费观看的视频| 日韩有码中文字幕| 日韩欧美 国产精品| 欧美日韩福利视频一区二区| 久久久久久免费高清国产稀缺| 男插女下体视频免费在线播放| 国产精品国产高清国产av| 国产亚洲精品久久久久5区| 久久 成人 亚洲| 亚洲国产欧美人成| 怎么达到女性高潮| 久久这里只有精品中国| 久久久久国产一级毛片高清牌| 国产成人系列免费观看| 久久这里只有精品19| 天天躁夜夜躁狠狠躁躁| 99精品久久久久人妻精品| 日本黄色视频三级网站网址| 久久天躁狠狠躁夜夜2o2o| 老司机在亚洲福利影院| 最近最新中文字幕大全免费视频| 少妇熟女aⅴ在线视频| 少妇裸体淫交视频免费看高清 | 国产精品久久久久久亚洲av鲁大| 国产成人一区二区三区免费视频网站| 国产爱豆传媒在线观看 | 人妻夜夜爽99麻豆av| 91国产中文字幕| 国产精品亚洲一级av第二区| 亚洲五月天丁香| 成人国产综合亚洲| 欧美zozozo另类| 亚洲一区中文字幕在线| 久久久精品欧美日韩精品| 韩国av一区二区三区四区| 久久久精品欧美日韩精品| 国产亚洲精品综合一区在线观看 | 欧美一级a爱片免费观看看 | 国产又黄又爽又无遮挡在线| 国产精品,欧美在线| 黄色丝袜av网址大全| 久久99热这里只有精品18| 亚洲天堂国产精品一区在线| 精品国产美女av久久久久小说| 国产精品永久免费网站| 色噜噜av男人的天堂激情| 国产久久久一区二区三区| 99久久精品热视频| 黄色丝袜av网址大全| 日韩欧美一区二区三区在线观看| 国产人伦9x9x在线观看| 国产蜜桃级精品一区二区三区| 午夜成年电影在线免费观看| 99久久精品热视频| 久久国产精品人妻蜜桃| 国产精品免费一区二区三区在线| 国产精品一区二区免费欧美| 黄片小视频在线播放| 精品久久久久久久久久免费视频| 男女床上黄色一级片免费看| 淫妇啪啪啪对白视频| 久久香蕉国产精品| 精品久久久久久久人妻蜜臀av| 亚洲无线在线观看| 91在线观看av| 国内毛片毛片毛片毛片毛片| 日韩精品中文字幕看吧| a级毛片在线看网站| 亚洲熟女毛片儿| 午夜免费激情av| 12—13女人毛片做爰片一| 国产精品美女特级片免费视频播放器 | 国产精品亚洲av一区麻豆| 人人妻人人澡欧美一区二区| 性欧美人与动物交配| 免费一级毛片在线播放高清视频| 亚洲熟妇熟女久久| 午夜日韩欧美国产| av视频在线观看入口| 国产乱人伦免费视频| 麻豆av在线久日| 在线观看www视频免费| 国产成人aa在线观看| 午夜福利在线在线| 久久久国产成人精品二区| 一区福利在线观看| 黄色a级毛片大全视频| 老汉色av国产亚洲站长工具| 国产1区2区3区精品| 日韩欧美免费精品| 欧美激情久久久久久爽电影| 亚洲精华国产精华精| 我的老师免费观看完整版| 欧美日韩福利视频一区二区| 久久久国产欧美日韩av|