亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        強(qiáng)磁性體ΔT異常計(jì)算的誤差分析研究

        2015-05-12 01:02:50袁曉雨姚長(zhǎng)利鄭元滿李澤林王君恒
        地球物理學(xué)報(bào) 2015年12期
        關(guān)鍵詞:磁化率磁化模量

        袁曉雨, 姚長(zhǎng)利, 鄭元滿, 李澤林, 王君恒

        地下信息探測(cè)技術(shù)與儀器教育部重點(diǎn)實(shí)驗(yàn)室和地質(zhì)過程與礦產(chǎn)資源國(guó)家重點(diǎn)實(shí)驗(yàn)室,中國(guó)地質(zhì)大學(xué)(北京), 北京 100083

        ?

        強(qiáng)磁性體ΔT異常計(jì)算的誤差分析研究

        袁曉雨, 姚長(zhǎng)利*, 鄭元滿, 李澤林, 王君恒

        地下信息探測(cè)技術(shù)與儀器教育部重點(diǎn)實(shí)驗(yàn)室和地質(zhì)過程與礦產(chǎn)資源國(guó)家重點(diǎn)實(shí)驗(yàn)室,中國(guó)地質(zhì)大學(xué)(北京), 北京 100083

        本文主要針對(duì)當(dāng)前磁法勘探中高精度處理解釋的需求,對(duì)強(qiáng)磁性體ΔT異常計(jì)算存在的誤差進(jìn)行分析研究.我們首先通過理論模型計(jì)算試驗(yàn),證明常規(guī)計(jì)算采用的投影關(guān)系的ΔT與實(shí)際測(cè)量的模量差ΔT之間的誤差E在磁異常幅值大時(shí)是明顯存在的,其影響不容忽視.其次,當(dāng)磁性強(qiáng)且剩磁存在時(shí),投影ΔT曲線及其誤差曲線在磁化方向與地磁場(chǎng)方向改變時(shí)具有一定的對(duì)稱性;地磁場(chǎng)T0、磁性體形態(tài)(如二度水平圓柱體模型的半徑r、柱體埋深R)和磁性參數(shù)(如磁化率κ)等參數(shù)確定的情況下,最大誤差值出現(xiàn)在磁性體正上方,且其大小與磁性參數(shù)(κ)和模型體規(guī)模(如r/R)之間皆是指數(shù)關(guān)系;另外,研究還發(fā)現(xiàn)ΔT的計(jì)算誤差曲線的一些其他規(guī)律特點(diǎn),如在各緯度帶上,ΔT計(jì)算誤差的最大值Emax曲線的極值主要分布在中緯度地區(qū);磁異常矢量Ta與地磁場(chǎng)T0的夾角θ逐漸變化時(shí),隨θ變化Emax曲線的極值分布在θ=90°~120°范圍內(nèi);當(dāng)磁異常幅值小于10000 nT時(shí),最大誤差近似為磁異常矢量垂直于地磁場(chǎng)方向的測(cè)點(diǎn)附近的誤差值;另外,磁性體(圓柱體為例)的半徑(即尺度)與埋深的比值r/R超過0.5,且磁化率超過0.1SI時(shí)誤差已達(dá)到3.9 nT,磁化率增大與對(duì)應(yīng)的Emax的值呈指數(shù)增長(zhǎng)特點(diǎn).因此,我們的研究表明,在強(qiáng)磁性體、磁異常幅值大的數(shù)據(jù)處理、反演及解釋時(shí),現(xiàn)有方法會(huì)產(chǎn)生較大的誤差,應(yīng)該基于嚴(yán)格的模量差ΔT,完善相應(yīng)的處理以及反演方法.

        總場(chǎng)異常ΔT; ΔT計(jì)算誤差; 磁異常模量Ta; 剩磁; 磁化率

        1 引言

        磁法是礦產(chǎn)資源勘探的重要手段,是尋找一些金屬礦類的主要方法之一.當(dāng)今的高精度磁力儀如質(zhì)子磁力儀、光泵磁力儀等,其靈敏度已分別達(dá)到0.01 nT和0.001 nT,因此,相應(yīng)的高精度處理轉(zhuǎn)換以及反演解釋逐步受到重視,例如,強(qiáng)磁性體的剩磁影響和退磁效應(yīng)近年來就受到普遍關(guān)注,相應(yīng)的研究得到不斷深化(Guo et al., 2001; Krahenbuhl and Li, 2007; Li et al., 2010),從而提高了處理解釋的質(zhì)量.

        我們關(guān)注的是一個(gè)更加基礎(chǔ)的問題,即ΔT異常近似計(jì)算的誤差問題.因?yàn)樵趯?shí)際的情況中,當(dāng)磁異常較弱時(shí),這個(gè)誤差通常非常小,所以在定性分析為主時(shí)這個(gè)問題一直被忽視,但當(dāng)磁性體磁性較強(qiáng)時(shí),所產(chǎn)生的磁異常幅值比較大,對(duì)后續(xù)定量計(jì)算和分析推斷可能會(huì)造成比較大的影響.所以,在當(dāng)今高精度勘探及處理解釋越來越受到重視的新階段,有必要對(duì)其影響進(jìn)行重新評(píng)估和深入分析,以期設(shè)法避免不必要的誤差,形成相應(yīng)的更高精度的處理解釋方法.

        現(xiàn)在勘探所用的主力儀器——質(zhì)子磁力儀(及光泵磁力儀)等直接測(cè)量地磁場(chǎng)總強(qiáng)度矢量T的幅值大小,將其與地磁場(chǎng)正常場(chǎng)T0的幅值相減以消除背景的影響,所獲得的模量之差稱為總場(chǎng)磁異常ΔT,基于研究的需要,這里我們稱其為模量差ΔT,簡(jiǎn)稱ΔT模量差.顯然,它是由地下磁性體產(chǎn)生的磁異常矢量Ta引起的,與其大小和方向變化密切相關(guān).航空磁測(cè)中Ta一般在幾個(gè)nT到幾千nT之間,一般不超過5000 nT(Blakely,1996;管志寧,2005),這時(shí)T0遠(yuǎn)大于Ta,ΔT就近似相當(dāng)于磁異常矢量Ta在地磁場(chǎng)T0方向上的投影分量.這樣,ΔT就和垂直分量異常Za以及水平分量異常Hax、Hay具有一樣的物理意義,都相當(dāng)于是磁異常矢量Ta的一個(gè)投影分量,我們稱其為投影分量ΔT,簡(jiǎn)稱ΔT投影.其重要意義是:由于其相當(dāng)于分量了,則磁異常ΔT投影與磁性體磁化強(qiáng)度變成簡(jiǎn)單的線性關(guān)系,與地磁場(chǎng)T0無直接關(guān)系,后續(xù)的分量互換等計(jì)算將變得可行,并且正演及反演計(jì)算也將得到大大簡(jiǎn)化.而采用嚴(yán)格的ΔT模量差磁異常形式,上述優(yōu)勢(shì)則不存在,即地磁場(chǎng)T0在計(jì)算表達(dá)式中顯式出現(xiàn),現(xiàn)有的分量互換方法將變得不可行.而實(shí)測(cè)磁異常由于其復(fù)雜性,其分析解釋離不開必要的定量計(jì)算分析,故處理轉(zhuǎn)換計(jì)算以及正演和反演計(jì)算逐步成為很重要的中間環(huán)節(jié).為此,當(dāng)強(qiáng)磁性體產(chǎn)生了大幅值磁異常時(shí),準(zhǔn)確的ΔT模量差與看作投影分量的ΔT投影之間的誤差到底有多大,對(duì)后續(xù)的處理以及反演會(huì)帶來多大的誤差和影響,正是本文關(guān)注和研究的內(nèi)容.

        2 總場(chǎng)異常ΔT計(jì)算誤差原理

        圖1a是ΔT及其他相關(guān)物理量關(guān)系的示意圖.圖中O點(diǎn)代表測(cè)量點(diǎn)(例如圖1b中的任一測(cè)點(diǎn)D點(diǎn)),該測(cè)點(diǎn)上地磁場(chǎng)正常場(chǎng)矢量為T0,地下的磁性體在該測(cè)點(diǎn)產(chǎn)生的磁異常矢量為Ta,測(cè)點(diǎn)O點(diǎn)處總磁場(chǎng)為這兩者的矢量合成量T,質(zhì)子磁力儀等磁場(chǎng)模量測(cè)量?jī)x器所測(cè)量的是T的大小即模量.顯然,如果沒有磁性體產(chǎn)生的磁異常Ta的影響,磁力儀測(cè)量的將只是背景地磁場(chǎng)T0的大小.所以,通常將實(shí)測(cè)磁場(chǎng)值減去背景地磁場(chǎng),獲得的是由磁性體引起的剩余磁場(chǎng),稱為磁異常.磁力儀測(cè)得的地磁場(chǎng)模量T減去背景地磁場(chǎng)T0,所得磁異常稱為總場(chǎng)異常ΔT.圖1b為示意顯示的地下一磁性體在斜磁化時(shí)產(chǎn)生的感應(yīng)磁場(chǎng)即磁異常Ta的空間分布情況,圖中的磁力線的分布直觀地表現(xiàn)了Ta的空間形態(tài)和方向變化,以及在地表其磁異常的矢量表達(dá).圖中A、B、C、D、E等為地面不同測(cè)點(diǎn),M為磁性體的磁化強(qiáng)度矢量,在地面不同測(cè)點(diǎn)處產(chǎn)生不同的磁異常矢量(紅色矢量線),與地磁場(chǎng)矢量T0合成,將作用于磁測(cè)儀器,產(chǎn)生觀測(cè)磁場(chǎng).

        圖1a所示的任意觀測(cè)點(diǎn)O處,地磁場(chǎng)總強(qiáng)度T、正常場(chǎng)T0,以及磁異常矢量Ta三者之間具有以下關(guān)系式:

        T=T0+Ta,

        (1)

        ΔT是T與T0的模量差,我們記作ΔT模量差,即

        ΔT模量差=|T|-|T0|=T-T0.

        (2)

        顯然,這里的ΔT模量差主要取決于Ta,但ΔT模量差并不等于模量Ta.基于(2)式及圖1a所示關(guān)系,可知

        ΔT模量差≤Ta.

        (3)

        另外前人的分析已指出,航空磁測(cè)中Ta一般不超過5000 nT(Blackely,1996;管志寧,2005),這時(shí),ΔT就近似為Ta在T0方向的投影,記作ΔT投影,即

        ΔT模量差≈Tacosθ=Tacos(Ta,T0)=ΔT投影.

        (4)

        結(jié)合圖1a及表達(dá)式(2)和(4)可以發(fā)現(xiàn),ΔT模量差和ΔT投影都主要由Ta控制,但又都與T0密切有關(guān).其中,ΔT模量差受T0的大小和方向雙重影響,而ΔT投影只受T0的方向影響(在該方向上的投影),與其大小無關(guān),這是兩者的關(guān)鍵區(qū)別.

        磁性體產(chǎn)生的磁異常Ta可以分解成三個(gè)分量Za、Hax和Hay,在正演計(jì)算中,通常也是通過先計(jì)算這三個(gè)分量,然后再將其合成獲得磁異常Ta的.現(xiàn)在,ΔT也是Ta的一個(gè)投影分量,將其看成是投影的一個(gè)重要優(yōu)點(diǎn)是,ΔT與Za等分量一樣,與磁性體磁性(如磁化強(qiáng)度)之間就是簡(jiǎn)單的線性關(guān)系,而與地磁場(chǎng)大小T0無關(guān),后續(xù)的磁異常處理轉(zhuǎn)換及正反演計(jì)算都能得到大大簡(jiǎn)化.例如,磁異常三分量之間可以互換,ΔT作為分量,也將同樣具有與三分量互換計(jì)算的理論基礎(chǔ);正演和反演中,如果只是磁性發(fā)生變化,也將是個(gè)線性問題.正因?yàn)槿绱耍L(zhǎng)期以來對(duì)于測(cè)量所得的ΔT模量差都直接將其當(dāng)作投影分量ΔT投影.另外,至今所有文獻(xiàn)中的處理轉(zhuǎn)換和正反演計(jì)算中,ΔT都是當(dāng)成ΔT投影對(duì)待的.

        圖1 地磁場(chǎng)及磁異常關(guān)系圖

        顯然,基于圖1所示原理圖不難發(fā)現(xiàn),由于地磁場(chǎng)T0的范圍約為30000~60000 nT,如果Ta很大,即如實(shí)際的磁鐵礦類等勘探中會(huì)遇到的數(shù)千納特(例如:俄羅斯的庫(kù)爾斯克鐵礦磁異常超過100000 nT,已遠(yuǎn)大于地磁場(chǎng)大?。┑拇女惓?,那么上述簡(jiǎn)化將引起明顯的誤差.這個(gè)誤差還和哪些因素直接相關(guān)?具有哪些特點(diǎn)?下面我們將通過必要的分析和模型實(shí)驗(yàn),加以研究評(píng)價(jià).

        由圖1所示關(guān)系容易得到,實(shí)測(cè)值ΔT模量差與ΔT投影即計(jì)算值之間存在的誤差E為

        E=ΔT模量差-ΔT投影,

        (5)

        根據(jù)矢量三角形余弦定理,可得

        (6)

        (7)

        結(jié)合前面(3)式可知

        ΔT模量差≥ΔT投影.

        (8)

        另外,分析誤差關(guān)系式(7)不難發(fā)現(xiàn),磁性體產(chǎn)生的磁異常幅值Ta大(甚至與背景地磁場(chǎng)T0在數(shù)量級(jí)上接近),誤差就會(huì)相當(dāng)大.外在因素對(duì)ΔT模量差的影響是比較復(fù)雜的,因?yàn)榇判泽w分布有限,其Ta的大小和方向在測(cè)線剖面上變化明顯(見圖1b),所以T(=Ta+T0)隨其相應(yīng)變化,這樣,ΔT模量差的影響會(huì)使誤差E隨測(cè)點(diǎn)變化而變化.當(dāng)Ta的值一定,ΔT模量差越?。ù藭r(shí)T與T0的值接近),誤差就會(huì)變大.

        為方便后面的定量對(duì)比分析,我們將ΔT模量差看作為實(shí)測(cè)數(shù)據(jù),ΔT投影看作計(jì)算數(shù)據(jù)(至今都是這么對(duì)待的),并且將計(jì)算數(shù)據(jù)相對(duì)于實(shí)測(cè)數(shù)據(jù)的均方誤差值與實(shí)測(cè)數(shù)據(jù)本身的標(biāo)準(zhǔn)差的比值定義為計(jì)算數(shù)據(jù)相對(duì)于實(shí)測(cè)數(shù)據(jù)的相對(duì)誤差,這也是ΔT投影相對(duì)于ΔT模量差的相對(duì)誤差.

        3 強(qiáng)磁性體總場(chǎng)異常ΔT投影計(jì)算誤差定量分析

        關(guān)于巖石磁性的研究,前人進(jìn)行過系統(tǒng)的巖石物理分析(Clark,1983;ClarkandEmerson,1991).其中,富含鐵的帶狀建造及含磁鐵礦的矽卡巖,其磁化率常見值可達(dá)2.0SI,最大可達(dá)4.0SI(ClarkandEmerson, 1991),而塊狀粗粒的(鈦)磁鐵礦和磁赤鐵礦的磁化率還要更高,其常見值在1.0SI到6.0SI之間, 此特征在國(guó)內(nèi)的強(qiáng)磁性鐵礦區(qū)都有體現(xiàn)(方華竹,1978; 朱永剛等,2006;張恒磊等,2011;朱朝吉等,2011;ZhangandHu,2012);此外,變質(zhì)巖中的條帶狀含鐵構(gòu)造的磁化率一般在0.01SI到3.0SI之間.我們的模型試驗(yàn)中磁性上限可據(jù)此確定.為簡(jiǎn)化問題,這里我們采用最簡(jiǎn)單的均勻磁化的二度水平圓柱體模型,分析ΔT的計(jì)算誤差的變化規(guī)律.

        由磁性體正演理論可知,磁性體產(chǎn)生的磁異常,其幅值主要受磁性體的規(guī)模、埋深,以及其磁化強(qiáng)度大小決定,而磁異常的變化特征(如剖面和平面分布特征)不僅與磁性體的幾何形態(tài)有關(guān),受磁化方向影響更大.為簡(jiǎn)單起見,我們首先選擇了最簡(jiǎn)單的垂直磁化條件進(jìn)行計(jì)算實(shí)驗(yàn).取觀測(cè)剖面長(zhǎng)度為200 m,測(cè)點(diǎn)間隔為1 m,圓柱體半徑為30 m,其中心埋深為40 m.地磁場(chǎng)強(qiáng)度設(shè)為50000 nT,磁傾角即為90°.這里地磁場(chǎng)為垂直向下,表明總場(chǎng)異常ΔT投影相當(dāng)于是垂直分量Za.為使ΔT計(jì)算誤差對(duì)比效果更加顯著,我們?cè)O(shè)模型的磁化率為3.0 SI,即圓柱體模型具有強(qiáng)磁性.通常強(qiáng)磁性體不僅具有較強(qiáng)的感應(yīng)磁化強(qiáng)度(感磁),而且具有較強(qiáng)的剩余磁化強(qiáng)度(剩磁),由于磁化方向?yàn)楦写排c剩磁的合成矢量的方向,如果與地磁場(chǎng)方向一致,即可看成無剩磁或剩磁方向與感磁方向一致;反之,如果磁化方向與地磁場(chǎng)方向不一致,則表明有剩磁的影響,且剩磁方向肯定與地磁場(chǎng)方向不同,從而改變了合成的磁化強(qiáng)度方向.但為簡(jiǎn)化分析過程,這里的磁化率可看成是產(chǎn)生總的磁化強(qiáng)度的等效磁化率,這只是便于分析而假設(shè)的,不影響結(jié)果和分析.該圓柱體模型的ΔT計(jì)算結(jié)果見圖2a所示.

        計(jì)算中所采取的方法是先計(jì)算二度圓柱體在剖面上各測(cè)點(diǎn)處的磁異常垂直分量Za和水平分量Ha,其合成得到總磁異常矢量Ta和其大小,即模量Ta,再由Ta與T0合成即得總的磁場(chǎng)矢量T,進(jìn)一步基于(2)式和(4)式即可得到ΔT的兩種計(jì)算結(jié)果:

        準(zhǔn)確的ΔT模量差和近似的ΔT投影.由圖2a可以看出:水平圓柱體的磁異常模量Ta對(duì)稱,極大值位置對(duì)應(yīng)柱體的中心點(diǎn),異常幅值已達(dá)42000 nT,接近地磁場(chǎng)大小了,屬于強(qiáng)磁異常.其次,ΔT模量差和ΔT投影在形態(tài)上一致,但顯然在數(shù)值上已出現(xiàn)明顯的誤差,從誤差值E曲線可以看出,在垂直磁化條件下,誤差最大值不在異常最大值處(異常最大值處誤差為0),而是在異常極值兩側(cè),誤差值已達(dá)5000 nT,達(dá)到異常幅值的18.5%,顯然這已經(jīng)完全不能忽略了,會(huì)直接影響后續(xù)的定量評(píng)價(jià).實(shí)際上基于理論公式可知,對(duì)于地磁場(chǎng)方向和磁化方向都垂直向下條件下的水平圓柱體,其磁異常ΔT投影與垂直分量Za完全一樣,其表達(dá)式由任意磁化的(9)式(Blakely,1996;管志寧,2005;李才明和李軍,2013)

        ΔT投影=

        (9)轉(zhuǎn)化為

        (10)

        而ΔT模量差的值是地磁場(chǎng)總強(qiáng)度T與地磁場(chǎng)T0的模量差值,由(2)式和(9)式可得具體表達(dá)式如下:

        -T0,

        (11)其中,R為圓柱體的中心埋深,μ0為真空中的磁導(dǎo)率,ms為單位長(zhǎng)度的磁矩,I0和I分別為地磁場(chǎng)傾角和磁化強(qiáng)度傾角.顯然,由式(10)可得到ΔT投影的兩個(gè)0值點(diǎn)的距離正好是圓柱體埋深的2倍,與由圖2a的ΔT投影曲線得到0值點(diǎn)的距離為80m完全一致,這就是簡(jiǎn)單反演方法中的特征點(diǎn)法.但是,從實(shí)測(cè)角度來看,對(duì)該模型實(shí)測(cè)的將是嚴(yán)格的ΔT模量差,而不是近似的ΔT投影,如果仍然取圖2a中ΔT模量差曲線的0值點(diǎn)的距離(該距離值約為98m),應(yīng)用式(10)推測(cè)的圓柱體中心埋深則為49m,與實(shí)際深度40m相比,誤差將達(dá)23%.顯然,此時(shí)再將ΔT模量差當(dāng)作近似的ΔT投影對(duì)待,在定量計(jì)算時(shí)會(huì)造成明顯的誤差,這是強(qiáng)磁異常定量計(jì)算及分析時(shí)需要注意的方面.

        我們?cè)傺a(bǔ)充一個(gè)反向磁化實(shí)驗(yàn),具體為保持地磁場(chǎng)不變,磁化強(qiáng)度方向改為反向,即磁化方向向上,與之對(duì)應(yīng)的ΔT投影和ΔT模量差曲線對(duì)比如圖2b所示.可以看出,對(duì)稱反向磁化后,兩次ΔT投影曲線變?yōu)檎?fù)鏡像關(guān)系,是嚴(yán)格對(duì)稱的,對(duì)稱軸為0值線,但ΔT模量差曲線則不存在嚴(yán)格對(duì)稱性這樣的特征,表現(xiàn)為只是正負(fù)鏡像,形態(tài)不嚴(yán)格對(duì)稱,鏡像軸抬高離開0值線.還可以看出,ΔT模量差在負(fù)值區(qū)比ΔT投影曲線“瘦”,在正值區(qū)比ΔT投影曲線“胖”,這是因?yàn)槭冀K存在ΔT模量差≥ΔT投影,即E總是大于等于0的緣故.

        由于ΔT模量差≥ΔT投影,如果我們?cè)趯?shí)際情況中將ΔT投影當(dāng)做ΔT模量差進(jìn)行后續(xù)的定量計(jì)算,如反演,其反演結(jié)果所對(duì)應(yīng)的場(chǎng)值和規(guī)模會(huì)偏大,這是不容忽視的,后面的模型計(jì)算實(shí)驗(yàn)也證明了這點(diǎn).

        上述分析的例子中,如果磁化方向與地磁場(chǎng)方向不一致,則表明肯定有剩磁的影響.保持上述其他條件及參數(shù)不變,只改變磁化方向及地磁場(chǎng)方向,對(duì)圓柱體模型算例進(jìn)行進(jìn)一步分析,則該水平圓柱體的磁異常各參量沿剖面將發(fā)生變化.

        圖3是磁化方向和地磁場(chǎng)方向都在變化的幾個(gè)算例,概括地展示了磁異常的變化特征.由圖3可見,雖然磁異常矢量Ta與磁化方向密切有關(guān),但其大小即模量Ta不受磁化方向的影響,這是二度體磁異常的特征,理論上已證明這點(diǎn)(Kontis and Young,1964).但對(duì)三度體而言,則不具有該特征,三度體的磁異常模量受磁化方向影響,只是與磁異常分量(ΔT投影、Za、Hax和Hay),相比較小,前人對(duì)此已作大量分析,這里也不再贅述.

        從式(9) 和式(10)可以看出,ΔT投影和ΔT模量差與磁化傾角I與地磁場(chǎng)傾角I0的和有關(guān),為簡(jiǎn)化起見,我們?cè)O(shè)S=I+I0.結(jié)合式(9)、(10)可以得出,當(dāng)S為定值時(shí),二度水平圓柱體模型在剖面上的ΔT投影曲線、ΔT模量差曲線以及誤差E的曲線都完全一致,如圖3(a、e、i)、圖3(b、f)或圖3(d、h);以上三條曲線在S和360°-S情況下是關(guān)于y軸對(duì)稱的;由于ΔT模量差≥ΔT模影,所以S和180°-S情況下ΔT投影曲線關(guān)于原點(diǎn)對(duì)稱,但是ΔT模量差曲線卻不具備這個(gè)對(duì)稱性的特點(diǎn),相應(yīng)的誤差E的曲線也不再具有對(duì)稱性;同理,S和180°+S情況下ΔT投影曲線關(guān)于原點(diǎn)對(duì)稱,但是ΔT模量差與其相對(duì)應(yīng)的誤差E曲線不再具有對(duì)稱性.

        S=0°時(shí),可有多種組合方式,如:磁化傾角I=0°,地磁場(chǎng)傾角I0=0°,即均為水平方向,如圖3c所示;或如:I=-90°,I0=90°,結(jié)果如前面圖2b中的磁化方向向上的狀態(tài),誤差E曲線關(guān)于y軸對(duì)稱分布,其極大值為11000 nT,剖面相對(duì)誤差為38.6%,誤差明顯.

        當(dāng)S=180°時(shí),也有多種組合方式,如:I=90°,I0=90°,結(jié)果如圖3g所示;或I=0°,I0=180°,誤差E曲線關(guān)于y軸對(duì)稱分布,其誤差的極大值在所有情況里面是最小的,為5673.3 nT,剖面相對(duì)誤差為18.8%.

        圖2 (a) 水平圓柱體磁異常各參量剖面(地磁場(chǎng)及磁化方向均垂直向下); (b) 地磁場(chǎng)T0垂直向下且磁化方向分別為向下和向上時(shí)的水平圓柱體模型的ΔT的模量差曲線與投影曲線對(duì)比

        圖3 二度水平圓柱體的ΔT異常及其誤差曲線

        當(dāng)S=90°以及S=270°時(shí),如圖3(a、e、i)所示,ΔT投影曲線本身關(guān)于原點(diǎn)對(duì)稱,此時(shí)誤差E是一條單峰曲線,最大值達(dá)到16820 nT,是此處Ta(41770 nT)的40.3%,剖面相對(duì)誤差為33.93%,影響顯然較大.另外注意,此時(shí)誤差最大值點(diǎn)位置并不在模型體中心點(diǎn)正上方,而是稍有偏移.

        通過上面的分析可知,S(=I+I0)一定時(shí),圓柱體對(duì)應(yīng)的ΔT投影、ΔT模量差以及誤差E曲線都是相同的,所以我們這里來探討在當(dāng)前磁參數(shù)下S在0°~360°范圍時(shí)剖面上誤差E的最大值隨S的變化規(guī)律,結(jié)果如圖4所示.可以看出,曲線關(guān)于S=180°對(duì)稱,在S=180°處出現(xiàn)一個(gè)極小值,且曲線的極大值出現(xiàn)在S=65°附近,為17798 nT,是此時(shí)異常幅值Ta(42187 nT)的42.2%.該最大誤差值位于模型體中心點(diǎn)正上方位置,此時(shí)整個(gè)剖面的相對(duì)誤差為37.3%.

        圖4 ΔT誤差最大值Emax隨S在0°~360°范圍內(nèi)的變化

        (12)

        此時(shí)ΔT模量差=0.

        從式(12)可以看出,誤差最大值與T0之間是線性關(guān)系,與κ和r/R之間皆是指數(shù)關(guān)系,可見磁化率κ對(duì)誤差值的影響呈指數(shù)遞增,一旦磁性體的磁化率較大,它造成的誤差影響我們更加不能忽視.

        4 ΔT計(jì)算誤差的幾個(gè)特點(diǎn)規(guī)律分析

        4.1 地表隨緯度變化ΔT計(jì)算誤差的特點(diǎn)

        在地球表面,不同的緯度對(duì)應(yīng)的磁傾角、磁偏角、總場(chǎng)大小等磁參數(shù)都不相同,所以,各個(gè)點(diǎn)上的ΔT的計(jì)算誤差E值也隨著這些磁參數(shù)的變化而改變.通過對(duì)東經(jīng)120°、緯度-90°~90°進(jìn)行一系列地磁參數(shù)的換算,得到各緯度所對(duì)應(yīng)的磁傾角、磁偏角、總場(chǎng)等數(shù)值.為評(píng)價(jià)地表各緯度帶上誤差E的變化情況,我們假設(shè)每個(gè)緯度取一個(gè)剖面進(jìn)行誤差E的計(jì)算,然后,分析每個(gè)剖面上E的最大值分布.繼續(xù)用之前的二度水平圓柱體模型,為觀察對(duì)比,等效磁化率分別取0.1 SI、1.0 SI 和3.0 SI.

        在東經(jīng)120° 經(jīng)線上,地理緯度-90°~90°對(duì)應(yīng)的磁傾角范圍為-85°~89°,地磁場(chǎng)總強(qiáng)度大小的變化范圍為39954~65960 nT.誤差E的最大值隨磁傾角變化的規(guī)律如圖5a所示,可以看出,曲線關(guān)于磁傾角0°是對(duì)稱關(guān)系,曲線在磁赤道和磁兩極處總是處于一個(gè)極小值的狀態(tài);并且不管磁化率κ如何變化,曲線極大值總是出現(xiàn)在中緯度范圍,整個(gè)曲線的幅值隨磁化率κ的增加而增加.可以看出,磁化率κ=0.1 SI時(shí),ΔT誤差最大值曲線和相對(duì)誤差曲線幅值都很小,在圖中近似一條直線,可以忽略.磁化率κ=1.0 SI時(shí),曲線在I=±50°處達(dá)到2000 nT的極值,是此時(shí)磁異常Ta的14.3%,且從圖5b中可以看到此時(shí)相對(duì)誤差整體在7%以上,最大值達(dá)到13.0%;隨磁化率增長(zhǎng)到κ=3.0 SI時(shí),誤差最大值曲線整體大于6000 nT,最大值達(dá)到16000 nT,是此時(shí)磁異常Ta的38.6%,而相對(duì)誤差也整體處于15%~40%之間,可見已不能忽略.所以,如果以ΔT的計(jì)算誤差和相對(duì)誤差5.0%為誤差衡量的一個(gè)標(biāo)準(zhǔn),那么磁化率超過1.0 SI時(shí)我們通常就要開始重視ΔT計(jì)算誤差了.

        磁化率不同,對(duì)應(yīng)的ΔT誤差最大值曲線和相對(duì)誤差曲線形態(tài)基本一致.磁化率越高,曲線幅值就越大,且集中在磁傾角±30°~±50°區(qū)域,對(duì)應(yīng)的緯度帶為-20°~-60°及20°~60°的中低緯度地區(qū),這說明中低緯度地區(qū)不僅ΔT誤差更大,相對(duì)誤差也更大.

        4.2T0、Ta與ΔT計(jì)算誤差關(guān)系

        假設(shè)在某測(cè)點(diǎn)上,地磁場(chǎng)T0與磁異常大小Ta均不變,但T0與Ta夾角θ發(fā)生改變時(shí),兩者關(guān)系可用圖6a表示.顯然,由公式(7)可知ΔT模量差=0時(shí),誤差為極大值:

        (13)

        此時(shí),圖6a中T0、T、Ta三者構(gòu)成等腰三角形,Ta與T0的角度θ為:

        (14)

        從式(7)和圖6b也可以看出,Ta變化對(duì)應(yīng)的其曲線幅值并不是線性增長(zhǎng)關(guān)系.我們?cè)O(shè)正常場(chǎng)強(qiáng)度取值為T0=50000 nT,磁異常Ta取值為Ta/T0=0.1,0.5及1.0,對(duì)應(yīng)的誤差E隨二者夾角θ變化如圖6b所示.可以看出,當(dāng)θ角為0°、180°時(shí),Ta與T0的方向在一條直線上,誤差E的值為0,且誤差曲線都關(guān)于θ=180°對(duì)稱.

        綜上,Ta/T0在0~1范圍時(shí),Emax出現(xiàn)的位置范圍是θ=90°~120°(和θ=240°~270°)之間,且Emax≤0.5Ta.

        根據(jù)誤差極值關(guān)系式(13),我們可以很容易得到一些指標(biāo)性信息,見表1.例如,現(xiàn)在的質(zhì)子磁力儀和光泵磁力儀的測(cè)量精度已分別達(dá)到0.01 nT和0.001 nT,則當(dāng)測(cè)點(diǎn)上磁異常幅值達(dá)到Ta=10 nT即可引起0.001 nT的最大誤差;高精度測(cè)量通常要求精度達(dá)到0.01 nT,則易得Ta=31.6 nT即可引起這樣的誤差;而100 nT的中等強(qiáng)度磁異常,可引起0.1 nT的誤差;磁異常達(dá)到1000 nT的幅值時(shí),引起的誤差E極大值可達(dá)10 nT,以高精度處理解釋的標(biāo)準(zhǔn)衡量,這個(gè)誤差已經(jīng)不能忽視了.過去認(rèn)為磁異常小于5000 nT的幅值時(shí),其誤差可忽略(Blackely,1996;管志寧,2005),而此時(shí)引起的誤差E極大值可達(dá)250 nT,這顯然表明,在強(qiáng)磁異常及高精度勘探階段,其誤差及其影響是需要考慮的了.

        表1 誤差最大值Emax與磁異常Ta的相互對(duì)應(yīng)關(guān)系

        圖5 東經(jīng)120°各緯度帶的ΔT誤差分布

        圖6 (a) Ta與T0夾角θ在0°~360°范圍內(nèi)變化示意圖; (b) Ta/T0=0.1/0.5/0.1時(shí)ΔT計(jì)算誤差隨二者夾角變化規(guī)律

        另外,我們進(jìn)一步比較分析Ta垂直于T0的特殊情況.此時(shí)ΔT投影≡0,所以ΔT的計(jì)算誤差E=ΔT模量差,如圖7a所示.由于此時(shí)Ta與T0方向不變,假設(shè)地磁場(chǎng)T0大小不變,分析此種情況下E隨Ta大小變化的情況.則由圖7a的圖形關(guān)系,利用下式非常容易計(jì)算該誤差:

        (15)

        例如,T0為50000 nT保持不變,Ta=5000 nT時(shí),E⊥=249.39 nT,非常接近由式(13)確定的最大誤差Emax=250 nT.當(dāng)磁異常Ta=10000 nT時(shí),E⊥=990.20 nT,也接近由式(13)確定的最大誤差Emax=1000 nT,所以,可以認(rèn)為,當(dāng)磁異常幅值小于10000 nT時(shí),最大誤差為磁異常矢量垂直于地磁場(chǎng)方向的測(cè)點(diǎn)附近,這對(duì)于我們基于磁異常剖面定性判斷最大誤差區(qū)域是有幫助的.當(dāng)Ta=50000 nT時(shí),此時(shí)誤差很大,為E⊥=20710.68 nT,而最大誤差Emax=25000 nT,此時(shí)兩者差異已很大,當(dāng)然磁異常如此之大的情況,在實(shí)際勘探中幾乎沒有,所以,可以不考慮.

        4.3 圓柱體ΔT誤差最大值隨r/R變化的特征

        從圖8中可以看出,磁化率在0.1SI時(shí),Emax曲線走勢(shì)平緩,r/R=0.5時(shí),Emax=3.91 nT;r/R=0.8時(shí),Emax=25.60 nT;r/R=1.0,即模型已出露地表時(shí),最大值Emax是62.5nT,高精度勘探的處理解釋中已不能忽略該誤差了.磁化率為1.0SI時(shí),r/R超過0.5之后,Emax曲線增長(zhǎng)趨勢(shì)明顯,r/R=0.8時(shí),Emax達(dá)到2560nT,這個(gè)誤差已經(jīng)非常大了.當(dāng)磁化率達(dá)到3.0SI時(shí),r/R超過0.3之后,Emax曲線增長(zhǎng)趨勢(shì)便非常明顯.

        綜上分析可以看出,ΔT計(jì)算誤差可能會(huì)對(duì)反演工作造成后續(xù)影響.用ΔT模量差的值來進(jìn)行反演計(jì)算造成的誤差是不可忽視的,或者說,采用嚴(yán)格的ΔT模量差值來進(jìn)行反演,其計(jì)算公式需要重新確定,不能使用基于ΔT投影值的方法.由于篇幅原因其反演結(jié)果對(duì)比在這里不再展開論述.

        圖8 磁化率κ變化時(shí),一條剖面上的誤差最大值隨r/R的變化關(guān)系

        圖7 Ta垂直T0時(shí)(a)ΔT計(jì)算誤差示意圖, (b)E⊥及Emax隨Ta的變化曲線

        5 結(jié)論

        基于高精度勘探的需要,我們對(duì)強(qiáng)磁性體ΔT異常計(jì)算存在的誤差進(jìn)行了分析研究和評(píng)估,取得以下認(rèn)識(shí).

        (1) ΔT異常存在的誤差來源于用近似的ΔT投影代替實(shí)測(cè)的ΔT模量差.之所以選擇ΔT投影,是因?yàn)槠溆?jì)算只與地磁場(chǎng)T0方向有關(guān)而與其大小無直接關(guān)系,因而ΔT投影與地磁場(chǎng)T0之間是良好的線性關(guān)系,且便于正反演及處理轉(zhuǎn)換中的計(jì)算;而ΔT模量差與地磁場(chǎng)T0的方向和大小均密切相關(guān),其與T0之間不具有線性關(guān)系特征,因而不便于一系列相應(yīng)的處理計(jì)算,所以,在弱磁異常及非高精度條件下,ΔT投影具有應(yīng)用上的便捷性.

        (2) 地磁場(chǎng)T0方向和大小在勘探測(cè)區(qū)范圍內(nèi)可以認(rèn)為是不變的,而磁性體產(chǎn)生的磁異常Ta的方向和大小則在測(cè)區(qū)內(nèi)不同測(cè)點(diǎn)之間具有明顯的變化;當(dāng)磁異常Ta幅值大時(shí),ΔT投影誤差會(huì)比較大,并且ΔT投影誤差與測(cè)點(diǎn)和場(chǎng)源的相對(duì)位置也有密切的關(guān)系,誤差比例最大的測(cè)點(diǎn)位于磁異常Ta接近垂直于地磁場(chǎng)T0的區(qū)域;從當(dāng)前磁力儀測(cè)量精度等指標(biāo)衡量,10 nT的磁異常能引起的最大誤差只有0.001 nT,100 nT的磁異常能引起的最大誤差也只有0.1 nT,而1000 nT的較強(qiáng)磁異常能引起的最大誤差就達(dá)到10 nT,5000 nT的強(qiáng)磁異常能引起的最大誤差可達(dá)250 nT,隨著磁異常的增大,誤差按近似指數(shù)關(guān)系迅速增大.

        (3) 應(yīng)用二維柱體模型進(jìn)行計(jì)算實(shí)驗(yàn),結(jié)果表明:在地表實(shí)際環(huán)境下,中緯度區(qū)帶的ΔT投影誤差明顯較大;在理論分析研究方面,二維圓柱體模型磁異常的一些屬性會(huì)發(fā)生變化,例如,僅改變磁化方向?yàn)榉聪?,兩次ΔT投影曲線為嚴(yán)格的正負(fù)鏡像關(guān)系,但ΔT模量差曲線則不具有嚴(yán)格對(duì)稱性特征;磁性體的尺度(即半徑)與埋深的比值r/R超過0.5且磁化率超過0.1 SI時(shí)誤差已經(jīng)達(dá)到3.9 nT.所以,強(qiáng)磁化磁性體產(chǎn)生了大幅度磁異常時(shí),其定量處理及解釋需要關(guān)注該誤差因素,具體情況視精度要求而定.

        綜上所述,在進(jìn)入高精度勘探精細(xì)處理解釋階段,ΔT計(jì)算誤差是不容忽視的,特別是強(qiáng)磁異常而言.

        本文只對(duì)該誤差進(jìn)行了分析研究,基于嚴(yán)格無誤差的ΔT模量差,進(jìn)行相應(yīng)的處理解釋,有必要繼續(xù)深化,開展此項(xiàng)研究工作.

        Blakely R J. 1996. Potential Theory in Gravity and Magnetic Applications. Cambridge: Cambridge University Press.

        Clark D A. 1983. Comments on magnetic petrophysics.ExplorationGeophysics, 14(2): 49-62.

        Clark D A, Emerson D W. 1991. Notes on rock magnetization characteristics in applied geophysical studies.ExplorationGeophysics, 22(3): 547-555.

        Fang H Z. 1978. Calculation of a magnetic anomaly by three-dimensional bodies of arbitrary shapes in the condition of high susceptibility (in Chinese): Acta Geol.Sinica.ActaGeol.Sinica, 52(1): 63-78.

        Guan Z N. 2005. Geomagnetic Field and Magnetic Exploration (in Chinese). Beijing: Geological Publish House.

        Guo W W, Dentith M C, Bird R T, et al. 2001. Systematic error analysis of demagnetization and implications for magnetic interpretation.Geophysics, 66(2): 562-570.

        Kontis A L, Young G A. 1964. Approximation of residual total-magnetic-intensity anomalies.Geophysics, 29(4): 623-627.

        Krahenbuhl R A, Li Y G. 2007. Influence of self-demagnetization effect on data interpretation in strongly magnetic environments.ASEGExtendedAbstracts, 713-717.

        Li C M, Li J. 2013. Gravity and Magnetic Prospecting Principle and Method (in Chinese). Beijing: Science Publishing House.

        Li Y G, Krahenbuhl R A. 2011. A General Framework for 3D Interpretation of Magnetic Data Affected by Remanence and Self-Demagnetization.∥12th International Congress of the Brazilian Geophysical Society.

        Zhang H L, Liu T Y, Zhu C J, et al. The effects of applying high-precision magnetic survey: a case study of the Galinge ore district in Qinghai province.GeophysicalandGeochemicalExploration(in Chinese), 35(1): 12-16.

        Zhang H L, Hu X Y. 2012. A stable downward continuation and its application for mineral is exploration.∥2012 SEG Annual Meeting. Society of Exploration Geophysicists.

        Zhu C J, Zhou Z W, Liu T Y, et al. 2011. Application of high-precision magnetic survey to prospecting: A case study in the Galinge Deposit of Qinghai Province.GeologyandExploration(in Chinese), 47(2:277-283).

        Zhu Y G, Yu C C. 2006. Analysis of magnetic properties of Daye iron deposit in Hubei province.ContributionstoGeologyandMineralResourcesResearch(in Chinese), 21(S1): 155-159.

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

        方華竹. 1978. 任意形狀強(qiáng)磁性三度體磁異常的計(jì)算. 地質(zhì)學(xué)報(bào), 52(1): 63-78.

        管志寧. 2005. 地磁場(chǎng)與磁力勘探. 北京:地質(zhì)出版社.

        李才明,李軍. 2013. 重磁勘探原理與方法. 北京:科學(xué)出版社.

        張恒磊, 劉天佑, 朱朝吉, 等. 2011. 高精度磁測(cè)找礦效果—以青海尕林格礦區(qū)為例. 物探與化探, 35(1): 12-16.

        朱朝吉,周肇武, 劉天佑等. 2011. 高精度磁測(cè)找礦效果: 以青海尕林格礦區(qū)為例. 地質(zhì)與勘探, 47(2): 277-283.

        朱永剛, 于長(zhǎng)春. 2006. 湖北省大冶鐵礦區(qū)內(nèi)礦石磁性特征分析. 地質(zhì)找礦論叢, 21(B10): 155-159.

        (本文編輯 胡素芳)

        Error analysis of calculation of total field anomaly due to highly magnetic bodies

        YUAN Xiao-Yu, YAO Chang-Li*, ZHENG Yuan-Man, LI Ze-Lin, WANG Jun-Heng

        KeyLaboratoryofGeo-detection(ChinaUniversityofGeosciences,Beijing),MinistryofEducation,Beijing100083,ChinaStateKeyLaboratoryofGeologicalProcessesandMineralResources,Beijing100083,China

        Currently, the magnetic prospecting has a great demand for high-precision processing and interpretation. In order to meet this requirement, we studied the approximate errorEbetween ΔTact(the total-field anomaly) and ΔTpro(the projection of anomalous field vector onto the direction of geomagnetic field). Generally speaking, the errorEis extremely small when the bodies have weakly magnetic susceptibilities. However, when the bodies have highly magnetic susceptibilities, the errorEmight be large. This will lead to significant effects on subsequent quantitative inference. Therefore, we investigated the errorEdue to highly magnetic bodies. In this paper, a systematic error analysis was made by using a 2-D elliptic cylinder model. We investigated the errorEthrough numerical experiments of magnetic anomalies produced by high-susceptibility bodies. Normally, for high susceptibilities, we found that the existed remanence significantly affected the errorE. The error analysis showed that the magnitude of ΔTactwas usually larger than that of ΔTpro. This implied that a theoretical anomaly computed without accounting for the errorEwould overestimate the anomaly associated with the body.we used the 2-D elliptic cylinder model to illustrate the importance of the error analysis when the bodies were highly magnetic. Firstly, we demonstrated through numerical experiments that the errorEwas obvious and should not be ignored. Secondly, we showed that the curves of ΔTproand the errorEhad a certain symmetry when the directions of magnetization and geomagnetic field changed. And we also displayed that theEmax(the maximum value of the errorE) appeared above the center of the magnetic body when the parameters such as the vector of geomagnetic fieldT0, the geometric form (e.g., the radiusrand the depthRof the 2D horizontal cylinder body) and the magnetic parameters (e.g., the susceptibilityκ) are determined, and that the relationship between theEmaxand magnetic parameters or size of the model (e.g.,r/R) is exponential. Then we discovered some other characteristics about the errorE. For instance, the curve ofEmaxwith respect to the latitude was symmetrical on both sides of magnetic equator, and the extremum of theEmaxcan always be found in the mid-latitudes. WhenT0was perpendicular toTa(the vector of the anomaly), the errorErised with the increase ofTa. Whenθ(the included angle betweenTaandT0) changed, theEmaxdue to the cylinder model is found inθ=90°~120°. Whenr/Ris larger than 0.5 and the susceptibility is larger than 0.1SI, theEmaxreaches to 3.9 nT. As a result of the numerical experiments, we concluded that when the bodies have highly magnetic susceptibilities, the errorEmay be great and will affect the subsequent magnetic processing and inversion. Therefore, the errorEcannot be ignored when the magnetic data are processed, inverted, and interpreted in highly magnetic environments.

        Total-field anomaly; Approximate error; Vector of magnetic anomaly; Remanence; Susceptibility

        國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目(2014AA06A613)和國(guó)家自然科學(xué)基金項(xiàng)目(41304104)資助.

        袁曉雨,主要從事重磁處理與反演研究.E-mail:792010739@qq.com

        *通信作者 姚長(zhǎng)利,主要從事重磁勘探理論與方法技術(shù)研究. E-mail:clyao@cugb.edu.cn

        10.6038/cjg20151235.

        10.6038/cjg20151235

        P318,P631

        2015-04-13,2015-10-09收修定稿

        袁曉雨, 姚長(zhǎng)利, 鄭元滿等. 2015. 強(qiáng)磁性體ΔT異常計(jì)算的誤差分析研究.地球物理學(xué)報(bào),58(12):4756-4765,

        Yuan X Y, Yao C L, Zheng Y M, et al. 2015. Error analysis of calculation of total field anomaly due to highly magnetic bodies.ChineseJ.Geophys. (in Chinese),58(12):4756-4765,doi:10.6038/cjg20151235.

        猜你喜歡
        磁化率磁化模量
        高勁度模量瀝青混合料在京臺(tái)高速車轍維修段的應(yīng)用
        室內(nèi)回彈模量和回彈再壓縮模量試驗(yàn)參數(shù)探討
        山西建筑(2020年11期)2020-06-04 00:09:48
        東北豐磁化炭基復(fù)合肥
        雙色球磁化炭基復(fù)合肥
        關(guān)于現(xiàn)行規(guī)范路基頂面回彈模量的理解和應(yīng)用
        上海公路(2018年4期)2018-03-21 05:57:24
        基于超拉普拉斯分布的磁化率重建算法
        基于磁化能量的鋰電池串模塊化均衡方法
        巖(礦)石標(biāo)本磁化率測(cè)定方法試驗(yàn)及認(rèn)識(shí)
        超強(qiáng)磁場(chǎng)下簡(jiǎn)并電子氣體的磁化
        溫度對(duì)不同初始狀態(tài)ising模型磁化強(qiáng)度和磁化率的影響
        河南科技(2014年18期)2014-02-27 14:15:23
        国产精品98视频全部国产| 少妇被粗大的猛烈进出69影院一| 欧美亚洲国产片在线播放| 亚洲日韩图片专区小说专区| 亚洲一区二区女优av| 99国产精品久久一区二区三区| 亚洲av日韩综合一区二区三区| 六月丁香婷婷色狠狠久久| 中文字幕大乳少妇| 日本韩国亚洲三级在线| 国产亚洲美女精品久久久2020| 欧美不卡视频一区发布| 国产自精品在线| 国产熟女白浆精品视频二| 台湾佬中文网站| 久久久久亚洲精品天堂| 吃下面吃胸在线看无码| 国产一级黄色片在线播放| 亚洲国产精品无码中文字| 亚洲av熟妇高潮30p| 亚洲va精品va国产va| 丰满人妻中文字幕一区三区| 在线成人一区二区| 日韩成人精品在线| 国产精品一区又黄又粗又猛又爽| 亚洲av高清一区二区三| 中年熟妇的大黑p| 久久精品国产亚洲AV无码不| av一区二区三区有码| 国产精品久久久久一区二区三区| 国产99久久亚洲综合精品| 麻豆国产AV网站| 亚洲女人的天堂网av| 成l人在线观看线路1| YW亚洲AV无码乱码在线观看| 久久亚洲av熟女国产| 亚洲av无码专区在线| 粗了大了 整进去好爽视频| 国产一区二区在线观看我不卡| 亚洲av区,一区二区三区色婷婷| 高潮迭起av乳颜射后入|