[發明專利]一種高密度電阻率成像方法

在審

• 申請號: 201610325150.9
• 申請日: 2016-05-17
• 公開（公告）號: CN105785476A
• 公開（公告）日: 2016-07-20
• 發明（設計）人: 盧德寶;王鳳;陳曉東;歐健
• 申請（專利權）人: 浙江水利水電學院
• 主分類號: G01V11/00
• 分類號: G01V11/00
• 代理公司: 北京科億知識產權代理事務所(普通合伙) 11350
• 代理人: 湯東鳳

1.一種高密度電阻率成像方法，其特征在于，包括以下步驟：(1).針對勘測區域，根據 地形特點以及高密度電阻率測量與人工折射地震測量的需要建立適宜的網格剖分方案； (2).根據網格劃分，建立人工折射地震數據的反演模型進行反演計算，從而得離散化的速 度模型；(3).根據速度與電阻率的經驗公式，將離散化的速度模型轉換成離散化的電阻率 模型；(4).將(3)中得到電阻率模型作為電阻率反演的約束初始模型進行反演計算，從而得 到真實電阻率分布；(5).對所得數據進行成像，輸出勘測區的高精度電阻率成像結果。

2.根據權利要求1所述的高密度電阻率成像方法，其特征在于，針對待測區域建立合理 的測量布極方式與計算網格剖分方案，在現場勘探時，采用相同的電極極距布置方式，對應 采集電阻率數據與人工折射地震數據，根據需要選擇二維測定或者三維測定，并利用開源 四面體網格生成程序Tetgen生成數據正演計算所需四面體或三角形剖分網格。

3.根據權利要求1所述的高密度電阻率成像方法，其特征在于，人工折射地震數據的處 理是通過利用Coppens的滑動視窗能量比法將測得的人工折射地震原始數據根進行時窗計 算，以得到折射地震走時數據，在上述得到的剖分網格的基礎上，利用開源地震數據處理程 序Seismicunix(SU)建立人工折射地震數據差分反演模型，進行反演計算，輸出離散化的 速度場模型。

4.根據權利要求1所述的高密度電阻率成像方法，其特征在于，速度與電阻率之間的關 系轉換，根據速度與電阻率之間的轉換經驗公式Log₁₀ρ＝mLog₁₀V_p+c，將速度模型轉換為電 阻率模型，其中ρ為電阻率，V_p為p波的速度，m、c為經驗常數，可根據研究區的巖性特征進行 確定。

5.根據權利要求1所述的高密度電阻率成像方法，其特征在于，代入先決條件的電阻率 反演計算，根據全局牛頓-高斯方程來建立最小化目標方程可以得到電阻率反演的擾動方 程式中s為感度矩陣，d為觀測數據，ρ為模 型的電阻率，G為正演算子，R_dd與R_mm分別表示觀測數據與模擬數據的協方差矩陣，ρ₀為初始 電阻率模型，Δρ_k為第k此迭代后的模擬電阻率，傳統的初始電阻率模型一般采用均值模 型，此處則采用上述轉換而來的電阻率模型，通過迭代計算，即可得到高精度的電阻率值。

6.根據權利要求1所述的高密度電阻率成像方法，其特征在于，利用矢量圖形描述語言 Asymptote將得到的電阻率值進行插值成像，即可輸出高質量的高密度電阻率成像矢量圖。