基于钻孔数据的地层力学参数有限元建模方法及终端设备

技术特征：
1.一种基于钻孔数据的地层力学参数有限元建模方法，其特征在于，包括以下步骤：s1、建立无地层参数的有限元模型；采集工程现场钻孔数据；所述工程现场钻孔数据包括各钻孔位置及编号、各钻孔内各土层的埋深及力学参数；s2、离散化所述现场钻孔数据，并对所述有限元模型赋予相应参数；s3、根据有限元模型单元间的相关作用关系，计算非钻孔位置单元的参数；s4、更新并分类非钻孔位置地层参数，并将更新后的非钻孔位置地层参数赋予有限元模型。2.根据权利要求1所述的基于钻孔数据的地层力学参数有限元建模方法，其特征在于，步骤s1中，建立无地层参数的有限元模型的具体实现过程包括：设置有限元模型的几何尺寸大小，完成有限元模型整体尺寸的几何建模，得到有限元整体几何模型；对有限元整体几何模型进行网格划分，并对划分后的单元和节点进行编号；导出所述单元和节点的对应关系以及各节点的坐标值。3.根据权利要求2所述的基于钻孔数据的地层力学参数有限元建模方法，其特征在于，所述单元和节点的对应关系存储于矩阵a中，所述矩阵a第一列为单元编号，其余列为节点编号；矩阵a每一行中第一个数字为单元编号，其余数字均为该单元所对应的节点编号。4.根据权利要求2所述的基于钻孔数据的地层力学参数有限元建模方法，其特征在于，所述各节点的坐标值存储于矩阵b中，矩阵b第一列为节点编号，其余列为坐标值；其中矩阵b每一行中第一个数字为节点编号，其余数字均为该节点编号对应的坐标值。5.根据权利要求1所述的基于钻孔数据的地层力学参数有限元建模方法，其特征在于，步骤s2中，对钻孔数据离散化，并对所述有限元模型赋予相应参数的具体实现过程包括：计算所述有限元模型中各单元中心点的坐标值；根据所述各单元中心点的坐标值，查找平面坐标等于各钻孔平面坐标的单元，若没有与某钻孔平面坐标相同的单元，则查找与该钻孔平面坐标距离最近的单元，将查找的所有单元组成矩阵e；对比钻孔内的土层埋深和单元的埋深坐标信息，将钻孔相应埋深的参数赋予所述矩阵e。6.根据权利要求5所述的基于钻孔数据的地层力学参数有限元建模方法，其特征在于，步骤s3的具体实现过程包括：将各单元中心点的坐标值存储于矩阵d中；将所述矩阵d与矩阵e做差，得到非钻孔位置单元的中心坐标集合f。7.根据权利要求1～6之一所述的基于钻孔数据的地层力学参数有限元建模方法，其特征在于，步骤s4的具体实现过程包括：将非钻孔位置第i个单元参数f
ij
与钻孔位置每个单元的参数相减，挑选出最小的差值，将该差值最小的钻孔位置单元参数作为非钻孔位置第i个单元参数，代替f
ij
；遍历完所有非钻孔位置单元，完成全部非钻孔位置单元参数的更新；更新后的每个非钻孔位置单元的参数都是钻孔位置单元参数的一种，完成单元参数分类，即实现地层分类；将更新好的单元参数按照单元编号赋予到有限元模型中，得到全地层参数的有限元模型；
其中，e
mj
为钻孔位置单元矩阵e中第m个单元的第j个力学参数；n为钻孔位置单元矩阵e的总个数；p
im
为计算非钻孔位置第i个单元参数时钻孔位置单元矩阵e中第m个单元参数所占权重，e中第m个单元参数所占权重，dx
im
、dy
im
分别为钻孔位置第m个单元和非钻孔位置第i个单元之间的水平向和竖向距离，dx
il
、dy
il
分别为钻孔位置第l个单元和非钻孔位置第i个单元之间的水平向和竖向距离，k为比例系数，n为钻孔位置单元矩阵e中的单元个数。8.根据权利要求7所述的基于钻孔数据的地层力学参数有限元建模方法，其特征在于，当有限元模型为三维模型时，p
im
通过下式计算：9.一种终端设备，包括存储器、处理器及存储在存储器上的计算机程序；其特征在于，所述处理器执行所述计算机程序，以实现权利要求1～8之一所述方法的步骤。10.一种计算机可读存储介质，其上存储有计算机程序/指令；其特征在于，所述计算机程序/指令被处理器执行时实现权利要求1～8之一所述方法的步骤。

技术总结
本发明公开了一种基于钻孔数据的地层力学参数有限元建模方法及终端设备，包括以下步骤：1)建立无地层参数的有限元模型；2)工程现场钻孔数据采集；3)钻孔数据离散化并对有限元模型赋予相应参数；4)根据单元间相互作用关系求解非钻孔位置单元参数；5)将非钻孔位置地层参数更新并分类，随后赋予到有限元模型。该发明避免了直接在有限元模型进行地层几何划分所带来的网格畸变等问题，可根据已有钻孔数据在网格均一的有限元模型中生成不同地层参数，以此反映不同地层起伏变化规律，从而实现带倾斜地层的有限元数值模拟计算。斜地层的有限元数值模拟计算。斜地层的有限元数值模拟计算。


技术研发人员：孙晓贺 施成华 雷明锋 郑可跃 王祖贤 张轩煜 康瑞强 贾朝军 彭铸
受保护的技术使用者：中南大学
技术研发日：2022.09.14
技术公布日：2022/12/16