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

    基于振幅—相位張量分解的大地電磁電流型畸變校正

    2022-12-09 07:12:16余才國(guó)湯井田皇祥宇徐錦通
    石油地球物理勘探 2022年6期
    關(guān)鍵詞:張量畸變振幅

    余才國(guó) 肖 曉*②③④ 湯井田②③④ 皇祥宇 徐錦通

    (①中南大學(xué)地球科學(xué)與信息物理學(xué)院,湖南長(zhǎng)沙 410083;②有色資源與地質(zhì)災(zāi)害探查湖南省重點(diǎn)實(shí)驗(yàn)室,湖南長(zhǎng)沙 410083;③中南大學(xué)有色金屬成礦預(yù)測(cè)與地質(zhì)環(huán)境監(jiān)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室,湖南長(zhǎng)沙 410083;④自然資源部覆蓋區(qū)深部資源勘查工程技術(shù)創(chuàng)新中心,安徽合肥 230001)

    0 引言

    大地電磁測(cè)深(MT)通過觀測(cè)地表相互正交的電磁場(chǎng)分量獲取地下電性結(jié)構(gòu)相關(guān)信息,具有不受高阻層屏蔽、對(duì)高導(dǎo)層分辨能力強(qiáng)、橫向電性差異分辨能力較強(qiáng)、勘探深度大、勘探費(fèi)用低、施工方便、資料處理和解釋技術(shù)成熟等優(yōu)勢(shì),廣泛應(yīng)用于區(qū)域地質(zhì)調(diào)查、資源勘探和工程勘探等領(lǐng)域[1-3]。近年來,隨著研究涉及的地質(zhì)構(gòu)造復(fù)雜程度明顯增加,對(duì)大地電磁測(cè)深資料的處理解釋要求亦更高[4]。當(dāng)近地表存在幾何尺寸遠(yuǎn)遠(yuǎn)小于電磁波趨膚深度的三維電性不均勻小異常體時(shí),觀測(cè)數(shù)據(jù)會(huì)受到電流型畸變的影響,若不消除這種影響,會(huì)影響后續(xù)的反演解釋[5-10]。

    阻抗張量分解方法是目前常用的大地電磁數(shù)據(jù)處理分析方法,可以去除局部畸變帶來的影響,獲得真實(shí)的區(qū)域阻抗、構(gòu)造方位角、構(gòu)造維性等信息。Swift[11]提出了經(jīng)典的阻抗張量分解方法,具體做法是:旋轉(zhuǎn)阻抗張量,尋找阻抗張量反對(duì)角元素最大時(shí)或者主對(duì)角元素最小時(shí)所對(duì)應(yīng)的角度,該角度即為所求的電性主軸方向。但是,此方法獲得的電性主軸方向存在90°的多解性,故需要依據(jù)實(shí)際地質(zhì)信息進(jìn)一步判斷。該方法為后續(xù)的阻抗分解方法奠定了堅(jiān)實(shí)的基礎(chǔ)。Berdichevsky等[12]研究了水平電性非均勻性引起的畸變問題,并將畸變劃分為感應(yīng)型畸變和電流型畸變。Bahr[13]根據(jù)局部畸變的觀測(cè)阻抗的同一列元素相位相等的特性進(jìn)行了畸變校正,提出了Bahr分解方法,但是對(duì)非嚴(yán)格意義上的三維/二維構(gòu)造效果不理想。Bahr[14]在此方法的基礎(chǔ)上提出了Bahr相位偏差法,即阻抗張量同一列元素的阻抗相位存在一個(gè)角度差,這個(gè)角度差稱為相位偏差角,并定義了與維性密切相關(guān)的表征參數(shù)。Groom等[15]提出了經(jīng)典的GB分解法,將局部畸變分為可確定部分(剪切因子和扭曲因子)和不可確定部分(各向異性因子和靜位移因子),然而GB分解得到的模型是一個(gè)欠定問題:8個(gè)已知參數(shù)(4個(gè)復(fù)數(shù)形式的觀測(cè)阻抗張量元素),9個(gè)未知參數(shù)(電性主軸角,4個(gè)實(shí)數(shù)形式的畸變矩陣元素,2個(gè)復(fù)數(shù)形式的區(qū)域阻抗張量元素)。需要特別強(qiáng)調(diào)的是,盡管GB分解方法得到了廣泛應(yīng)用,取得了很好的效果,但實(shí)際應(yīng)用中往往忽略了其應(yīng)用前提,即地下介質(zhì)為二維或近似二維的電性結(jié)構(gòu)[16]。Caldwell等[17]提出了相位張量的概念,基于相位張量不受電場(chǎng)畸變的影響這一特點(diǎn),從觀測(cè)阻抗中提取區(qū)域阻抗的構(gòu)造電性主軸方向以及二維偏度等表征構(gòu)造屬性的參數(shù),是目前最常用的維性分析方法之一。Bibby等[18]分析了阻抗張量畸變與相位張量的關(guān)系,將相位張量用于局部畸變校正。蔡軍濤等[19]提出一種共軛阻抗變換法(CCZ),通過對(duì)大地電磁阻抗張量做共軛阻抗變換,求取畸變因子和最佳電性主軸方位角,該方法具有與相位張量分解相似的優(yōu)點(diǎn),即不受局部電場(chǎng)畸變的影響,且無(wú)需假設(shè)區(qū)域維性。尹曜田等[5]將遺傳算法引入GB分解的求解過程,通過實(shí)際資料處理證明了該算法的有效性。謝成良等[20]使用相位張量分解獲得電性主軸方位角,將旋轉(zhuǎn)后的區(qū)域阻抗作為GB分解區(qū)域阻抗元素的初值,并將各向異性參數(shù)引入目標(biāo)函數(shù),采用共軛梯度法進(jìn)行模型最優(yōu)化求解,提高了GB算法的穩(wěn)定性及可靠性。

    但前人工作大多集中于二維區(qū)域構(gòu)造。對(duì)于三維區(qū)域結(jié)構(gòu),一些學(xué)者已經(jīng)開展了一些研究。Utada等[21]利用觀測(cè)阻抗張量和觀測(cè)地磁傳遞函數(shù)(即傾子)確定畸變參數(shù),壓制電流型畸變,獲得三維區(qū)域電性構(gòu)造。Garcia等[22]假設(shè)相鄰測(cè)點(diǎn)的區(qū)域電性結(jié)構(gòu)變化不大但受畸變的影響不同,提出了一種類似于GB分解且適用于三維區(qū)域的畸變校正方法,使用牛頓法求解三維區(qū)域結(jié)構(gòu)中的畸變參數(shù)。Bibby等[18]利用不受電流型畸變影響的相位張量獲得地下維度信息,進(jìn)而采用低維度數(shù)據(jù)求解畸變參數(shù)。這些方法只適用于特定條件,對(duì)于更一般的三維/三維模型不適用。Neukirch等[23]提出了與相位張量互補(bǔ)的振幅張量,利用“相位張量和振幅張量反映同樣的地下介質(zhì)信息,且所有電流型畸變信息都保存在振幅張量中”這一假說,采用遺傳算法實(shí)現(xiàn)了畸變校正,該方法的實(shí)測(cè)效果及理論效果均較好[24]。

    本文借鑒Neukirch等[24]的研究成果,基于阻抗振幅張量與相位張量的聯(lián)合分解,采用改進(jìn)粒子群算法求解畸變參數(shù),對(duì)三維模型正演數(shù)據(jù)和安徽省霍山地區(qū)中生代火山巖盆地實(shí)測(cè)MT數(shù)據(jù)進(jìn)行處理。

    1 大地電磁振幅相位張量分解基本理論

    1.1 GB分解理論

    對(duì)于三維/二維電性結(jié)構(gòu)模型,Groom等[15]只考慮電流型畸變影響,不考慮磁場(chǎng)感應(yīng)型畸變,將區(qū)域構(gòu)造主軸坐標(biāo)系順時(shí)針旋轉(zhuǎn)角度θ到測(cè)量坐標(biāo)系,則觀測(cè)阻抗Zm為

    Zm=R(θ)CZ2DRT(θ)

    =R(θ)gTSAZ2DRT(θ)

    (1)

    式中:R(θ)為單位旋轉(zhuǎn)矩陣;C為畸變矩陣;Z2D為區(qū)域阻抗矩陣;g表示增益因子;T表示扭曲角矩陣;S表示剪切角矩陣;A表示各向異性因子矩陣?;兙仃嘋可表示為

    (2)

    式中:扭曲因子t=tanφt,其中φt∈[-90°,90°]為扭曲角;剪切因子e=tanφe,其中φe∈[-45°,45°]為剪切角;各向異性因子s=tanφs,其中φs∈[-45°,45°]為各向異性因子角。

    GB分解把畸變矩陣分解為可確定部分(扭曲角矩陣T、剪切角矩陣S)和不確定部分(各向異性因子矩陣A與增益因子g的乘積),并將不確定部分并入?yún)^(qū)域阻抗矩陣Z2D,得到

    (3)

    式中ZTE和ZTM分別表示二維區(qū)域構(gòu)造下TE模式和TM模式的阻抗張量元素。此時(shí),式(1)變?yōu)?/p>

    (4)

    上式中右邊8個(gè)參數(shù)(觀測(cè)阻抗的實(shí)部、虛部)為已知,7個(gè)參數(shù)(畸變參數(shù)t、e、θ和區(qū)域阻抗反對(duì)角元素的實(shí)部、虛部)為未知,利用最優(yōu)化求解即可獲得其解??梢钥闯?,對(duì)于三維介質(zhì),阻抗張量的區(qū)域阻抗對(duì)角元素中也存在未知數(shù),未知數(shù)為11個(gè),故傳統(tǒng)的GB分解無(wú)法適用于三維介質(zhì)的情況。使用GB分解的前提是區(qū)域構(gòu)造為二維,對(duì)于三維區(qū)域構(gòu)造不再適用,因?yàn)槿S區(qū)域構(gòu)造增加了4個(gè)未知數(shù),此時(shí)的未知數(shù)個(gè)數(shù)大于方程個(gè)數(shù),傳統(tǒng)的阻抗分解方法失效,需要探索新的方法求解畸變參數(shù)。

    1.2 相位張量

    大地電磁阻抗張量Z為一個(gè)復(fù)數(shù)張量,既含有電阻率信息,也含有相位信息,可以表示為Z=X+iY,這里X和Y分別表示阻抗張量的實(shí)部和虛部。根據(jù)Caldwell等[17]對(duì)相位張量的定義,相位張量Φ的表達(dá)式為

    (5)

    式中:Xpq、Ypq分別表示阻抗張量元素的實(shí)部和虛部,其中p=1、2,q=1、2;det(X)=X11X22-X12X21;Φpq表示相位張量Φ的元素。

    只考慮電場(chǎng)畸變影響時(shí),畸變矩陣C是一個(gè)與頻率無(wú)關(guān)的實(shí)數(shù)矩陣,觀測(cè)阻抗Zm可以寫為

    Zm=CZ=C(X+iY)

    (6)

    Xm=CX

    (7)

    Ym=CY

    (8)

    根據(jù)式(6),若能求得畸變矩陣C,即可計(jì)算得到未畸變的阻抗張量Z

    Z=C-1Zm

    (9)

    根據(jù)相位張量的定義,此時(shí)的觀測(cè)相位張量Φm可根據(jù)下式求得

    (10)

    由上式可以看出,相位張量不受局部電場(chǎng)畸變的影響,因此可廣泛應(yīng)用于區(qū)域構(gòu)造走向的研究和維度分析。相位張量概念的提出對(duì)于阻抗張量分解具有里程碑式的意義,為大地電磁數(shù)據(jù)分析提供了新的思路和方法。但同時(shí)也意味著采用相位張量無(wú)法判斷測(cè)點(diǎn)受畸變的程度,因而只采用不含振幅信息的相位張量進(jìn)行反演,會(huì)很大程度增加反演的非唯一性,導(dǎo)致反演不穩(wěn)定。

    1.3 振幅張量

    通過相位張量的定義發(fā)現(xiàn),相位張量是阻抗相位的混合,只包括阻抗的相位信息,僅是阻抗張量信息的一半,這對(duì)于測(cè)點(diǎn)本就稀疏、數(shù)據(jù)不足的大地電磁來說,是一項(xiàng)非常大的損失。Neukirch等[23]為了提取振幅信息,創(chuàng)造性地提出了類似于相位張量的振幅張量。

    根據(jù)Neukirch等[23]的研究成果,一個(gè)復(fù)數(shù)M可用歐拉公式展開為一個(gè)實(shí)數(shù)L和一個(gè)角度β的函數(shù)

    M=Lei β=L(cosβ+isinβ)

    (11)

    按照式(11)的形式,可將單個(gè)復(fù)數(shù)擴(kuò)展到二階復(fù)數(shù)張量的形式,則阻抗張量Z可以表示為

    Z=PeiΦ=P[cos(Φ)+isin(Φ)]

    (12)

    式中:為了書寫方便,cos(·)和sin(·)表示對(duì)一個(gè)矩陣的各元素分別求余弦和正弦;P為振幅張量。由于采用張量表示形式,下面給出振幅張量的具體求取方法。

    張量下的約束條件為

    cos(Φ)[cos(Φ)]T+sin(Φ)[sin(Φ)]T=I

    (13)

    式中I為二階單位矩陣。通過式(13)和式(5)可得

    (14)

    (15)

    eiΦ=cos(Φ)+isin(Φ)

    (16)

    將式(16)代入式(12),便可得到振幅張量P的唯一表示形式

    (17)

    2 振幅—相位張量分解最優(yōu)化畸變校正算法

    2.1 目標(biāo)函數(shù)

    Booker[25]提出了一種相位張量參數(shù)化分解方法,把相位張量的奇異值分解與區(qū)域構(gòu)造走向、維度等地質(zhì)信息聯(lián)系起來,這種分解方法同樣適用于同為實(shí)數(shù)矩陣的振幅張量。振幅張量和相位張量都是地下介質(zhì)、構(gòu)造信息的反映,在某種程度上其表現(xiàn)形式相似,而且相位張量不受電流型畸變的影響,所有的畸變信息都隱含在振幅張量里,通過尋找與相位張量相似性最高的振幅張量,便可實(shí)現(xiàn)大地電磁的畸變校正。其中,畸變矩陣C為待求參數(shù),采用式(9)對(duì)觀測(cè)張量阻抗Zm進(jìn)行校正,根據(jù)式(17)獲得振幅張量P,根據(jù)式(5)獲得相位張量Φ,并對(duì)這個(gè)張量進(jìn)行分解。求解畸變參數(shù)C時(shí),振幅張量P的分解參數(shù)與相位張量Φ的分解參數(shù)最相似。因此,根據(jù)以上思路構(gòu)建了關(guān)于畸變參數(shù)的求解策略,將畸變求解問題轉(zhuǎn)化為一個(gè)最優(yōu)化問題。

    相位張量矩陣使用奇異值分解可表示為[25]

    (18)

    式中VΦ=R(-θΦ)、ΣΦ和WΦ=R(-ψΦ-θΦ)是奇異值分解后的三個(gè)矩陣,其中旋轉(zhuǎn)矩陣R的表達(dá)式為

    (19)

    θΦ=arccos[Tr(VΦ)/2]表示笛卡爾坐標(biāo)與矩陣坐標(biāo)之間的角度,這里Tr(·)表示求矩陣的跡,在二維情況下,θΦ表示區(qū)域構(gòu)造走向;φ1和φ2是對(duì)角矩陣ΣΦ的奇異值;偏離角ψΦ可表示為

    --------------------

    (20)

    --------------------

    同樣地,對(duì)振幅張量進(jìn)行上述分解,可以得到與相位張量分解參數(shù)相對(duì)應(yīng)的各個(gè)分解參數(shù)。

    為了使相位張量與振幅張量的奇異值具有可比性,Neukirch等[23]定義了相位相對(duì)各向異性參數(shù)ηφ和振幅相對(duì)各向異性參數(shù)ηP,其表達(dá)式分別為

    (21)

    (22)

    式中P1和P2分別是振幅張量奇異值分解的兩個(gè)奇異值。

    假設(shè)某一測(cè)點(diǎn)有N個(gè)觀測(cè)頻率,且其電流型畸變與頻率無(wú)關(guān)??上惹蟮眠@N個(gè)頻率的振幅張量和相位張量,然后進(jìn)行上述分解,得到各頻點(diǎn)相應(yīng)的偏離角ψ、相位相對(duì)各向異性參數(shù)ηφ和振幅相對(duì)各向異性參數(shù)ηP,所有頻率同一分解參數(shù)皆可組成1×N階矩陣。為了對(duì)振幅張量和相位張量的分解參數(shù)進(jìn)行對(duì)比,定義目標(biāo)函數(shù)

    (23)

    通過目標(biāo)函數(shù)值的大小可以衡量相位張量與振幅張量的相似程度,從而通過求目標(biāo)函數(shù)的最小值尋找與相位張量相似程度最高的振幅張量,進(jìn)而得到未畸變阻抗張量,實(shí)現(xiàn)畸變校正。

    2.2 粒子群優(yōu)化算法

    為了求取目標(biāo)函數(shù)(式(23))的最小值,可以使用簡(jiǎn)單有效的粒子群優(yōu)化算法,將式(23)的值作為粒子適應(yīng)度,此問題只有三個(gè)未知數(shù),使用粒子群算法能夠有效求解。

    粒子群算法(Particle Swarm Optimization,PSO)是由Kennedy等[26]提出的一種基于群體智能進(jìn)化的計(jì)算方法。在給定粒子個(gè)數(shù)后,隨機(jī)初始化粒子的位置和速度,然后通過迭代方法優(yōu)化問題。每一次迭代中,粒子通過跟蹤兩個(gè)“極值”更新自身的位置和速度:一個(gè)是粒子本身所找到的最優(yōu)解,即個(gè)體極值Pbest,另一個(gè)是整個(gè)粒子群中所有粒子在歷代搜索過程中所達(dá)到的最優(yōu)解,即全局極值Gbest。每個(gè)粒子都不斷地改變自身在解空間中的速度,當(dāng)找到這兩個(gè)最優(yōu)解后,可進(jìn)一步實(shí)現(xiàn)PSO中最重要的“加速”過程,以盡可能地朝Pbest和Gbest所指向的區(qū)域“飛”去。

    假設(shè)對(duì)于M個(gè)粒子,在一個(gè)D維空間進(jìn)行搜索,粒子i的信息可用兩個(gè)D維向量表示,即位置向量xi和速度向量vi

    xi=(xi1,xi2,…,xiD)T

    (24)

    vi=(vi1,vi2,…,viD)T

    (25)

    通過測(cè)試發(fā)現(xiàn),使用標(biāo)準(zhǔn)PSO算法容易陷入局部極小,不能很好地求解畸變參數(shù)。為了提高種群多樣性,避免產(chǎn)生早熟收斂現(xiàn)象,本文對(duì)Clerc[27]提出的具有收縮因子的PSO算法進(jìn)行了改進(jìn)。找到Pbest和Gbest后,粒子即可根據(jù)下式更新速度和位置

    vid=k[vid+c1r1(Pbest,d-xid)+

    c2r2(Gbest,d-xid)]

    (26)

    xid=xid+vid

    i=1,2,…,M;d=1,2,…,D

    (27)

    c1=c2=2.0+rand(0,0.5)

    (28)

    通過改進(jìn)學(xué)習(xí)因子,收縮因子k在每一次迭代中也會(huì)相應(yīng)改變,這樣既可以增強(qiáng)粒子的活性,加強(qiáng)算法的全局探索能力,也可以避免因選取的收縮因子過大而錯(cuò)過最優(yōu)值。

    數(shù)值實(shí)驗(yàn)結(jié)果表明,當(dāng)粒子群體數(shù)量M為32時(shí),速度限制為[-3,3],針對(duì)此場(chǎng)景問題表現(xiàn)出了良好的效果[28]。

    2.3 計(jì)算流程

    增益因子g的求解非常困難,到目前為止并沒有一個(gè)特別好的求解策略,因此本文將增益因子g并入?yún)^(qū)域阻抗,只考慮扭曲因子t、剪切因子e和各向異性因子s。利用本文方法進(jìn)行畸變校正的流程見圖1,具體求解過程如下。

    圖1 本文方法畸變校正流程

    (1)給定32個(gè)初始粒子,隨機(jī)初始化每個(gè)粒子的空間位置x(φt,φe,φs)和速度v,空間位置的三維坐標(biāo)與求解的畸變參數(shù)扭曲角φt、剪切角φe和各向異性因子角φs一一對(duì)應(yīng),需同時(shí)給定畸變參數(shù)的取值范圍;

    (2)利用生成粒子的位置信息生成畸變矩陣的逆C-1,利用式(9)校正阻抗張量Zm,然后利用式(5)和式(17)計(jì)算振幅張量P和相位張量Φ;

    (3)對(duì)振幅張量P和相位張量Φ進(jìn)行參數(shù)化分解,根據(jù)分解得到的參數(shù)計(jì)算式(23)目標(biāo)函數(shù)的值F(C),將得到的目標(biāo)函數(shù)值作為粒子的適應(yīng)度;

    (4)根據(jù)步驟(3)計(jì)算得到的值不斷更新粒子的位置x(φt,φe,φs)和速度v,直到達(dá)到預(yù)設(shè)條件,迭代達(dá)到最優(yōu),迭代終止。

    本文迭代終止條件設(shè)置為固定的迭代次數(shù)。經(jīng)過數(shù)值實(shí)驗(yàn),迭代次數(shù)設(shè)置為1000可以滿足此問題的求解。由于求解參數(shù)只有3個(gè),故無(wú)需考慮求解速度。

    3 模型實(shí)驗(yàn)和實(shí)例分析

    3.1 理論模型

    為了驗(yàn)證方法的可靠性,設(shè)計(jì)了一組三維模型,對(duì)三維阻抗張量進(jìn)行畸變校正。圖2所示為三維電阻率模型及測(cè)點(diǎn)分布圖。在100Ω·m的均勻半空間中分布著三個(gè)電阻率異常體,其電阻率分別為1000、5000、10Ω·m。點(diǎn)距為200m,共25個(gè)測(cè)點(diǎn),測(cè)點(diǎn)site01~site25從左到右測(cè)點(diǎn)號(hào)依次增大,計(jì)算頻率為0.01~1000Hz,19個(gè)頻點(diǎn)對(duì)數(shù)等間隔分布。使用MT3D軟件進(jìn)行正演計(jì)算[29]。

    圖2 三維理論模型及測(cè)點(diǎn)示意圖

    對(duì)正演得到的數(shù)據(jù)施加局部電場(chǎng)畸變影響,即人為添加畸變[30-31]。選取8個(gè)點(diǎn)位分布比較典型的測(cè)點(diǎn)(site01、site05、site08、site10、site13、site15、site19、site23),隨機(jī)生成了各測(cè)點(diǎn)的畸變參數(shù)。使用本文方法對(duì)其中4個(gè)測(cè)點(diǎn)(site10、site13、site15和site19)畸變數(shù)據(jù)進(jìn)行校正,結(jié)果如圖3所示。

    由圖3可見,施加局部電場(chǎng)畸變對(duì)相位的影響明顯弱于對(duì)視電阻率的影響,且視電阻率曲線首支出現(xiàn)不同程度的分離;校正后的視電阻率曲線與未施加畸變的視電阻率曲線幾乎重合,可見很好地恢復(fù)了區(qū)域電性結(jié)構(gòu),證明此方法對(duì)于恢復(fù)除靜位移因子外的其他三個(gè)畸變參數(shù)(t、e和s)是有效的。所選取測(cè)點(diǎn)隨機(jī)生成的畸變參數(shù)及恢復(fù)的參數(shù)如表1所示??梢姡齭ite19受畸變嚴(yán)重的點(diǎn)恢復(fù)效果有所欠缺外,其他3個(gè)測(cè)點(diǎn)各畸變參數(shù)相對(duì)誤差都小于3%,對(duì)于人為添加畸變數(shù)據(jù)利用本文方法恢復(fù)效果良好。

    表1 三維理論模型畸變參數(shù)表

    圖3 三維理論模型測(cè)點(diǎn)site10(a)、site13(b)、site15(c)、site19(d)畸變數(shù)據(jù)校正前、后視電阻率(上)和相位(下)對(duì)比

    3.2 實(shí)測(cè)數(shù)據(jù)

    安徽霍山—舒城盆地位于金寨斷裂南側(cè),屬山間拗陷盆地(圖4),呈東西向,基底除見有少量盧鎮(zhèn)關(guān)巖群、佛子嶺巖群變質(zhì)巖外,多為中、晚侏羅世三尖鋪組、鳳凰臺(tái)組山間磨拉石建造,盆地充填物中有毛坦廠組火山巖[32]。對(duì)霍山地區(qū)中生代火山巖盆地的1條測(cè)線上的15個(gè)MT測(cè)點(diǎn)(圖4中紅框區(qū)域)實(shí)測(cè)數(shù)據(jù)應(yīng)用本文方法進(jìn)行校正。

    圖4 安徽霍山地區(qū)MT測(cè)點(diǎn)分布圖

    MT數(shù)據(jù)的采集使用加拿大鳳凰公司生產(chǎn)的MTU-5A系統(tǒng),該系統(tǒng)配備MTC-30磁傳感器,實(shí)際測(cè)量中最低頻率為0.00055Hz,最高頻率為320Hz,共40個(gè)頻點(diǎn)。

    由于測(cè)區(qū)存在強(qiáng)烈的電磁干擾,需要對(duì)MT數(shù)據(jù)進(jìn)行精細(xì)的整理和處理。首先對(duì)工區(qū)內(nèi)MT01測(cè)線上的15個(gè)MT測(cè)點(diǎn)都進(jìn)行了時(shí)間域及頻率域數(shù)據(jù)預(yù)處理,對(duì)干擾嚴(yán)重的點(diǎn)進(jìn)行了功率譜挑選、飛點(diǎn)剔除、曲線擬合等處理,然后使用相位張量方法對(duì)這15個(gè)測(cè)點(diǎn)進(jìn)行了維性分析,圖5所示為該測(cè)線相位張量分析圖[33]。

    一般地,當(dāng)二維偏離度小于3°時(shí),表明地下構(gòu)造為二維或準(zhǔn)二維,反之則趨向于三維。分析MT01測(cè)線相位張量圖,發(fā)現(xiàn)該剖面10Hz以下頻段的二維偏離度高于其閾值3°,即表現(xiàn)出明顯的三維性,因而傳統(tǒng)的阻抗分解方法不適用于該剖面。在大地電磁勘探中,由于第四系沉積物的存在,一般表現(xiàn)為一維電性結(jié)構(gòu),理論上其阻抗張量滿足一維地電模型條件下的特征,即TE和TM模式視電阻率曲線的首支重合,實(shí)際上這也是判斷資料是否存在畸變的一個(gè)重要標(biāo)準(zhǔn)[19],故可以通過對(duì)比校正前、后視電阻率曲線的首支是否重合評(píng)價(jià)校正結(jié)果。

    以測(cè)點(diǎn)site01、site02、site07、site15為例分析校正效果。圖6為這4個(gè)測(cè)點(diǎn)經(jīng)本文方法校正前、后視電阻率和相位的對(duì)比??梢姡U暗臏y(cè)點(diǎn)site01和測(cè)點(diǎn)site02的XY和YX模式視電阻率偏離較?。恍U笫字Ы踔睾?,測(cè)點(diǎn)site01視電阻率曲線在10Hz附近分離,測(cè)點(diǎn)site02視電阻率曲線在50Hz附近分離,這與圖5的相位張量分析結(jié)果一致,在相應(yīng)的重合頻段呈現(xiàn)近圓形,分離的頻段相位張量橢圓率增大。測(cè)點(diǎn)site07同樣具有這種特征。測(cè)點(diǎn)site15的原始視電阻率曲線首支分離程度較大,校正后首支重合,這在一定程度上佐證了本文方法的有效性。

    圖5 安徽霍山地區(qū)MT01測(cè)線相位張量分析結(jié)果

    圖6 測(cè)點(diǎn)site01(a)、site02(b)、site07(c)、site15(d)經(jīng)本文方法校正前、后視電阻率(上)、相位(下)曲線對(duì)比

    4 結(jié)論與討論

    傳統(tǒng)阻抗分解算法一般是對(duì)區(qū)域構(gòu)造進(jìn)行一維或者二維假設(shè),然后選擇相應(yīng)的方法進(jìn)行校正,對(duì)于三維區(qū)域構(gòu)造,傳統(tǒng)的阻抗分解方法失效。在前人的研究基礎(chǔ)上,本文提出基于大地電磁振幅—相位張量分解進(jìn)行大地電磁阻抗畸變校正,無(wú)需進(jìn)行區(qū)域構(gòu)造維性假設(shè),對(duì)三維構(gòu)造MT數(shù)據(jù)的畸變校正適用。此外,本文方法無(wú)需將各向異性參數(shù)并入?yún)^(qū)域阻抗,通過改進(jìn)粒子群算法實(shí)現(xiàn)了畸變參數(shù)的可靠求取。理論模型和實(shí)測(cè)數(shù)據(jù)計(jì)算結(jié)果分析表明本文方法能夠有效進(jìn)行大地電磁阻抗張量畸變校正。

    本文算法存在以下不足:畸變校正只針對(duì)單一測(cè)點(diǎn),若實(shí)際觀測(cè)數(shù)據(jù)受噪聲影響嚴(yán)重,則校正結(jié)果不夠理想,有待在后續(xù)研究中改進(jìn),進(jìn)一步討論多站點(diǎn)、其他條件約束下的畸變校正方法。

    猜你喜歡
    張量畸變振幅
    偶數(shù)階張量core逆的性質(zhì)和應(yīng)用
    四元數(shù)張量方程A*NX=B 的通解
    在Lightroom中校正鏡頭與透視畸變
    十大漲跌幅、換手、振幅、資金流向
    擴(kuò)散張量成像MRI 在CO中毒后遲發(fā)腦病中的應(yīng)用
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    輻射誘導(dǎo)染色體畸變的快速FISH方法的建立
    工程中張量概念的思考
    河南科技(2014年19期)2014-02-27 14:15:33
    《癌變·畸變·突變》第六屆編委會(huì)第2次會(huì)議紀(jì)要
    久久精品国产鲁丝片午夜精品| 亚洲婷婷狠狠爱综合网| 国产精品久久久久成人av| 国产淫片久久久久久久久| 国产日韩欧美视频二区| 中文字幕人妻熟人妻熟丝袜美| 亚洲av成人精品一区久久| 色5月婷婷丁香| 黑丝袜美女国产一区| 国内少妇人妻偷人精品xxx网站| 欧美日韩视频精品一区| 汤姆久久久久久久影院中文字幕| 亚洲欧洲国产日韩| 日韩av免费高清视频| 日本av手机在线免费观看| 久久午夜综合久久蜜桃| 热99国产精品久久久久久7| a级毛色黄片| 丝瓜视频免费看黄片| 欧美精品高潮呻吟av久久| a级片在线免费高清观看视频| 色94色欧美一区二区| 成人亚洲欧美一区二区av| 亚洲天堂av无毛| 一级二级三级毛片免费看| 狂野欧美白嫩少妇大欣赏| 亚洲第一av免费看| 久久人妻熟女aⅴ| 国产成人免费观看mmmm| 欧美 亚洲 国产 日韩一| av线在线观看网站| 精品国产露脸久久av麻豆| 视频区图区小说| 最近中文字幕高清免费大全6| 国产日韩欧美视频二区| 人妻制服诱惑在线中文字幕| 亚洲av成人精品一二三区| 日本爱情动作片www.在线观看| 青春草国产在线视频| 国产精品国产av在线观看| 边亲边吃奶的免费视频| a级毛片免费高清观看在线播放| 亚洲欧美中文字幕日韩二区| 精品一品国产午夜福利视频| 天天操日日干夜夜撸| 日韩中文字幕视频在线看片| 人妻人人澡人人爽人人| 日韩精品免费视频一区二区三区 | 99久久综合免费| 亚洲国产成人一精品久久久| 黑丝袜美女国产一区| 菩萨蛮人人尽说江南好唐韦庄| 国产成人a∨麻豆精品| 中文字幕免费在线视频6| 国产探花极品一区二区| 亚洲天堂av无毛| 久久久久久久精品精品| 内地一区二区视频在线| 欧美日韩综合久久久久久| 亚洲美女视频黄频| 美女福利国产在线| 蜜桃久久精品国产亚洲av| 人妻制服诱惑在线中文字幕| 美女国产视频在线观看| 男人爽女人下面视频在线观看| 国产成人精品一,二区| 全区人妻精品视频| 免费av中文字幕在线| 水蜜桃什么品种好| 欧美人与善性xxx| 街头女战士在线观看网站| 亚洲国产精品一区三区| 欧美日韩亚洲高清精品| 波野结衣二区三区在线| 在现免费观看毛片| 在线观看免费视频网站a站| 美女大奶头黄色视频| 亚洲欧洲日产国产| 在线观看www视频免费| 国产精品99久久久久久久久| 亚洲美女黄色视频免费看| 天美传媒精品一区二区| 高清在线视频一区二区三区| 欧美精品一区二区大全| 久久综合国产亚洲精品| 亚洲美女搞黄在线观看| 亚洲欧美一区二区三区黑人 | 国产毛片在线视频| 哪个播放器可以免费观看大片| 久久久久久久精品精品| 久久久国产欧美日韩av| 天天躁夜夜躁狠狠久久av| 激情五月婷婷亚洲| 欧美3d第一页| 免费av不卡在线播放| 亚洲国产精品专区欧美| 中文乱码字字幕精品一区二区三区| 国产欧美日韩精品一区二区| 国产真实伦视频高清在线观看| 午夜免费鲁丝| 极品人妻少妇av视频| 亚洲国产av新网站| 国产精品久久久久久精品电影小说| 秋霞伦理黄片| 成人黄色视频免费在线看| 国产av国产精品国产| 嘟嘟电影网在线观看| 九草在线视频观看| 在线观看三级黄色| 青春草国产在线视频| 久久精品国产亚洲av涩爱| 观看av在线不卡| av卡一久久| 黑人巨大精品欧美一区二区蜜桃 | 91成人精品电影| 不卡视频在线观看欧美| 又大又黄又爽视频免费| 人人妻人人添人人爽欧美一区卜| 乱人伦中国视频| 国产熟女欧美一区二区| 多毛熟女@视频| 在线观看www视频免费| 国产欧美亚洲国产| 啦啦啦在线观看免费高清www| 亚洲欧美中文字幕日韩二区| 亚洲欧美日韩东京热| 69精品国产乱码久久久| 在现免费观看毛片| 久久精品国产亚洲av涩爱| 亚洲美女视频黄频| 乱系列少妇在线播放| 久久av网站| 国产熟女午夜一区二区三区 | 极品教师在线视频| 大陆偷拍与自拍| 日韩中文字幕视频在线看片| 亚洲成人一二三区av| 大香蕉久久网| 久久久国产一区二区| 免费在线观看成人毛片| 免费av不卡在线播放| 亚洲国产欧美在线一区| 啦啦啦啦在线视频资源| 成人亚洲欧美一区二区av| 高清av免费在线| 99热6这里只有精品| 久久鲁丝午夜福利片| 草草在线视频免费看| 婷婷色综合www| 日韩精品免费视频一区二区三区 | 亚洲,欧美,日韩| 久久精品国产自在天天线| 久久午夜福利片| 一级爰片在线观看| 欧美日韩视频精品一区| 老司机亚洲免费影院| 少妇的逼水好多| 国产黄频视频在线观看| 精品久久久精品久久久| 蜜桃在线观看..| 哪个播放器可以免费观看大片| av国产久精品久网站免费入址| 国产免费又黄又爽又色| 婷婷色综合大香蕉| 简卡轻食公司| 国产精品.久久久| 丰满人妻一区二区三区视频av| 国产精品一区www在线观看| 在线精品无人区一区二区三| 又粗又硬又长又爽又黄的视频| 韩国高清视频一区二区三区| 黑丝袜美女国产一区| √禁漫天堂资源中文www| 婷婷色av中文字幕| 日本黄大片高清| 搡老乐熟女国产| 久久久久久人妻| 永久网站在线| av在线播放精品| 丝袜喷水一区| 成人黄色视频免费在线看| 一本久久精品| 国产一区二区三区综合在线观看 | 精品国产乱码久久久久久小说| 丰满少妇做爰视频| 日日摸夜夜添夜夜添av毛片| 亚洲欧美成人综合另类久久久| 国内精品宾馆在线| 精品酒店卫生间| 在线观看三级黄色| 亚洲精品乱码久久久v下载方式| 欧美精品一区二区大全| 国产免费视频播放在线视频| 18+在线观看网站| 日韩精品有码人妻一区| 亚洲精品成人av观看孕妇| 欧美激情极品国产一区二区三区 | av在线观看视频网站免费| 久久国产亚洲av麻豆专区| 在线观看一区二区三区激情| 看免费成人av毛片| 一边亲一边摸免费视频| 纯流量卡能插随身wifi吗| 大话2 男鬼变身卡| 国产精品久久久久久av不卡| 97超碰精品成人国产| 亚洲美女视频黄频| 另类精品久久| 日韩av不卡免费在线播放| 久久韩国三级中文字幕| h视频一区二区三区| 亚洲精品亚洲一区二区| 嫩草影院入口| 人妻一区二区av| 最近最新中文字幕免费大全7| 亚洲欧美日韩另类电影网站| 欧美区成人在线视频| 久久午夜综合久久蜜桃| 在线观看一区二区三区激情| 国产欧美日韩综合在线一区二区 | 尾随美女入室| 国产男人的电影天堂91| 亚洲精品乱码久久久v下载方式| 两个人的视频大全免费| 好男人视频免费观看在线| 国产成人午夜福利电影在线观看| 各种免费的搞黄视频| 啦啦啦啦在线视频资源| 欧美xxⅹ黑人| 卡戴珊不雅视频在线播放| 男女无遮挡免费网站观看| 一区在线观看完整版| 久久久久网色| 色网站视频免费| 一级毛片久久久久久久久女| 久久久精品免费免费高清| 99热国产这里只有精品6| 水蜜桃什么品种好| 欧美xxⅹ黑人| 国产免费一区二区三区四区乱码| 亚洲国产精品999| 中国三级夫妇交换| 久久久久久久大尺度免费视频| 男的添女的下面高潮视频| 亚洲成人手机| 精品人妻熟女毛片av久久网站| xxx大片免费视频| 在线亚洲精品国产二区图片欧美 | 欧美少妇被猛烈插入视频| 又爽又黄a免费视频| 免费高清在线观看视频在线观看| 亚洲国产色片| 精品熟女少妇av免费看| 久久6这里有精品| 久久这里有精品视频免费| 成人国产麻豆网| 人妻系列 视频| 国产精品久久久久久久电影| 卡戴珊不雅视频在线播放| 国产又色又爽无遮挡免| 九色成人免费人妻av| 老熟女久久久| 婷婷色麻豆天堂久久| 99久久精品国产国产毛片| 99久久精品一区二区三区| 国产女主播在线喷水免费视频网站| 91在线精品国自产拍蜜月| 中文字幕免费在线视频6| 国产精品一区二区在线观看99| 成人影院久久| 久久婷婷青草| 午夜激情久久久久久久| 久久久a久久爽久久v久久| 哪个播放器可以免费观看大片| 亚洲av免费高清在线观看| 国产精品一区二区三区四区免费观看| 久久久久人妻精品一区果冻| 中文字幕久久专区| 成人国产av品久久久| videos熟女内射| av免费观看日本| 亚洲人与动物交配视频| 成人毛片a级毛片在线播放| 亚洲成人手机| 三级经典国产精品| 最近中文字幕高清免费大全6| 国产欧美日韩精品一区二区| 又粗又硬又长又爽又黄的视频| 久久精品国产鲁丝片午夜精品| 成人黄色视频免费在线看| 51国产日韩欧美| 国产熟女欧美一区二区| 免费观看无遮挡的男女| 日韩,欧美,国产一区二区三区| 国国产精品蜜臀av免费| 国产在线一区二区三区精| 亚洲欧美精品专区久久| 韩国av在线不卡| 高清黄色对白视频在线免费看 | 日韩欧美 国产精品| 午夜老司机福利剧场| 九草在线视频观看| 国产av国产精品国产| 国产伦理片在线播放av一区| 大片电影免费在线观看免费| 丁香六月天网| 99久久中文字幕三级久久日本| 在线亚洲精品国产二区图片欧美 | 亚洲精品久久午夜乱码| a级片在线免费高清观看视频| 国产男女内射视频| 大又大粗又爽又黄少妇毛片口| 99久久精品热视频| 高清毛片免费看| 哪个播放器可以免费观看大片| 秋霞在线观看毛片| 天天躁夜夜躁狠狠久久av| 在现免费观看毛片| 全区人妻精品视频| 午夜91福利影院| 看免费成人av毛片| 黄色毛片三级朝国网站 | 插阴视频在线观看视频| 女性生殖器流出的白浆| 少妇猛男粗大的猛烈进出视频| 精品少妇黑人巨大在线播放| 久久ye,这里只有精品| 日韩视频在线欧美| 亚洲内射少妇av| 亚洲第一av免费看| 中文精品一卡2卡3卡4更新| 老女人水多毛片| 国内揄拍国产精品人妻在线| 黑人巨大精品欧美一区二区蜜桃 | 汤姆久久久久久久影院中文字幕| 99热6这里只有精品| 久久热精品热| 欧美精品亚洲一区二区| 久久久亚洲精品成人影院| 久久久久网色| 国产永久视频网站| 水蜜桃什么品种好| 亚洲av中文av极速乱| 简卡轻食公司| 欧美精品人与动牲交sv欧美| 最近2019中文字幕mv第一页| 日韩人妻高清精品专区| 亚洲欧美一区二区三区国产| 亚洲av国产av综合av卡| 秋霞在线观看毛片| 国产精品成人在线| 另类精品久久| 最黄视频免费看| 亚洲成人一二三区av| 一区二区三区精品91| 97精品久久久久久久久久精品| 高清黄色对白视频在线免费看 | 国产精品偷伦视频观看了| 在线免费观看不下载黄p国产| 久久热精品热| 国产日韩欧美视频二区| 如日韩欧美国产精品一区二区三区 | 久久久国产一区二区| 亚洲第一av免费看| 一级a做视频免费观看| 成人特级av手机在线观看| 国产淫语在线视频| 多毛熟女@视频| a级一级毛片免费在线观看| 国产精品三级大全| 人妻 亚洲 视频| 亚洲内射少妇av| 精品一品国产午夜福利视频| 日韩强制内射视频| av视频免费观看在线观看| 欧美国产精品一级二级三级 | 91久久精品国产一区二区三区| 午夜免费观看性视频| 国产亚洲午夜精品一区二区久久| 久久久欧美国产精品| 欧美少妇被猛烈插入视频| 日韩一区二区视频免费看| √禁漫天堂资源中文www| 国产伦精品一区二区三区四那| 国产精品欧美亚洲77777| 久久这里有精品视频免费| 夫妻午夜视频| 久久亚洲国产成人精品v| 2018国产大陆天天弄谢| 国产精品人妻久久久久久| 91久久精品国产一区二区三区| 午夜免费男女啪啪视频观看| 亚洲情色 制服丝袜| 亚洲无线观看免费| 久久99热这里只频精品6学生| 麻豆成人av视频| 免费大片18禁| 九草在线视频观看| 亚洲国产最新在线播放| 高清黄色对白视频在线免费看 | 欧美日韩精品成人综合77777| 欧美日韩亚洲高清精品| 少妇人妻久久综合中文| 日本91视频免费播放| 最近最新中文字幕免费大全7| 国产精品麻豆人妻色哟哟久久| 免费不卡的大黄色大毛片视频在线观看| 久久久久久久亚洲中文字幕| 日本猛色少妇xxxxx猛交久久| 赤兔流量卡办理| 亚洲av电影在线观看一区二区三区| 午夜福利,免费看| 亚洲精品久久午夜乱码| 老司机影院毛片| 最近手机中文字幕大全| 三级国产精品欧美在线观看| 国产成人aa在线观看| 亚洲精品亚洲一区二区| 国产一级毛片在线| 亚洲欧洲国产日韩| 一级爰片在线观看| 精品一区在线观看国产| 99re6热这里在线精品视频| 黄色毛片三级朝国网站 | 国产黄色视频一区二区在线观看| 2022亚洲国产成人精品| 亚洲精品456在线播放app| 欧美bdsm另类| 高清在线视频一区二区三区| 久久人人爽人人爽人人片va| 蜜桃在线观看..| 黄色怎么调成土黄色| 国产成人精品福利久久| 高清欧美精品videossex| 亚洲精品aⅴ在线观看| 亚洲人与动物交配视频| 免费观看a级毛片全部| 午夜精品国产一区二区电影| 美女国产视频在线观看| 久久婷婷青草| 观看美女的网站| 欧美激情极品国产一区二区三区 | 日本wwww免费看| 乱码一卡2卡4卡精品| 日韩中文字幕视频在线看片| 久久99精品国语久久久| 特大巨黑吊av在线直播| 我的女老师完整版在线观看| 美女中出高潮动态图| 精品一区在线观看国产| 久久久久国产网址| 午夜福利影视在线免费观看| 久久热精品热| 免费看日本二区| 日韩制服骚丝袜av| 精品久久久噜噜| 成人特级av手机在线观看| 亚洲,欧美,日韩| 国产成人a∨麻豆精品| 亚洲电影在线观看av| 久久久久精品性色| 亚洲美女黄色视频免费看| 久久久久久久大尺度免费视频| 久久久久久久久大av| 成年女人在线观看亚洲视频| 免费大片18禁| 最近中文字幕2019免费版| 热99国产精品久久久久久7| 亚洲精品视频女| 极品教师在线视频| 各种免费的搞黄视频| 亚洲人成网站在线观看播放| 高清不卡的av网站| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩av久久| 一区二区三区精品91| 自线自在国产av| 亚洲av福利一区| 人人妻人人看人人澡| 啦啦啦视频在线资源免费观看| 夜夜看夜夜爽夜夜摸| 久久精品夜色国产| 免费播放大片免费观看视频在线观看| 成人综合一区亚洲| 欧美性感艳星| 久久久久久久久久人人人人人人| 成年人午夜在线观看视频| 国产黄频视频在线观看| av黄色大香蕉| √禁漫天堂资源中文www| 国产探花极品一区二区| 我要看黄色一级片免费的| 这个男人来自地球电影免费观看 | kizo精华| 美女内射精品一级片tv| 国产永久视频网站| 亚洲精品国产av蜜桃| 中国三级夫妇交换| 熟女av电影| 六月丁香七月| 日韩av不卡免费在线播放| 中文字幕亚洲精品专区| 午夜激情久久久久久久| 国产在线视频一区二区| 一本一本综合久久| 久久国产乱子免费精品| 午夜影院在线不卡| 美女xxoo啪啪120秒动态图| 久久久久久久久久人人人人人人| 97精品久久久久久久久久精品| 久久精品国产自在天天线| 国产免费又黄又爽又色| 久久久久久久精品精品| 亚洲第一区二区三区不卡| 免费在线观看成人毛片| 美女国产视频在线观看| 日韩在线高清观看一区二区三区| 男男h啪啪无遮挡| 久久99一区二区三区| 成人二区视频| 大码成人一级视频| 十八禁网站网址无遮挡 | 国产精品国产三级专区第一集| 久久精品国产亚洲av天美| 国产精品欧美亚洲77777| 大香蕉久久网| 欧美精品亚洲一区二区| 亚洲欧美成人综合另类久久久| 一级毛片aaaaaa免费看小| av国产精品久久久久影院| 亚洲国产精品一区二区三区在线| 久久99热6这里只有精品| 一级毛片 在线播放| 另类亚洲欧美激情| 99九九在线精品视频 | 香蕉精品网在线| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品999| 亚洲综合精品二区| 日本与韩国留学比较| 精品亚洲乱码少妇综合久久| 丰满乱子伦码专区| 老司机影院毛片| 韩国av在线不卡| 久久久久久伊人网av| a级片在线免费高清观看视频| 午夜老司机福利剧场| 国产亚洲5aaaaa淫片| 久久精品国产a三级三级三级| 精品一区在线观看国产| 大陆偷拍与自拍| 建设人人有责人人尽责人人享有的| 少妇人妻久久综合中文| 日韩成人伦理影院| 国产亚洲精品久久久com| 欧美日韩视频精品一区| 男女无遮挡免费网站观看| 日本av手机在线免费观看| 国产精品一区二区性色av| 三级国产精品片| 亚洲精品乱久久久久久| 日韩在线高清观看一区二区三区| 日韩免费高清中文字幕av| 欧美日韩在线观看h| 在线观看人妻少妇| 精品少妇久久久久久888优播| 国产成人精品无人区| 欧美日韩亚洲高清精品| 欧美国产精品一级二级三级 | 亚洲精品日本国产第一区| 26uuu在线亚洲综合色| 日韩视频在线欧美| 国精品久久久久久国模美| 午夜福利网站1000一区二区三区| 亚洲人与动物交配视频| 日韩制服骚丝袜av| 国产毛片在线视频| 国产精品嫩草影院av在线观看| 久久精品国产自在天天线| 欧美变态另类bdsm刘玥| 亚洲成人手机| 日本与韩国留学比较| 国产色婷婷99| 少妇猛男粗大的猛烈进出视频| 九九在线视频观看精品| 狂野欧美白嫩少妇大欣赏| 精品久久久久久电影网| 日韩一区二区视频免费看| 好男人视频免费观看在线| 久久99蜜桃精品久久| 中文字幕久久专区| 亚洲成人手机| 国产精品福利在线免费观看| av在线老鸭窝| 国产黄片视频在线免费观看| 日韩人妻高清精品专区| 中文字幕免费在线视频6| 欧美少妇被猛烈插入视频| 国产男女超爽视频在线观看| a级毛色黄片| 又粗又硬又长又爽又黄的视频| 国产黄片美女视频| 水蜜桃什么品种好| 日本vs欧美在线观看视频 | 毛片一级片免费看久久久久| 男女边摸边吃奶| 亚洲图色成人| 在线观看美女被高潮喷水网站| 亚洲精品aⅴ在线观看| 欧美日韩亚洲高清精品| 亚洲精品国产色婷婷电影| 欧美xxⅹ黑人| 久久精品国产亚洲av天美| 天堂8中文在线网| 国内少妇人妻偷人精品xxx网站|