• <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不卡| 99久国产av精品国产电影| 丰满饥渴人妻一区二区三| 国产精品麻豆人妻色哟哟久久| 久久人人爽人人爽人人片va| 9色porny在线观看| 王馨瑶露胸无遮挡在线观看| 18禁观看日本| 岛国毛片在线播放| 69精品国产乱码久久久| 成人影院久久| 国产老妇伦熟女老妇高清| 乱人伦中国视频| 中国美白少妇内射xxxbb| 91久久精品电影网| 大码成人一级视频| 成人无遮挡网站| 中文字幕免费在线视频6| 亚洲,欧美,日韩| 2022亚洲国产成人精品| 亚洲综合色网址| 国产精品久久久久成人av| 十八禁高潮呻吟视频| 国产白丝娇喘喷水9色精品| 男人操女人黄网站| 视频中文字幕在线观看| 丰满饥渴人妻一区二区三| 在线观看免费视频网站a站| 日韩精品免费视频一区二区三区 | 秋霞伦理黄片| 免费人成在线观看视频色| 老女人水多毛片| 大又大粗又爽又黄少妇毛片口| 新久久久久国产一级毛片| 久久精品国产a三级三级三级| 免费黄网站久久成人精品| 在线免费观看不下载黄p国产| 在线观看www视频免费| 久久精品国产a三级三级三级| 大又大粗又爽又黄少妇毛片口| 亚洲欧美成人精品一区二区| 亚洲性久久影院| 亚洲精品视频女| 色哟哟·www| 下体分泌物呈黄色| 精品视频人人做人人爽| 22中文网久久字幕| 激情五月婷婷亚洲| 亚洲五月色婷婷综合| 两个人免费观看高清视频| 精品人妻偷拍中文字幕| 一级毛片电影观看| 一级黄片播放器| 亚洲图色成人| 黑丝袜美女国产一区| 女性生殖器流出的白浆| 日韩av不卡免费在线播放| 最新中文字幕久久久久| 久久综合国产亚洲精品| 国产精品99久久99久久久不卡 | 18+在线观看网站| 99国产综合亚洲精品| 国产黄频视频在线观看| 久久久久久久久大av| xxxhd国产人妻xxx| 一区二区三区四区激情视频| 久久99热这里只频精品6学生| 精品少妇内射三级| 国产日韩欧美视频二区| 狠狠精品人妻久久久久久综合| 日本91视频免费播放| 精品人妻在线不人妻| 日韩熟女老妇一区二区性免费视频| 男女啪啪激烈高潮av片| 精品国产乱码久久久久久小说| 亚洲精品色激情综合| 精品少妇久久久久久888优播| 我要看黄色一级片免费的| 最近2019中文字幕mv第一页| 亚洲国产精品国产精品| 日韩成人伦理影院| 三级国产精品欧美在线观看| 日韩中字成人| 久久午夜综合久久蜜桃| 久久午夜福利片| 麻豆乱淫一区二区| 亚洲精品456在线播放app| 日本欧美视频一区| 成人免费观看视频高清| 在线天堂最新版资源| freevideosex欧美| 99热网站在线观看| 啦啦啦中文免费视频观看日本| 亚洲欧美日韩另类电影网站| 国产一区有黄有色的免费视频| 日韩一区二区视频免费看| 亚洲国产精品999| 免费高清在线观看日韩| 2022亚洲国产成人精品| 桃花免费在线播放| 一本色道久久久久久精品综合| 免费av中文字幕在线| 国产精品久久久久久av不卡| 国产精品一二三区在线看| 又黄又爽又刺激的免费视频.| 午夜视频国产福利| 九草在线视频观看| 国产女主播在线喷水免费视频网站| 91国产中文字幕| 18禁裸乳无遮挡动漫免费视频| 日韩成人av中文字幕在线观看| 久久韩国三级中文字幕| 久久 成人 亚洲| 9色porny在线观看| 欧美国产精品一级二级三级| 黑人高潮一二区| 成人18禁高潮啪啪吃奶动态图 | 国产高清三级在线| 亚洲成色77777| 欧美日韩在线观看h| 91久久精品国产一区二区三区| 人妻人人澡人人爽人人| 内地一区二区视频在线| 性高湖久久久久久久久免费观看| 中文天堂在线官网| 女人精品久久久久毛片| 在线观看免费高清a一片| 国产精品一二三区在线看| 亚洲欧美成人精品一区二区| 2022亚洲国产成人精品| av又黄又爽大尺度在线免费看| 日日摸夜夜添夜夜爱| 婷婷色av中文字幕| 中国美白少妇内射xxxbb| 欧美最新免费一区二区三区| 久久国产精品大桥未久av| 五月开心婷婷网| 国产成人精品福利久久| 夫妻性生交免费视频一级片| 插逼视频在线观看| 日本色播在线视频| 午夜激情久久久久久久| 9色porny在线观看| 国产黄片视频在线免费观看| 丰满少妇做爰视频| 久久99热6这里只有精品| 欧美变态另类bdsm刘玥| 高清黄色对白视频在线免费看| 亚洲精品乱久久久久久| 国产精品熟女久久久久浪| 国产乱人偷精品视频| 国产在线免费精品| 观看av在线不卡| 一区二区三区四区激情视频| 婷婷色av中文字幕| 久久av网站| 蜜臀久久99精品久久宅男| 国产亚洲午夜精品一区二区久久| 夜夜骑夜夜射夜夜干| 精品99又大又爽又粗少妇毛片| 两个人免费观看高清视频| 狠狠婷婷综合久久久久久88av| 亚洲精华国产精华液的使用体验| 国产熟女欧美一区二区| 久久久国产一区二区| 久久精品久久精品一区二区三区| 欧美成人精品欧美一级黄| 欧美+日韩+精品| 亚洲激情五月婷婷啪啪| 狂野欧美激情性xxxx在线观看| 亚洲婷婷狠狠爱综合网| 女的被弄到高潮叫床怎么办| 麻豆成人av视频| a级毛片在线看网站| 亚洲国产欧美在线一区| 大片免费播放器 马上看| 久久人人爽人人片av| 两个人免费观看高清视频| 国产免费福利视频在线观看| 哪个播放器可以免费观看大片| 99九九线精品视频在线观看视频| 亚洲精品日本国产第一区| 久久精品国产鲁丝片午夜精品| 99re6热这里在线精品视频| 最新的欧美精品一区二区| 色哟哟·www| 2021少妇久久久久久久久久久| 亚洲精品乱码久久久v下载方式| 波野结衣二区三区在线| 亚洲美女视频黄频| 国产精品久久久久久av不卡| 国产精品一二三区在线看| 赤兔流量卡办理| 亚洲av日韩在线播放| 国产精品秋霞免费鲁丝片| 成人影院久久| 少妇人妻精品综合一区二区| 久久久久久人妻| 九草在线视频观看| 日本色播在线视频| 91久久精品国产一区二区三区| 夜夜骑夜夜射夜夜干| 亚洲精品中文字幕在线视频| 一级爰片在线观看| 久久精品人人爽人人爽视色| 交换朋友夫妻互换小说| 少妇人妻 视频| 一级片'在线观看视频| 久久国产亚洲av麻豆专区| 青春草亚洲视频在线观看| 色哟哟·www| 丝袜美足系列| 日韩一区二区三区影片| 欧美成人午夜免费资源| 中国美白少妇内射xxxbb| 狠狠婷婷综合久久久久久88av| 国产欧美日韩综合在线一区二区| 黄片播放在线免费| 晚上一个人看的免费电影| 日韩电影二区| 大香蕉久久网| 亚洲欧洲国产日韩| 日韩成人av中文字幕在线观看| 另类精品久久| 久久婷婷青草| 成人毛片60女人毛片免费| 亚洲精品日韩av片在线观看| a级毛片免费高清观看在线播放| 大香蕉97超碰在线| 一个人免费看片子| 国产熟女欧美一区二区| 日韩中字成人| 色吧在线观看| 久久久久人妻精品一区果冻| 国产精品熟女久久久久浪| 亚洲国产毛片av蜜桃av| 亚洲第一区二区三区不卡| 精品一区二区免费观看| 久久 成人 亚洲| 亚洲国产成人一精品久久久| 黄色配什么色好看| 少妇被粗大猛烈的视频| 日本爱情动作片www.在线观看| 超碰97精品在线观看| 国产成人freesex在线| 视频在线观看一区二区三区| 国产亚洲最大av| 国产亚洲一区二区精品| 欧美变态另类bdsm刘玥| 久久人妻熟女aⅴ| 久久久久视频综合| 能在线免费看毛片的网站| 黑人巨大精品欧美一区二区蜜桃 | 久热这里只有精品99| 国产永久视频网站| 亚洲综合色网址| av不卡在线播放| 国产日韩欧美在线精品| 男的添女的下面高潮视频| 亚洲av福利一区| 国内精品宾馆在线| 久久久a久久爽久久v久久| av有码第一页| 久久99精品国语久久久| 中文字幕最新亚洲高清| 人妻系列 视频| 欧美变态另类bdsm刘玥| 久久热精品热| 欧美另类一区| 99热全是精品| 18+在线观看网站| 亚洲精品,欧美精品| 久久综合国产亚洲精品| 亚洲人成网站在线播| 欧美老熟妇乱子伦牲交| 国产永久视频网站| 制服人妻中文乱码| 丰满饥渴人妻一区二区三| 99国产精品免费福利视频| av播播在线观看一区| 国产成人a∨麻豆精品| 国产精品久久久久久精品古装| 91精品三级在线观看| 欧美人与性动交α欧美精品济南到 | 啦啦啦中文免费视频观看日本| 国产乱来视频区| 伊人亚洲综合成人网| 久久精品人人爽人人爽视色| 老女人水多毛片| 九草在线视频观看| 日韩av在线免费看完整版不卡| 成人国产av品久久久| 免费高清在线观看日韩| 日韩三级伦理在线观看| 久久精品夜色国产| 日韩成人伦理影院| 人人妻人人澡人人看| 老司机影院成人| 成人午夜精彩视频在线观看| 五月伊人婷婷丁香| xxxhd国产人妻xxx| 亚洲av在线观看美女高潮| 免费观看av网站的网址| 日韩伦理黄色片| 丰满乱子伦码专区| av电影中文网址| 999精品在线视频| 中国美白少妇内射xxxbb| 国产极品天堂在线| 国产精品人妻久久久影院| 欧美xxxx性猛交bbbb| 尾随美女入室| 午夜福利网站1000一区二区三区| 成人国产av品久久久| 特大巨黑吊av在线直播| 国产熟女欧美一区二区| 久久人人爽av亚洲精品天堂| 欧美成人午夜免费资源| 成年美女黄网站色视频大全免费 | av视频免费观看在线观看| 欧美日韩成人在线一区二区| xxxhd国产人妻xxx| 欧美精品高潮呻吟av久久| 男女边摸边吃奶| 有码 亚洲区| 91精品国产九色| 插逼视频在线观看| 一级毛片电影观看| 国产亚洲一区二区精品| 大香蕉久久网| 18禁裸乳无遮挡动漫免费视频| 97精品久久久久久久久久精品| 日韩熟女老妇一区二区性免费视频| 18禁在线播放成人免费| 国产午夜精品一二区理论片| 欧美另类一区| 国产精品国产av在线观看| 日韩av免费高清视频| 精品久久久噜噜| 97超视频在线观看视频| 午夜影院在线不卡| 亚洲色图 男人天堂 中文字幕 | 成人二区视频| 看免费成人av毛片| 亚洲国产av影院在线观看| 少妇猛男粗大的猛烈进出视频| 新久久久久国产一级毛片| 日韩av不卡免费在线播放| 人人妻人人添人人爽欧美一区卜| av黄色大香蕉| 亚洲国产av新网站| 晚上一个人看的免费电影| 成人毛片a级毛片在线播放| 又粗又硬又长又爽又黄的视频| 国产精品三级大全| 婷婷成人精品国产| 99国产综合亚洲精品| 国产男女内射视频| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久久久免| 久久99一区二区三区| 久久97久久精品| 国产精品国产av在线观看| 爱豆传媒免费全集在线观看| 狠狠精品人妻久久久久久综合| 交换朋友夫妻互换小说| 午夜福利视频精品| kizo精华| 国产欧美另类精品又又久久亚洲欧美| 香蕉精品网在线| 久久热精品热| 最后的刺客免费高清国语| 欧美+日韩+精品| 国产男人的电影天堂91| 超碰97精品在线观看| 男女国产视频网站| 亚洲精品日本国产第一区| 日日摸夜夜添夜夜添av毛片| 七月丁香在线播放| 久久久精品区二区三区| 亚洲国产精品一区二区三区在线| xxxhd国产人妻xxx| 人妻 亚洲 视频| 桃花免费在线播放| 91精品国产九色| 精品人妻偷拍中文字幕| 日韩欧美精品免费久久| 精品少妇黑人巨大在线播放| 一本—道久久a久久精品蜜桃钙片| 欧美丝袜亚洲另类| 欧美精品一区二区大全| 国模一区二区三区四区视频| 欧美+日韩+精品| 亚洲内射少妇av| 高清在线视频一区二区三区| 国产午夜精品久久久久久一区二区三区| 色婷婷久久久亚洲欧美| 亚洲av二区三区四区| 一级二级三级毛片免费看| 日韩视频在线欧美| 日韩不卡一区二区三区视频在线| 国产无遮挡羞羞视频在线观看| 永久免费av网站大全| videossex国产| 亚洲欧美成人精品一区二区| 国产精品人妻久久久影院| 国产精品蜜桃在线观看| 国产免费一级a男人的天堂| 日日摸夜夜添夜夜爱| a级片在线免费高清观看视频| 精品国产乱码久久久久久小说| 国产精品三级大全| 日韩 亚洲 欧美在线| 高清午夜精品一区二区三区| 91精品国产国语对白视频| 18+在线观看网站| 国产亚洲欧美精品永久| 精品人妻一区二区三区麻豆| 日韩中字成人| 久久精品国产亚洲网站| 成人亚洲欧美一区二区av| 国产精品一区二区在线不卡| 亚洲少妇的诱惑av| av不卡在线播放| av有码第一页| 美女cb高潮喷水在线观看| 两个人的视频大全免费| 韩国av在线不卡| 91aial.com中文字幕在线观看| 中文字幕制服av| 美女国产视频在线观看| 久久狼人影院| 亚洲情色 制服丝袜| 免费看光身美女| 好男人视频免费观看在线| 狠狠婷婷综合久久久久久88av| 国产精品久久久久久精品电影小说| 亚洲精品日韩在线中文字幕| 亚洲av综合色区一区| 亚洲欧美成人综合另类久久久| 91成人精品电影| 少妇的逼好多水| a级片在线免费高清观看视频| 久久青草综合色| 黑丝袜美女国产一区| 女性被躁到高潮视频| 精品人妻偷拍中文字幕| 国产精品久久久久久久久免| 中国三级夫妇交换| 国产成人午夜福利电影在线观看| 日日爽夜夜爽网站| 亚洲欧美一区二区三区黑人 | 男女国产视频网站| 九草在线视频观看| 99久国产av精品国产电影| 亚洲精品视频女| 亚洲国产欧美日韩在线播放| 99热网站在线观看| 国产av一区二区精品久久| 黄色怎么调成土黄色| 老司机影院毛片| 曰老女人黄片| 18禁动态无遮挡网站| 一区二区三区精品91| 自线自在国产av| 777米奇影视久久| 亚洲欧美成人精品一区二区| 毛片一级片免费看久久久久| 三级国产精品欧美在线观看| 韩国高清视频一区二区三区| 色婷婷av一区二区三区视频| 久久这里有精品视频免费| 国产成人av激情在线播放 | 观看美女的网站| 色网站视频免费| 水蜜桃什么品种好| 久久精品国产自在天天线| 欧美精品高潮呻吟av久久| 狂野欧美激情性xxxx在线观看| 精品酒店卫生间| 亚洲av免费高清在线观看| 蜜臀久久99精品久久宅男| 精品人妻熟女毛片av久久网站| 国产免费视频播放在线视频| 成人毛片60女人毛片免费| 欧美丝袜亚洲另类| 最近最新中文字幕免费大全7| 在线观看免费视频网站a站| 天美传媒精品一区二区| 狠狠婷婷综合久久久久久88av| 亚洲五月色婷婷综合| 日韩精品有码人妻一区| 成年女人在线观看亚洲视频| 一区二区三区四区激情视频| 久久99一区二区三区| 99热这里只有是精品在线观看| 国产国语露脸激情在线看| 欧美一级a爱片免费观看看| 搡老乐熟女国产| 女人久久www免费人成看片| 亚洲精品日本国产第一区| 国产一区二区三区av在线| 18禁动态无遮挡网站| 少妇被粗大猛烈的视频| 中文字幕制服av| 亚洲精品久久成人aⅴ小说 | 国产亚洲午夜精品一区二区久久| 亚洲av成人精品一区久久| 少妇人妻久久综合中文| 纯流量卡能插随身wifi吗| 午夜av观看不卡| 欧美 日韩 精品 国产| 久久精品国产自在天天线| 午夜精品国产一区二区电影| 国产 一区精品| 丝袜在线中文字幕| 最黄视频免费看| 精品99又大又爽又粗少妇毛片| 男女边摸边吃奶| 国产爽快片一区二区三区| 精品久久久久久久久av| 国产精品.久久久| 在线 av 中文字幕| 五月伊人婷婷丁香| 久久婷婷青草| av卡一久久| 亚洲精品亚洲一区二区| 欧美成人午夜免费资源| 男女国产视频网站| 大陆偷拍与自拍| 三级国产精品片| 色吧在线观看| 99国产精品免费福利视频| 2021少妇久久久久久久久久久| 欧美日韩在线观看h| 午夜影院在线不卡| 啦啦啦中文免费视频观看日本| 人妻夜夜爽99麻豆av| 欧美日韩亚洲高清精品| 高清av免费在线| 黄色配什么色好看| 日韩 亚洲 欧美在线| 国产精品秋霞免费鲁丝片| 另类亚洲欧美激情| 午夜老司机福利剧场| 国产精品偷伦视频观看了| 国产日韩欧美亚洲二区| 国产片内射在线| 国产精品不卡视频一区二区| 日本爱情动作片www.在线观看| 日韩一区二区视频免费看| 最近中文字幕2019免费版| 美女福利国产在线| 国产高清不卡午夜福利| av电影中文网址| 三级国产精品欧美在线观看| 久久久精品94久久精品| 久久狼人影院| 黑人猛操日本美女一级片| 少妇的逼水好多| 亚洲国产精品一区三区| 热99久久久久精品小说推荐| 成人漫画全彩无遮挡| www.av在线官网国产| 大片免费播放器 马上看| 亚洲欧美一区二区三区黑人 | 我的老师免费观看完整版| 精品人妻熟女av久视频| av播播在线观看一区| 三上悠亚av全集在线观看| 亚洲av男天堂| h视频一区二区三区| 黄片播放在线免费| 伊人久久精品亚洲午夜| 精品一区二区三区视频在线| 一区二区av电影网| 亚洲无线观看免费| 国产片特级美女逼逼视频| 日韩精品有码人妻一区| 亚洲欧美成人精品一区二区| 久久国内精品自在自线图片| 精品国产乱码久久久久久小说| 亚洲人成网站在线观看播放| 婷婷色综合www| 精品久久久久久久久亚洲| 考比视频在线观看| 午夜福利影视在线免费观看| 观看av在线不卡| 欧美精品高潮呻吟av久久| 免费av中文字幕在线| 国产高清三级在线| 国产熟女欧美一区二区| 成人无遮挡网站| 欧美精品一区二区大全| 婷婷成人精品国产| 亚洲av.av天堂| 秋霞伦理黄片| 人人澡人人妻人| 女性生殖器流出的白浆| 最后的刺客免费高清国语| 我的老师免费观看完整版| 一本一本综合久久| 大香蕉久久成人网| 九九爱精品视频在线观看| 99久久中文字幕三级久久日本| 久久婷婷青草| 成人漫画全彩无遮挡| 精品酒店卫生间| 免费观看的影片在线观看| 男女边摸边吃奶| 欧美成人精品欧美一级黄| 亚洲伊人久久精品综合| 日韩精品免费视频一区二区三区 | 日产精品乱码卡一卡2卡三|