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

    軸向各向異性巷-孔瞬變電磁三分量響應(yīng)特征研究

    2022-08-09 06:51:38郭建磊高小偉侯彥威
    煤田地質(zhì)與勘探 2022年7期
    關(guān)鍵詞:全空間電導(dǎo)率電磁

    郭建磊,高小偉,侯彥威

    (中煤科工集團西安研究院有限公司,陜西 西安 710077)

    煤礦開采條件復(fù)雜,安全生產(chǎn)問題突出,其中,在掘進工作面發(fā)生重特大水災(zāi)事故占比達51.16%[1],因此,提高超前預(yù)報的準(zhǔn)確性對保障煤礦安全生產(chǎn)至關(guān)重要。煤礦超前探測中一般采用鉆探和物探兩種方法,其中,鉆探方法結(jié)果可靠,但施工周期長、費用高,對巷道掘進生產(chǎn)影響較大。物探方法中,礦井瞬變電磁法[2]因無損、快速、信息豐富且對含水體反應(yīng)靈敏等獨特優(yōu)點而被廣泛應(yīng)用,但受巷道復(fù)雜環(huán)境影響較大。巷-孔瞬變電磁超前探測方法[3]有效結(jié)合了礦井瞬變電磁法和超前鉆孔的優(yōu)勢,該方法將發(fā)射線框布置在掘進工作面處,接收探頭放入鉆孔中,按間隔逐點采集三分量響應(yīng)。孫懷鳳等[4]通過物理模擬試驗證明可以通過巷-孔瞬變電磁垂直分量信號在幅值和方向反號時間上存在的差異判斷掘進工作面前方是否存在異常構(gòu)造。陳丁等[5]采用積分方程法以煤層底板下方存在低阻體為模型研究了全空間條件下巷-孔瞬變電磁場的響應(yīng)特征,結(jié)果表明,水平分量響應(yīng)的延續(xù)時間優(yōu)于垂向分量,垂向分量響應(yīng)的幅值強于水平分量,在實際觀測中不僅要觀測垂向分量,也要觀測水平分量。范濤[6-7]利用鉆孔瞬變電磁垂直分量實現(xiàn)二維擬地震反演,并進行了三維數(shù)值模擬、水槽物理模擬和井下現(xiàn)場模擬試驗;后續(xù)利用趨勢面提取技術(shù)提取了水平分量純異常場,通過異常場形態(tài)組合和幅值關(guān)系實現(xiàn)異常體中心方位角定位,并在此基礎(chǔ)上通過聚類算法對水平分量異常形態(tài)進行自動分類,完成異常體空間角度定位[8]。該方法達到“一孔多用”的目的,有效克服了電磁和人文干擾,具有觀測點距離異常體近、采集數(shù)據(jù)信噪比高等優(yōu)點,在全國多個煤礦獲得了廣泛應(yīng)用[9-10]。

    巷-孔瞬變電磁法在資料處理解釋過程中將地層假設(shè)為各向同性介質(zhì),但隨著電磁勘探方法的發(fā)展,越來越多的學(xué)者認識到電導(dǎo)率各向異性在地下構(gòu)造中普遍分布[11-13],并對介質(zhì)電導(dǎo)率各向異性特性與矢量電磁場關(guān)系做了大量研究。當(dāng)前,關(guān)于頻率域電磁法,如大地電磁法[14-18]、海洋可控源電磁法[19-21]、直流電法[22-24]、井間電磁法[25-26]和航空電磁法[27-28]等方法的各向異性正反演研究均取得了較為豐富的成果。關(guān)于時間域電磁法的各向異性研究主要針對半空間介質(zhì)和電導(dǎo)率各向異性對垂直分量的影響。周建美等[29]基于有限體積法實現(xiàn)雙軸各向異性瞬變電磁三維正演,劉亞軍等[30]基于有限體積法實現(xiàn)任意各向異性瞬變電磁三維正演。在全空間介質(zhì)方面,程久龍等[31]基于時域有限差分法研究了不同主軸各向異性介質(zhì)對礦井瞬變電磁場垂直分量的影響特征,解釋了探測結(jié)果存在偏差的原因,并結(jié)合實例進行說明。上述研究并未涉及電導(dǎo)率各向異性對三分量響應(yīng)的影響。

    目前,巷-孔瞬變電磁法從硬件到軟件均實現(xiàn)了三分量信號的采集與解釋,通過將地層假設(shè)為各向同性介質(zhì),依據(jù)水平分量形態(tài)組合及幅值關(guān)系判定孔旁異常體的分布象限,依據(jù)垂直分量反演計算異常體相關(guān)參數(shù),最終實現(xiàn)巷-孔瞬變電磁立體成像解釋。由于資料解釋受制于對介質(zhì)電導(dǎo)率各向同性的假設(shè),導(dǎo)致在存在明顯電導(dǎo)率各向異性的勘探區(qū)會產(chǎn)生較大的解釋誤差,因此,非常有必要研究全空間條件下電導(dǎo)率各向異性對巷-孔瞬變電磁響應(yīng)的影響方式及程度,為進一步提高超前精度提供指導(dǎo)依據(jù)。本文研究采用時域有限差分算法(Finite Difference Time Domain,F(xiàn)DTD)[32]實現(xiàn)各向異性瞬變電磁三維正演,通過引入電導(dǎo)率各向異性張量構(gòu)建控制方程,以差分代替微分對控制方程進行離散,空間離散采用后向差分,時間離散采用中心差分,以電流密度的形式加入各向異性Maxwell 方程組的安培環(huán)路定理實現(xiàn)任意電流源的加載,進而推導(dǎo)出有源和無源區(qū)域電磁場迭代表達式。最后,構(gòu)建全空間模型、層狀模型和三維塊狀模型進行正演并分析電導(dǎo)率軸向各向異性對巷-孔瞬變電磁三分量響應(yīng)的影響程度與方式。

    1 正演理論基礎(chǔ)

    1.1 控制方程

    電導(dǎo)率各向異性、無源媒質(zhì)中的Maxwell 方程組為:

    式中:E為電場強度,V/m;H為磁場強度,A/m;為張量電導(dǎo)率,S/m;μ為磁導(dǎo)率,Wb/(A·m);ε為介電常數(shù),F(xiàn)/m;t為時間,s。

    地球物理低頻電磁勘探中一般忽略位移電流,為了滿足三維時域有限差分計算要求,因此,加入虛擬介電常數(shù)構(gòu)成顯式的時間迭代格式,式(1b)變?yōu)椋?/p>

    式中:γ為虛擬介電常數(shù)。

    采用經(jīng)典Yee 氏網(wǎng)格對計算區(qū)域進行剖分,電場分量位于棱邊中心點,磁場分量位于面中心點,每一個電場(磁場)分量均被相應(yīng)的4 個磁場(電場)分量包圍,電場分量和磁場分量在時間軸上交替采樣。

    采用非均勻網(wǎng)格形式[32]進行模型剖分,將式(2)在t(n+1/2)時刻進行空間和時間離散,以差分代替微分,其中,空間上采用后向差分、時間上采用中心差分,無源區(qū)域得到如下軸向各向異性的電場強度迭代公式:

    在有源媒質(zhì)區(qū)域,式(2)必須包含源電流密度項,將式(2)修改為:

    式中:Js為源電流密度。

    巷-孔瞬變電磁法采用回線源進行發(fā)射,坐標(biāo)系滿足“右手定則”,假定,巷道方向為z方向,上下為x方向,左右為y方向,發(fā)射源將被放置在xoy平面,因此,式(6)中僅存在x、y方向上的源電流密度項(Jsx、Jsy)??紤]低頻近似后,在直角坐標(biāo)系中將式(6)展開得到下式,按照差分格式正常離散,可以得到有源區(qū)域網(wǎng)格電場迭代公式。磁場迭代公式不包含電導(dǎo)率項,因此,各向異性條件下磁場迭代公式與各向同性條件下磁場迭代公式相同[32]。

    1.2 算法驗證

    構(gòu)建均勻全空間各向同性模型進行正演模擬,通過將數(shù)值解與解析解進行對比驗證算法精度。全空間介質(zhì)電導(dǎo)率為0.01 S/m,采用笛卡爾坐標(biāo)系,回線源中心位于坐標(biāo)系原點,邊長為3 m×3 m,電流為1 A,發(fā)射源上升沿和下降沿均為1 μs,持續(xù)時間為10 ms,二次場采樣時間為10 ms。最小網(wǎng)格尺寸為1 m,相鄰網(wǎng)格放大系數(shù)為1.05,剖分網(wǎng)格數(shù)為221×221×200,接收點位于(0.5 m,0.5 m,0 m),在接收點處接收?Bx/?t、?By/?t、?Bz/?t響應(yīng)。

    數(shù)值解與解析解對比結(jié)果如圖1 所示。由圖1 可知,數(shù)值解和解析解的三分量響應(yīng)曲線基本重合,當(dāng)時間達到0.003 ms 時誤差控制在5%以內(nèi),其中,?Bx/?t和?By/?t的相對誤差稍大于?Bz/?t的相對誤差,由于接收點位于坐標(biāo)系對角線上,故?Bx/?t與?By/?t響應(yīng)曲線重合且相對誤差一致。故本文算法滿足計算精度要求。

    圖1 數(shù)值解與解析解對比結(jié)果Fig.1 Comparison of numerical solutions and analytical solutions

    2 各向異性巷-孔瞬變電磁響應(yīng)特征

    2.1 各向異性全空間模型

    巷-孔瞬變電磁法勘探環(huán)境為全空間介質(zhì),因此,建立如圖2 所示模型研究全空間介質(zhì)電導(dǎo)率各向異性對巷-孔瞬變電磁三分量響應(yīng)的影響方式與程度。如圖2 所示,采用上節(jié)的坐標(biāo)系系統(tǒng),巷道內(nèi)空氣電導(dǎo)率為0.000 1 S/m,空腔范圍為(-3~3 m,-3~3 m,0~-∞ m)。發(fā)射源平貼放置于巷道工作面且中心位于坐標(biāo)原點,測點位于超前鉆孔中,鉆孔孔口位于(0.5 m,0.5 m,0 m),鉆孔平行于z軸,長度范圍為0~80 m,按照1 m 的點間距從孔底至孔口逐點采集三分量響應(yīng)數(shù)據(jù)。模型a為各向同性,電導(dǎo)率為σx=σy=σz=0.01 S/m;模型b 為x軸各向異性,電導(dǎo)率為σx=0.1 S/m、σy=σz=0.01 S/m;模型c 為y軸各向異性,電導(dǎo)率為σy=0.1 S/m、σx=σz=0.01 S/m;模型d 為z軸各向異性,電導(dǎo)率為σz=0.1 S/m、σx=σy=0.01 S/m。模型剖分與發(fā)射源參數(shù)與1.2 節(jié)一致。

    圖2 全空間各向異性巷-孔瞬變電磁超前探測模型Fig.2 Full-space anisotropic tunnel-hole transient electromagnetic advance detection model

    圖3 為4 個模型的三分量響應(yīng)多測道圖,時間道的時間為0.06~0.21 ms。由圖3 可知,4 個模型的?Bx/?t、?By/?t多測道圖形態(tài)呈反“C”型、?Bz/?t多測道圖形態(tài)呈半“C”型;模型a 的?Bx/?t、?By/?t多測道圖形態(tài)和幅值一致,最大幅值位于25 m;?Bz/?t多測道圖最大幅值位于0 m,隨深度增加逐漸衰減;與模型a 相比,模型b 的?Bx/?t多測道圖形態(tài)和幅值在發(fā)射源處發(fā)生微弱變化,?By/?t多測道圖幅值整體增大且最大幅值位于10 m、至40 m 趨于0,?Bz/?t多測道圖幅值整體增大,至40 m 趨于0;與模型a 相比,模型c 的?Bx/?t多測道圖幅值整體增大且最大幅值位于10 m、至40 m 趨于0,?By/?t多測道圖形態(tài)和幅值在發(fā)射源處發(fā)生微弱變化;?Bz/?t多測道圖幅值整體增大、至40 m處趨于0。模型d 的三分量響應(yīng)多測道圖特征與模型a 的基本一致。

    圖3 三分量響應(yīng)多測道圖Fig.3 Three-component responses multi-channel diagram

    將4 個模型的0.06 ms 時刻的三分量響應(yīng)曲線進行對比(圖4)發(fā)現(xiàn),x軸各向異性主要影響?By/?t和?Bz/?t響應(yīng),y軸各向異性主要影響?Bx/?t和?Bz/?t響應(yīng),z軸各向異性對三分量響應(yīng)基本不產(chǎn)生影響。根據(jù)“煙圈效應(yīng)”,回線源產(chǎn)生的感應(yīng)電流主要在水平方向,三分量響應(yīng)受到水平方向電導(dǎo)率的影響較大,而受到垂直方向的電導(dǎo)率影響很小;x方向感應(yīng)電流垂直穿過y軸電導(dǎo)率,y方向感應(yīng)電流垂直穿過x軸電導(dǎo)率,因此,x軸電導(dǎo)率主要影響?Bx/?t響應(yīng),y軸電導(dǎo)率主要影響?By/?t響應(yīng)。

    圖4 三分量響應(yīng)曲線對比(t=0.06 ms)Fig.4 Comparison of three-component responses (t=0.06 ms)

    2.2 各向異性層狀模型

    對于水平鉆孔而言,結(jié)合煤炭固有的地層特征,巷-孔瞬變電磁法勘探環(huán)境基本處于層狀介質(zhì)。當(dāng)回線源位置固定后,地下不同深度地層受到一次場激發(fā)的程度不同,地層電導(dǎo)率各向異性對三分量響應(yīng)的影響程度也不一樣。

    以各向同性全空間模型為基礎(chǔ),分別將距回線源不同距離的地層設(shè)置為各向異性(圖5)。如圖5 所示,模型a 為上部地層(40~∞ m)電導(dǎo)率分別為各向同性和各向異性、模型b 為中間地層(-40~40 m)電導(dǎo)率分別為各向同性和各向異性、模型c 為下部地層(-40~-∞ m)電導(dǎo)率分別為各向同性和各向異性。當(dāng)?shù)貙訛楦飨蛲詴r,電導(dǎo)率σx=σy=σz=0.01 S/m;當(dāng)x軸為各向異性時,電導(dǎo)率σx=0.1 S/m、σy=σz=0.01 S/m;當(dāng)y軸為各向異性時,電導(dǎo)率σy=0.1 S/m、σx=σz=0.01 S/m;當(dāng)z軸為各向異性時,電導(dǎo)率σz=0.1 S/m、σx=σy=0.01 S/m。模型剖分與發(fā)射源參數(shù)與1.2 節(jié)一致。

    圖5 層狀模型Fig.5 Schematic diagram of the layered models

    對比分析模型a-模型c 的不同電導(dǎo)率參數(shù)下0.06 ms 時刻的三分量響應(yīng)曲線(圖6)。由圖6a 可知,對于回線源上部的地層,y軸各向異性對?Bx/?t響應(yīng)影響明顯,且使?Bx/?t幅值由負變正,對?By/?t、?Bz/?t響應(yīng)幾乎沒有影響;同時,x軸、z軸各向異性與各向同性的三分量響應(yīng)曲線幾乎重合。由圖6c 可知,y軸各向異性對?Bx/?t響應(yīng)影響明顯,與模型a 的?Bx/?t響應(yīng)反映不同的是,模型c 的?Bx/?t響應(yīng)幅值正負不變且幅值增大。因此,通過水平分量形態(tài)和幅值的相互關(guān)系可以判斷各向異性介質(zhì)所在方位及電導(dǎo)率主軸方向,這主要是因為模型a 的地層位于回線源上方,模型c的地層位于回線源下方,巷-孔瞬變電磁三分量響應(yīng)具有方向性。由圖6b 可知,x軸、z軸各向異性和各向同性的?Bx/?t響應(yīng)曲線重合且偏離y軸各向異性的?Bx/?t響應(yīng)曲線;y軸、z軸各向異性和各向同性的?By/?t響應(yīng)曲線重合且偏離x軸各向異性的?By/?t響應(yīng)曲線;x軸和y軸各向異性的?Bz/?t響應(yīng)曲線重合且與z軸各向異性和各向同性的?Bz/?t響應(yīng)曲線偏離較大,說明距離發(fā)射源較近地層的電導(dǎo)率各向異性對三分量響應(yīng)的影響較大,也再次說明巷-孔瞬變電磁三分量響應(yīng)主要受水平方向電導(dǎo)率的影響。

    圖6 三分量響應(yīng)曲線對比(t=0.06 ms)Fig.6 Comparison of three-component responses (t=0.06 ms)

    2.3 各向異性塊狀模型

    煤炭資源廣泛分布于沉積巖中,沉積巖地區(qū)由于層理發(fā)育導(dǎo)致地下介質(zhì)電阻率隨電流方向發(fā)生變化[33];煤炭開采產(chǎn)生不同形態(tài)、高度和密度的采空裂隙和裂縫,裂隙和裂縫充水后形成裂隙含水體[34],上述情況均會導(dǎo)致地層介質(zhì)具有電導(dǎo)率各向異性特性。因此,本節(jié)主要研究各向同性地層中含有各向異性塊狀模型和各向異性地層中含有各向同性塊體模型。

    1)各向同性地層中含有各向異性塊狀模型

    建立如圖7 所示的模型,異常體1 的規(guī)模為20 m×20 m×20 m,異常體2 的規(guī)模為30 m×20 m×20 m,異常體3 的規(guī)模為20 m×30 m×20 m,3 個異常體的中心均位于(20 m,20 m,20 m)。將3 個異常體分別設(shè)置為各向同性和各向異性,為各向同性時電導(dǎo)率為σx=σy=σz=0.1 S/m;為x軸各向 異性時電導(dǎo)率 為σx=1 S/m、σy=σz=0.1 S/m;為y軸各向異性時電導(dǎo)率為σy=1 S/m、σx=σz=0.1 S/m;為z軸各向異性時電導(dǎo)率為σz=1 S/m、σx=σy=0.1 S/m。模型剖分與發(fā)射源參數(shù)與1.2 節(jié)一致。

    圖7 塊狀異常體模型俯視圖Fig.7 Top view of block abnormal model

    將3 個異常體的0.06 ms 時刻的三分量響應(yīng)曲線進行對比(圖8)。由圖8 可知,當(dāng)3 個異常體均為各向同性時,?Bx/?t響應(yīng)大小為異常體3>異常體2>異常體1,說明y方向邊長改變影響?Bx/?t響應(yīng);?By/?t響應(yīng)的大小為異常體2>異常體3>異常體1,說明改變x方向邊長改變影響?By/?t響應(yīng);?Bz/?t響應(yīng)大小為異常體2=異常體3 >異常體1。當(dāng)3 個異常體均為z軸各向異性時,其三分量響應(yīng)形態(tài)、規(guī)律與各向同性時相同,但其水平分量響應(yīng)幅度大于各向同性;當(dāng)異常體均為x軸和y軸各向異性時,其?Bx/?t和?By/?t響應(yīng)形態(tài)、規(guī)律與各向同性的相同,不同之處在于,當(dāng)3 個異常體均為x軸各向異性時,異常體2 的?Bz/?t響應(yīng)>異常體3,當(dāng)3 個異常體均為y軸各向異性時,異常體3 的?Bz/?t響應(yīng)>異常體2。

    圖8 三分量響應(yīng)曲線對比(t=0.06 ms)Fig.8 Comparison of three-component responses (t=0.06 ms)

    將水平分量響應(yīng)曲線進行進一步對比(圖9)。當(dāng)3 個異常體均為x軸各向異性時,異常體1 的?By/?t響應(yīng)>?Bx/?t響應(yīng),異常體3 的?Bx/?t(?By/?t)響應(yīng)>異常體2 的?By/?t(?Bx/?t)響應(yīng)。當(dāng)3 個異常體均為y軸各向異性時,異常體1 的?Bx/?t響應(yīng)>?By/?t響應(yīng),異常體2 的?By/?t(?Bx/?t)響應(yīng)>異常體3 的?Bx/?t(?By/?t)響應(yīng)。當(dāng)3 個異常體均為z軸各向異性時,異常體1 的?Bx/?t響應(yīng)和?By/?t響應(yīng)曲線重合,異常體2 的?By/?t(?Bx/?t)響應(yīng)與異常體3 的?Bx/?t(?By/?t)響應(yīng)曲線重合,且異常體2 的?By/?t(異常體3 的?Bx/?t)響應(yīng)>異常體2 的?Bx/?t(異常體3 的?By/?t)響應(yīng)。上述情況說明異常體邊長不同對水平分量響應(yīng)的影響強弱不同,與模型長軸平行方向的電導(dǎo)率改變對三分量響應(yīng)的影響較大。

    圖9 水平分量響應(yīng)對比(t=0.06 ms)Fig.9 Comparison diagram of horizontal component responses (t=0.06 ms)

    2)各向異性地層中含有各向同性塊體模型

    基于圖7 中的異常體1 模型,將異常體1 設(shè)置為各向同性,電導(dǎo)率為σx=σy=σz=0.1 S/m,將全空間介質(zhì)分別設(shè)置為各向同性和各向異性,為各向同性時電導(dǎo)率為σx=σy=σz=0.01 S/m;為x軸各向異性時電導(dǎo)率為σx=0.1 S/m、σy=σz=0.01 S/m;為y軸各向異性時電導(dǎo)率為σy=0.1 S/m、σx=σz=0.01 S/m;為z軸各向異性時電導(dǎo)率為σz=0.1 S/m、σx=σy=0.01 S/m。

    將全空間介質(zhì)為各向同性和各向異性時0.06 ms時刻的三分量響應(yīng)曲線進行對比(圖10)。由圖10 可知,全空間介質(zhì)為各向同性與z軸各向異性的三分量響應(yīng)幾乎重合,?Bx/?t和?By/?t響應(yīng)曲線在距離發(fā)射源40 m 后出現(xiàn)偏離;改變?nèi)臻g介質(zhì)x軸和y軸電導(dǎo)率對異常體的三分量響應(yīng)均影響較大,使?Bx/?t和?By/?t異常幅值由正變負,且改變y軸電導(dǎo)率對?Bx/?t響應(yīng)的影響>改變x軸電導(dǎo)率、改變x軸電導(dǎo)率對?By/?t響應(yīng)的影響>改變y軸電導(dǎo)率,改變x軸和y軸電導(dǎo)率對?Bz/?t響應(yīng)的影響一樣,這與上述小節(jié)的規(guī)律一致。全空間介質(zhì)為各向異性的三分量異常響應(yīng)大于各向同性的三分量異常響應(yīng),說明異常體產(chǎn)生的異常響應(yīng)被淹沒在全空間介質(zhì)電導(dǎo)率各向異性產(chǎn)生的異常響應(yīng)中,如果此時不考慮介質(zhì)的各向異性特性會對巷-孔瞬變電磁解釋造成較大誤差。

    圖10 三分量響應(yīng)對比(t=0.06 ms)Fig.10 Comparison of three-component responses (t=0.06 ms)

    3 結(jié)論

    a.介質(zhì)距回線源的距離不同,其各向異性特性對巷-孔瞬變電磁三分量響應(yīng)的影響強弱和方式不同,可以通過其三分量響應(yīng)的形態(tài)和幅值的相互關(guān)系辨別各向異性介質(zhì)所在方位及主軸電導(dǎo)率方向,同時由于三維異常體的邊長不同對其三分量響應(yīng)的影響強弱不同,進一步增加了各向異性條件下三分量解釋的復(fù)雜性。

    b.回線源產(chǎn)生的感應(yīng)電流主要是水平方向,并且x方向感應(yīng)電流垂直穿過y軸電導(dǎo)率,y方向感應(yīng)電流垂直穿過x軸電導(dǎo)率,因此,z軸各向異性對三分量響應(yīng)基本沒有影響,x軸和y軸各向異性對巷-孔瞬變電磁三分量響應(yīng)有較大影響,同時?Bx/?t響應(yīng)主要受y軸電導(dǎo)率影響,?By/?t響應(yīng)主要受x軸電導(dǎo)率影響。

    c.當(dāng)異常體所處介質(zhì)呈電導(dǎo)率各向異性時,異常體產(chǎn)生的異常響應(yīng)被淹沒在全空間介質(zhì)電導(dǎo)率各向異性產(chǎn)生的異常響應(yīng)中,介質(zhì)的電導(dǎo)率各向異性特征在巷-孔瞬變電磁法解釋過程中不可忽略,因此,本文研究內(nèi)容為進一步研究各向異性提供一定的基礎(chǔ),也為后續(xù)進行各向異性反演提供參考。

    d.未來的工作包括引入更復(fù)雜的電導(dǎo)率張量、進行亞網(wǎng)格剖分、添加吸收邊界條件等,實現(xiàn)任意各向異性的不規(guī)則復(fù)雜模型三維快速正演,為實際生產(chǎn)等提供重要的參考意義。

    猜你喜歡
    全空間電導(dǎo)率電磁
    三維多孔電磁復(fù)合支架構(gòu)建與理化表征
    基于比較測量法的冷卻循環(huán)水系統(tǒng)電導(dǎo)率檢測儀研究
    低溫脅迫葡萄新梢電導(dǎo)率和LT50值的研究
    掌握基礎(chǔ)知識 不懼電磁偏轉(zhuǎn)
    全空間瞬變電磁場低阻層屏蔽效應(yīng)數(shù)值模擬研究
    中國煤炭(2016年1期)2016-05-17 06:11:33
    多層圓線圈在電磁計算中全空間磁感應(yīng)強度B的分布
    全空間上一類半線性雙調(diào)和方程正解的衰減
    高電導(dǎo)率改性聚苯胺的合成新工藝
    電磁換向閥應(yīng)用探討
    河南科技(2014年16期)2014-02-27 14:13:21
    瞬變電磁法在煤礦防治水中的應(yīng)用
    河南科技(2014年8期)2014-02-27 14:07:41
    一级毛片 在线播放| 久久精品综合一区二区三区| 啦啦啦中文免费视频观看日本| 下体分泌物呈黄色| 精品久久久久久久人妻蜜臀av| 日本免费在线观看一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产乱人偷精品视频| 亚洲国产日韩一区二区| 久久99蜜桃精品久久| 2018国产大陆天天弄谢| 熟女人妻精品中文字幕| 国产精品国产av在线观看| 美女脱内裤让男人舔精品视频| 22中文网久久字幕| 亚洲欧美日韩无卡精品| 街头女战士在线观看网站| 人妻夜夜爽99麻豆av| 国产成人精品久久久久久| 中文天堂在线官网| 免费高清在线观看视频在线观看| 久久国产乱子免费精品| eeuss影院久久| 国产精品一及| 国产亚洲午夜精品一区二区久久 | 欧美日本视频| 久久99精品国语久久久| 日韩成人av中文字幕在线观看| 秋霞在线观看毛片| 国产伦理片在线播放av一区| 亚洲一区二区三区欧美精品 | 亚洲,欧美,日韩| 日韩av在线免费看完整版不卡| 国产黄a三级三级三级人| 国产高清有码在线观看视频| 婷婷色av中文字幕| 搞女人的毛片| 精华霜和精华液先用哪个| 欧美性感艳星| 国产精品一二三区在线看| 在线播放无遮挡| 国产黄色免费在线视频| 69人妻影院| 人妻系列 视频| 日韩电影二区| 91午夜精品亚洲一区二区三区| 久久人人爽人人爽人人片va| 亚洲av免费高清在线观看| 国产高清三级在线| 久久久久久久久久人人人人人人| 少妇人妻 视频| 国产高清有码在线观看视频| 嘟嘟电影网在线观看| 国产探花极品一区二区| 久久精品久久久久久噜噜老黄| av在线播放精品| 欧美成人精品欧美一级黄| 内射极品少妇av片p| 黄片无遮挡物在线观看| 最近中文字幕高清免费大全6| 欧美丝袜亚洲另类| 欧美日本视频| 亚洲欧美成人综合另类久久久| 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| www.av在线官网国产| 精华霜和精华液先用哪个| 国产av国产精品国产| 嘟嘟电影网在线观看| 亚洲最大成人中文| 亚洲综合色惰| 五月天丁香电影| 精品视频人人做人人爽| 一区二区三区精品91| 国产精品精品国产色婷婷| 韩国高清视频一区二区三区| 欧美日本视频| 少妇人妻 视频| 欧美xxxx黑人xx丫x性爽| 中文在线观看免费www的网站| 在线a可以看的网站| 99久久精品热视频| 亚洲最大成人手机在线| 18+在线观看网站| 神马国产精品三级电影在线观看| 日韩强制内射视频| 亚洲精品456在线播放app| 看免费成人av毛片| 99热国产这里只有精品6| 一个人观看的视频www高清免费观看| 九九爱精品视频在线观看| 国产黄片美女视频| 大片电影免费在线观看免费| 日韩 亚洲 欧美在线| 免费av观看视频| 免费看a级黄色片| 亚洲av一区综合| 超碰97精品在线观看| 国产欧美亚洲国产| 各种免费的搞黄视频| 下体分泌物呈黄色| av在线老鸭窝| 大香蕉久久网| 女人久久www免费人成看片| 欧美精品一区二区大全| 97人妻精品一区二区三区麻豆| 国产成人精品婷婷| av在线播放精品| 亚洲成人久久爱视频| 丝袜喷水一区| 色综合色国产| 国产亚洲5aaaaa淫片| 2021天堂中文幕一二区在线观| 91aial.com中文字幕在线观看| 久久久久久九九精品二区国产| 免费播放大片免费观看视频在线观看| xxx大片免费视频| 嫩草影院新地址| 国产老妇伦熟女老妇高清| 熟妇人妻不卡中文字幕| 青春草视频在线免费观看| 九草在线视频观看| 午夜日本视频在线| 免费电影在线观看免费观看| 欧美性感艳星| 精品一区二区三区视频在线| 亚洲av不卡在线观看| 小蜜桃在线观看免费完整版高清| 国产av码专区亚洲av| 亚洲av福利一区| 国产中年淑女户外野战色| 亚洲性久久影院| 视频区图区小说| 亚洲精品国产av成人精品| 国产欧美日韩精品一区二区| 日韩av在线免费看完整版不卡| 色视频在线一区二区三区| 高清毛片免费看| 永久网站在线| 校园人妻丝袜中文字幕| 亚洲在线观看片| 日韩av不卡免费在线播放| 亚洲av电影在线观看一区二区三区 | 丝袜美腿在线中文| 黄色视频在线播放观看不卡| 亚洲国产精品国产精品| 国产伦在线观看视频一区| 国产免费一级a男人的天堂| 亚洲精品日韩av片在线观看| 又粗又硬又长又爽又黄的视频| 嫩草影院入口| 成人二区视频| 97超碰精品成人国产| 1000部很黄的大片| 国产精品久久久久久精品电影| 黄色怎么调成土黄色| 蜜桃久久精品国产亚洲av| 肉色欧美久久久久久久蜜桃 | 亚洲欧美日韩东京热| 日日撸夜夜添| 热re99久久精品国产66热6| 久久影院123| 亚洲国产成人一精品久久久| 精品熟女少妇av免费看| 黄色一级大片看看| av.在线天堂| 视频区图区小说| 久久人人爽av亚洲精品天堂 | 少妇熟女欧美另类| 国产探花在线观看一区二区| 久久久久久久久久人人人人人人| 免费观看a级毛片全部| 国产亚洲最大av| 亚洲精品日本国产第一区| 99久久人妻综合| 插逼视频在线观看| 一级毛片久久久久久久久女| 高清毛片免费看| 国产高潮美女av| 少妇裸体淫交视频免费看高清| av线在线观看网站| 国产黄频视频在线观看| 一级毛片我不卡| 国产在线男女| 色吧在线观看| 免费看日本二区| 熟女电影av网| 偷拍熟女少妇极品色| 极品少妇高潮喷水抽搐| 街头女战士在线观看网站| 午夜老司机福利剧场| 国产午夜精品一二区理论片| 久久精品久久久久久噜噜老黄| 五月天丁香电影| 十八禁网站网址无遮挡 | 婷婷色综合大香蕉| 少妇人妻精品综合一区二区| 免费看光身美女| xxx大片免费视频| 免费看日本二区| 亚洲av日韩在线播放| 嫩草影院新地址| 午夜精品一区二区三区免费看| 欧美亚洲 丝袜 人妻 在线| videossex国产| 少妇的逼水好多| 精品国产三级普通话版| 国产乱来视频区| av天堂中文字幕网| 少妇裸体淫交视频免费看高清| 观看免费一级毛片| 精品国产三级普通话版| 欧美3d第一页| 日韩人妻高清精品专区| 插阴视频在线观看视频| 伊人久久精品亚洲午夜| 综合色av麻豆| 涩涩av久久男人的天堂| 精华霜和精华液先用哪个| 好男人在线观看高清免费视频| 内射极品少妇av片p| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久国产精品人妻一区二区| 亚洲色图av天堂| 97在线视频观看| 久久精品国产亚洲av天美| 亚洲av国产av综合av卡| 国产老妇伦熟女老妇高清| 欧美激情久久久久久爽电影| 国产男女内射视频| 亚洲天堂国产精品一区在线| 日日啪夜夜撸| 下体分泌物呈黄色| 免费观看av网站的网址| 国产免费视频播放在线视频| 国产精品国产av在线观看| 国产成人91sexporn| 久久久久性生活片| 国产午夜福利久久久久久| av国产精品久久久久影院| av专区在线播放| 免费播放大片免费观看视频在线观看| 狂野欧美激情性bbbbbb| 亚洲欧美精品自产自拍| 成人无遮挡网站| 特级一级黄色大片| 激情 狠狠 欧美| 国产精品麻豆人妻色哟哟久久| 深夜a级毛片| 天天躁日日操中文字幕| 极品少妇高潮喷水抽搐| 亚洲国产高清在线一区二区三| 国产精品国产三级国产av玫瑰| 美女内射精品一级片tv| 亚洲人成网站在线播| 欧美三级亚洲精品| 高清av免费在线| 免费观看性生交大片5| 在线观看三级黄色| 2021少妇久久久久久久久久久| 国产高清三级在线| 内射极品少妇av片p| 精品国产露脸久久av麻豆| 欧美亚洲 丝袜 人妻 在线| 99久久精品热视频| 亚洲综合精品二区| 中文字幕免费在线视频6| 国产精品.久久久| 精品一区二区免费观看| 一区二区三区乱码不卡18| 国产高清三级在线| 成人毛片60女人毛片免费| 夫妻午夜视频| 国产 一区 欧美 日韩| 亚洲精品自拍成人| 18禁在线播放成人免费| 嘟嘟电影网在线观看| 欧美激情久久久久久爽电影| 亚洲欧洲国产日韩| 国产免费一区二区三区四区乱码| 国产色爽女视频免费观看| 大话2 男鬼变身卡| 97在线视频观看| 色视频www国产| 久久鲁丝午夜福利片| 搡女人真爽免费视频火全软件| videossex国产| 国产成人a区在线观看| 中文天堂在线官网| 中国国产av一级| 国产精品一二三区在线看| 成年版毛片免费区| 国产精品av视频在线免费观看| videossex国产| 亚洲av成人精品一区久久| 久久久午夜欧美精品| 午夜亚洲福利在线播放| 久久久久久久国产电影| 国产精品久久久久久久电影| 亚洲av福利一区| 免费黄网站久久成人精品| 波野结衣二区三区在线| 99久久中文字幕三级久久日本| 亚洲精品日本国产第一区| 欧美高清成人免费视频www| 久久国产乱子免费精品| 能在线免费看毛片的网站| 大又大粗又爽又黄少妇毛片口| 午夜福利在线在线| 色播亚洲综合网| av在线老鸭窝| 午夜福利在线在线| 精品久久久久久久末码| 国产色婷婷99| 午夜福利在线在线| 午夜福利视频精品| 18禁动态无遮挡网站| 欧美成人午夜免费资源| 97在线人人人人妻| 国产亚洲5aaaaa淫片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国语对白做爰xxxⅹ性视频网站| 91精品一卡2卡3卡4卡| 亚洲美女视频黄频| 人妻少妇偷人精品九色| 国产成人精品福利久久| 一边亲一边摸免费视频| 国产亚洲一区二区精品| 国内精品美女久久久久久| 国产成人免费无遮挡视频| 青春草视频在线免费观看| 在线观看一区二区三区| 在线免费十八禁| 伊人久久精品亚洲午夜| 国产成人精品婷婷| 免费在线观看成人毛片| 搞女人的毛片| 国产成人免费观看mmmm| 乱系列少妇在线播放| 在线精品无人区一区二区三 | 亚洲综合色惰| 一本久久精品| 亚洲丝袜综合中文字幕| 91aial.com中文字幕在线观看| 午夜福利视频1000在线观看| 又黄又爽又刺激的免费视频.| 天天躁夜夜躁狠狠久久av| 卡戴珊不雅视频在线播放| 免费人成在线观看视频色| 亚洲精品色激情综合| 老师上课跳d突然被开到最大视频| 九九久久精品国产亚洲av麻豆| 中文字幕免费在线视频6| 欧美日韩精品成人综合77777| 性色av一级| 久久久午夜欧美精品| 亚洲成人久久爱视频| 国产爱豆传媒在线观看| 亚洲欧美日韩另类电影网站 | eeuss影院久久| 超碰97精品在线观看| 黄片无遮挡物在线观看| 日本黄色片子视频| 哪个播放器可以免费观看大片| 天堂俺去俺来也www色官网| 蜜臀久久99精品久久宅男| 久久久久性生活片| 国产一区二区三区av在线| 婷婷色av中文字幕| 国产一区二区在线观看日韩| 在线亚洲精品国产二区图片欧美 | 午夜精品一区二区三区免费看| 各种免费的搞黄视频| 真实男女啪啪啪动态图| 亚洲av福利一区| 成人综合一区亚洲| 日本猛色少妇xxxxx猛交久久| 可以在线观看毛片的网站| 精品久久久久久久久av| 免费av观看视频| 听说在线观看完整版免费高清| 日日摸夜夜添夜夜添av毛片| 在线观看美女被高潮喷水网站| 国产伦理片在线播放av一区| 亚洲av一区综合| 亚洲在久久综合| 午夜免费鲁丝| 成人高潮视频无遮挡免费网站| 久久久久久久精品精品| 亚洲色图av天堂| 国产免费一区二区三区四区乱码| 在线天堂最新版资源| 亚洲精品国产av成人精品| 国产午夜福利久久久久久| 好男人在线观看高清免费视频| 日韩欧美 国产精品| 亚洲国产精品成人久久小说| 久久精品人妻少妇| 啦啦啦啦在线视频资源| 在线观看一区二区三区| 国产综合懂色| 有码 亚洲区| 亚洲人成网站在线观看播放| 激情 狠狠 欧美| 欧美激情国产日韩精品一区| 嫩草影院新地址| 亚洲,欧美,日韩| 麻豆国产97在线/欧美| 男人添女人高潮全过程视频| 亚洲国产精品专区欧美| 久久精品国产亚洲av天美| 七月丁香在线播放| 成年女人在线观看亚洲视频 | 国产精品99久久久久久久久| 嫩草影院新地址| 久久精品国产亚洲av涩爱| 久久久精品欧美日韩精品| 成人综合一区亚洲| 国产色婷婷99| 亚洲精品日韩在线中文字幕| av在线播放精品| 国产亚洲av片在线观看秒播厂| 亚洲精品,欧美精品| 精品人妻熟女av久视频| 黄色日韩在线| 麻豆国产97在线/欧美| 成人二区视频| 男的添女的下面高潮视频| 美女内射精品一级片tv| 人人妻人人看人人澡| 亚洲av免费在线观看| 热re99久久精品国产66热6| 亚洲电影在线观看av| 亚洲精品久久午夜乱码| 黄色欧美视频在线观看| 日韩不卡一区二区三区视频在线| 国产视频内射| 纵有疾风起免费观看全集完整版| 久久韩国三级中文字幕| 熟女av电影| 一边亲一边摸免费视频| 久久精品国产亚洲av天美| 欧美日韩在线观看h| 97精品久久久久久久久久精品| 熟女电影av网| 国产av国产精品国产| 午夜视频国产福利| 久久精品久久精品一区二区三区| 免费av不卡在线播放| 高清欧美精品videossex| 男的添女的下面高潮视频| 99热这里只有是精品在线观看| 欧美日韩精品成人综合77777| 伊人久久国产一区二区| 麻豆国产97在线/欧美| 色5月婷婷丁香| 丝袜美腿在线中文| 天堂中文最新版在线下载 | 欧美老熟妇乱子伦牲交| 久久韩国三级中文字幕| 国产免费一级a男人的天堂| 韩国av在线不卡| 欧美bdsm另类| 日本一本二区三区精品| 久久久精品94久久精品| 岛国毛片在线播放| 美女cb高潮喷水在线观看| 亚洲激情五月婷婷啪啪| 成人午夜精彩视频在线观看| 搡老乐熟女国产| 欧美激情久久久久久爽电影| 毛片女人毛片| 婷婷色麻豆天堂久久| 在线亚洲精品国产二区图片欧美 | 免费看不卡的av| 一区二区三区四区激情视频| 纵有疾风起免费观看全集完整版| 韩国av在线不卡| 欧美bdsm另类| 欧美一级a爱片免费观看看| a级毛片免费高清观看在线播放| 国产精品99久久久久久久久| 精品久久久精品久久久| 国产午夜福利久久久久久| 亚洲精品成人av观看孕妇| 亚洲成人中文字幕在线播放| 国产高清三级在线| 国产精品精品国产色婷婷| 一个人看的www免费观看视频| 国产老妇伦熟女老妇高清| 亚洲精品,欧美精品| 日本三级黄在线观看| 2018国产大陆天天弄谢| 国产探花极品一区二区| 高清欧美精品videossex| freevideosex欧美| 别揉我奶头 嗯啊视频| 国产男人的电影天堂91| 亚洲性久久影院| 国模一区二区三区四区视频| 国产爱豆传媒在线观看| 中文字幕亚洲精品专区| 又粗又硬又长又爽又黄的视频| 国产伦在线观看视频一区| 97热精品久久久久久| 久久精品人妻少妇| 亚洲av不卡在线观看| 大码成人一级视频| 91狼人影院| 午夜老司机福利剧场| 亚洲一区二区三区欧美精品 | xxx大片免费视频| av在线亚洲专区| 性色avwww在线观看| 亚洲精品视频女| 欧美日韩亚洲高清精品| 各种免费的搞黄视频| 日韩成人伦理影院| 欧美日韩国产mv在线观看视频 | 在线观看国产h片| 欧美最新免费一区二区三区| 国产精品99久久99久久久不卡 | 免费人成在线观看视频色| 久久久久久国产a免费观看| 99久久精品一区二区三区| 亚洲国产欧美人成| 国产综合懂色| 国产成人精品婷婷| 久久久精品免费免费高清| 视频区图区小说| 晚上一个人看的免费电影| 国产爱豆传媒在线观看| 亚洲欧美清纯卡通| 两个人的视频大全免费| 美女cb高潮喷水在线观看| 亚洲欧美精品自产自拍| 好男人视频免费观看在线| 少妇高潮的动态图| 一二三四中文在线观看免费高清| 人人妻人人澡人人爽人人夜夜| 精品久久久精品久久久| 天堂中文最新版在线下载 | 亚洲av中文字字幕乱码综合| 狂野欧美激情性xxxx在线观看| 看黄色毛片网站| 国产高清国产精品国产三级 | 夫妻性生交免费视频一级片| 能在线免费看毛片的网站| 成年女人看的毛片在线观看| 久久久久精品性色| 日韩人妻高清精品专区| 人妻系列 视频| av免费观看日本| 久久精品国产亚洲av天美| 高清日韩中文字幕在线| 国产成人精品一,二区| 有码 亚洲区| 欧美高清成人免费视频www| 久久综合国产亚洲精品| 国产男女超爽视频在线观看| 建设人人有责人人尽责人人享有的 | 神马国产精品三级电影在线观看| 中文字幕人妻熟人妻熟丝袜美| 我要看日韩黄色一级片| 国产av码专区亚洲av| 成人鲁丝片一二三区免费| 国产精品人妻久久久影院| 国产欧美亚洲国产| 亚洲欧美一区二区三区国产| 久久久久久伊人网av| 国产伦在线观看视频一区| 一本一本综合久久| 在线观看一区二区三区| 免费看不卡的av| 亚洲欧美清纯卡通| 亚洲精品成人久久久久久| 日韩电影二区| 美女内射精品一级片tv| 2021天堂中文幕一二区在线观| 亚洲人成网站高清观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品乱码久久久久久按摩| 中文字幕人妻熟人妻熟丝袜美| 在线观看美女被高潮喷水网站| 国产女主播在线喷水免费视频网站| 蜜桃久久精品国产亚洲av| 久久久精品欧美日韩精品| 亚洲第一区二区三区不卡| 国产欧美日韩一区二区三区在线 | 亚洲成色77777| 国产精品一二三区在线看| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美人与善性xxx| 国产精品嫩草影院av在线观看| 视频区图区小说| 中国美白少妇内射xxxbb| 日本欧美国产在线视频| 精品国产三级普通话版| 日韩成人伦理影院| 国产老妇伦熟女老妇高清| 草草在线视频免费看| 欧美日本视频| 国产片特级美女逼逼视频| 精品久久久久久电影网| 色5月婷婷丁香| 久久久久精品性色| 亚洲av男天堂| 国产欧美亚洲国产| 免费在线观看成人毛片| 中文精品一卡2卡3卡4更新| 亚洲经典国产精华液单| 啦啦啦啦在线视频资源| 男女边吃奶边做爰视频| 99精国产麻豆久久婷婷| 特大巨黑吊av在线直播|