一种基于电阻率的水渗流岩体电各向异性数值模拟方法

1.本发明涉及地质监测设备技术领域，特别涉及一种基于电阻率的水渗流岩体电各向异性数值模拟方法。


背景技术：

2.大坝、水库等储水体的渗漏过程、冰雪融化过程、灌溉过程、雨水在地面的下渗过程、海水入侵过程、存储生活垃圾和核废料的储体渗漏过程、煤矿的突水过程、浆液进入注浆体的过程以及各种水文探测中示踪剂在水中的流动过程均可引起周围介质的电性变化，水渗的本质是采动围岩中水渗流的失稳，由于采动或其他岩土工程活动，常常会造成岩体开裂甚至破碎，而破碎岩体裂隙中的渗透率要远比孔隙渗透率高，这就是渗流引起的重大灾害事故往往多发生在破碎岩体中的原因。
3.单纯只靠实地监测电阻率的方式，需要对多个布局点进行人工监测，不仅耗费了大量的人力物力，且只通过电阻率来监测水渗流岩体的典型变化不够全面，导致对电性变化的测量产生误差，不能准确得到水渗流岩体的典型变化，导致不能准确地预测重大灾害事故的发生。


技术实现要素：

4.本发明提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，可以对水渗流岩体的结构进行预测，及时对水渗流岩体采取应对措施，避免灾害发生。
5.一种基于电阻率的水渗流岩体电各向异性数值模拟方法，包括：
6.s1：对水渗流岩体进行岩体勘探，得到水渗流岩体模型，并基于所述水渗流岩体模型设定电阻率监测点，得到在电阻率监测点的电阻率；
7.s2：基于所述电阻率监测点和电阻率，确定所述水渗流岩体模型的电各向异性取值；
8.s3：基于所述水渗流岩体模型，确定各个水渗流岩体的结构特征，利用数值模拟的方法确定电各向异性和结构特征之间的关系。
9.优选的，s1中，对水渗流岩体进行岩体勘探，得到水渗流岩体模型，包括：
10.利用自动化探测设备对所述水渗流岩体进行三维点云数据获取和平面图像获取；
11.对所述三维点云数据得到所述水渗流岩体的结构布局和尺寸，基于所述二维图像确定所述水渗流岩体的材质特征；
12.基于所述结构布局、尺寸和材质特征，建立所述水渗流岩体模型。
13.优选的，s1中，基于所述水渗流岩体模型设定电阻率监测点，得到在电阻率监测点的电阻率，包括：
14.基于所述水渗流岩体模型的岩体密度和岩体孔隙率，确定在所述水渗流岩体模型的电阻率监测点；
15.从实体的水渗流岩体获取所述电阻率监测点对应的实际点，将电极插入泥地中，
并接收电极的信号，基于所述信号，确定电阻率。
16.优选的，s2中，基于所述电阻率监测点和电阻率，确定所述水渗流岩体模型的电各向异性取值，包括：
17.基于电阻率监测点之间的相对位置特征和电阻率监测点位置所处的岩体结构特征，确定每个电阻率监测点与其他电阻率监测点的关联系数；
18.基于所述关联系数和电阻率，确定所述渗流岩体模型的电各向异性取值。
19.优选的，s3中，基于所述水渗流岩体模型，确定各个水渗流岩体的结构特征，包括：
20.将所述水渗流岩体模型按照网格划分为多个单元，按照预设分析方法对每个单元进行分析，得到第一单元特征；
21.根据所述第一单元特征，对所述多个单元进行合并，得到整体单元，并根据所述第一单元特征，得到整体单元的第二单元特征；
22.对所述第二单元特征进行融合处理，得到水渗流岩体的结构特征。
23.优选的，s3中，利用数值模拟的方法确定电各向异性和结构特征之间的关系包括：
24.基于结构属性，对所述结构特征进行一级标记和二级标记，基于电各向异性取值对所述结构特征进行结果标记，并根据标记结果，得到结果标记与标记组合的对应关系，其中，所述标记组合包括某一种一级标记和二级标记，所述对应关系为结果标记对应标记组合，对应方式为一对多；
25.基于所述电各向异性、结构特征，在所述水渗流岩体模型的基础上进行标记建立静态模型，并基于所述对应关系，赋予所述二级标记的动态调整功能，建立动态模型；
26.对动态模型进行再次数值模拟，确定电各向异性和结构特征之间的关系。
27.优选的，根据标记结果，得到结果标记与标记组合的对应关系，包括：
28.基于结构属性，对所述结构特征进行特征分析，得到每个结构属性下的属性子特征，并按照结构属性的等级划分规则，将所述属性子特征进行等级划分，得到属性子特征的等级；
29.根据所述属性子特征，对结构特征进行一级标记，根据所述等级，对结构特征进行二级标记，基于所述电各向异性取值对所述结构特征进行结果标记；
30.将所述一级标记和二级标记作为输入，所述结果标记作为分类标准，获取在相同结果标记对应的一级标记和二级标记的组合，得到结果标记与标记组合的对应关系。
31.优选的，对动态模型进行再次数值模拟，确定电各向异性和结构特征之间的关系，包括：
32.在所述动态模型中按照一级标记，对所述二级标记进行单一性调整，基于所述对应关系，对所述动态模型进行数值模拟，输出单一性调整对应的单一动态结果标记，改变所述一级标记的单一性调整属性，直到得到所有一级标记对应的单一动态结果标记；
33.基于所述所有一级标记及其对应的单一动态结果标记之间的标记关系，对所有一级标记中至少两个一级标记同时进行调整的情况进行再次数值模拟，得到综合结果标记，并以综合结果标记对应的电各项异性数值为横坐标，二级标记为对应的属性子特征取值为纵坐标，建立得到关系曲线图。
34.优选的，基于所述水渗流岩体模型的岩体密度和岩体孔隙率，确定在所述水渗流岩体模型的电阻率监测点，包括：
35.按照所述水渗流岩体模型的岩体密度，对所述水渗流岩体模型进行划分，得到部分水渗流岩体，并基于部分水渗流岩体的岩体密度和面积，设定在所述部分水渗流岩体的监测点个数；
36.对所述部分水渗流岩体再次进行更大精度的岩体密度检测，并根据检测结果，确定在所述部分水渗流岩体的监测区域，其中，所述监测区域对应所述监测点个数；
37.获取所述监测区域中岩体孔隙率最高的位置，作为电阻率监测点。
38.优选的，获取所述监测区域中岩体孔隙率最高的位置，作为电阻率监测点后，还包括：
39.确定所述位置在所述水渗流岩体模型中坐标点，获取所述坐标点在实地的水渗流岩体的实地位置；
40.确定所述实地位置的地质信息，并判断所述地质信息是否满足电阻率的监测环境要求；
41.若是，确定所述位置为最终的电阻率监测点；
42.否则，选取与所述位置最近且满足所述监测环境要求的临近位置作为最终的电阻率监测点。
43.本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
44.下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
45.附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：
46.图1为本发明实施例中一种基于电阻率的水渗流岩体电各向异性数值模拟方法的流程图；
47.图2为本发明实施例中获取水渗流岩体模型的流程图；
48.图3为本发明实施例中确定各个水渗流岩体的结构特征的流程图。
具体实施方式
49.以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
50.实施例1
51.一种基于电阻率的水渗流岩体电各向异性数值模拟方法，如图1所示，包括：
52.s1:对水渗流岩体进行岩体勘探，得到水渗流岩体模型，并基于所述水渗流岩体模型设定电阻率监测点，得到在电阻率监测点的电阻率；
53.s2:基于所述电阻率监测点和电阻率，确定所述水渗流岩体模型的电各向异性取值；
54.s3:基于所述水渗流岩体模型，确定各个水渗流岩体的结构特征，利用数值模拟的方法确定电各向异性和结构特征之间的关系。
55.在该实施例中所述电各向异性取值用于表示此水渗流岩体的电学性能沿不同方向的差异，差异越大，对应的取值越大。
56.在该实施例中，对水渗流岩体进行岩体勘探包括对岩体材质、岩体见空隙、岩体外形等进行勘探。
57.在该实施例中，所述电阻率监测点根据水渗流岩体模型中水渗流岩体的结构布局决定。
58.在该实施例中，所述结构特征包括材料、尺寸、空隙等的特征。
59.在该实施例中，利用数值模拟的方法确定电各向异性和结构特征之间的关系后，可以明确电各向异性的取值对结构特征的影响关系，基于所述影响关系，在确定电各向异性的取值后，可以对水渗流岩体的结构进行预测，根据预测结果及时对水渗流岩体采取应对措施，避免灾害发生。
60.上述设计方案的有益效果是：通过根据对岩体的勘探，得到水渗流岩体模型，通过水渗流岩体模型设置电阻率监测点，得到在电阻率监测点的电阻率，然后，基于所述电阻率监测点和电阻率，确定所述水渗流岩体模型的电各向异性取值，利用数值模拟的方法对水渗流岩体模型的结构特征和电各向异性取值进行分析，确定水渗流岩体模型的结构特征和电各向异性取值之间的关联度，可以在在确定电各向异性的取值后，可以对水渗流岩体的结构进行预测，根据预测结果及时对水渗流岩体采取应对措施，避免灾害发生。
61.实施例2
62.基于实施例1的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，如图2所示，s1中，对水渗流岩体进行岩体勘探，得到水渗流岩体模型，包括：
63.s101：利用自动化探测设备对所述水渗流岩体进行三维点云数据获取和平面图像获取；
64.s102：对所述三维点云数据得到所述水渗流岩体的结构布局和尺寸，基于所述二维图像确定所述水渗流岩体的材质特征；
65.s103：基于所述结构布局、尺寸和材质特征，建立所述水渗流岩体模型。
66.在该实施例中，利用自动化探测设备对所述水渗流岩体进行三维点云数据为对水渗流岩体进行三维将官扫描，得到三维点云数据。
67.在该实施例中，基于所述二维图像确定所述水渗流岩体的材质特征具体为对图像的像素进行分析，来确定水渗流岩体的像素特征，根据像素特征和材质特征之间的对应关系，来确定材质特征。
68.上述设计方案的有益效果是：通过利用自动化探测设备对所述水渗流岩体进行三维点云数据获取和平面图像获取得到岩体的特征来建立水渗流岩体模型，保证了水渗流岩体模型的准确性，为数值模拟提供基础。
69.实施例3
70.基于实施例1的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，s1中，基于所述水渗流岩体模型设定电阻率监测点，得到在电阻率监测点的电阻率，包括：
71.基于所述水渗流岩体模型的岩体密度和岩体孔隙率，确定在所述水渗流岩体模型
的电阻率监测点；
72.从实体的水渗流岩体获取所述电阻率监测点对应的实际点，将电极插入泥地中，并接收电极的信号，基于所述信号，确定电阻率。
73.在该实施例中，所述岩体密度越大的位置设置的电阻率监测点较多，岩体孔隙率较大的位置容易发送灾害，所以设置的电阻率监测点也相应较多。
74.上述设计方案的有益效果是：通过根据水渗流岩体模型的岩体密度和岩体孔隙率，确定在所述水渗流岩体模型的电阻率监测点，使得电阻率监测点的布局合理，从而使得计算得到的电各向异性的取值更加准确，可以更加准确的了解水渗流岩体的电性变化，提高后续数值模拟结果的准确性，可以对水渗流岩体的结构进行预测，根据预测结果及时对水渗流岩体采取应对措施，避免灾害发生。
75.实施例4
76.基于实施例1的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，如图3所示，s2中，基于所述电阻率监测点和电阻率，确定所述水渗流岩体模型的电各向异性取值，包括：
77.基于电阻率监测点之间的相对位置特征和电阻率监测点位置所处的岩体结构特征，确定每个电阻率监测点与其他电阻率监测点的关联系数；
78.所述关联系数的计算公式如下：
[0079][0080]
其中，t
(a，b)
表示当前电阻率监测点和其他电阻率监测点之间的关联系数，s
(a，b)
表示所述当前电阻率监测点和其他电阻率监测点之间的距离，s0表示参考距离，n表示所述当前电阻率监测点或其他电阻率监测点的结构属性个数，σi表示第i个结构属性对电阻率的影响系数，取值为(0.2，0.8)，ai表示所述当前电阻率监测点的第i个结构属性的属性特征值，取值为(0，1)，bi表示其他电阻率监测点的第i个结构属性的属性特征值，取值为(0，1)；
[0081]
基于所述关联系数和电阻率，确定所述渗流岩体模型的电各向异性取值；
[0082]
所述当前电阻率监测点的电各向异性取值的计算公式如下：
[0083][0084]
其中，ka表示所述当前电阻率监测点的电各向异性取值，m表示他电阻率监测点的个数，t
(a，j)
表示当前电阻率监测点和第j个其他电阻率监测点之间的关联系数，fa表示所述当前电阻率监测点的电阻率，fj表示第j个其他电阻率监测点的电阻率。
[0085]
在该实施例中，所述电阻率监测点的属性特征值根据属性特征预先设定。
[0086]
在该实施例中，所述关联系数越大，对应的电各向异性取值的取值越小。
[0087]
在该实施例中，所有电阻率监测点的电各向异性取值共同组成所述渗流岩体模型的电各向异性取值。
[0088]
在该实施例中所述电各向异性取值用于表示此水渗流岩体的电学性能沿不同方向的差异，差异越大，对应的取值越大。
[0089]
上述设计方案的有益效果是：通过对所述电阻率监测点和电阻率，结合阻率监测点之间的相对位置特征和电阻率监测点位置所处的岩体结构特征，确定所述水渗流岩体模型的电各向异性取值，保证得到电各向异性取值的准确性，为数值模拟提供基础。
[0090]
实施例5
[0091]
基于实施例1的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，如图3所示，s3中，基于所述水渗流岩体模型，确定各个水渗流岩体的结构特征，包括：
[0092]
s301：将所述水渗流岩体模型按照网格划分为多个单元，按照预设分析方法对每个单元进行分析，得到第一单元特征；
[0093]
s302：根据所述第一单元特征，对所述多个单元进行合并，得到整体单元，并根据所述第一单元特征，得到整体单元的第二单元特征；
[0094]
s303：对所述第二单元特征进行融合处理，得到水渗流岩体的结构特征。
[0095]
在该实施例中，将所述水渗流岩体模型按照网格划分为多个单元中每个单元的大小相同。
[0096]
在该实施例中，所述预设分析方法为对渗流岩体模型的每个单元进行扫描，得到单元的外形和结构，并根据扫描图像确定每个单元的材质特征，组成所述第一单元特征，通过先将所述水渗流岩体模型按照网格划分为多个单元进行细化，对得到每个细小的单元进行分析得到第一单元特征，保证了对渗流岩体模型结构分析的精确性，从而得到的特征更加准确。
[0097]
在该实施例中，根据所述第一单元特征，对所述多个单元进行合并，得到整体单元具体为：
[0098]
比较各个单元的第一单元特征的关联度，若关联度大于预设关联度，将对应的单元进行合并，最终得到整体单元，根据第一单元特征，对所述多个单元进行合并，得到整体单元，使得整体单元作为一个整体的水渗流岩体结构，便于了对水渗流岩体模型整体的分析，保证获取结构特征的准确度。
[0099]
在该实施例中，对所述第二单元特征进行融合处理集体为根据所述第二单元特征的关联度，若关联度大于预设关联度，将对应的第二单元特征进行融合，得到整体特征，否则，保留原来的第二单元特征，这样可以简明水渗流岩体的结构特征的描述，同时，又保证水渗流岩体的结构特征可以更好地表示实际的水渗流岩体的结构。
[0100]
上述设计方案的有益效果是：通过将所述水渗流岩体模型按照网格划分为多个单元，按照预设分析方法对每个单元进行分析，得到第一单元特征，通过先将所述水渗流岩体模型按照网格划分为多个单元进行细化，对得到每个细小的单元进行分析得到第一单元特征，保证了对渗流岩体模型结构分析的精确性，从而得到的特征更加准确，根据第一单元特征，对所述多个单元进行合并，得到整体单元，使得整体单元作为一个整体的水渗流岩体结构，便于了对水渗流岩体模型整体的分析，保证获取结构特征的准确度，对所述第二单元特征进行融合处理，可以简明水渗流岩体的结构特征的描述，同时，又保证水渗流岩体的结构特征可以更好地表示实际的水渗流岩体的结构。
[0101]
实施例6
[0102]
基于实施例1的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向
异性数值模拟方法，s3中，利用数值模拟的方法确定电各向异性和结构特征之间的关系包括：
[0103]
基于结构属性，对所述结构特征进行一级标记和二级标记，基于电各向异性取值对所述结构特征进行结果标记，并根据标记结果，得到结果标记与标记组合的对应关系，其中，所述标记组合包括某一种一级标记和二级标记，所述对应关系为结果标记对应标记组合，对应方式为一对多；
[0104]
基于所述电各向异性、结构特征，在所述水渗流岩体模型的基础上进行标记建立静态模型，并基于所述对应关系，赋予所述二级标记的动态调整功能，建立动态模型；
[0105]
对动态模型进行再次数值模拟，确定电各向异性和结构特征之间的关系。
[0106]
在该实施例中，对所述结构特征进行一级标记和二级标记包括基于结构属性，对所述结构特征进行特征分析，得到每个结构属性下的属性子特征，并按照结构属性的等级划分规则，将所述属性子特征进行等级划分，得到属性子特征的等级；根据所述属性子特征，对结构特征进行一级标记，根据所述等级，对结构特征进行二级标记。
[0107]
在该实施例中，对动态模型进行再次数值模拟，确定电各向异性和结构特征之间的关系包括在所述动态模型中按照一级标记，对所述二级标记进行单一性调整，基于所述对应关系，对所述动态模型进行数值模拟，输出单一性调整对应的单一动态结果标记，改变所述一级标记的单一性调整属性，直到得到所有一级标记对应的单一动态结果标记；基于所述所有一级标记及其对应的单一动态结果标记之间的标记关系，对所有一级标记中至少两个一级标记同时进行调整的情况进行再次数值模拟，得到综合结果标记，并以综合结果标记对应的电各项异性数值为横坐标，二级标记为对应的属性子特征取值为纵坐标，建立得到关系曲线图。
[0108]
在该实施例中，不同的结构特征可能导致相同的电各项异性，所以所述对应关系为结果标记对应标记组合，对应方式为一对多。
[0109]
在该实施例中，所述结构属性包括材质、空隙、表面等。
[0110]
在该实施例中，所述一级标记为属性标记，二级标记为在一级标记下对属性的特征数值标记。
[0111]
上述设计方案的有益效果是：首先根据结构属性和电各向异性，明确结构属性与电各向异性之间的对应关系，由于又实际得到的数值对应关系中数值的取值个数有限，只能得到具体的数值对应关系，不能准确得到随机一个取值所对应的关系，因此，根据基于所述电各向异性、结构特征，在所述水渗流岩体模型的基础上进行标记建立静态模型，并基于所述对应关系，赋予所述二级标记的动态调整功能，建立动态模型，进行数值模拟，可明确所有特征取值对应的电各项异性取值，使得不仅可以对现有的水渗流岩体结构进行分析，还可以根据电各向异性的取值后，对水渗流岩体的结构进行预测，根据预测结果及时对水渗流岩体采取应对措施，避免灾害发生。
[0112]
实施例7
[0113]
基于实施例6的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，根据标记结果，得到结果标记与标记组合的对应关系，包括：
[0114]
基于结构属性，对所述结构特征进行特征分析，得到每个结构属性下的属性子特征，并按照结构属性的等级划分规则，将所述属性子特征进行等级划分，得到属性子特征的
等级；
[0115]
根据所述属性子特征，对结构特征进行一级标记，根据所述等级，对结构特征进行二级标记，基于所述电各向异性取值对所述结构特征进行结果标记；
[0116]
将所述一级标记和二级标记作为输入，所述结果标记作为分类标准，获取在相同结果标记对应的一级标记和二级标记的组合，得到结果标记与标记组合的对应关系。
[0117]
在该实施例中，所述属性子特征的等级按照取值进行划分，例如属性子特征取值在第一范围内的属性子特征作为第一等级，属性子特征取值在第二范围内的属性子特征作为第二等级。
[0118]
在该实施例中，所述标记组合包括某一种一级标记和二级标记，所述对应关系为结果标记对应标记组合，对应方式为一对多。
[0119]
上述设计方案的有益效果是：通过根据结构属性，对所述结构特征进行两次标记，明确了属性既属性下的特征取值两个方面对电各项异性的对应关系，为数值模拟提供参考基础。
[0120]
实施例8
[0121]
基于实施例6的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，对动态模型进行再次数值模拟，确定电各向异性和结构特征之间的关系，包括：
[0122]
在所述动态模型中按照一级标记，对所述二级标记进行单一性调整，基于所述对应关系，对所述动态模型进行数值模拟，输出单一性调整对应的单一动态结果标记，改变所述一级标记的单一性调整属性，直到得到所有一级标记对应的单一动态结果标记；
[0123]
基于所述所有一级标记及其对应的单一动态结果标记之间的标记关系，对所有一级标记中至少两个一级标记同时进行调整的情况进行再次数值模拟，得到综合结果标记，并以综合结果标记对应的电各项异性数值为横坐标，二级标记为对应的属性子特征取值为纵坐标，建立得到关系曲线图。
[0124]
在该实施例中，所述结果标记与二级标记的对应方式为一对多，即横坐标的取值有多个纵坐标取值共同决定。
[0125]
在该实施例中，在所述动态模型中按照一级标记，对所述二级标记进行单一性调整，设定所述动态模型的拟合精度，并基于所述对应关系，对所述动态模型进行数值模拟，输出单一性调整对应的满足所述拟合精度的单一动态结果标记。
[0126]
在该实施例中，根据所述关系曲线图，即可确定电各向异性和结构特征之间的关系。
[0127]
在该实施例中，首先，在动态模型的其他特征属性对应的特征取值不变的前提下，对某一特征属性的特征取值进行改变，得到此属性对应的单一动态结果标记，然后在根据单一动态结果标记，对至少两个属性的特征取值进行改变，得到综合结果标记，并以综合结果标记对应的电各项异性数值为横坐标，二级标记为对应的属性子特征取值为纵坐标，建立得到关系曲线图，明确了在多个属性的特征取值综合改变的情况下，依然能够预测得到对应的电各向异性的取值，可以水渗流岩体的结构进行预测，根据预测结果及时对水渗流岩体采取应对措施，避免灾害发生。
[0128]
实施例9
[0129]
基于实施例3的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，基于所述水渗流岩体模型的岩体密度和岩体孔隙率，确定在所述水渗流岩体模型的电阻率监测点，包括：
[0130]
按照所述水渗流岩体模型的岩体密度，对所述水渗流岩体模型进行划分，得到部分水渗流岩体，并基于部分水渗流岩体的岩体密度和面积，设定在所述部分水渗流岩体的监测点个数；
[0131]
对所述部分水渗流岩体再次进行更大精度的岩体密度检测，并根据检测结果，确定在所述部分水渗流岩体的监测区域，其中，所述监测区域对应所述监测点个数；
[0132]
获取所述监测区域中岩体孔隙率最高的位置，作为电阻率监测点。
[0133]
上述设计方案的有益效果是：通过根据水渗流岩体模型的岩体密度和岩体孔隙率，确定在所述水渗流岩体模型的电阻率监测点，使得电阻率监测点的布局合理，从而使得计算得到的电各向异性的取值更加准确，可以更加准确的了解水渗流岩体的电性变化，提高后续数值模拟结果的准确性，可以对水渗流岩体的结构进行预测，根据预测结果及时对水渗流岩体采取应对措施，避免灾害发生。
[0134]
实施例10
[0135]
基于实施例9的基础上，本发明实施例提供一种基于电阻率的水渗流岩体电各向异性数值模拟方法，获取所述监测区域中岩体孔隙率最高的位置，作为电阻率监测点后，还包括：
[0136]
确定所述位置在所述水渗流岩体模型中坐标点，获取所述坐标点在实地的水渗流岩体的实地位置；
[0137]
确定所述实地位置的地质信息，并判断所述地质信息是否满足电阻率的监测环境要求；
[0138]
若是，确定所述位置为最终的电阻率监测点；
[0139]
否则，选取与所述位置最近且满足所述监测环境要求的临近位置作为最终的电阻率监测点。
[0140]
在该实施例中，由于监测电阻率要电极插入泥地中，为了更好地进行电极地监测，要保证电极插入泥地的稳定，因此泥地的选择尤为重要。
[0141]
上述设计方案的有益效果是：通过对实地的水渗流岩体的实地位置的检测，选择检测效果好的位置为最终的电阻率监测点，便于更好地进行电阻率地监测。
[0142]
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。