(19)中华 人民共和国 国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202110408595.4 (22)申请日 2021.04.16 (65)同一申请的已公布的文献号 申请公布号 CN 113516347 A (43)申请公布日 2021.10.19 (73)专利权人 江苏省地质调查研究院 地址 210049 江苏省南京市珠 江路700号 (72)发明人 尚通晓　郝社锋　关艺晓　朱首峰　 盛君　姜素　许书刚　姜国庆　 张大莲　 (74)专利代理 机构 北京德崇智捷知识产权代理 有限公司 1 1467 代理人 杨楠 (51)Int.Cl. G06Q 10/06(2012.01) G06Q 10/04(2012.01) G06F 30/27(2020.01) G06K 9/62(2022.01) G06Q 50/06(2012.01)G06T 17/00(2006.01) (56)对比文件 CN 211652186 U,2020.10.09 CN 106203729 A,2016.12.07 CN 111965720 A,2020.1 1.20 CN 103487843 A,2014.01.01 US 2016195 631 A1,2016.07.07 US 2016334531 A1,2016.1 1.17 田少兵 等. “直流电测深法在沿海平原区地 质填图中的应用—以苏 北盆地连云港灌云地区 为例”. 《物探与化探》 .2019,第43卷(第4期),第 783-793页. Asfahani J 等.“Geoelectrical and hydrochemical i nvestigations for characterizi ng the salt water i ntrusion in the Kha-nas ser Valley”. 《Acta Geophysica》 .2013,第62卷(第2期),第42 2-444 页. 审查员 田瑜 (54)发明名称 基于地球物理与水文数据融合的地下水水 质评估方法 (57)摘要 本发明公开了一种基于地球物理与水文数 据融合的地下水水质评估方法， 包括以下步骤： 用电阻率测深法获取换算电阻率数据， 并获取钻 井的水文地质分层数据、 电测井数据和矿化度分 析数据； 利用电测井数据对换算电阻率数据进行 校正； 基于水文地质分层数据构建多层含水层结 构三维模型， 并利用校正换算电阻率数据对模型 的空间分辨率进行优化； 根据校正换算电阻率数 据确定咸淡水分界面； 根据钻井的矿化度分析数 据和钻井旁的平均校正换算电阻率， 确定各含水 层的矿化度 ‑电阻率相关数学模型， 由各含水层 的平均校正换算电阻率反演 矿化度， 建立各含水 层矿化度空间分布。 本发明将地球物理数据和水文数据融合， 可以低成本获得高分辨率地下水水 质空间分布。 权利要求书2页 说明书6页 附图6页 CN 113516347 B 2022.01.14 CN 113516347 B 1.基于地球物理与水文数据融合的地下 水水质评估方法， 其特 征在于， 包括以下步骤： 步骤A、 用电阻率测深法获取研究区域的换算电阻率数据， 并获取研究区域中一组钻井 的水文数据， 包括水文地质分层数据、 电测井数据和矿化度分析 数据； 步骤B、 利用所述电测井数据对换算电阻率数据进行校正， 得到研究区域的校正换算电 阻率数据； 步骤C、 基于所述水文地质分层数据构建研究区域的多层含水层结构三维模型， 并利用 校正换算电阻率数据对所述多层含水层结构三维模型的空间分辨 率进行优化； 步骤D、 根据校正换算电阻率数据确定研究区域的咸淡水分界面； 步骤E、 根据钻井的矿化度分析数据和钻井旁的平均 校正换算电阻率， 确定各含水层的 矿化度‑电阻率相关数学模型， 并根据所述数学模 型， 由各含 水层的平均校正换算电阻率反 演矿化度， 建立各含水层矿化度空间分布。 2.如权利要求1所述地下水水质评估方法， 其特征在于， 根据换算电阻率数据对研究区 域进行电性分区， 并按照电性分区布置所述钻井， 以保证 每个电性分区都有至少一个钻井。 3.如权利要求1所述地下水水质评估方法， 其特征在于， 按照以下方法对换算电阻率数 据进行校正： 分别找出每个钻井的电测井数据中电阻率 ‑深度曲线以及钻井旁电阻率测深 法测点的换算电阻率 ‑深度曲线中电阻率变化最大的前N个拐点， 按照深度从小到大的排序 方式， 用电测井数据中电阻率 ‑深度曲线中第i个拐点的深度与钻井旁电阻率测深法测点的 换算电阻率 ‑深度曲线中第i个拐点的深度之 间的比值作为钻井旁测点的换算电阻率 ‑深度 曲线中第i ‑1个拐点至第i个拐点间深度数据所对应的深度校正系数， N为大于等于1的整 数， i＝1,2, …,N； 使用插值或拟合的方法获得研 究区域中其它测点换算电阻率 ‑深度曲线 的深度校正系 数； 最后用各测点的深度校正系 数对各测点的换算电阻率 ‑深度曲线中的相 应深度数据进行 校正， 得到研究区域的校正换算电阻率数据。 4.如权利要求3所述 地下水水质评估方法， 其特 征在于， N 为1或2。 5.如权利要求1所述地下水水质评估方法， 其特征在于， 具体使用以下方法构建研究区 域的多层含水层结构三维模型： 利用所述水文地质分层数据， 通过插值或拟合计算出各含 水层、 隔水层以及基岩界面， 从而构建多层含水层结构三维模型。 6.如权利要求1所述地下水水质评估方法， 其特征在于， 所述利用校正换算电阻率数据 对所述多层含水层结构三维模型 的空间分辨率进行优化， 具体如下： 对于所述多层含水层 结构三维模型中校正换算电阻率数据存在明显电性差异的相 邻含水层， 结合校正换算电阻 率数据构建空间分辨 率更高的所述相邻含水层之间的含水层界面。 7.如权利要求1所述地下水水质评估方法， 其特征在于， 所述换算电阻率通过以下公式 计算： 权　利　要　求　书 1/2 页 2 CN 113516347 B 2式中， ρ(i)为第 i层换算电阻率， ρS(i)为第i个极距实测视电阻率， h(i)为第i层厚度， r (i)为第i个AB/2极距， K(i)为第i层反射系数， 是双对 数坐标下电测深曲线在第i ‑1到第i个 极距的斜 率， K′(i)为斜率校正后的反射系数。权　利　要　求　书 2/2 页 3 CN 113516347 B 3