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

    多源頻率域地空系統(tǒng)三維電磁響應分析

    2021-04-07 01:53:44張繼鋒劉寄仁馮兵董巖鄭一安
    地球物理學報 2021年4期
    關鍵詞:多源測線全域

    張繼鋒, 劉寄仁, 馮兵*, 董巖, 鄭一安

    1 長安大學地質工程與測繪學院, 西安 710054 2 中南大學地球科學與信息物理學院, 長沙 410083

    0 引言

    傳統(tǒng)的地面電磁法(包括可控源音頻大地電磁法、瞬變電磁法等)是資源探查的重要手段,特別是在礦產資源及環(huán)境工程領域發(fā)揮了重要的作用(底青云等, 2019; Xue et al.,2020).隨著“十三五”國家“四深”發(fā)展戰(zhàn)略的提出,地球深部礦產資源的精細化探測提上日程,迫切需要發(fā)展新的技術方法.多道瞬變電磁法、廣域電磁法、短偏移距瞬變電磁以及地空電磁法等一批新的探測技術應運而生(陳衛(wèi)營等,2016;薛國強等,2020;Ji et al., 2016).地空電磁法融合了地面電磁法大功率發(fā)射和無人機空中快速非接觸式采集的雙重優(yōu)點,能夠進入地形復雜區(qū)域開展大深度資源勘探,近年來成為地球物理電磁法研究的熱點(Allah et al., 2013, 2014; 李貅等, 2015;張瑩瑩和李貅,2017;Xue et al., 2018;Wu et al., 2019).該方法與無人機飛行平臺相結合,有利于在中小區(qū)域開展大深度地下結構精細探測.相對于地面和航空電磁法,更加經濟、安全和便捷,因此具有廣闊的市場前景和應用價值.

    地空電磁法的發(fā)展最早可追溯到20世紀70年代,加拿大TURAIR系統(tǒng)研制成功并用于安大略湖北部金屬硫化礦的探測(Pemberton et al., 1970; Bosschart and Seigel, 1972).該系統(tǒng)采用地面大回線作為發(fā)射源,在近空區(qū)域測量特定頻率下兩個水平線圈或正交線圈的磁場比值以及相位差,發(fā)射功率可達15 kW,探測深度可達到200 m.90年代后期,一系列新的地空儀器設備相繼出現,包括澳大利亞的FLAIRTEM系統(tǒng),加拿大的TerraAir系統(tǒng)和日本GREATEM系統(tǒng)f621(Elliott,1998; Mogi et al., 1998 ).這幾種設備采集方式都是觀測一次場關斷后瞬變電磁響應,解釋仍然采用地面回線源的解釋方法,只是前兩種系統(tǒng)仍然采用大回線發(fā)射源激勵,GREATEM系統(tǒng)采用地面電性源激勵方式.進入21世紀,Smith等(2001)對地面、地空和航空瞬變電磁系統(tǒng)進行了對比研究,系統(tǒng)分析了三種方法的噪聲水平和探測效果.結果表明:地空電磁系統(tǒng)的探測深度和數據質量均優(yōu)于航空電磁系統(tǒng),僅次于地面電磁系統(tǒng).此后,電性源地空系統(tǒng)GREATEM先后用于火山區(qū)和海岸帶地電結構探測,深度可達800 m,成為一種重要的深部電磁勘探方法(Mogi et al., 2009; Ito et al., 2011, 2014).德國BGR研究所(Nittinger et al.,2017)也推出了一套頻率域地空電磁探測系統(tǒng),并成功應用于金屬礦探測中.目前地空電磁系統(tǒng)的研究主要以單輻射源為主,大量的研究集中在儀器研制和信號的去噪方面(嵇艷鞠等,2013;Wang et al., 2013;Qu et al., 2017; Li et al., 2017; 劉富波等,2017;Yuan et al., 2019).由于地空電磁法接收信號干擾較大,單源發(fā)射功率有限,探測深度和測量范圍都會受到限制,因此,迫切需要發(fā)展地空多源電磁勘探系統(tǒng)和方法.

    多源勘探概念最早應用在直流電阻率法中,Bibby(1986)提出多源測量模式并引入電阻率張量,該方法克服了視電阻率對目標體位置和方向的依賴,減少了單源測量結果的“假異?!爆F象,證明多源測量能夠獲得更可靠的信息.Caldwell和Bibby(1998)在長偏移距瞬變電磁(LOTEM)勘探中提出了多源瞬時電阻率張量,保留了直流電阻率張量的許多有用性質,特別是在地下三維電阻率分布情況下,許多解釋技術可以應用于時間域電磁數據.90年代初,阻抗張量和傾子矢量被引入到地面可控源電磁測量中(Li and Pedersen,1991),該方法一般采用兩個相互垂直的激勵源,分別從不同的方向激勵,兩個源單獨工作,能夠獲取更加豐富的地電信息,解釋結果更加可信,后來進一步得到發(fā)展(Caldwell et al., 2002).席振銖等(2011)研究了正交磁偶源電磁場分布規(guī)律,實現了可控源高頻大地電磁張量測量;王顯祥等(2014)等研究了“L”型源的各分量表達式,設計了一種新的信號發(fā)射模式,實現了全張角范圍內測量,顯示了張量可控源方法的優(yōu)勢.

    近年來,多源地空電磁系統(tǒng)的研究也逐漸受到了重視,李貅等(2015)提出了多輻射場源的全域視電阻率定義方式,不僅提供了一種定性直觀的解釋方法,也為三維多源地空瞬變電磁的發(fā)展提供了理論依據.Zhou等(2018a,b)年提出頻率域多源電磁探測方法,主要從儀器的研制、激勵信號的同步以及多源之間的相互影響等方面進行了研究,表明多源地空電磁勘探具有廣闊的發(fā)展前景.因此,需要進一步從三維角度對多源電磁響應進行模擬和分析.目前地球物理電磁法三維數值模擬主要集中在單激勵源或被動源三維模擬方面(李建慧等, 2016; Ren et al., 2013, 2014; 殷長春等, 2017),多源頻率域地空三維電磁響應特征研究比較少,特別是全域視電阻率分布特征和規(guī)律.

    本文采用多源三維有限元數值模擬方法,從多源輻射角度出發(fā),設置了雙源、三源以及四源輻射方式,建立了低阻體、高阻體模型以及兩個低阻體模型,分別從總場、二次場以及視電阻率參數等方面對不同發(fā)射源模式進行分析,總結每種源的電磁響應和全域視電阻率特征,為實際選擇多源輻射方式進行地空電磁勘探提供理論基礎.

    1 多源頻率域地空系統(tǒng)

    地空頻率域電磁探測方法是采用地面布設發(fā)射源,空中測量垂直磁場響應的工作模式.該方法由吉林大學教授林君提出,其儀器(GAFEM)已經研制成功,采用偽隨機發(fā)射波形,在一些老礦區(qū)進行了野外試驗,取得了良好效果.單輻射源能量有限,信號強度弱,加之受運動噪聲的影響,信噪比和探測深度都會受到一定的影響.多場源激勵可改變空間磁場分布,提高信號強度,擴大觀測區(qū)域,提高信噪比.通過多場源不同角度和不同方向激勵,可加強發(fā)射源與目標體的耦合程度,改善對目標體的分辨率(圖1).頻率域地空電磁法測量空中磁場垂直分量,不再測量電場分量,場的特征和解釋方法完全不同,地面的處理解釋方法不再適用.因此,必須基于空中磁場分量,研究三維地電模型的電磁場空間分布特征和新的解釋方法.

    圖1 頻率域多源地空電磁采集方式示意圖Fig.1 Schematic diagram showing data acquisition for multiple-source frequency domain grounded-airborne electromagnetics

    2 多源頻率域地空有限元分析

    2.1 基于總場的公式

    在準靜態(tài)極限條件下的電場E和磁場H所滿足的頻率域方程為(湯井田和何繼善,2005):

    (1)

    (2)

    其中假設電磁波的諧變因子為eiω t,式中μ為磁導率,本文μ選取真空磁導率μ0,σ為電導率,ω為圓頻率,Js為傳導電流密度.

    由此便可導出電場雙旋度方程:

    (3)

    (4)

    其中,V為研究區(qū)域的體積,在進行節(jié)點有限元求解時,該方程不滿足散度條件,產生的偽解直接影響計算精度.為了提高計算精度,可在方程中增加一個罰項來削弱偽解產生的影響(金建銘,1998):

    (5)

    2.2 有限元分析

    本文采用非結構化網格剖分方法,利用TetGen生成高質量的四面體網格,將源和異常體附近場值變化劇烈的區(qū)域進行網格加密,網格邊界大于10倍的趨膚深度,以保證場值在邊界上已衰減為零,如圖2所示.

    圖2 非結構化網格剖分Fig.2 Unstructured grid

    圖3 四面體單元Fig.3 Tetrahedron element

    (6)

    若將整個研究區(qū)域剖分為n個單元,離散化后的變分方程為:

    (7)

    那么地空多源正演問題轉化為求解該變分方程的極值問題,該方程可看成多元函數,多元函數取極值應當滿足下列方程(徐世浙,1994):

    (8)

    其中i為節(jié)點編號,e為單元編號.

    由此可得到所有單元所滿足的矩陣方程:

    (9)

    按照單元節(jié)點的編號,可將每個單元生成的子矩陣組合成總體剛度矩陣,得到最終的離散矩陣方程:

    (C+iωM)·E=-iωb.

    (10)

    由于產生的剛度矩陣是大型的對稱稀疏矩陣,直接存儲將會浪費大量存儲空間,因此本文采用雙向鏈表組裝矩陣,按CSR方式存儲.由于多頻可控源正演模擬的頻點之間是相互獨立的,解方程將會浪費大量的時間,嚴重影響運算效率,而經前人研究發(fā)現,可以通過模型降階方法,將一個復雜的系統(tǒng)簡化,只要經過一次模型降階,便可實現多頻點的計算,提高運算速率(蔣耀林,2010;B?rner et al., 2008;周建美等,2018).

    2.3 模型降階方法

    對于矩陣方程(10),我們可以做如下變形:

    (M-1C+iωI)·E=-iωM-1b,

    (11)

    令A=M-1C,X=M-1b,gω(A)=-iω(A+iωI)-1,解得E(ω)=-iω(A+iωI)-1X=gω(A)X.那么求解E(ω)的關鍵便是求解gω(A),因此我們采用Krylov子空間(張繼鋒等,2009, 2020)投影算法,可降低矩陣運算的維數,大大提高運算效率.

    模型降階算法實際上就是尋找一個有效逼近rm(A),使得E≈rm(A)成立,求解rm(A)等價于求解矩陣A在m維Krylov子空間Qm+1=qm(A)-1κm+1(A,X)上的投影(周建美等,2018).由于有限元離散矩陣的特點,矩陣M并非單位矩陣,為了避免因傳統(tǒng)正交算法造成運算的復雜性,我們引入一種M正交運算(蔣耀林,2010; B?rner, et al., 2015),進而得到列正交向量Vm+1,下面給出模型降階后的最終表達式:

    E=‖X‖MVm+1gω(Am+1)e1,

    (12)

    2.4 垂直磁場的求取

    求解方程(10)可得到整個區(qū)域所有節(jié)點的電場三分量值,而在任意單元e中,任意點(x0,y0,z0)處的電場大小可通過插值基函數求得:

    (13)

    在地空系統(tǒng)中,電偶源在地面發(fā)射電流,線圈在空中接收磁場,因此本文通過研究z方向的磁場值Hz來研究大地的電性特征.由式(1)可得:

    (14)

    其中相應場在某一點的導數值可通過式(13)代入式(14)求得.在實際數據處理和資料解釋中,為了能夠直觀的反映地下介質的電性特征,我們需要將磁場轉化為視電阻率.張瑩瑩等(2015)定義了時間域瞬變電磁場的全域視電阻率,我們將其引入到頻率域,利用泰勒級數展開,求得視電阻率的近似表達式,然后通過迭代算法得到全域視電阻率.

    3 多源三維電磁響應分析

    3.1 有限元解和解析解的比較

    圖4a、b分別給出了有限元解和解析解的垂直磁場分量以及相對誤差.可以看出,三維有限元計算結果同一維正演結果吻合的很好,磁場最大誤差為1.8%左右.

    為了進一步驗證程序的正確性,我們設計了一個兩層模型,第一層電阻率為200 Ωm,厚度為100 m,第二層電阻率為50 Ωm.第一個源沿x軸正方向放置,電偶源的中心位置為(0 m,-4000 m,0 m),第二個源沿x軸負方向放置,電偶源的中心坐標為(0 m,2000 m,0 m),設置測點位置為(1600 m,-2000 m,-100 m).如圖5所示,磁場相對誤差最大不超過3%,說明基于模型降階的有限元計算結果滿足精度要求.

    3.2 三維電磁響應分析

    3.2.1 低阻體模型

    (1)單源和多源比較

    為了說明頻率域多源地空電磁方法的優(yōu)勢,設計一個低阻體模型,針對垂直磁場分量,比較幾種多源組合模式的電磁場分布.低阻體模型如圖6所示,異常體電阻率為10 Ωm,大小為500 m×500 m×300 m,異常體中心坐標為(0 m,0 m,250 m), 圍巖電阻率為100 Ωm.圖6c為非結構化網格剖分圖.為便于對比多源組合方式,假設所有的電性源偶極矩為1000 Am,飛行高度30 m.

    圖4 均勻半空間3D-Forward、1D-Forward的磁場Hz及對應的誤差曲線(a) 磁場Hz; (b) 相對誤差曲線.Fig.4 Magnetic field Hz and relative error for homogenous half space from 3D- and 1D-forward modeling(a) Magnetic field Hz; (b) Relative error.

    圖5 兩層模型3D-Forward與1D-Forward的磁場Hz及對應的誤差曲線(a) 磁場Hz; (b) 相對誤差.Fig.5 Magnetic field Hz and relative error of two-layer model from 3D- and 1D-forward modeling(a) Magnetic field Hz; (b) Relative error.

    圖7為多個源的位置示意圖,我們分析幾種組合方式:①單源,①②組合,①③組合,①②③組合,①②③④組合,①⑤⑥組合.各個源的中心坐標已在圖7中給出.

    我們首先給出256 Hz各種組合源作用下的均勻半空間電磁響應.圖8為各種組合源作用下的垂直磁場切片圖,可以看出,靠近源的一側磁場強,而遠離源的一側磁場逐漸衰減,同時我們也注意到,這幾種組合源都將使得背景場值增大,信號強度增強,這是多源場之間的相互疊加作用導致的.圖9給出了中心測線全域視電阻率測深曲線及相對誤差,可以看出,全域視電阻率基本接近于背景電阻率100 Ωm,誤差不超過2.5%,完全滿足計算精度要求,也說明本文全域視電阻率算法的正確性.

    圖10給出了單源和多源作用下,含異常體的總場幅值分布切片圖,圖11為相應的二次場分布圖,圖12是單源和多源情況下異常幅度的分布圖,圖13給出了單源及多源作用下的全域視電阻率分布圖.從中可以看出,若在單源作用下測量垂直磁場時,在靠近源的一側和遠離源的一側由于在界面上電阻率突變,導致界面產生積累電荷,磁場響應增強,異常體的中心與異常響應的極值并不對應,在靠近源的一側邊界將會使得場值增大,而在遠離源的一側邊界將會使場值減小,如圖10a和圖11a所示,而磁場增大將會使得全域視電阻率增大,因而在靠近源的一側邊界處出現出高阻極大值,而在遠離源的一側邊界處出現了低阻極小值,如圖13a所示.可見,對于三維地電模型的電磁響應分布比較復雜,異常體的位置和磁場及全域視電阻率形態(tài)并不一定吻合,需要進行仔細分辨.磁場分布不僅與異常體的形態(tài)有關,也與異常體與源的相對位置有很大關系.

    圖6 低阻體模型(a) 剖面圖; (b) 平面圖; (c) 非結構化網格剖分圖.紅色代表低阻體,綠色代表背景.Fig.6 Model of a conductive body(a) Cross-section; (b) Plane view; (c) Unstructured grid.Red shows the conductive body and green shows background.

    分析①②組合源和①②③組合源的電磁響應,發(fā)現總場和二次場的幅度都在一定程度上有所增強,如圖10b、圖10d、圖11b、圖11d所示.這也是多源的優(yōu)勢之一,可增強觀測區(qū)域的信號強度.為了更好的反映異常響應分布,需要分析其異常響應的相對幅值,如圖12b、d所示,我們發(fā)現在邊緣處的極大極小值的相對幅度在減小,這是不同方向上的源產生的場在邊界區(qū)域進行疊加的結果,導致異常場的分布形態(tài)發(fā)生了變化,而圖13b、d視電阻率分布圖也說明了這一點,在邊緣處的電阻率極大值在減小,而極小值在增大.圖10c、圖11c是①③組合的總場幅值和二次場大小,這與單源情況下相比,也提高了信號強度,而從相對幅度分布圖12c及全域視電阻率圖13c來看,在異常體的中心出現了低阻圈閉,在邊緣部位出現了高阻響應.這種源的組合方式在異常形態(tài)分布上優(yōu)于前面分析的三種方式,能更好的反映異常的中心位置,削弱了異常體邊界處的假極值,但是在源的方向上,出現了拉長的低阻,而在垂直源的方向上出現了高阻.

    圖8 256 Hz單源和多源的垂直磁場分布圖(a) ①; (b) ①②; (c) ①③; (d) ①②③; (e) ①②③④; (f) ①⑤⑥.Fig.8 Vertical magnetic field of single and multiple sources at 256 Hz

    圖9 中心測點全域視電阻率及相對誤差(a) 全域視電阻率; (b) 相對誤差.Fig.9 Full-domain apparent resistivity and relative error curves of central measuring point(a) Full-domain apparent resistivity; (b) Relative error.

    圖10 256 Hz單源和多源的總場幅值分布圖(a) ①; (b) ①②; (c) ①③; (d) ①②③; (e) ①②③④; (f) ①⑤⑥.Fig.10 Total magnetic field Hs at 256 Hz of single and multiple sources

    圖11 256 Hz單源和多源的二次場響應分布圖(a) ①; (b) ①②; (c) ①③; (d) ①②③; (e) ①②③④; (f) ①⑤⑥.Fig.11 Secondary magnetic field at 256 Hz of single and multiple sources

    圖12 256 Hz單源和多源的異常相對幅度分布圖(a) ①; (b) ①②; (c) ①③; (d) ①②③; (e) ①②③④; (f) ①⑤⑥.Fig.12 Relative percentage of anomaly at 256 Hz of single and multiple sources

    圖13 256 Hz單源和多源的全域視電阻率分布圖(a) ①; (b) ①②; (c) ①③; (d) ①②③; (e) ①②③④; (f) ①⑤⑥.Fig.13 Full-domain apparent resistivity at 256 Hz of single and multiple sources

    圖14 單源和多源的主軸剖面的廣域視電阻率圖(a)、(g) ①; (b)、(h) ①②; (c)、(i) ①③; (d)、(j) ①②③; (e)、(k) ①②③④; (f)、(l) ①⑤⑥.Fig.14 Full-domain apparent resistivity of single-source and multiple sources spindle sections

    圖15 256 Hz單源和多源全域視電阻率、總場幅值和二次場響應(a)—(c) ①; (d)—(f) ①②③④; (g)—(i) ①⑤⑥.Fig.15 Full-domain apparent resistivity, total field and secondary filed responses at 256 Hz of single and multiple sources

    圖16 單源和多源的主軸剖面的廣域視電阻率圖(a)、(b) ①; (c)、(d) ①②③④; (e)、(f) ①⑤⑥.Fig.16 Full-domain apparent resistivity of spindle sections of single source and multiple sources

    從①③組合方式來看,當源關于異常體的中心出現對稱分布時,對異常體的分辨較好,因而本文還設計了①②③④組合和①⑤⑥組合,而①⑤⑥組合三個源相互呈120°對稱,從總場幅值和二次場分布可看出,在測網中心出現了低阻圈閉,且①②③④組合的總場極小值比其余組合方式的極小值更大,說明①②③④組合在增強信號強度方面是最強的,而在二次場分布方面,由于場源和異常體的對稱性,這兩種方式基本消除了異常體邊界的極值,突出了異常中心的電性特征,而異常相對幅度和全域視電阻率分布更好的說明了這一點,在沒有異常體的區(qū)域,受到異常體的影響很小,幾乎接近于背景值,而在異常體中心形成了很明顯的低阻異常圈閉,因而這兩種方式都能夠提高磁場的信號強度,同時使全域視電阻率分布形態(tài)更接近于真實的地電分布.若不考慮其他因素,①②③④組合是相對較優(yōu)的多源組合方式,若考慮野外施工效率,也可選擇①⑤⑥多源組合方式.

    接下來我們分析每個組合方式的主測線剖面圖,包括x=0測線和y=0測線,如圖14所示,圖14a—f表示單源和多源組合的x=0測線剖面,圖14g—l表示單源和多源組合的y=0測線剖面.從中可看出,單源激勵下,兩條測線僅平行于源的測線出現了低阻圈閉,另一條垂直測線靠近源的一側呈現高阻,遠離源的一側呈現低阻,這和上述切片圖是一致的(圖14a、 g).組合源①②的兩條測線均沒有在異常體位置出現低阻圈閉.對于①③組合,在切片圖中已展現出較寬的低阻圈閉,而兩條中心測線的剖面圖也能較好的體現異常體的位置,但平行源的測線低阻圈閉橫向范圍寬(圖14i),而其x=0測線在邊界處呈現了高阻響應,中心低阻圈閉較小(圖14c).①②③組合的切片圖形與單源形式是類似的,因而其剖面圖的形態(tài)也類似,只是出現低阻圈閉的測線相反.①②③④組合和①⑤⑥組合的兩條中心測線都出現了明顯的低阻圈閉,而且與異常體的位置對應的較好,周圍圍巖呈現出背景電阻率,這兩種組合更容易識別異常體的分布.

    3.2.2 高阻體模型

    同樣,我們也分析了高阻體的地空三維電磁響應分布,圖15a—c為單源作用下的高阻異常體響應,可以看出在靠近源的一側表現出了低阻性質,而在遠離源的一側表現出了高阻性質,這與低阻體的響應是相反的.從①⑤⑥組合和①②③④組合的響應來看,異常響應呈現高阻圈閉,消除了邊緣的假極值.而圖16給出了相應的主軸剖面,從中看出,單源情況下的x=0測線剖面和y=0測線剖面的形態(tài)和低阻體單源情況下的剖面形態(tài)是一致的,只是高阻和低阻發(fā)生了倒轉,如圖16a、b所示,也可以看出,y=0測線的高阻圈閉形態(tài)出現誤差,這是因為地空響應對高阻體不敏感造成的,產生的異常響應和圍巖的電阻率較為接近,因而對高阻體的靈敏度不如低阻體.從①②③④組合和①⑤⑥組合的中心測線剖面可以看出,每種組合的兩條中心測線都出現了高阻圈閉,雖然異常幅度較低,但是能夠很好的反映異常體的位置.

    3.2.3 兩個低阻體模型

    建立兩個大小相同的低阻體,大小均為400 m×500 m×300 m,電阻率都為10 Ωm,關于y軸對稱,左邊異常體中心為(-400 m,0 m,250 m),右側異常體中心為(400 m,0 m,250 m),背景電阻率為100 Ωm,圖17a為模型剖面圖,圖17b為模型平面圖,圖17c為源的位置示意圖,圖17d為非結構化網格剖分圖.

    共設計了三種場源的組合方式,在單源①的作用下,靠近源的一側呈現兩個了高阻,而遠離源的一側呈現了兩個低阻圈閉,且對稱分布,如圖18a—c所示,這與前文的結果是吻合的.而①⑤⑥組合和①②③④組合的響應也呈對稱分布,但是較單源而言,低阻異常的位置同實際位置更加接近,兩側呈對稱分布的高阻異常也已大致趨于背景值,而異常響應在兩個低阻體相互靠近的一側邊緣出現了低阻極小值,雖然和異常體并沒有完全吻合,但比單源作用下的響應形態(tài)更簡單,更能反映異常體的電性分布特征.

    為了更好的比較,我們給出了不同頻率下兩條中心測線(x=0和y=0)的垂直磁場變化曲線,如在圖19, 可以發(fā)現,在y=0測線上,多源激發(fā)下垂直磁場都出現了兩個低值凹陷,但明顯四源和三源情況下的凹陷位置更接近于異常體的中心,而單源激發(fā)下,在中間位置磁場響應卻增大了,與實際模型電性不一致(圖19e).從剖面圖整體上可以看出,多源的信號強度明顯強于單源的信號,這對于實際測量非常有利,在相同噪聲背景下能夠提高信噪比,如圖19a、c、e所示.而對于x=0測線,四源和三源在中心位置出現了低阻凹陷,單源在中心位置兩側則出現了一個突起和一個凹陷,這個與平面圖的異常特征一致,多源激勵不僅可以增強信號的強度,而且可以改變目標體的異常響應特征,如圖19b、d、f所示.

    為了方便說明,在圖20中給出了512 Hz的中心測線(x=0和y=0)全域視電阻率變化曲線,對于y=0測線,四源和三源均出現了兩個低阻凹陷,且剛好對應異常體的中心位置,而單源在異常體邊緣處出現了極小值,中間位置則出現了極大值(圖20a).而對于x=0測線,四源和三源情況下均出現了低阻凹陷,這是因為雖然測線沒有經過異常體,但旁側兩個低阻異常體在多源激勵下出現視電阻率低值,而單源激勵下兩側分別出現了較強的極大極小值分布,不利于異常體電性特征及平面位置的定性分析(圖20b).

    4 結論

    本文采用基于電場雙旋度方程的非結構化節(jié)點有限元方法,成功實現了多源地空電磁響應特征的模擬和分析.通過Krylov子空間模型降階技術大大降低了大型稀疏矩陣的維度,然后直接求逆矩陣便可得到傳遞函數,進而求得多頻電磁方程的解.采用均勻半空間和層狀模型,通過和一維解析解進行對比,驗證了本文算法的正確性和有效性.通過泰勒級數展開取線性主部的迭代算法計算多源激勵的全域視電阻率,并將其應用于三維地下介質的定性分析解釋.三維地電模型的多源電磁響應分析表明:對于單個常體,一般單源激勵會在靠近源的一側出現高視電阻率,遠離源的一側出現低視電阻率;而對稱的三源或四源激勵模式可削弱場在邊界上的響應,其全域視電阻率特征相對簡單,和地下介質電性分布較吻合,這對于尋找淺地表獨立的目標體,如未爆炸的炮彈、地雷等是有利的.對于鄰近兩個低阻異常體,多源激勵模式不僅能夠增加信號的強度,而且能夠提高橫向分辨率.實際地下介質的電性變化非常復雜,多源激勵的電磁響應及全域視電阻率可為定性解釋提供參考,精細化的解釋將進一步依賴于地空三維電磁反演方法的發(fā)展.

    圖17 模型示意圖(a) 剖面圖; (b) 平面圖; (c) 源的加載方式; (d) 網格剖分. 紅色代表低阻體,綠色代表背景.Fig.17 Model of two conductive bodies(a) Cross-section; (b) Plane view; (c) Loading manner of source; (d) Unstructured grid. Red shows the conductive body and green shows background.

    圖18 256 Hz兩個低阻體的單源和多源全域視電阻率、總場幅值和二次場響應(a)—(c)①; (d)—(f)①⑤⑥; (g)—(i)①②③④.Fig.18 Full-domain apparent resistivity, total and secondary vertical magnetic field at 256 Hz of single and multiple sources

    附錄A

    (A1)

    (A2)

    (A3)

    (A4)

    (A5)

    圖19 不同頻率兩條中心測線的垂直磁場變化曲線圖(a)、(b) ①②③④; (c)、(d) ①⑤⑥; (e)、(f) ①.Fig.19 Vertical magnetic field curves of different frequencies in two central survey lines

    圖20 512 Hz下中心測線的全域視電阻率曲線(a) y=0; (b) x=0.Fig.20 Full-domain apparent resistivity curves at 512 Hz

    猜你喜歡
    多源測線全域
    極地海洋多波束測量測線布設系統(tǒng)設計及實現
    基于動態(tài)規(guī)劃的多波束測線布設模型
    ETC拓展應用場景下的多源異構交易系統(tǒng)
    用一粒米撬動全域旅游
    基于GNSS硬件在環(huán)的多源融合定位高逼真仿真方法
    “全域人人游”火爆周寧
    海峽姐妹(2017年9期)2017-11-06 08:39:37
    謀全域 圖四時 大連金普新區(qū)的全域“旅游+”
    遼寧經濟(2017年5期)2017-07-12 09:39:33
    全域旅游向更廣更深發(fā)展
    一種利用點特征和互信息的多源遙感影像配準方法
    遙感信息(2015年3期)2015-12-13 07:26:54
    空間位置信息的多源POI數據融合
    精品人妻视频免费看| 12—13女人毛片做爰片一| 熟妇人妻久久中文字幕3abv| 免费观看精品视频网站| 国产色婷婷99| 国产 一区 欧美 日韩| 精品无人区乱码1区二区| 欧美成人一区二区免费高清观看| 18禁在线播放成人免费| АⅤ资源中文在线天堂| 他把我摸到了高潮在线观看| 欧美日本亚洲视频在线播放| 日韩人妻高清精品专区| av.在线天堂| 国产精品99久久久久久久久| 久久久久国产精品人妻aⅴ院| 精品国内亚洲2022精品成人| 欧美另类亚洲清纯唯美| 少妇人妻精品综合一区二区 | 国产伦人伦偷精品视频| 少妇丰满av| 全区人妻精品视频| 小蜜桃在线观看免费完整版高清| 91精品国产九色| 美女cb高潮喷水在线观看| 九九热线精品视视频播放| 自拍偷自拍亚洲精品老妇| 成人av一区二区三区在线看| 少妇的逼水好多| 99久久成人亚洲精品观看| 久久久久性生活片| 简卡轻食公司| 欧美黑人巨大hd| 少妇裸体淫交视频免费看高清| 国产av不卡久久| 亚洲精品成人久久久久久| 成人毛片a级毛片在线播放| 国语自产精品视频在线第100页| 免费av观看视频| 99久久成人亚洲精品观看| 女的被弄到高潮叫床怎么办 | 免费无遮挡裸体视频| 国产亚洲精品综合一区在线观看| 乱码一卡2卡4卡精品| 极品教师在线视频| 午夜激情福利司机影院| 国产亚洲av嫩草精品影院| 国产成人福利小说| 一a级毛片在线观看| 人妻夜夜爽99麻豆av| 成人高潮视频无遮挡免费网站| 精品久久久久久久久久免费视频| 麻豆av噜噜一区二区三区| 国产亚洲av嫩草精品影院| 日本免费a在线| 一区二区三区免费毛片| 露出奶头的视频| 97超视频在线观看视频| 亚洲精品成人久久久久久| 亚洲美女黄片视频| 99久久久亚洲精品蜜臀av| 在线天堂最新版资源| 一个人看的www免费观看视频| 国产精品亚洲美女久久久| 亚洲av中文字字幕乱码综合| 日本一本二区三区精品| 亚洲不卡免费看| 国产亚洲欧美98| 麻豆成人午夜福利视频| 成年人黄色毛片网站| 中国美女看黄片| 99国产精品一区二区蜜桃av| 少妇裸体淫交视频免费看高清| a级毛片a级免费在线| 亚洲精品色激情综合| 久久99热这里只有精品18| 久久热精品热| 亚洲成人免费电影在线观看| 国产午夜精品论理片| 亚洲三级黄色毛片| 欧美高清性xxxxhd video| 99热6这里只有精品| 春色校园在线视频观看| 亚洲七黄色美女视频| 国内精品宾馆在线| 午夜老司机福利剧场| 两个人视频免费观看高清| 在线观看美女被高潮喷水网站| 亚洲精品国产成人久久av| АⅤ资源中文在线天堂| 国产男人的电影天堂91| 国产精品国产三级国产av玫瑰| 久久精品国产亚洲av香蕉五月| 国产主播在线观看一区二区| 国产精品久久久久久亚洲av鲁大| 亚洲美女黄片视频| 午夜免费成人在线视频| 久久欧美精品欧美久久欧美| 22中文网久久字幕| 国产精品一区二区性色av| 日韩中文字幕欧美一区二区| 免费av不卡在线播放| 亚洲av五月六月丁香网| 日韩精品中文字幕看吧| 又黄又爽又免费观看的视频| 内射极品少妇av片p| 中文字幕av在线有码专区| 日本一本二区三区精品| 免费看日本二区| 人妻少妇偷人精品九色| 午夜免费激情av| 狠狠狠狠99中文字幕| a级毛片免费高清观看在线播放| 久久中文看片网| 国产免费av片在线观看野外av| 黄片wwwwww| 校园人妻丝袜中文字幕| 一个人免费在线观看电影| 亚洲av一区综合| 亚洲av一区综合| 久久6这里有精品| 欧美zozozo另类| 欧美日韩乱码在线| 免费大片18禁| 一级毛片久久久久久久久女| 校园人妻丝袜中文字幕| 日本a在线网址| 成人午夜高清在线视频| 欧美日韩乱码在线| 国产高清有码在线观看视频| 嫩草影视91久久| 日本色播在线视频| 久久中文看片网| 18禁在线播放成人免费| 久久天躁狠狠躁夜夜2o2o| 好男人在线观看高清免费视频| 18禁黄网站禁片午夜丰满| 亚洲一区高清亚洲精品| 亚洲图色成人| 欧美一级a爱片免费观看看| 精品久久久久久久人妻蜜臀av| 久久久久久久久久成人| 免费av不卡在线播放| 99久久久亚洲精品蜜臀av| 日本与韩国留学比较| 搡女人真爽免费视频火全软件 | 欧美精品啪啪一区二区三区| av福利片在线观看| 欧美一区二区精品小视频在线| 午夜影院日韩av| 国产免费av片在线观看野外av| 国产精品三级大全| 在线观看美女被高潮喷水网站| 中国美女看黄片| 狂野欧美白嫩少妇大欣赏| 成年女人毛片免费观看观看9| 久久精品国产鲁丝片午夜精品 | 亚洲自偷自拍三级| 69人妻影院| 在线观看午夜福利视频| 亚洲国产精品成人综合色| 嫩草影院入口| 给我免费播放毛片高清在线观看| 亚洲四区av| 午夜老司机福利剧场| 亚洲av一区综合| a级毛片免费高清观看在线播放| 国产精品国产三级国产av玫瑰| 成年女人永久免费观看视频| 偷拍熟女少妇极品色| 哪里可以看免费的av片| 国产亚洲91精品色在线| 成年女人永久免费观看视频| 国产精品日韩av在线免费观看| 国国产精品蜜臀av免费| 欧美性猛交黑人性爽| 日韩强制内射视频| 午夜亚洲福利在线播放| 久久久久久久精品吃奶| 欧美xxxx性猛交bbbb| 免费在线观看成人毛片| 中出人妻视频一区二区| 精品久久久久久久久久久久久| 99热精品在线国产| 欧美绝顶高潮抽搐喷水| 亚洲人成网站在线播放欧美日韩| 久久久久久九九精品二区国产| 国产av在哪里看| 国产高清激情床上av| 国产高清视频在线观看网站| 乱人视频在线观看| 夜夜爽天天搞| 熟女人妻精品中文字幕| 色视频www国产| 久久精品国产亚洲av天美| 联通29元200g的流量卡| 国产精品不卡视频一区二区| 日韩,欧美,国产一区二区三区 | 亚洲欧美日韩东京热| 国产老妇女一区| 中文字幕av成人在线电影| 九九在线视频观看精品| 亚洲成人久久性| 亚洲精品在线观看二区| 国产精品无大码| 婷婷六月久久综合丁香| 天堂动漫精品| or卡值多少钱| 成人午夜高清在线视频| a在线观看视频网站| 国产人妻一区二区三区在| 男人狂女人下面高潮的视频| 国产大屁股一区二区在线视频| 欧美潮喷喷水| 又黄又爽又免费观看的视频| 国产极品精品免费视频能看的| 色在线成人网| 日日撸夜夜添| av在线天堂中文字幕| 欧美xxxx黑人xx丫x性爽| 国产成人aa在线观看| 校园春色视频在线观看| 久久久久久久亚洲中文字幕| 国产精品不卡视频一区二区| 中文亚洲av片在线观看爽| 自拍偷自拍亚洲精品老妇| 久久中文看片网| 听说在线观看完整版免费高清| 国模一区二区三区四区视频| 热99re8久久精品国产| 亚洲图色成人| 欧美日本视频| 国产精品久久久久久久电影| 2021天堂中文幕一二区在线观| 午夜福利欧美成人| 在线免费观看不下载黄p国产 | 老熟妇仑乱视频hdxx| 亚洲av中文av极速乱 | 在线看三级毛片| 免费黄网站久久成人精品| 联通29元200g的流量卡| 国产免费男女视频| 亚洲中文字幕日韩| 免费看a级黄色片| 亚洲综合色惰| 中文字幕熟女人妻在线| 禁无遮挡网站| 亚洲av不卡在线观看| 亚洲中文日韩欧美视频| 亚洲中文字幕一区二区三区有码在线看| 国产高清激情床上av| 欧美激情在线99| 91精品国产九色| 久久久久免费精品人妻一区二区| 91午夜精品亚洲一区二区三区 | 尤物成人国产欧美一区二区三区| 久久精品国产亚洲av涩爱 | 99热这里只有是精品在线观看| 中文字幕久久专区| 99视频精品全部免费 在线| avwww免费| а√天堂www在线а√下载| 久久久久性生活片| 国产精品野战在线观看| 不卡一级毛片| 久久中文看片网| 欧美成人一区二区免费高清观看| 国产精品久久久久久久久免| 日韩在线高清观看一区二区三区 | 美女高潮喷水抽搐中文字幕| 欧美黑人巨大hd| 特大巨黑吊av在线直播| 不卡视频在线观看欧美| 精华霜和精华液先用哪个| av视频在线观看入口| 99热这里只有精品一区| 亚洲成人免费电影在线观看| 88av欧美| 男人和女人高潮做爰伦理| 久久久久精品国产欧美久久久| 两个人的视频大全免费| 成人特级av手机在线观看| 国产精品久久电影中文字幕| 精品一区二区三区视频在线观看免费| 国产亚洲91精品色在线| 成年免费大片在线观看| 亚洲色图av天堂| 精品乱码久久久久久99久播| 成人av在线播放网站| 国产精品人妻久久久久久| av国产免费在线观看| 性插视频无遮挡在线免费观看| 最近最新中文字幕大全电影3| 亚洲成人久久性| 在线观看美女被高潮喷水网站| 国产av麻豆久久久久久久| 99热只有精品国产| 老师上课跳d突然被开到最大视频| 久久精品国产自在天天线| 亚洲人成伊人成综合网2020| 国产精品久久久久久亚洲av鲁大| 热99re8久久精品国产| 亚洲中文字幕一区二区三区有码在线看| 男女视频在线观看网站免费| 成年版毛片免费区| 天天躁日日操中文字幕| 国产精品98久久久久久宅男小说| 天天躁日日操中文字幕| 一个人看的www免费观看视频| 国产黄色小视频在线观看| 99精品在免费线老司机午夜| АⅤ资源中文在线天堂| 啦啦啦韩国在线观看视频| 色吧在线观看| 精品人妻一区二区三区麻豆 | 精品国产三级普通话版| 国产精品久久电影中文字幕| 国产高清激情床上av| 老司机福利观看| 亚洲精品一区av在线观看| 国产美女午夜福利| 97碰自拍视频| 在线看三级毛片| 亚洲 国产 在线| 成人av在线播放网站| 欧美日韩黄片免| 日韩欧美精品免费久久| 亚洲中文字幕日韩| 亚洲最大成人手机在线| 动漫黄色视频在线观看| 免费观看的影片在线观看| 一级黄色大片毛片| 午夜福利在线观看吧| 精品福利观看| 嫩草影院入口| 婷婷精品国产亚洲av| 丝袜美腿在线中文| 婷婷精品国产亚洲av| 91久久精品国产一区二区成人| 国产精品乱码一区二三区的特点| 亚洲欧美日韩卡通动漫| 伦理电影大哥的女人| 岛国在线免费视频观看| 伦理电影大哥的女人| 欧美激情在线99| 欧美精品啪啪一区二区三区| 欧美在线一区亚洲| h日本视频在线播放| 免费人成视频x8x8入口观看| 亚洲av五月六月丁香网| videossex国产| 免费无遮挡裸体视频| 丝袜美腿在线中文| 在线观看66精品国产| 夜夜爽天天搞| 99视频精品全部免费 在线| 久久久久国内视频| 国产极品精品免费视频能看的| 桃红色精品国产亚洲av| 免费观看在线日韩| 亚洲国产色片| 日韩人妻高清精品专区| 麻豆精品久久久久久蜜桃| 精品人妻1区二区| 欧美xxxx性猛交bbbb| 大型黄色视频在线免费观看| 麻豆精品久久久久久蜜桃| 婷婷精品国产亚洲av在线| 欧美xxxx性猛交bbbb| 又粗又爽又猛毛片免费看| 亚洲人与动物交配视频| 人人妻,人人澡人人爽秒播| 在线看三级毛片| 人人妻人人澡欧美一区二区| 老司机午夜福利在线观看视频| 日本色播在线视频| 日韩欧美在线二视频| 欧美bdsm另类| 天堂影院成人在线观看| 波多野结衣高清作品| 免费在线观看日本一区| 淫秽高清视频在线观看| 内射极品少妇av片p| 国产黄色小视频在线观看| 日韩,欧美,国产一区二区三区 | 999久久久精品免费观看国产| 日本欧美国产在线视频| 午夜福利高清视频| 桃色一区二区三区在线观看| 亚洲精品日韩av片在线观看| 国产精品久久电影中文字幕| 亚洲成a人片在线一区二区| 日本五十路高清| 啦啦啦啦在线视频资源| 最新中文字幕久久久久| 最近最新中文字幕大全电影3| 亚洲精品色激情综合| 成人综合一区亚洲| 狂野欧美激情性xxxx在线观看| 中国美白少妇内射xxxbb| 观看美女的网站| 久久久久国产精品人妻aⅴ院| 丝袜美腿在线中文| 亚洲中文字幕日韩| 久99久视频精品免费| 嫩草影视91久久| 成年免费大片在线观看| 国产一区二区在线观看日韩| 午夜久久久久精精品| 国产精品人妻久久久影院| 99精品在免费线老司机午夜| 国产精品精品国产色婷婷| 老司机深夜福利视频在线观看| 校园人妻丝袜中文字幕| 变态另类成人亚洲欧美熟女| 真人一进一出gif抽搐免费| 午夜福利高清视频| 又爽又黄无遮挡网站| 少妇被粗大猛烈的视频| 最近最新中文字幕大全电影3| 亚洲av成人av| 韩国av一区二区三区四区| 搡老妇女老女人老熟妇| 小说图片视频综合网站| 嫁个100分男人电影在线观看| 亚洲中文字幕一区二区三区有码在线看| 简卡轻食公司| 亚洲成人精品中文字幕电影| 别揉我奶头 嗯啊视频| 久久久久久久久久黄片| 精品久久久久久久人妻蜜臀av| 国产精品精品国产色婷婷| 国产精品乱码一区二三区的特点| 一边摸一边抽搐一进一小说| 老熟妇仑乱视频hdxx| 97超级碰碰碰精品色视频在线观看| 国产 一区精品| 亚洲图色成人| 国产av不卡久久| 欧美日本视频| 嫩草影院入口| 12—13女人毛片做爰片一| 亚洲欧美激情综合另类| 日本 av在线| 观看免费一级毛片| 琪琪午夜伦伦电影理论片6080| av在线亚洲专区| 免费搜索国产男女视频| 999久久久精品免费观看国产| 国产麻豆成人av免费视频| 国产精品无大码| 别揉我奶头 嗯啊视频| 国产高清视频在线观看网站| 欧美日本视频| 亚洲精华国产精华精| 久久久色成人| ponron亚洲| 日本黄色片子视频| 国产蜜桃级精品一区二区三区| 亚洲在线自拍视频| 99九九线精品视频在线观看视频| 精品福利观看| 国产真实乱freesex| 亚洲乱码一区二区免费版| 精品久久久久久成人av| 免费高清视频大片| 欧美+日韩+精品| 一区二区三区免费毛片| 精品人妻偷拍中文字幕| 如何舔出高潮| 亚洲精品色激情综合| 日韩欧美在线乱码| 性插视频无遮挡在线免费观看| 三级国产精品欧美在线观看| 免费看av在线观看网站| 国产私拍福利视频在线观看| 午夜免费男女啪啪视频观看 | 欧美日韩亚洲国产一区二区在线观看| 在线免费观看的www视频| 日日摸夜夜添夜夜添小说| 一夜夜www| 亚洲成人免费电影在线观看| 亚洲欧美激情综合另类| 99久久精品热视频| 天堂av国产一区二区熟女人妻| 搡老妇女老女人老熟妇| 国产午夜福利久久久久久| 亚洲欧美激情综合另类| 精品一区二区免费观看| 成人精品一区二区免费| 男女视频在线观看网站免费| 色综合色国产| 色播亚洲综合网| 国产乱人伦免费视频| 久久天躁狠狠躁夜夜2o2o| 在线观看免费视频日本深夜| 久久久久久久精品吃奶| 黄色配什么色好看| 天美传媒精品一区二区| 亚洲性久久影院| 亚洲最大成人av| 日日摸夜夜添夜夜添小说| 亚洲内射少妇av| 国产av一区在线观看免费| 国产一区二区在线观看日韩| 51国产日韩欧美| 中文字幕av成人在线电影| 在线免费十八禁| 亚洲一区二区三区色噜噜| 亚洲最大成人手机在线| 欧美最新免费一区二区三区| 夜夜爽天天搞| 成年版毛片免费区| 色播亚洲综合网| 成人高潮视频无遮挡免费网站| 最近视频中文字幕2019在线8| 国产精品人妻久久久久久| 久久人人爽人人爽人人片va| 国产在视频线在精品| 久久久久久伊人网av| 亚洲美女搞黄在线观看 | 国产精品爽爽va在线观看网站| 网址你懂的国产日韩在线| 久久午夜亚洲精品久久| 美女免费视频网站| 亚洲在线观看片| 国内毛片毛片毛片毛片毛片| 动漫黄色视频在线观看| 成年女人毛片免费观看观看9| 国产精品人妻久久久影院| 熟妇人妻久久中文字幕3abv| 男人和女人高潮做爰伦理| 深夜a级毛片| 天堂网av新在线| 日本a在线网址| 永久网站在线| 精品日产1卡2卡| 精品久久久久久久久av| 白带黄色成豆腐渣| 午夜激情欧美在线| 香蕉av资源在线| av在线亚洲专区| 给我免费播放毛片高清在线观看| 日本色播在线视频| 国产精品电影一区二区三区| 夜夜爽天天搞| 人人妻,人人澡人人爽秒播| 中文在线观看免费www的网站| 欧美性猛交╳xxx乱大交人| 欧美区成人在线视频| 成年免费大片在线观看| 国产精品人妻久久久影院| 亚洲最大成人中文| 乱系列少妇在线播放| 久久精品国产亚洲av香蕉五月| 日韩欧美在线乱码| 男插女下体视频免费在线播放| 日韩欧美精品v在线| 国产色婷婷99| 啪啪无遮挡十八禁网站| 国产精品久久久久久亚洲av鲁大| 国产 一区 欧美 日韩| 一个人免费在线观看电影| 亚洲成人久久性| 国产精品久久电影中文字幕| 亚洲最大成人av| 亚洲国产高清在线一区二区三| 久久精品国产99精品国产亚洲性色| 亚洲18禁久久av| 黄色配什么色好看| 热99re8久久精品国产| 国产精品一及| 亚洲va日本ⅴa欧美va伊人久久| 精品一区二区三区视频在线| 少妇丰满av| 欧美高清成人免费视频www| 日韩欧美精品免费久久| 亚洲国产精品sss在线观看| 免费人成在线观看视频色| 最近最新免费中文字幕在线| 国产精品福利在线免费观看| 黄色欧美视频在线观看| 69av精品久久久久久| 久久欧美精品欧美久久欧美| 又黄又爽又刺激的免费视频.| 日日干狠狠操夜夜爽| 91午夜精品亚洲一区二区三区 | 亚洲狠狠婷婷综合久久图片| 日本-黄色视频高清免费观看| 色av中文字幕| 久久久精品欧美日韩精品| 亚洲中文字幕一区二区三区有码在线看| 91精品国产九色| 亚洲av日韩精品久久久久久密| 麻豆成人av在线观看| 最新在线观看一区二区三区| 十八禁网站免费在线| 婷婷精品国产亚洲av| 啦啦啦啦在线视频资源| 有码 亚洲区| 亚洲人成网站在线播放欧美日韩| 91精品国产九色| 淫秽高清视频在线观看| 国产极品精品免费视频能看的| 欧美另类亚洲清纯唯美| av在线天堂中文字幕| 国产高潮美女av| 直男gayav资源| 日韩欧美免费精品| 日本免费一区二区三区高清不卡| 午夜久久久久精精品| 亚洲欧美日韩卡通动漫| 国产av在哪里看| 成人一区二区视频在线观看| 99riav亚洲国产免费|