基于联合性升序边界元法的弹性层状介质裂缝问题分析.pdfVIP

第45卷第6期湖北科技学院学报Vol.45,No.6

2025年12月JournalofHubeiUniversityofScienceandTechnologyDec.2025

引文格式：周嘉麒，龙恭博，徐朝阳，等.基于联合性升序边界元法的弹性层状介质裂缝问题分析[J].湖北科技学院学报，2025，

45(6):108-115.

基于联合性升序边界元法的弹性层状

介质裂缝问题分析

周嘉麒，龙恭博*，徐朝阳,姚华楠，张皖琦，王启明，张证法

（武汉工程大学机电工程学院，湖北武汉430205）

摘要：本文提出一种联合性升序边界元法（JBEM）用于解决弹性层状介质中的裂缝问题，计算裂

缝的张开位移，以克服传统直接边界元法（DBEM）求解效率低和位移不连续法（DDM)适用性差的

局限性。JBEM在求解过程中使用DBEM量化表征层状介质各层界面对裂缝张开位移的影响，同时

使用DDM从层状介质底层开始，利用各层界面连续条件升序至顶层，逐层计算裂缝张开位移。

JBEM最终得到的数学方程形式较简洁，只包含需要求解的裂缝张开位移这一未知量，不包含各层

界面上的未知量。数值算例表明,JBEM比无法解决一般多层弹性介质裂缝问题的DDM适用性更

好。同时，随着层数增加，涉及的元素数量越多，JBEM的效率优势越明显，且计算效率比DBEM快。

关键词：边界元；零阶元素；裂缝问题；弹性介质；数值算例

中图分类号:TE357.1文献标识码：A文章编号：2095-4654(2025)06-0108-08

开放科学（资源服务）标识码（OSID）：

水力压裂是有效开采我国储量丰富的非常规油无法使用边界元法，限制了该方法的适用性。因此，

使用边界元法解决弹性层状介质裂缝问题的关键就

气资源的主要技术[1-2]。在压裂过程中,水力裂缝受

压裂液压力作用在弹性层状岩石介质中张开。因此，是找到基础解析解，同时精确、高效地建立边界元方

水力压裂数值模拟的核心要素之一就是在已知压力程，从而求解裂缝张开位移[9-10]

作用下计算弹性层状介质中的裂缝张开位移[341]。边界元法有两个分支：位移不连续法（displace-

在众多计算裂缝张开位移的数值方法中，边界mentdiscontinuitymethod，DDM）和直接边界元法

元法得到了广泛地使用。相比于有限元法、有限差其中,

分法这些需要对研究对象整体进行离散化处理的方DDM利用裂缝张开位移和不连续位移之间的物理

法,边界元法只需对研究对象的边界进行离散化处相似性,以常数不连续位移的应力场和位移场为基

理。边界元法通过离散化得到的元素个数明显少于础解析,在解决裂缝问题时，无需考虑界面元素，因

有限元法和有限差分法，因此提高了计算效率[5-6]。此它是计算裂缝张开位移最高效的数值方法[。

另一方面,边界元法是一种以基础解析解为基础的但是，该方法只适用于均匀介质或两均匀介质[12-3]。

数值方法，所以其计算结果比有限元法等其他数值当介质层数大于2时，无法找到对应的基础解析解，

方法更精确[7-1]。然而，当基础解析解不存在时，就也就无法使用该方法。另一方面,DBEM以弹性力

*收稿日期：2024-10－09

基金项目：武汉工程大学第十五届研究生教育创新基金项目（CX2023215）；湖北省教育厅重点项目“水力压裂过程中岩

石渗透率的分形研究”（。

通信作者：龙恭博，邮箱：longgongbo@163.com

2025年第6期周嘉麒，龙恭博，徐朝阳，姚华楠，张皖琦，王启明，张证法：基于联合性升序边界元法的弹性层状介质裂缝问题分析

学中的经典开尔文解为基础解析解，能够解决多层一、裂缝问题模型的推导

弹性层状介质中的裂缝问题[13]。但是，在使用