[發明專利]一種基于電阻率勘探法的地下水脆弱性評價方法

有效

• 申請號: 201610512180.0
• 申請日: 2016-06-27
• 公開（公告）號: CN105938204B
• 公開（公告）日: 2018-06-15
• 發明（設計）人: 盧德寶;王鳳;王蕙;陳曉東;歐劍
• 申請（專利權）人: 浙江水利水電學院
• 主分類號: G01V3/00
• 分類號: G01V3/00;G01N27/04
• 代理公司: 北京科億知識產權代理事務所(普通合伙) 11350
• 代理人: 湯東鳳

1.一種基于電阻率勘探法的地下水脆弱性評價方法,其特征在于，包括以下步驟：步驟一，根據待評價區的地形地貌特點，制定一維及二維電阻率測定方案，采集待評價區的一維及二維視電阻率數據；步驟二，對采集到的數據進行反演計算，得到待評價區的真實電阻率分布及二維電阻率成像結果；步驟三，根據以上結果，提取待評價區脆弱性評價所需關鍵參數；步驟四，依據提取的參數，根據IEC計算公式，計算出待評價區的IEC值；步驟五，根據制定的IEC分級標準及其輻射范圍編制地下水脆弱等級圖，輸出脆弱性評價結果；

其中在所述的步驟一中，針對待評價區的特點，建立電阻率測量方案，分別使用一維電阻率、二維電阻率成像測量設備進行待評價區的視電阻率測量，獲得視電阻率數據；

在所述的步驟二中，使用全局牛頓-高斯方程法對一維、二維電阻率數據進行反演計算，具體方程如下：

公式(1)中S為感度矩陣，S_κ為第k次迭代后的感度矩陣，第k次迭代后的感度矩陣轉置，d為觀測數據，ρ為模型的電阻率，G為正演算子，R_dd與R_mm分別表示觀測數據與模擬數據的協方差矩陣，ρ₀為初始電阻率模型，ρ_κ為第k次迭代后的電阻率，Δρ_k為第k次迭代后的電阻率變化量；

在所述的步驟三中，根據以上反演計算的結果，得到電阻率分布與二維電阻率成像結果，分別提取出測量點含水層上方第i層沉積層電阻率ρ_i、含水層上方第i層沉積層厚度h_i，根據公式(2)

在所述的步驟三和四中，計算出待評價區內所有測量點的含水層上方每一層IEC值，作為該測量點的評價參數；

在所述的步驟五中，根據待評價區的實際需求，依據IEC值將待評價區分成脆弱區、較脆弱區、較安全區、安全區，結合每一測量點的輻射范圍，即可得到待評價區的地下水脆弱性評價結果空間分布圖。