程 勃，丁彥禮，徐志鋒

（桂林理工大學(xué)　a.廣西隱伏金屬礦產(chǎn)勘查重點(diǎn)實(shí)驗(yàn)室；b.地球科學(xué)學(xué)院，廣西桂林 541004）

統(tǒng)計(jì)學(xué)建模法在面積性高密度電阻率數(shù)據(jù)反演中的應(yīng)用

程 勃，丁彥禮，徐志鋒

（桂林理工大學(xué)a.廣西隱伏金屬礦產(chǎn)勘查重點(diǎn)實(shí)驗(yàn)室；b.地球科學(xué)學(xué)院，廣西桂林 541004）

在面積性的勘查工作中，為了能夠通過高密度電阻率數(shù)據(jù)反演獲得詳細(xì)的地層結(jié)構(gòu)和巖性信息，采用了基于統(tǒng)計(jì)學(xué)方法建立初始模型的反演算法。該方法充分地利用高密度實(shí)測數(shù)據(jù)建立初始模型，對單個剖面實(shí)施遺傳算法反演，將其反演結(jié)果作為初始模型用于相鄰剖面數(shù)據(jù)的反演。實(shí)例表明，該方法大幅度降低了初次擬合誤差，減少了反演計(jì)算量，使反演結(jié)果較好地反映出測區(qū)范圍內(nèi)的地層結(jié)構(gòu)特征和巖性分布情況。

電阻率法；二維反演；統(tǒng)計(jì)學(xué)建模；天然氣管道

0 引 言

高密度電阻率法廣泛應(yīng)用于工程地質(zhì)勘查工作中，主要解決分層及探測地質(zhì)構(gòu)造等問題。由于高密度電阻率法同時具有電阻率剖面法和測深法的特點(diǎn)，在工程場地勘查中往往進(jìn)行包含多條平行測線的面積性測量，人們期待高密度電阻率數(shù)據(jù)反演能夠給出盡可能詳盡的三維地下結(jié)構(gòu)特征［1-3］。

由于三維反演計(jì)算時間長，還不能廣泛用于實(shí)際工作［4-6］。對于多條測線的面積性高密度電阻率測量，還是以二維反演為主。在常規(guī)二維反演解釋過程中，每條測線是獨(dú)立進(jìn)行解釋的，反演圖像是電阻率連續(xù)變化的等值線圖，是地下介質(zhì)的電性成像，沒有明確的地層結(jié)構(gòu)信息，需要數(shù)據(jù)處理人員根據(jù)經(jīng)驗(yàn)分析解釋出地層結(jié)構(gòu)、分層和電阻率信息［7-9］。有時，測線間的距離不大，解釋軟件給出的不同測線地電斷面電性差異卻比較大，增大了解釋人員綜合分析推斷的難度［10-12］。

統(tǒng)計(jì)學(xué)建模反演方法［1-3］以電阻率測深數(shù)據(jù)的統(tǒng)計(jì)學(xué)特征建立二維地層結(jié)構(gòu)初始模型，這種模型具有地下的地層結(jié)構(gòu)和電性參數(shù)。利用改進(jìn)的遺傳算法修改地層結(jié)構(gòu)和層參數(shù)，并在修改模型的反演過程中充分考慮先驗(yàn)信息。這樣反演的最終結(jié)果具有明確的地層結(jié)構(gòu)和具體的層參數(shù)信息，并可以直接與地質(zhì)構(gòu)造、巖性等信息對應(yīng)。在面積性測量時，測區(qū)內(nèi)包含多條距離較小且相鄰的剖面，根據(jù)常理和經(jīng)驗(yàn)可以假設(shè)，相鄰剖面的地層結(jié)構(gòu)和巖性有相似之處。因此，通過處理一個剖面資料獲得的反演地層結(jié)構(gòu)和層參數(shù)可以推廣到相鄰剖面上。

本文嘗試在面積性勘查的高密度電法數(shù)據(jù)處理中使用統(tǒng)計(jì)學(xué)建模方法反演，利用原始數(shù)據(jù)的統(tǒng)計(jì)學(xué)信息并結(jié)合地質(zhì)資料，建立最貼近實(shí)際的初始模型。用多剖面統(tǒng)計(jì)學(xué)特征值判斷地層結(jié)構(gòu)、獲得層參數(shù)并建立二維地質(zhì)模型。對某一剖面進(jìn)行二維反演，完成單個剖面的反演后，使用反演結(jié)果作為相鄰剖面的反演初始模型，對其他剖面反演。這樣既提高了反演速度，又控制了區(qū)域內(nèi)的模型參數(shù)差異。這一處理方法與常規(guī)方法相比的優(yōu)點(diǎn)是：反演結(jié)果有明確的地層結(jié)構(gòu)和具體的層參數(shù)信息，可以直接與地質(zhì)構(gòu)造、巖性等信息對應(yīng)，在測區(qū)內(nèi)地質(zhì)條件變化不大、相鄰測線高密度電阻率異常特征相似的情況下，將某一剖面的反演結(jié)果用于相鄰測線，既提高了反演速度，又可以使測區(qū)的地層結(jié)構(gòu)和層參數(shù)具有統(tǒng)一性，利于進(jìn)一步的地質(zhì)解釋。值得注意的是，當(dāng)測區(qū)內(nèi)地質(zhì)條件變化大、電阻率異常特征差異大時，需要解釋人員先將測線分區(qū)，然后再進(jìn)行反演解釋。

本文用統(tǒng)計(jì)學(xué)建模方法對某天然氣管道基礎(chǔ)勘探測區(qū)3條測線的高密度電阻率數(shù)據(jù)進(jìn)行反演。反演出的測區(qū)勘探深度范圍內(nèi)的地層結(jié)構(gòu)特征和巖性分布狀態(tài)，將高密度電法反演解釋結(jié)果直接與地層結(jié)構(gòu)聯(lián)系起來，提高了資料解釋的準(zhǔn)確性，增強(qiáng)高密度電法勘探的準(zhǔn)確性和實(shí)用性，得到了較好的地質(zhì)效果。

1　統(tǒng)計(jì)學(xué)建模反演方法和步驟

1.1統(tǒng)計(jì)學(xué)建模流程

統(tǒng)計(jì)學(xué)方法建模流程如圖1所示，主要步驟分為：（1）處理實(shí)測數(shù)據(jù)并判斷地層結(jié)構(gòu)類型（圖1主干部分）；（2）確定地質(zhì)體空間位置（圖1右側(cè)部分）；（3）加入先驗(yàn)信息（圖1左側(cè)部分）；（4）建立初始模型。

1.2統(tǒng)計(jì)權(quán)重的分配

在圖1中，當(dāng)模型類別被判斷為含橫向變化（如存在斷裂）時，需要確認(rèn)電性橫向變化的位置及兩側(cè)的地層結(jié)構(gòu)和層參數(shù)。為此，統(tǒng)計(jì)電阻率測深數(shù)據(jù)各極距視電阻率的水平方向梯度。水平梯度的極值位置信息可以輔助完成初始模型結(jié)構(gòu)類別判斷，并獲得斷裂或巖性接觸面的初始位置［1-2］。數(shù)值模擬結(jié)果和實(shí)測資料證明，淺部的局部異常體同樣會造成水平梯度出現(xiàn)極值，可能造成對主要異常體位置判斷失誤。因此在水平梯度統(tǒng)計(jì)過程中根據(jù)探測目的設(shè)置統(tǒng)計(jì)權(quán)重，例如，在探測斷裂、巖性接觸帶或深部異常體時，對深部信息賦予較大的統(tǒng)計(jì)權(quán)重。

圖1　統(tǒng)計(jì)學(xué)建模流程Fig.1 Process of statistic modeling

以低阻斷裂為例，斷裂左邊界確定函數(shù)為

斷裂右邊界確定函數(shù)為

其中：Δρi，j/Δi為視電阻率在水平方向上的變化率，亦為權(quán)重因子；i為測點(diǎn)號；j為極距號。P（x）的作用是求x在測線上的位置。

如果是巖性接觸帶類模型，則視接觸面兩側(cè)的電阻率大小關(guān)系采用式（1）和式（2）中的一個。

1.3測區(qū)內(nèi)電阻率參數(shù)的確定

對地層中含有橫向電性變化的情況下，初次統(tǒng)計(jì)的層參數(shù)數(shù)學(xué)期望只能作為判斷的信息，不能作為初始模型層參數(shù)初值。利用式（1）、（2）獲得各剖面的橫向電性變化位置后，分別統(tǒng)計(jì)橫向電性變化位置（斷裂或接觸面）兩側(cè)測點(diǎn)的一維反演層參數(shù)、測深曲線類型、同類型曲線的數(shù)量及所占總曲線數(shù)量的份額。由此，以占較大份額的曲線類型作為某一側(cè)的初始地層電性結(jié)構(gòu)。

1.4先驗(yàn)信息的加入

測區(qū)的先驗(yàn)性信息，包括地質(zhì)構(gòu)造、地層巖性及地球物理信息，是地球物理勘探資料解釋中必須參考的信息。結(jié)合先驗(yàn)信息判斷地層結(jié)構(gòu)，能降低局部和淺部異常體對整體地層結(jié)構(gòu)判斷的干擾。

在以尋找斷裂為目標(biāo)的面積性勘查工作中，測線布置往往垂直斷裂走向，因此測區(qū)內(nèi)間隔不太遠(yuǎn)且平行的測線很可能位于同種類型的地層結(jié)構(gòu)之上，而且很多模型參數(shù)都是相互近似的。測區(qū)內(nèi)的鉆孔資料提供了地層結(jié)構(gòu)和巖性信息，加入這些信息，使初始模型更接近實(shí)際的地電結(jié)構(gòu)。在反演過程中，由先驗(yàn)信息設(shè)置的參數(shù)不參與遺傳變異，子模型數(shù)量減少，對提高反演速度也有一定的作用。

1.5修正的遺傳算法修改層參數(shù)

建立初始模型之后，采用有限元方法計(jì)算模型的視電阻率，根據(jù)模型計(jì)算值和實(shí)測值的差異，采用修正的遺傳算法修改模型，反演出地層結(jié)構(gòu)參數(shù)和層參數(shù)［1-3］。

遺傳算法是非線性反演中的常用方法，遺傳算法遺傳、變異樣本并正演計(jì)算，需要的計(jì)算時間與模型數(shù)量呈正比。利用統(tǒng)計(jì)學(xué)建模方法已經(jīng)設(shè)定的初始模型，在繁殖子模型的過程中，限制了參數(shù)的變異方式：地層或地質(zhì)體電阻率的變異，每次變異的幅度為當(dāng)前值的10%；空間位置類參數(shù)，包括斷裂的左右邊界、接觸面的位置、地層的厚度等這些參數(shù)每次改變量為一個單元網(wǎng)格（有限單元法正演中的剖分網(wǎng)格）。

在修改參數(shù)的過程中，首先調(diào)整電阻率參數(shù)，在初始模型的基礎(chǔ)上于10%的范圍內(nèi)搜索更接近于實(shí)際模型的層參數(shù)，然后固定電阻率參數(shù)，根據(jù)局部誤差修改各點(diǎn)的層厚度。這樣可以節(jié)省計(jì)算時間，又不至于大幅度修改初始模型。

這樣的修改參數(shù)方式相當(dāng)于搜索法，從初始模型開始，繁殖出參數(shù)變異的子模型，每個子模型每次僅有1個參數(shù)發(fā)生變異。因?yàn)槊總€參數(shù)有兩種變異方式（增加/減少，提高/降低，左移/右移），所以每次繁殖的子模型數(shù)量是參數(shù)個數(shù)的2倍。對這些模型進(jìn)行有限元數(shù)值模擬并計(jì)算擬合誤差。誤差最小的子模型存活并成為新的繁殖起點(diǎn)，開始新一輪的繁殖和篩選。這樣循環(huán)若干次，直至擬合誤差無法降低時反演結(jié)束，最終存活的模型為反演結(jié)果。這種修正的遺傳算法，在修改初始模型參數(shù)的同時，不修改初始模型的地層類型，減少了計(jì)算量。

當(dāng)需要解釋的資料包含多個相鄰剖面時，先采用單剖面反演，反演結(jié)果作為相鄰剖面的反演初始模型，這樣做的好處是合理的降低了反演計(jì)算量，保持相鄰剖面結(jié)構(gòu)統(tǒng)一以及主要參數(shù)接近，便于在多剖面反演結(jié)束時對測區(qū)的三維結(jié)構(gòu)進(jìn)行推斷。

2　應(yīng)用實(shí)例

廣西天然氣輸氣管道工程跨越數(shù)百千米，沿途地質(zhì)條件復(fù)雜，斷裂或巖性接觸帶附近是需要重點(diǎn)勘探的區(qū)域。天然氣管道穿越柳州市雒容鎮(zhèn)附近的洛清江。洛清江西岸測線范圍內(nèi)，大部分為第四系沖洪積（）粘土、粉質(zhì)粘土，測線小號點(diǎn)約40～50 m附近有小范圍砂巖出露。測線大號點(diǎn)約180 m處，鉆探取心的情況為：第四系沖洪積）粘土、粉質(zhì)粘土層；含礫粘土、卵石層；下部為石炭系中統(tǒng)（C2）白云巖、灰?guī)r，風(fēng)化嚴(yán)重。要求利用高密度電法勘探查明斷裂（砂巖和白云巖的接觸位置）及斷裂兩側(cè)地層分布等情況。

2.1測線布置及高密度電阻率數(shù)據(jù)的采集方法

沿設(shè)計(jì)管道軸線在地表的投影以及左右兩側(cè)各10 m布置了3條高密度電法測線（a1、a2、a3），3條測線長度均為840 m。測量儀器為WDJD-2型多功能數(shù)字高密度電阻率儀。根據(jù)勘查要求，采用聯(lián)合三極裝置，即在一條測線上同時進(jìn)行AMN和MN-B測量。采用的測量電極距為10 m，點(diǎn)距10 m，測量層數(shù)為8層，對應(yīng)供電極距AO為15～85 m。

2.2實(shí)測數(shù)據(jù)的初步統(tǒng)計(jì)分析判斷地層結(jié)構(gòu)

根據(jù)地質(zhì)及鉆孔資料可知，斷裂位于測線40～180 m范圍內(nèi)，因此斷裂定位解釋主要在該范圍內(nèi)進(jìn)行。使用實(shí)測的2組三極裝置數(shù)據(jù)組合出對稱四極測深數(shù)據(jù)，3條0～300 m剖面范圍內(nèi)的視電阻率斷面等值線如圖2所示?？梢钥闯?，3條剖面的視電阻率變化規(guī)律基本相似，在150 m附近有明顯的電性差異，小號點(diǎn)一側(cè)視電阻率較低，在50～150 Ωm范圍內(nèi)，大號點(diǎn)一側(cè)在150 Ωm以后，呈現(xiàn)出起伏層狀特征。

3條剖面的視電阻率特征基本一致，但使用常規(guī)的反演軟件解釋出的電性參數(shù)差異較大，無法根據(jù)解釋結(jié)果統(tǒng)一給出斷裂位置和各層深度等信息。因此，采用統(tǒng)計(jì)學(xué)建模方法進(jìn)行反演，處理的步驟如下：

按照圖1統(tǒng)計(jì)學(xué)建模流程圖的步驟，首先根據(jù)鉆孔資料和視電阻率分布情況，將探測深度范圍內(nèi)的地層確定為3層，對a1、a2、a3剖面上各測深點(diǎn)的數(shù)據(jù)進(jìn)行一維反演，獲得各點(diǎn)的電阻率和地層厚度參數(shù)，對這些層參數(shù)分剖面統(tǒng)計(jì)數(shù)學(xué)期望和方差。其中，a1剖面統(tǒng)計(jì)結(jié)果如表1所示。

圖2　實(shí)測視電阻率斷面等值線圖Fig.2 　Contour map of measured apparent resistivity section

表　1 a1剖面一維反演參數(shù)統(tǒng)計(jì)結(jié)果Table 1 Layer parameters statistic result of Profile a1

對照文獻(xiàn)［1］中的判斷規(guī)則，3個剖面的統(tǒng)計(jì)數(shù)據(jù)都支持判斷地層結(jié)構(gòu)為直立異常體（斷裂）。文獻(xiàn)［1］中的地層結(jié)構(gòu)判斷方法是根據(jù)各種地層結(jié)構(gòu)模型的正演數(shù)據(jù)統(tǒng)計(jì)歸納后得出的統(tǒng)計(jì)規(guī)律。

作出測區(qū)地層結(jié)構(gòu)為直立異常體（斷裂）初始判斷后，參考地質(zhì)資料修正地層結(jié)構(gòu)。從測區(qū)地質(zhì)資料可知：測線小號點(diǎn)一側(cè)有砂巖出露，大號點(diǎn)一側(cè)基巖為灰?guī)r，判斷測區(qū)內(nèi)地質(zhì)結(jié)構(gòu)為有覆蓋層的接觸面結(jié)構(gòu)，接觸面兩側(cè)的巖層電性不同。

2.3各剖面斷裂（巖性接觸面）位置的確定方法

為了確定各剖面斷裂（巖性接觸面）位置，按照圖1統(tǒng)計(jì)學(xué)建模流程圖的步驟對實(shí)測視電阻率數(shù)據(jù)求各極距的水平梯度，用1.2節(jié)中求斷裂邊界的式（2）確定接觸面位置，這一方法的實(shí)質(zhì)是利用不同極距的視電阻率的水平梯度變化確定電性變化的邊界。為確定深部的電性變化位置，以大極距水平梯度信息為主要依據(jù)。a1剖面的水平梯度曲線如圖3a所示，包括極距為75和85 m的水平梯度曲線，水平梯度的極大值點(diǎn)在剖面的150 m處。圖3b是極距為75和85 m的視電阻率曲線，可以看出，以剖面中間為界，兩側(cè)有明顯的電性差異。其他兩條剖面（a2、a3）水平梯度和視電阻率曲線與a1剖面類似。由此，初步確定巖性接觸面在a1剖面的150 m處。

圖3　a1剖面的水平梯度和視電阻率曲線Fig.3 Horizontal grade and apparent resistivit curves of Profile a1

2.4初始模型層參數(shù)的確定

按照統(tǒng)計(jì)學(xué)建模流程圖（圖1）的步驟，確定了巖性接觸面的位置后，對a1剖面150 m（斷裂位置）兩側(cè)重新分別統(tǒng)計(jì)電阻率測深曲線類型及層參數(shù)的數(shù)學(xué)期望，各統(tǒng)計(jì)參數(shù)如表2所示，左側(cè)表示小于150 m一側(cè)，右側(cè)表示大于150 m一側(cè)。從表2可知，巖性接觸面的左側(cè)，電阻率測深曲線的類型主要為H型，右側(cè)的曲線類型主要為K型。

表2　a1剖面接觸面兩側(cè)曲線類型統(tǒng)計(jì)結(jié)果Table 2 Curve type and layer parameters’expect for each side of Profile a1

2.5使用統(tǒng)計(jì)學(xué)參數(shù)建模

綜合2.2、2.3和2.4節(jié)的統(tǒng)計(jì)學(xué)方法確定的巖性接觸面位置、斷層邊界兩側(cè)層參數(shù)等信息建立的a1剖面的反演初始模型見圖4，初始模型中各區(qū)塊的電阻率如模型中圖例所示。從a1剖面的初始模型可以直觀地看到測線范圍內(nèi)的地層結(jié)構(gòu)，電性分布的大致情況。

圖　4 a1剖面初始模型Fig.4 Initial model of Profile a1

3　修正的遺傳算法反演效果

使用圖4的初始模型進(jìn)行反演擬合，初次擬合誤差在100%以上。按照1.5節(jié)中的修正的遺傳算法反演步驟，分別修改巖性接觸帶邊界、電阻率和層厚度參數(shù)，最終獲得的反演結(jié)果斷面如圖5所示。a1測線的擬合誤差為22.4%。與初始模型相比，地層結(jié)構(gòu)沒有變化，反演改變了電阻率和界面深度。

由于a2與a1相鄰且間距只有10 m，a2剖面的反演用a1的反演結(jié)果作為初始模型，這樣選擇初始模型的好處是大幅度降低了計(jì)算量，也省去了圖1中統(tǒng)計(jì)學(xué)建模的各個步驟，初次計(jì)算擬合誤差僅為25.4%。反演完成時a2剖面的擬合誤差為21.9%，從初始模型到反演結(jié)束所需計(jì)算量大大少于a1剖面。

使用a2剖面的反演結(jié)果作為a3剖面的反演初始模型，初次計(jì)算擬合誤差為26.7%。反演完成時，a3剖面的擬合誤差為21.6%。3條剖面的反演結(jié)果如圖5所示，各剖面均有明確清晰的地層結(jié)構(gòu)和層參數(shù)信息，各剖面的電阻率變化較小，從反演結(jié)果可以看到下層巖性接觸面從a1到a3剖面逐漸向大號點(diǎn)移動。

為了證明統(tǒng)計(jì)學(xué)建模方法的反演效果，以圖5各剖面反演結(jié)果為模型，用有限元方法計(jì)算了視電阻率斷面等值線圖，結(jié)果如圖6所示。與圖2實(shí)測數(shù)據(jù)的斷面等值線圖對比，兩者的等值線圖像的主要特征基本相似。正演計(jì)算與實(shí)測數(shù)據(jù)的平均擬合誤差為22.0%。

圖5的反演成果給出了各剖面巖性接觸帶（即通過測區(qū)的斷層）的位置、斷層兩側(cè)的地層分層信息，及各層的厚度變化情況。與鉆孔資料對比，層厚度基本準(zhǔn)確。綜合各剖面反演結(jié)果可以看出，接觸面位置在a1剖面的140 m、a2剖面的145 m、a3剖面的150 m處，從而推斷出接觸面的走向。

圖5　各剖面反演結(jié)果Fig.5 　Inversion model for each profile

圖6　各剖面反演結(jié)果的正演計(jì)算結(jié)果Fig.6 Stimulation result of inversion model for each profile

4　結(jié) 論

采用統(tǒng)計(jì)學(xué)建立初始模型的反演方法可以有效完成對地層結(jié)構(gòu)的判斷，并利用水平梯度確定接觸帶位置。將地質(zhì)資料作為先驗(yàn)信息參與到反演過程中，減少了遺傳算法的計(jì)算量。

應(yīng)用實(shí)例表明：應(yīng)用統(tǒng)計(jì)學(xué)、水平梯度計(jì)算和先驗(yàn)信息融合的方法盡可能利用了高密度電法實(shí)測數(shù)據(jù)和地質(zhì)信息，建立了最大限度接近真實(shí)的初始模型。初始模型建立后，使用遺傳算法配合有限元正演模擬調(diào)整模型參數(shù)，完成了非線性反演。在反演結(jié)果中，地層結(jié)構(gòu)被明確顯示出來，地層厚度、電阻率和分界面位置都有明確的數(shù)值。

在面積性高密度電法勘探資料解釋中，當(dāng)測區(qū)內(nèi)地質(zhì)條件變化不大、相鄰測線高密度電阻率異常特征相似時，有必要利用好相鄰測線的反演結(jié)果。將某一測線的反演結(jié)果用于全測區(qū)，既提高了反演速度，又可以使測區(qū)的地層結(jié)構(gòu)和層參數(shù)具有統(tǒng)一性，利于進(jìn)一步的地質(zhì)解釋，提高反演的速度和質(zhì)量。

電法勘探中還有其他測深類方法，例如激發(fā)極化測深、頻率測深及瞬變電磁測深等和電阻率測深有類似數(shù)據(jù)結(jié)構(gòu)的方法，也可以用統(tǒng)計(jì)學(xué)建模方法建立初始模型并進(jìn)行反演，預(yù)期可以獲得較好的地質(zhì)效果。由于影響巖土層電性變化的因素很多，自然界中地下電性結(jié)構(gòu)非常復(fù)雜。如何處理復(fù)雜的地層結(jié)構(gòu)以及含有各種干擾因素的野外實(shí)測數(shù)據(jù)，統(tǒng)計(jì)學(xué)建模法還需要進(jìn)一步探索研究和完善。

［1］程勃，底青云.基于遺傳算法和統(tǒng)計(jì)學(xué)的電阻率測深二維反演研究［J］.地球物理學(xué)進(jìn)展，2012，27（2）：788-795.

［2］程勃，底青云.基于偏導(dǎo)數(shù)和統(tǒng)計(jì)學(xué)方法的電阻率測深二維反演［J］.石油地球物理勘探，2012，47（6）：1006-1013.

［3］程勃，底青云.復(fù)雜地電結(jié)構(gòu)條件下統(tǒng)計(jì)學(xué)建模法電阻率測深二維反演 ［J］.地球物理學(xué)報(bào)，2014，57（3）：961-967.

［4］阮百堯，村上裕，徐世浙.電阻率/激發(fā)極化率數(shù)據(jù)的二維反演程序［J］.物探化探計(jì)算技術(shù)，1999，21（2）：116 -125.

［5］蘇朱劉，胡文寶，嚴(yán)良俊.電阻率和極化率測深法的正演修正法反演［J］.石油物探，2005，44（2）：194-198.

［6］何門貴，溫永輝.高密度電阻率法二維反演在工程勘探中的應(yīng)用［J］.物探與化探，2002，26（2）：156-159.

［7］底青云，倪大來，王若，等.高密度電阻率成像［J］.地球物理學(xué)進(jìn)展，2003，18（2）：323-326.

［8］趙改善，譯.求解非線性最優(yōu)化問題的遺傳算法［J］.地球物理學(xué)進(jìn)展，1992，7（1）：90-97.

［9］劉斌，聶利超，李術(shù)才，等.三維電阻率空間結(jié)構(gòu)約束反演成像方法［J］.巖石力學(xué)與工程學(xué)報(bào)，2012，31（11）：2258-2268.

［10］阮百堯.電阻率/激發(fā)極化法測深數(shù)據(jù)的一維最優(yōu)化反演方法［J］.桂林工學(xué)院學(xué)報(bào)，1999，19（4）：321-325.

［11］阮百堯.電阻率測深解釋中的一種新的反演方法 ［J］.桂林冶金地質(zhì)學(xué)院學(xué)報(bào)，1994，14（1）：80-85.

［12］呂玉增，阮百堯.高密度電法二維反演軟件設(shè)計(jì)與實(shí)現(xiàn)［J］.桂林工學(xué)院學(xué)報(bào)，2004，24（4）：417-421.

Statistic modeling 2D inversion application in parallel resistivity survey lines

CHENG Bo，DING Yan-li，XU Zhi-feng
（a.Guangxi Key Laboratory of Hidden Metallic Ore Deposit Exploration；b.College of Earth Sciences，Guilin University of Technology，Guilin 541004，China）

For the key geophysical information，statistic modeling inversion is applied in a survey area with a few parallel survey lines.This method makes full use of the field data to get initial model，and complete the inversion with genetic algorithm.One line's inversion result became the initial model of the next line.The example proves that this method reduces the fitting error effectively.The inversion results of lines show the geophysical characteristic in the survey area.

resistivity；2D inversion；statistic modeling；gas pipe

P631.322

A

1674-9057（2016）03-0452-06

10.3969/j.issn.1674-9057.2016.03.006

2015-04-28

國家自然科學(xué)基金項(xiàng)目（41404116）；廣西自然科學(xué)基金項(xiàng)目 （2014GXNSFBA118232；2014GXNSFAA118305）；廣西隱伏金屬礦產(chǎn)勘查重點(diǎn)實(shí)驗(yàn)室系統(tǒng)研究課題 （15-140-27-06）

程 勃 （1982—），男，博士，講師，研究方向：地球物理正反演，chengbo@glut.edu.cn。

引文格式：程勃，丁彥禮，徐志鋒.統(tǒng)計(jì)學(xué)建模法在面積性高密度電阻率數(shù)據(jù)反演中的應(yīng)用［J］.桂林理工大學(xué)學(xué)報(bào)，2016，36（3）：452-457.