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

    基于姿態(tài)變化的航空頻率域電磁法儀器偏置實時校正方法研究

    2016-11-08 03:00:01李光渠曉東陳潔黃玲方廣有
    地球物理學(xué)報 2016年10期

    李光, 渠曉東, 陳潔, 黃玲, 方廣有

    1 中國科學(xué)院電子學(xué)研究所, 北京 100190 2 中國科學(xué)院電磁輻射與探測技術(shù)重點實驗室, 北京 100190 3 中國科學(xué)院大學(xué), 北京 100049

    ?

    基于姿態(tài)變化的航空頻率域電磁法儀器偏置實時校正方法研究

    李光1,2,3, 渠曉東1,2,3, 陳潔1,2, 黃玲1,2, 方廣有1,2

    1 中國科學(xué)院電子學(xué)研究所, 北京100190 2 中國科學(xué)院電磁輻射與探測技術(shù)重點實驗室, 北京100190 3 中國科學(xué)院大學(xué), 北京100049

    磁偶極子的航空頻率域電磁法儀器在飛行測量的過程中由于儀器偏置的存在,且儀器偏置會隨著外部氣壓、溫度等環(huán)境因素以及收發(fā)線圈晃動的影響而呈現(xiàn)非線性變化,使得觀測數(shù)據(jù)出現(xiàn)誤差,因此需要對儀器偏置進(jìn)行校正.而傳統(tǒng)的在測線飛行前后將儀器抬至高空的“零場值”標(biāo)定方法具有成本高、受測區(qū)環(huán)境影響大以及采用線性插值獲取測線飛行過程中儀器偏置的精度低等缺點.本文根據(jù)儀器偏置與儀器姿態(tài)角變化無關(guān)的特性,通過測得儀器的姿態(tài)角信息,在滿足重疊偶極子模型的條件下,實現(xiàn)對儀器偏置的高精度實時校正.模型仿真結(jié)果表明,在30 m常規(guī)飛行高度下,該方法實時測得的儀器偏置精度接近于110 m高空測得的精度;校正后儀器偏置的絕對誤差與理論二次場的比值即相對誤差小于5%,滿足反演大地電導(dǎo)率的精度要求.該方法不僅減少了飛行的工作量,降低了飛行成本和飛行難度,而且可更加精確地獲得測線飛行過程中儀器偏置的非線性變化值,提高航測數(shù)據(jù)的觀測精度.

    航空頻率域電磁法; 實時; 偏置; 姿態(tài)角; 重疊偶極子模型

    1 引言

    航空頻率域電磁法以直升機(jī)、固定翼等作為飛行平臺,可實現(xiàn)快速、經(jīng)濟(jì)、大面積的探測,目前已被廣泛應(yīng)用于近地表的地質(zhì)環(huán)境普查、礦物勘探等領(lǐng)域.其中采用磁偶極子作為發(fā)射與接收裝置的頻率域電磁法的基本工作原理是:通過發(fā)射線圈產(chǎn)生單頻或多頻的一次磁場信號(即一次場);與大地或地下良導(dǎo)體作用產(chǎn)生感應(yīng)渦流,進(jìn)而產(chǎn)生相同頻率的二次磁場信號(即二次場);系統(tǒng)通過接收線圈接收該二次場信號,可達(dá)到探測大地電性參數(shù)分布的目的.航空頻率域電磁法儀器中較為常見的三種收發(fā)線圈架構(gòu)分別為水平共面(Horizontal Coplanar,即HCP)、垂直共面(Vertical Coplanar ,即VCP)、垂直共軸(Vertical Coaxial,即VCA).目前商業(yè)化的航空頻率域電磁法儀器有加拿大Fugro公司的DIGHEM系統(tǒng)(Fraser,1979,1986)、RESOLVE系統(tǒng)(Cain,2004)以及美國Geophex公司的GEM-2A系統(tǒng)(Won et al.,2003)等.

    依據(jù)航空頻率域電磁法儀器的工作原理,在沒有外部二次場輸入、收發(fā)線圈位置固定的情況下,由接收線圈接收到的一次場可以理論算出,進(jìn)而可以通過數(shù)據(jù)處理或者Bucking線圈(Won et al.,2003)消除掉接收線圈接收到的一次場值,實現(xiàn)系統(tǒng)輸出為零.但由于收發(fā)線圈的位置偏差、接收線圈特性以及后續(xù)調(diào)理電路等因素的影響使得系統(tǒng)的輸出不為零即為儀器偏置(Valleau,2000).一方面,該偏置的存在以及其所具有的隨氣壓、溫度、濕度等外部環(huán)境變化及收發(fā)線圈晃動影響而呈現(xiàn)的非線性變化特性,會導(dǎo)致相鄰測線數(shù)據(jù)之間的不連續(xù),出現(xiàn)“窗簾效應(yīng)”的條帶狀零漂誤差(Huang et al.,1999;Green,2003;李文杰,2007);無法實現(xiàn)對姿態(tài)誤差的高精度校正(Yin et al.,2004;Fitterman et al.,2004);且反演得到的層狀大地電導(dǎo)率、磁化率以及各層厚度等參數(shù)也會存在誤差,以至于無法實現(xiàn)對海水深度及冰層厚度等的高精度探測(Haas et al.,2009;Reid and Bishop,2004),因此需要對儀器偏置進(jìn)行校正.另一方面,與時間域電磁法不同,頻率域電磁法的一次場一直處于發(fā)射狀態(tài),且接收線圈接收到的二次場會受到外部環(huán)境如未知大地或良導(dǎo)體響應(yīng)的影響,使得儀器偏置的校正實現(xiàn)起來較為困難.

    對于地面頻率域電磁法儀器而言,目前較為常用的儀器偏置校正方法主要有以下三種:在同一測量點處利用不同高度的測量值聯(lián)合反演出儀器偏置(Sasaki et al.,2008;Minsley et al.,2014);將儀器抬至一定的高度直接測量得到儀器偏置(Mitsuhata and Imasato,2009);采用標(biāo)準(zhǔn)的電阻率法儀器如直流電阻率法儀器實現(xiàn)對儀器偏置的校正(Deszcz et al.,1998;Minsley et al.,2012).航空頻率域電磁法儀器偏置的校正可以采用相類似的校正方法,其中以將儀器抬至高空(如300 m)中直接測量儀器偏置的校正方法最為簡單、方便、也較為常用.因此航空頻率域電磁法常采用在測線飛行前后將儀器抬至高空中實現(xiàn)對儀器偏置的測量(Valleau,2000;Cain,2004),而飛行測量過程中的儀器偏置則可以通過線性插值方式獲得,然后再通過調(diào)平的方式進(jìn)一步消除“窗簾效應(yīng)”的影響(Huang and Fraser,1999;Green,2003;李文杰,2007).若儀器抬至的高度不夠高,則不能完全消除未知大地及良導(dǎo)體的影響,使得測得的儀器偏置精度不高.另外上述傳統(tǒng)的高空“零場值”標(biāo)定方法還具有飛行難度大,飛行成本高,受外部環(huán)境(如測區(qū)海拔,溫度,氣壓,濕度等)影響大,線性插值校正精度低等缺點.

    已知航空頻率域電磁法儀器在測線飛行的過程中,其姿態(tài)角發(fā)生變化是不可避免的.而儀器姿態(tài)角的變化,會對測量的二次場產(chǎn)生影響,需要對姿態(tài)角變化進(jìn)行測量從而實現(xiàn)對二次場的姿態(tài)校正(Yin et al.,2004;Fraser and Hodges,2007;嵇艷鞠等,2010;王琦等,2013;曲昕馨等,2014).另一方面,由磁通門傳感器偏置的校正算法(Bracken et al.,2005;Pang et al.,2013)可知,磁通門的偏置不受其姿態(tài)角變化的影響,而頻率域電磁法儀器與之類似.本文利用偏置與姿態(tài)角變化無關(guān)的特性以及三種常見線圈架構(gòu)測得的二次場之間滿足的等式關(guān)系,在航空頻率域電磁法儀器測線飛行的過程中,通過測量儀器姿態(tài)角的變化信息實現(xiàn)對儀器偏置的實時校正.

    2 層狀大地正演模型及重疊偶極子模型

    航空頻率域電磁法儀器的飛行平臺以直升機(jī)最為常見,將電磁法儀懸掛于直升機(jī)的下方,如圖1所示.收發(fā)線圈沿x軸向分布,發(fā)射線圈的坐標(biāo)為(0,0,-h),接收線圈的坐標(biāo)為(xr,yr,zr),其中h為頻率域電磁法儀器距離地面的高度.

    圖1 航空頻率域電磁法儀器探測層狀大地示意圖Fig.1 Schematicdiagram of airborne frequency domain electromagnetic instrument detecting the layered earth

    航空頻率域電磁法儀器的收發(fā)線圈分別取xyz軸向,可以得到九種線圈架構(gòu)(Tlbll et al.,2007).在層狀大地模型下,可相應(yīng)地正演出九種線圈架構(gòu)下二次場的張量矩陣HS(Ward and Hohmann,1988;Yin et al.,2004):

    (1)

    (2)

    (3)

    (4)

    對于航空頻率域電磁法儀器而言,其飛行高度h遠(yuǎn)大于收發(fā)距ρ,滿足重疊偶極子模型(WestandMacnae,1991;YinandFraser,2004;Fittermanetal., 2004;FraserandHodges,2007),此時HCP、VCP以及VCA三種常見線圈架構(gòu)下的二次場具有以下等式關(guān)系成立:

    (5)

    3 儀器偏置校正算法

    航空頻率域電磁法儀器在飛行探測的過程中其姿態(tài)角會發(fā)生變化,圖2所示為單個線圈姿態(tài)角的變化,線圈繞x軸旋轉(zhuǎn),產(chǎn)生搖擺角α;繞y軸旋轉(zhuǎn),產(chǎn)生俯仰角β;繞z軸旋轉(zhuǎn),產(chǎn)生偏航角γ.認(rèn)為線圈在其姿態(tài)角未發(fā)生變化時處于參考坐標(biāo)系(x,y,z)下,而發(fā)生姿態(tài)角旋轉(zhuǎn)后,則處于線圈坐標(biāo)系(x′,y′,z′)下.

    圖2 線圈姿態(tài)旋轉(zhuǎn)示意圖(a) 搖擺旋轉(zhuǎn); (b) 俯仰旋轉(zhuǎn); (c) 偏航旋轉(zhuǎn).Fig.2 The coil attitude rotation schematic(a) Swing rotation; (b) Pitch rotation; (c) Yaw rotation.

    在層狀大地模型下,偏航角γ的變化對于二次場的測量沒有影響,而搖擺角α,俯仰角β的變化則對二次場的測量影響較大,因此可以僅考慮α及β角變化的影響.由參考坐標(biāo)系旋轉(zhuǎn)為線圈坐標(biāo)系,需要相應(yīng)的姿態(tài)角旋轉(zhuǎn)矩陣R.首先給出搖擺角α的旋轉(zhuǎn)矩陣,即式(6),俯仰角β的旋轉(zhuǎn)矩陣,即式(7),進(jìn)而得到由搖擺角α,俯仰角β均發(fā)生變化時的旋轉(zhuǎn)矩陣R,即式(8),且滿足RT=R-1.

    (6)

    (7)

    (8)

    由于收發(fā)線圈之間為硬連接,具有相同的姿態(tài)角變化.考慮二次場張量矩陣HS,收發(fā)線圈姿態(tài)角旋轉(zhuǎn)矩陣R以及儀器偏置offset的影響,可得接收線圈接收到的二次場張量矩陣:

    (9)

    對于以z軸發(fā)射,z軸接收的HCP線圈架構(gòu)而言,其接收線圈測得的二次場為

    (10)

    將式(5)代入,可進(jìn)一步得到:

    (12)  已知三種常見線圈架構(gòu)下的二次場滿足成立,則有Δexx=-Δeyy成立,可進(jìn)一步得到HCP線圈架構(gòu)在其姿態(tài)角發(fā)生變化前后任意兩次測得的二次場值及:

    (13)

    (14)

    由式(11)可知,忽略式(13),(14)中的誤差項并聯(lián)立,可求得儀器偏置offsetm:

    (15)

    由實測儀器偏置offsetm與理論儀器偏置offset做差,可得到儀器偏置的絕對誤差Δoffset:

    (16)

    由式(16)可知,當(dāng)α1=α2=α?xí)r,有Δoffset=2Δexxcos2α成立.而當(dāng)β1=β2=β時,有Δoffset=-2Δexx/cos2β成立.即在一個參數(shù)(如α角)固定,另一個參數(shù)(如β角)取不同的值時,求出的儀器偏置值僅與固定的參數(shù)有關(guān),而與變化的參數(shù)無關(guān).需要指出的是,在采用式(15)求取儀器偏置時,應(yīng)避免分母為零的情況發(fā)生,而推導(dǎo)得到的offsetm需要與HCP線圈架構(gòu)下的一次場做比值才能得到儀器偏置的ppm值(×10-6).

    4 模型仿真

    以Fugro公司的RESOLVE系統(tǒng)參數(shù)作為本文的模型仿真參數(shù),其采用的HCP線圈架構(gòu)收發(fā)距為7.9m,飛行高度為30m(Cain,2004).而RESOLVE系統(tǒng)中VCA線圈架構(gòu)儀器偏置的求取與之類似(僅受俯仰角β影響),這里不再敘述.文中首先采用均勻大地模型分析了其他參數(shù)對儀器偏置校正精度的影響,并于文章最后采用三種常見的大地模型驗證了本文儀器偏置校正方法的可行性.

    4.1三種線圈架構(gòu)下的二次場

    利用VCA、VCP以及HCP線圈架構(gòu)下求取的二次場實部、虛部分別與相應(yīng)一次場的比值,可以得到相應(yīng)線圈架構(gòu)下同相分量(In-phase)與正交分量(Quadrature)的ppm值(×10-6)(Won et al.,2003),如圖3所示.

    由圖3可知,對于儀器發(fā)射頻率與大地電導(dǎo)率的乘積(即頻率電導(dǎo)率乘積fσ)較小時測得的二次場或VCA線圈架構(gòu)測得的二次場而言,因其同相分量與正交分量值較小,儀器偏置的存在將嚴(yán)重限制大地電導(dǎo)率的反演精度.由于儀器的飛行高度h遠(yuǎn)大于其收發(fā)距ρ,滿足重疊偶極子模型,此時VCA線圈架構(gòu)得到的同相分量與正交分量為負(fù)值,且大小分別約為VCP線圈架構(gòu)下同相分量與正交分量的1/2,HCP線圈架構(gòu)下同相分量與正交分量的1/4(Frischknecht et al.,1991;Yin et al.,2004).

    為進(jìn)一步驗證式(5)的正確性,式(12)中Δexx,Δeyy的實部、虛部分別與ΔHxx,ΔHyy的實部、虛部做比值,可分別求得Δexx,Δeyy實部、虛部的相對誤差.在收發(fā)距固定不變的情況下,儀器高度在30~100 m之間變化,求得Δexx,Δeyy實部、虛部的相對誤差,如圖4所示.

    由圖4可知,Δexx,Δeyy實部、虛部的相對誤差受儀器高度以及頻率電導(dǎo)率乘積變化的影響,且均隨著高度的增大而減小,隨著頻率電導(dǎo)率乘積的增大而增大,幅度分別可以達(dá)到2.4%、4.1%.因此在滿足重疊偶極子模型及5%的精度要求下,三種常見線圈架構(gòu)下的二次場值滿足關(guān)系等式即公式(5)成立.

    4.2偏置校正的誤差分析

    由式(15)可知,單獨由α角或β角變化,或α角、β角同時變化均可以實現(xiàn)對儀器偏置誤差的測量.為進(jìn)一步說明α角、β角的不同取值,對求取儀器偏置的影響.分析α角在-1~1 rad之間任意變化,β角取0、0.1、0.2、0.3 rad四個不同角度值時,β角的不同取值對求取儀器偏置測量的影響;類似地,分析β角由-1~1 rad之間任意變化,α角取0、0.1、0.2、0.3 rad四個不同角度值時,α角的不同取值對求取儀器偏置的影響.儀器偏置的大小在0.1~1000 ppm之間變化,系統(tǒng)采用的發(fā)射頻率為10 kHz,而均勻大地模型的電阻率為50 Ωm即頻率電導(dǎo)率乘積為200 Hz S/m.本文方法校正后儀器偏置的絕對誤差,如圖5、圖6所示.

    圖3 三種常用線圈架構(gòu)下(a)同相分量與(b)正交分量隨頻率電導(dǎo)率乘積的變化Fig.3 The (a) In-phase and (b) Quadrature of three common coil configurations with the change of frequency-conductivity products

    由圖5、圖6可知,本文方法校正后儀器偏置的絕對誤差與偏置的大小無關(guān),而與儀器所處的姿態(tài)角α、β有關(guān).當(dāng)采用α角單獨變化求取儀器偏置時,β角越大,偏置實部與虛部的絕對誤差越大,β角取0.3 rad時,偏置實部與虛部的絕對誤差最大,分別為4.09 ppm,11.03 ppm;而當(dāng)采用β角單獨變化求取儀器偏置時,α角越大,偏置實部與虛部的絕對誤差就越小,α角取零時,儀器偏置實部與虛部的絕對誤差最大,分別為-3.74 ppm,-10.06 ppm.相較于α角單獨變化求取儀器偏置,β角單獨變化求取儀器偏置的絕對誤差較小,這與式(16)的推導(dǎo)結(jié)果相吻合.

    圖4 Δexx,Δeyy的相對誤差:(a)Δexx實部;(b)Δexx虛部;(c)Δeyy實部;(d)Δeyy虛部.Fig.4 The relative errors of Δexx,Δeyy in (a), (c) real parts and (b), (d) imaginary parts respectively

    圖5 α角變化,β角取不同的值時,測得不同偏置(a)實部與(b)虛部的絕對誤差Fig.5 The measured absolute errors of different bias (a) real parts and (b) imaginary parts with the change of α and specific value of β

    由圖7、圖8可知,本文方法校正后頻率域電磁法儀器偏置實部與虛部的相對誤差均隨著頻率電導(dǎo)率乘積的增大而增大并逐漸趨于平穩(wěn).當(dāng)采用α角單獨變化求取儀器偏置時,β角越大,偏置實部與虛部的相對誤差越大,分別可以達(dá)到2.84%,4.82%;而當(dāng)采用β角單獨變化求取儀器偏置時,α角越小,偏置實部與虛部的相對誤差也越大,分別可以達(dá)到-2.59%,-4.40%,均滿足反演大地電導(dǎo)率的精度要求.同樣地,相較于α角單獨變化求取儀器偏置,采用β角單獨變化求取儀器偏置的相對誤差較小.

    圖6 β角變化,α角取不同的值時,測得不同偏置(a)實部與(b)虛部的絕對誤差Fig.6 The measured absolute errors of different bias (a) real parts and (b) imaginary parts with the change of β and specific value of α

    圖7 α角單獨變化,β角取不同的值時,不同頻率電導(dǎo)率乘積下儀器偏置(a)實部與(b)虛部的相對誤差Fig.7 Relative errors of bias (a) real parts and (b) imaginary parts with the change of frequency-conductivity products, α and specific values of β

    圖8 β角單獨變化,α角取不同的值時,不同頻率電導(dǎo)率乘積下儀器偏置(a)實部與(b)虛部的相對誤差Fig.8 Relative errors of bias (a) real parts and (b) imaginary parts with the change of frequency-conductivity products, β and specific values of α

    4.3與高空偏置校正的比較

    圖9所示為頻率電導(dǎo)率乘積在20 Hz S/m~2 MHz S/m之間變化時,儀器在30 m常規(guī)飛行高度下,通過測量姿態(tài)角變化(β角取零值,α角單獨變化)求得的儀器偏置值、將儀器抬至110 m高空測得的儀器偏置值與理論偏置值100 ppm的比較.從整個頻率電導(dǎo)率乘積范圍來看,在30 m處本文方法測得的儀器偏置精度接近于110 m高空測得的儀器偏置精度.具體地,在頻率電導(dǎo)率乘積小于34.03 kHz S/m時,本文方法測得偏置實部的精度較高,在頻率電導(dǎo)率乘積大于34.03 kHz S/m時,采用110 m高空測量偏置實部的精度較高;而在頻率電導(dǎo)率乘積小于506 Hz S/m時,本文方法測得偏置虛部的精度較高,在頻率電導(dǎo)率乘積大于506 Hz S/m時,采用110 m高空測量偏置虛部的精度較高.由式(16)、圖4、圖5、圖6的分析及高空偏置校正方法的原理可知,兩種方法測得偏置的絕對誤差值均與理論偏置大小無關(guān),而與儀器所處的高度有關(guān),且偏置精度均隨著高度的增大而提高,經(jīng)仿真可知儀器在60 m處采用本文方法測得的儀器偏置精度接近于340 m高空測得的儀器偏置精度.

    由于高空“零場值”偏置校正方法無法在測線飛行過程中實施,進(jìn)而無法測得期間儀器偏置的非線性變化情況,若采用線性插值方式,偏置校正的精度會很低.與之不同,本文的方法可以實時測得儀器偏置的非線性變化,偏置校正具有較高的精度.在初始偏置實部與虛部均為100 ppm,頻率電導(dǎo)率乘積為200 Hz S/m的情況下,若在測線飛行的20 min內(nèi)儀器偏置受外部環(huán)境或收發(fā)線圈晃動影響而呈現(xiàn)幅度為10 ppm的非線性變化(以正弦變化為例),如圖10所示.可知在測線飛行的過程中不斷地將儀器抬至190 m高空進(jìn)行偏置校正,盡管可以反映出偏置的非線性變化,但在現(xiàn)實測量中不具備可操作性;在測線飛行前后采用190 m的高空標(biāo)定方法測得的儀器偏置值相同,其實部與虛部分別為113.1 ppm、103.4 ppm,若采用線性插值的方式,該偏置值恒定不變,無法反映儀器偏置在飛行測量過程中的變化情況.而通過姿態(tài)角變化測得的儀器偏置則能夠反映出這一變化,從而可以實現(xiàn)儀器偏置的高精度實時校正.

    圖9 理論偏置為100 ppm,30 m高度姿態(tài)變化與110 m高空測量偏置(a)實部與(b)虛部的對比Fig.9 Comparison of bias (a) real parts and (b) imaginary parts measured by attitude angle change at a height of 30 m and measured in a high altitude of 110 m, while theoretical bias is 100 ppm

    圖10 理論偏置正弦波變化時,30 m高度姿態(tài)變化與190 m高空測量的偏置(a)實部與(b)虛部對比Fig.10 Comparison of bias (a) real parts and (b) imaginary parts measured by attitude angle change at a height of 30 m and measured in a high altitude of 190 m, while theoretical bias varies in sine wave

    為進(jìn)一步驗證本文基于姿態(tài)變化求取儀器偏置的可行性.在發(fā)射頻率為10 kHz的情況下,以下考慮三種常見的大地模型,在俯仰角β取不同值時,搖擺角α單獨變化測得儀器偏置的相對誤差,如圖11所示.其中K型的大地模型采用厚度分別為50 m、20 m,而電阻率分別為5、500 Ωm以及50 Ωm的三層大地模型;H型大地模型采用厚度分別為50、20 m,而電阻率分別為500、5 Ωm以及50 Ωm的三層大地模型;均勻大地模型電阻率為50 Ωm.可知,校正后三種常見大地模型的偏置相對誤差均隨著俯仰角的增大而增大,且均小于1.88%,驗證了本文方法的可行性.

    圖11 在三種常見大地模型下,測得的儀器偏置(a)實部與(b)虛部隨俯仰角β變化的相對誤差Fig.11 The measured relative errors of bias (a) real parts and (b) imaginary parts with the change of pitch β in three common earth models

    5 結(jié)論

    基于磁偶極子的航空頻率域電磁法儀器在飛行測量的過程中由于儀器偏置的存在,且儀器偏置會隨著外部環(huán)境以及收發(fā)線圈晃動的影響而發(fā)生變化,使得觀測數(shù)據(jù)出現(xiàn)誤差,因此需要對儀器偏置進(jìn)行校正.本文基于儀器偏置與姿態(tài)角變化無關(guān)的特性,通過測量姿態(tài)角變化信息,在滿足重疊偶極子模型的條件下,實現(xiàn)航空頻率域電磁法儀器測線飛行過程中儀器偏置的實時高精度校正.本文分別求取了三種常用線圈架構(gòu)下的二次場,并驗證了在5%的精度要求條件下,滿足相應(yīng)的等式關(guān)系即公式(5)成立;再利用該等式關(guān)系,推導(dǎo)出公式(15),進(jìn)而可以利用測得的儀器姿態(tài)角變化信息實現(xiàn)儀器偏置的測量,其絕對誤差與儀器偏置的大小無關(guān),而與儀器所處的姿態(tài)角、頻率電導(dǎo)率乘積以及高度有關(guān).文章最后采用三種常見的大地模型進(jìn)一步驗證了本文方法的可行性.針對HCP線圈架構(gòu)下儀器偏置的求取,可知具有如下規(guī)律成立:

    (1) 由搖擺角α單獨變化求取儀器偏置時,儀器所處俯仰角β的值越大,求得儀器偏置的絕對誤差與相對誤差就越大;由俯仰角β單獨變化求取儀器偏置時,儀器所處搖擺角α的值越小,求得儀器偏置的絕對誤差與相對誤差則越大,且采用俯仰角β單獨變化求取儀器偏置的絕對誤差與相對誤差均較??;

    (2) 當(dāng)頻率電導(dǎo)率乘積為200 Hz S/m時,校正后儀器偏置實部與虛部的絕對誤差分別小于4.09 ppm,11.03 ppm;而頻率電導(dǎo)率乘積在20 Hz S/m~2 MHz S/m之間變化時,校正后儀器偏置實部與虛部的絕對誤差與所測量二次場的比值即相對誤差均小于5%,滿足反演大地電導(dǎo)率的精度要求.

    (3) 頻率電導(dǎo)率乘積在20 Hz S/m~2 MHz S/m間變化時,儀器在30 m常規(guī)飛行高度下,通過姿態(tài)角變化實時測得的儀器偏置精度接近于110 m高空中測得的偏置精度,且能夠反映出飛行測量過程中儀器偏置的非線性變化,實現(xiàn)更高精度的實時偏置校正.

    采用本文的儀器偏置校正方法及Yin等提出的姿態(tài)校正算法(Yin,2004;Fitterman and Yin,2004)可以對三種常用線圈架構(gòu)的航空頻率域電磁法儀器,在飛行探測的過程中儀器偏置及姿態(tài)角變化所引起二次場誤差的聯(lián)合校正,也可以進(jìn)一步提高姿態(tài)校正的精度.通過姿態(tài)角變化校正儀器偏置的算法簡單,不僅減少了飛行的工作量,降低了飛行成本和飛行難度,而且還可以更加精確地實時獲得飛行測量過程中儀器偏置的非線性變化值,提高航測數(shù)據(jù)的觀測精度.

    Bracken R E, Smith D V, Brown P J, et al. 2005. Calibrating a tensor magnetic gradiometer using spin data. Scientific Investigations Report 2005-5045. Virginia State: US Geological Survey. Cain M J. 2004. Resolve survey for U.S. geological survey east poplar oil fields, Montana. Fugro Airborne Surveys Report #04034. Mississauga, Ontario: Fugro Airborne Surveys.

    Deszcz-Pan M, Fitterman D V, Labson V F. 1998. Reduction of inversion errors in helicopter EM data using auxiliary information.ExplorationGeophysics, 29(1-2): 142-146. Fitterman D V, Yin C C. 2004. Effect of bird maneuver on frequency-domain helicopter EM response.Geophysics, 69(5): 1203-1215.

    Fraser D C. 1979. The multicoil II airborne electromagnetic system.Geophysics, 44(8): 1367-1394.

    Fraser D C. 1986. Dighem resistivity techniques in airborne electromagnetic mapping. Airborne resistivity mapping. ∥Palacky G J ed. Airborne Resistivity Mapping. Canada: Geological Survey of Canada, 86-22. Fraser D C, Hodges G. 2007. Induction-response functions for frequency-domain electromagnetic mapping system for airborne and ground configurations.Geophysics, 72(2): F35-F44.

    Frischknecht F C, Labson V F, Spies B R, et al. 1991. Profiling methods using small sources. ∥Nabighian M N ed. Electromagnetic Methods in Applied Geophysics. Tulsa: Society of Exploration Geophysicists, 105-270. Green A. 2003. Correcting drift errors in HEM data. ∥ASEG 16th Geophysical Conference and Exhibition. ASEG Extended Abstracts, 2:1.

    Haas C, Lobach J, Hendricks S, et al. 2009. Helicopter-borne measurements of sea ice thickness, using a small and lightweight, digital EM system.JournalofAppliedGeophysics, 67(3): 234-241.Huang H P, Fraser D C. 1999. Airborne resistivity data leveling.Geophysics, 64(2): 378-385.

    Ji Y J, Lin J, Guan S S, et al. 2010. Theoretical study of concentric loop coils attitude correction in helicopter-borne TEM.ChineseJ.Geophys. (in Chinese), 53(1): 171-176, doi: 10.3969/j.issn.0001-5733.2010.01.019.

    Li W J. 2007. An automatic 2-D data leveling technique for frequency domain airborne electromagnetic survey.JournalofChengduUniversityofTechnology(Science&TechnologyEdition) (in Chinese), 34(4): 447-451.

    Minsley B J, Smith B D, Hammack R, et al. 2012. Calibration and filtering strategies for frequency domain electromagnetic data.JournalofAppliedGeophysics, 80: 56-66.

    Minsley B J, Kass M A, Hodges G, et al. 2014. Multielevation calibration of frequency-domain electromagnetic data.Geophysics, 79(5): E201-E216.

    Mitsuhata Y, ImasatoT. 2009. On-site bias noise correction in multi-frequency slingram-type electromagnetic induction measurements.JournalofEnvironmental&EngineeringGeophysics, 14(4): 179-188.Pang H F, Li J, Chen D X, et al. 2013. Calibration of three-axis fluxgate magnetometers with nonlinear least square method.Measurement, 46(4): 1600-1606.

    Qu X X, Li T L, Wang F. 2014. A Study on the effect of bird's attitude on helicopter EM data and the correction method.ChineseJ.Geophys. (in Chinese), 57(4): 1310-1321, doi: 10.6038/cjg20140428.

    Reid J E, Bishop J. 2004. Post-processing calibration of frequency-domain electromagnetic data for sea-ice thickness measurements.ExplorationGeophysics, 35(4): 283-287.Sasaki Y, Son J S, Kim C, et al. 2008. Resistivity and offset error estimations for the small-loop electromagnetic method.Geophysics, 73(3): F91-F95.Tlbll R J, Christensen N B. 2007. Sensitivity functions of frequency-domain magnetic dipole-dipole systems.Geophysics, 72(2): F45-F56.

    Valleau N C. 2000. HEM data processing-a practical overview.ExplorationGeophysics, 31(4): 584-594.

    Wang Q, Lin J, Yu S B, et al. 2013. Study on influence and correction of coil attitude and bird swing for the fixed-wing time-domain electromagnetic system.ChineseJ.Geophys. (in Chinese), 56(11): 3741-3750, doi: 10.6038/cjg20131115.Ward S H, Hohmann G W. 1988. Electromagnetic theory for geophysical applications. ∥Nabighian M N ed. Electromagnetic Methods in Applied Geophysics. Tulsa: Society of Exploration Geophysicists, 131-311.

    West G F, Macnae J C. 1991. Physics of the electromagnetic induction exploration method. ∥Nabighian M N ed. Electromagnetic Methods in Applied Geophysics-Applications, Part A and Part B. Tulsa: Society of Exploration Geophysicists, 5-45.Won I J, Oren A, Funak F. 2003. GEM-2A: a programmable broadband helicopter-towed electromagnetic sensor.Geophysics, 68(6): 1888-1895.Yin C C, Fraser D C. 2004. Attitude corrections of helicopter EM data using a superposed dipole model.Geophysics, 69(2): 431-439.

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

    李文杰. 2007. 用于頻率域航空電磁數(shù)據(jù)的二維自動調(diào)平. 成都理工大學(xué)學(xué)報(自然科學(xué)版), 34(4): 447-451.

    嵇艷鞠, 林君, 關(guān)珊珊等. 2010. 直升機(jī)航空TEM中心回線線圈姿態(tài)校正的理論研究. 地球物理學(xué)報, 53(1): 171-176, doi: 10.3969/j.issn.0001-5733.2010.01.019.

    曲昕馨, 李桐林, 王飛. 2014. 直升機(jī)吊艙姿態(tài)變化對電磁場測量的影響規(guī)律及其校正方法研究. 地球物理學(xué)報, 57(4): 1310-1321, doi: 10.6038/cjg20140428.

    王琦, 林君, 于生寶等. 2013. 固定翼航空電磁系統(tǒng)的線圈姿態(tài)及吊艙擺動影響研究與校正. 地球物理學(xué)報, 56(11): 3741-3750, doi: 10.6038/cjg20131115.

    (本文編輯汪海英)

    Bias real-time correction of airborne frequency domain electromagnetic system based on attitude change

    LI Guang1,2,3, QU Xiao-Dong1,2,3, CHEN Jie1,2, HUANG Ling1,2, FANG Guang-You1,2

    1InstituteofElectronics,ChineseAcademyofSciences,Beijing100190,China2KeyLaboratoryofElectromagneticRadiationandSensingTechnology,ChineseAcademyofSciences,Beijing100190,China3UniversityofChineseAcademyofSciences,Beijing100049,China

    There exists bias for airborne frequency domain electromagnetic instruments which are based on magnetic dipole in the course of flight measurement, and the bias will present nonlinear changes with the change of external air pressure, temperature and other environmental factors and the sway of the transmitter and receiver coils, leading to errors in observation data, so it is necessary to calibrate the bias. The traditional method, “zero field” calibration method, which moves the instrument to a high altitude before and after the flight measurement has a high cost, can be greatly affected by the zone environment, and has low precision by using linear interpolation during the flight. Based on the feature of not varying with the change of attitude angles, the bias of airborne frequency domain electromagnetic instruments can get real-time correction with high-precision during flight measurement by measuring the attitude angles, while meeting the condition of superposed dipole model. The model simulation shows that at the routine flight height of 30m, the precision of bias measured by this method is close to the precision of bias measured at a height of 110 m; the relative error of bias absolute error over the secondary field is less than 5% after correction, meeting the accuracy requirements of ground conductivity inversion. This approach not only reduces the workload and cost of flight, but also can more accurately obtain the nonlinear changes during flight measurement, improving the accuracy of aerial observation data.

    Airborne frequency domain electromagnetic; Real-time; Bias; Attitude angle; Superposed dipole model

    10.6038/cjg20161035.

    國家重大科研裝備研制項目“深部資源探測核心裝備研發(fā)”子項目“航空瞬變電磁勘探儀”(ZDYZ2012-1-03-05)資助.

    李光,男,1988年生,博士研究生,主要從事頻率域電磁法儀器設(shè)計及理論研究. E-mail:liguang706@163.com

    10.6038/cjg20161035

    P631

    2015-09-18,2016-07-20收修定稿

    李光, 渠曉東, 陳潔等. 2016. 基于姿態(tài)變化的航空頻率域電磁法儀器偏置實時校正方法研究. 地球物理學(xué)報,59(10):3917-3926,

    Li G, Qu X D, Chen J, et al. 2016. Bias real-time correction of airborne frequency domain electromagnetic system based on attitude change.ChineseJ.Geophys. (in Chinese),59(10):3917-3926,doi:10.6038/cjg20161035.

    淫妇啪啪啪对白视频| 五月伊人婷婷丁香| 久久久精品大字幕| 久久国产精品人妻蜜桃| 嫩草影院精品99| 一个人看的www免费观看视频| 99热6这里只有精品| 欧美午夜高清在线| 不卡av一区二区三区| 日本一本二区三区精品| 观看美女的网站| 午夜a级毛片| 2021天堂中文幕一二区在线观| 精品国产超薄肉色丝袜足j| 日本熟妇午夜| 看黄色毛片网站| 九色国产91popny在线| 亚洲av电影不卡..在线观看| xxxwww97欧美| 美女cb高潮喷水在线观看 | 午夜福利成人在线免费观看| x7x7x7水蜜桃| 最近在线观看免费完整版| 99久久国产精品久久久| 日韩欧美一区二区三区在线观看| 国产精品久久久久久亚洲av鲁大| 欧美成狂野欧美在线观看| 久久精品国产亚洲av香蕉五月| 欧美日韩国产亚洲二区| 国产精品亚洲av一区麻豆| 亚洲熟妇中文字幕五十中出| 日本免费一区二区三区高清不卡| 国产伦一二天堂av在线观看| 免费观看精品视频网站| 欧美中文综合在线视频| 一进一出抽搐动态| 成年人黄色毛片网站| 桃红色精品国产亚洲av| 麻豆av在线久日| 久久久久九九精品影院| 亚洲精品久久国产高清桃花| 国产精品av久久久久免费| 国产一区二区激情短视频| 欧美午夜高清在线| 亚洲精品色激情综合| 男人舔奶头视频| 黑人巨大精品欧美一区二区mp4| 日韩精品中文字幕看吧| 久久人妻av系列| 老司机午夜福利在线观看视频| 十八禁人妻一区二区| 一个人免费在线观看的高清视频| 国产精品,欧美在线| 伊人久久大香线蕉亚洲五| 手机成人av网站| 欧美日韩国产亚洲二区| 色尼玛亚洲综合影院| www国产在线视频色| 99国产精品一区二区蜜桃av| 99热只有精品国产| 91字幕亚洲| 欧美日韩精品网址| 丁香六月欧美| 午夜成年电影在线免费观看| xxxwww97欧美| a级毛片在线看网站| 久久久国产精品麻豆| 日韩欧美一区二区三区在线观看| 最好的美女福利视频网| 日韩中文字幕欧美一区二区| 日韩精品中文字幕看吧| 日日摸夜夜添夜夜添小说| 精品99又大又爽又粗少妇毛片 | 男人舔女人下体高潮全视频| 免费在线观看成人毛片| 日日干狠狠操夜夜爽| 变态另类丝袜制服| 亚洲国产看品久久| 欧美+亚洲+日韩+国产| 三级男女做爰猛烈吃奶摸视频| 一级毛片高清免费大全| 特大巨黑吊av在线直播| 久9热在线精品视频| 亚洲国产精品成人综合色| 国产午夜精品久久久久久| 亚洲美女视频黄频| 亚洲av成人不卡在线观看播放网| 悠悠久久av| 又粗又爽又猛毛片免费看| 激情在线观看视频在线高清| 亚洲欧美日韩高清在线视频| 两性夫妻黄色片| 国产精品自产拍在线观看55亚洲| 国产又色又爽无遮挡免费看| 伊人久久大香线蕉亚洲五| 中文亚洲av片在线观看爽| 亚洲专区字幕在线| 99热这里只有是精品50| 一卡2卡三卡四卡精品乱码亚洲| 一级作爱视频免费观看| 亚洲五月天丁香| 一个人观看的视频www高清免费观看 | 窝窝影院91人妻| 欧美极品一区二区三区四区| 亚洲一区高清亚洲精品| 三级毛片av免费| 国产亚洲av嫩草精品影院| 一夜夜www| 久久久久久大精品| 亚洲av美国av| 亚洲18禁久久av| 一二三四社区在线视频社区8| 国产精品香港三级国产av潘金莲| 久久精品人妻少妇| 国产精品久久久久久亚洲av鲁大| 白带黄色成豆腐渣| 国产熟女xx| 国产av麻豆久久久久久久| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产99精品国产亚洲性色| 熟女人妻精品中文字幕| 青草久久国产| 午夜精品在线福利| 嫩草影院入口| 成人无遮挡网站| 成人午夜高清在线视频| 欧美激情在线99| 舔av片在线| 成人特级av手机在线观看| 丰满的人妻完整版| 精品国产三级普通话版| 人人妻,人人澡人人爽秒播| 一进一出抽搐gif免费好疼| 黄色 视频免费看| 国产一区二区三区视频了| 亚洲色图av天堂| 国产乱人伦免费视频| 免费看光身美女| 国产精品电影一区二区三区| 国产人伦9x9x在线观看| 精品久久久久久久久久久久久| 亚洲av电影不卡..在线观看| 色综合欧美亚洲国产小说| 美女高潮的动态| 久9热在线精品视频| 亚洲成人久久爱视频| 亚洲成av人片免费观看| www.www免费av| 欧美xxxx黑人xx丫x性爽| 99国产精品一区二区蜜桃av| 亚洲av第一区精品v没综合| 我要搜黄色片| 叶爱在线成人免费视频播放| 嫩草影院精品99| av在线天堂中文字幕| 舔av片在线| 久久久久久大精品| 成在线人永久免费视频| 国产精品日韩av在线免费观看| 国产精品一及| 亚洲美女视频黄频| 国产精品久久久久久精品电影| 中文字幕人成人乱码亚洲影| 人人妻人人澡欧美一区二区| av在线蜜桃| 国产一级毛片七仙女欲春2| 成人特级黄色片久久久久久久| 免费看a级黄色片| 亚洲av成人不卡在线观看播放网| 亚洲精品在线美女| 99久久久亚洲精品蜜臀av| 一区二区三区高清视频在线| 18禁观看日本| 伊人久久大香线蕉亚洲五| 人人妻人人澡欧美一区二区| 日本五十路高清| 国产野战对白在线观看| 999久久久国产精品视频| 国产精品综合久久久久久久免费| 老熟妇仑乱视频hdxx| 日韩免费av在线播放| 久久精品人妻少妇| av中文乱码字幕在线| 午夜福利欧美成人| 美女高潮喷水抽搐中文字幕| 丰满的人妻完整版| 男女做爰动态图高潮gif福利片| 少妇丰满av| 国产成年人精品一区二区| 亚洲熟妇中文字幕五十中出| 国产午夜精品久久久久久| 亚洲av成人av| av在线天堂中文字幕| 久久精品综合一区二区三区| 国产亚洲欧美在线一区二区| 亚洲性夜色夜夜综合| 久久久精品欧美日韩精品| 亚洲av成人av| 成在线人永久免费视频| а√天堂www在线а√下载| 国产黄色小视频在线观看| 亚洲精华国产精华精| 成人无遮挡网站| 亚洲人与动物交配视频| 99久久国产精品久久久| 欧美一级毛片孕妇| 国产人伦9x9x在线观看| xxx96com| 亚洲国产日韩欧美精品在线观看 | 亚洲人成网站高清观看| 日本三级黄在线观看| 制服人妻中文乱码| 国产在线精品亚洲第一网站| 国产激情偷乱视频一区二区| 精品久久久久久久人妻蜜臀av| 日日摸夜夜添夜夜添小说| 免费在线观看视频国产中文字幕亚洲| 亚洲国产精品合色在线| 国产精品一区二区精品视频观看| 国产午夜福利久久久久久| 成人国产一区最新在线观看| 免费看a级黄色片| 丁香欧美五月| av中文乱码字幕在线| 无遮挡黄片免费观看| 亚洲国产精品999在线| 亚洲美女视频黄频| 久久精品影院6| 色综合站精品国产| 夜夜爽天天搞| 国产精华一区二区三区| 两人在一起打扑克的视频| 99精品在免费线老司机午夜| 黄色日韩在线| 91av网一区二区| 亚洲国产欧洲综合997久久,| 国产成人影院久久av| 国内精品久久久久久久电影| 久久久久精品国产欧美久久久| 男女那种视频在线观看| 国产精品电影一区二区三区| 99久久精品一区二区三区| 久久久精品大字幕| 日本一本二区三区精品| 男女那种视频在线观看| 国产激情偷乱视频一区二区| 亚洲国产高清在线一区二区三| av国产免费在线观看| 人人妻人人澡欧美一区二区| 色综合亚洲欧美另类图片| 欧美日韩中文字幕国产精品一区二区三区| 91在线精品国自产拍蜜月 | 免费看美女性在线毛片视频| 一本久久中文字幕| 精品久久久久久成人av| 看片在线看免费视频| 日本黄色视频三级网站网址| 久久亚洲精品不卡| 嫩草影院精品99| 久99久视频精品免费| 午夜免费激情av| 欧美另类亚洲清纯唯美| 亚洲av片天天在线观看| 亚洲专区中文字幕在线| 久久香蕉国产精品| 在线观看日韩欧美| 国产av麻豆久久久久久久| 精品国产超薄肉色丝袜足j| 久久中文看片网| 又大又爽又粗| 国产高清videossex| 色老头精品视频在线观看| 午夜福利18| 久久亚洲精品不卡| 男女之事视频高清在线观看| 人妻久久中文字幕网| 亚洲人成网站在线播放欧美日韩| 精品乱码久久久久久99久播| 日韩欧美免费精品| 日韩欧美免费精品| 午夜免费激情av| 婷婷丁香在线五月| 亚洲国产欧美网| 人人妻人人澡欧美一区二区| 日日干狠狠操夜夜爽| 欧美极品一区二区三区四区| 亚洲国产日韩欧美精品在线观看 | 免费看a级黄色片| 国产精品日韩av在线免费观看| 成人无遮挡网站| 亚洲欧美激情综合另类| 91在线观看av| 人人妻人人澡欧美一区二区| 亚洲国产欧美人成| 国产美女午夜福利| 一级毛片高清免费大全| 亚洲av成人不卡在线观看播放网| www.999成人在线观看| 久久中文字幕一级| 亚洲精品乱码久久久v下载方式 | 在线观看66精品国产| 久久久久久国产a免费观看| 亚洲av电影不卡..在线观看| 色哟哟哟哟哟哟| 日日干狠狠操夜夜爽| 免费在线观看日本一区| 天天躁狠狠躁夜夜躁狠狠躁| 日日干狠狠操夜夜爽| 人妻丰满熟妇av一区二区三区| 好男人电影高清在线观看| 亚洲黑人精品在线| xxx96com| 999久久久国产精品视频| 午夜亚洲福利在线播放| 成人18禁在线播放| 久久久久久九九精品二区国产| 我要搜黄色片| 国产人伦9x9x在线观看| 久久精品人妻少妇| 国产精品久久久久久人妻精品电影| 欧美大码av| 国产av在哪里看| 一边摸一边抽搐一进一小说| 一级毛片高清免费大全| 最新中文字幕久久久久 | 一区二区三区国产精品乱码| 美女午夜性视频免费| 香蕉国产在线看| 不卡av一区二区三区| 亚洲精品美女久久久久99蜜臀| 久久久国产成人精品二区| 国产极品精品免费视频能看的| 偷拍熟女少妇极品色| 久久午夜综合久久蜜桃| 久久99热这里只有精品18| 熟女少妇亚洲综合色aaa.| 成人三级黄色视频| 欧美色欧美亚洲另类二区| 国产三级中文精品| 国产成人av教育| 999久久久国产精品视频| 亚洲国产色片| 国产高潮美女av| 两性夫妻黄色片| 一区二区三区激情视频| 日韩高清综合在线| 精品一区二区三区视频在线观看免费| 日本 欧美在线| 欧美日韩瑟瑟在线播放| а√天堂www在线а√下载| 特级一级黄色大片| 岛国在线免费视频观看| 超碰成人久久| 色播亚洲综合网| 欧美精品啪啪一区二区三区| 啪啪无遮挡十八禁网站| 在线a可以看的网站| 99久久无色码亚洲精品果冻| 又紧又爽又黄一区二区| 美女cb高潮喷水在线观看 | 精品一区二区三区四区五区乱码| 免费看a级黄色片| 97超级碰碰碰精品色视频在线观看| 欧美又色又爽又黄视频| 在线免费观看不下载黄p国产 | 国产激情久久老熟女| 男女下面进入的视频免费午夜| 欧美日韩亚洲国产一区二区在线观看| av欧美777| 色av中文字幕| 久久久成人免费电影| 国产成人啪精品午夜网站| 亚洲美女视频黄频| 观看免费一级毛片| 天堂√8在线中文| 99热精品在线国产| 一个人免费在线观看的高清视频| 一本综合久久免费| 69av精品久久久久久| 国产精品一区二区三区四区免费观看 | 久久香蕉国产精品| 在线国产一区二区在线| 在线观看美女被高潮喷水网站 | 99久久综合精品五月天人人| 一a级毛片在线观看| 欧美极品一区二区三区四区| 色综合亚洲欧美另类图片| 级片在线观看| 亚洲精品美女久久久久99蜜臀| 精品国内亚洲2022精品成人| 久久久久久国产a免费观看| 久久久久久久久中文| 最近视频中文字幕2019在线8| 亚洲无线观看免费| 亚洲欧美一区二区三区黑人| 亚洲精品一卡2卡三卡4卡5卡| 免费看十八禁软件| 麻豆久久精品国产亚洲av| 成人性生交大片免费视频hd| 午夜福利在线在线| 亚洲第一欧美日韩一区二区三区| 国产高清激情床上av| 国产av一区在线观看免费| 色吧在线观看| 午夜精品一区二区三区免费看| 精品国产美女av久久久久小说| av黄色大香蕉| 亚洲精品美女久久久久99蜜臀| 日日夜夜操网爽| 九九热线精品视视频播放| 中文资源天堂在线| 成人午夜高清在线视频| 亚洲人成网站在线播放欧美日韩| 久久久国产成人免费| 欧美+亚洲+日韩+国产| 两人在一起打扑克的视频| 美女大奶头视频| 欧美成人免费av一区二区三区| 热99在线观看视频| 精品一区二区三区av网在线观看| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩国产亚洲二区| 男女视频在线观看网站免费| av天堂中文字幕网| 久久久国产成人免费| 国产三级中文精品| 欧美一级毛片孕妇| 欧美午夜高清在线| 国产欧美日韩精品一区二区| 香蕉丝袜av| 亚洲最大成人中文| 十八禁人妻一区二区| 两个人的视频大全免费| 夜夜夜夜夜久久久久| 99re在线观看精品视频| 久久人妻av系列| 国产一区二区在线av高清观看| 黄色丝袜av网址大全| 中文字幕最新亚洲高清| 欧美在线黄色| 精品久久蜜臀av无| www日本黄色视频网| 亚洲人成电影免费在线| 身体一侧抽搐| 国产精华一区二区三区| 国产毛片a区久久久久| 成人午夜高清在线视频| 熟妇人妻久久中文字幕3abv| 一卡2卡三卡四卡精品乱码亚洲| 欧美国产日韩亚洲一区| 97人妻精品一区二区三区麻豆| 久99久视频精品免费| 国产一区二区三区视频了| 亚洲国产欧美一区二区综合| 国内精品一区二区在线观看| 欧美一区二区国产精品久久精品| xxxwww97欧美| 小说图片视频综合网站| 搡老妇女老女人老熟妇| 麻豆一二三区av精品| 村上凉子中文字幕在线| 亚洲人成电影免费在线| 老汉色∧v一级毛片| 禁无遮挡网站| 成人特级黄色片久久久久久久| 午夜精品久久久久久毛片777| 日本五十路高清| 精品久久久久久成人av| 亚洲精品美女久久久久99蜜臀| 9191精品国产免费久久| 99久久无色码亚洲精品果冻| 亚洲成人久久性| 美女高潮的动态| 久久99热这里只有精品18| 脱女人内裤的视频| 搡老妇女老女人老熟妇| 日本黄色视频三级网站网址| aaaaa片日本免费| av欧美777| 午夜精品一区二区三区免费看| 中文字幕熟女人妻在线| 国产精品99久久99久久久不卡| 久久九九热精品免费| 国产精品影院久久| 美女扒开内裤让男人捅视频| 亚洲18禁久久av| 欧美乱码精品一区二区三区| 欧美成人一区二区免费高清观看 | 久久久久性生活片| 久久精品aⅴ一区二区三区四区| 夜夜躁狠狠躁天天躁| 俄罗斯特黄特色一大片| 久久国产精品影院| 在线免费观看不下载黄p国产 | 怎么达到女性高潮| 国产黄色小视频在线观看| 久久香蕉精品热| 99热6这里只有精品| 日本 av在线| 啦啦啦韩国在线观看视频| 国产伦人伦偷精品视频| 精品一区二区三区视频在线观看免费| 日本黄色片子视频| 久久久成人免费电影| 天堂av国产一区二区熟女人妻| 亚洲av熟女| 亚洲一区二区三区色噜噜| 欧美3d第一页| 美女高潮喷水抽搐中文字幕| 成人永久免费在线观看视频| 高清毛片免费观看视频网站| 国产男靠女视频免费网站| 久久久久精品国产欧美久久久| 视频区欧美日本亚洲| 国产黄片美女视频| 午夜精品在线福利| 97人妻精品一区二区三区麻豆| 国产在线精品亚洲第一网站| 免费观看人在逋| 精品国产亚洲在线| 国产成人欧美在线观看| 中文字幕人妻丝袜一区二区| 老司机深夜福利视频在线观看| 青草久久国产| 国产精品久久久av美女十八| 久久久久久久久免费视频了| 一夜夜www| 999精品在线视频| 欧美一级毛片孕妇| 亚洲熟妇熟女久久| 老熟妇乱子伦视频在线观看| 18禁黄网站禁片午夜丰满| 久久亚洲精品不卡| 看黄色毛片网站| 18美女黄网站色大片免费观看| 日韩人妻高清精品专区| 日本在线视频免费播放| 亚洲av电影在线进入| 欧美日韩综合久久久久久 | 99久久久亚洲精品蜜臀av| 性欧美人与动物交配| 88av欧美| 欧美精品啪啪一区二区三区| 在线观看日韩欧美| 久久精品aⅴ一区二区三区四区| 99国产精品99久久久久| 欧美xxxx黑人xx丫x性爽| 午夜免费成人在线视频| 久久人妻av系列| 一个人免费在线观看电影 | 夜夜夜夜夜久久久久| 丰满人妻一区二区三区视频av | 中文亚洲av片在线观看爽| 最新在线观看一区二区三区| 丰满人妻一区二区三区视频av | 亚洲熟女毛片儿| av片东京热男人的天堂| 亚洲精品色激情综合| 亚洲va日本ⅴa欧美va伊人久久| 麻豆国产av国片精品| 午夜a级毛片| 色综合婷婷激情| 桃红色精品国产亚洲av| 一级a爱片免费观看的视频| 脱女人内裤的视频| 久久亚洲真实| 亚洲真实伦在线观看| 国产精品一区二区免费欧美| 国产高清有码在线观看视频| 91字幕亚洲| 中文字幕最新亚洲高清| 国产伦人伦偷精品视频| 国产精品综合久久久久久久免费| 日本在线视频免费播放| 99国产综合亚洲精品| 国产一区二区三区视频了| 欧美另类亚洲清纯唯美| or卡值多少钱| 精品一区二区三区视频在线观看免费| 色综合欧美亚洲国产小说| 欧美一区二区国产精品久久精品| 国产精品亚洲一级av第二区| 少妇的丰满在线观看| 麻豆国产97在线/欧美| 变态另类成人亚洲欧美熟女| 亚洲一区高清亚洲精品| 床上黄色一级片| 一级黄色大片毛片| 国产单亲对白刺激| 女生性感内裤真人,穿戴方法视频| 黄色日韩在线| 女人被狂操c到高潮| 国产午夜精品久久久久久| 麻豆成人av在线观看| 免费看日本二区| 色综合站精品国产| 日本黄色视频三级网站网址| 欧美日韩亚洲国产一区二区在线观看| 淫秽高清视频在线观看| 国产一级毛片七仙女欲春2| www.精华液| 国产精品女同一区二区软件 | 亚洲欧美精品综合久久99| 欧美日韩中文字幕国产精品一区二区三区| 亚洲人成网站在线播放欧美日韩| 国产v大片淫在线免费观看| 欧美性猛交黑人性爽| 午夜福利在线观看免费完整高清在 | 一夜夜www| 欧美+亚洲+日韩+国产| 久久久久久久久免费视频了| 亚洲黑人精品在线| 日本与韩国留学比较| 99精品在免费线老司机午夜| 丝袜人妻中文字幕|