具体实施方式

以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。

图2是本发明一种实施方式提供的水力压裂激发Kraukls波模拟系统的系统结构图。如图2所示，本发明实施方式提供一种水力压裂激发Kraukls 波模拟系统，所述系统包括：采集单元10，用于获取人工裂缝的裂缝模型；处理单元20，用于：根据预设剖分规则进行裂缝模型离散化处理，获得网格单元模型；扫描所有网格单元，筛选出与裂缝开口相关的网格单元，并确定与裂缝开口相关的网格单元的属性参数；根据与裂缝开口相关的网格单元的属性参数和预设位移载荷计算公式逐一计算各与裂缝开口相关的网格单元对应裂缝开口处的位移载荷；根据位移载荷计算结果对所有网格单元进行位移载荷装配，并将装配完成的位移载荷加到常规有限元的位移项上，获得最终有限元方程；输出单元30，用于求解所述有限元方程，获得并输出模拟结果。

图1是本发明一种实施方式提供的水力压裂激发Kraukls波模拟方法的方法流程图。如图1所示，本发明实施方式提供一种水力压裂激发Kraukls 波模拟方法，所述方法包括：

步骤S10：获取所述人工裂缝的裂缝模型，并根据预设剖分规则进行所述裂缝模型离散化处理，获得网格单元模型。

具体的，采集单元10进行水力压裂后产生的人工裂缝模型采集，包括人工裂缝的空间形状和投影形状。因为裂缝的空间展布存在着复杂性，无法直接用形状模型进行后续模拟，所以需要使用网格单元模型，进行人工裂缝形状拟合，以保证后续移动载荷计算和装配的简略性，提高Kraukls波模拟效率。优选的，将人工裂缝模型中离散的点进行三角网格剖分，将形状复杂且无规律的裂缝形状拟合成多个相接的三角形，便于后续对裂缝进行形状分析。判断一个三角网格化产生的结果是否为最优，一般认为，三角网格化后得到的三角形中细长和尖锐的三角形越少越好，而类似与等边三角形和锐角三角形越多越好。而Delaunay等参三角网格单元模型因为其具有空圆特性和最大化最小角特性，使得避免了狭长三角形的产生，所以优选在本发明中选用Delaunay等参三角网格单元模型进行人工裂缝原始形状拟合。这种网格的尺寸可以随着拟合区域的尺寸灵活变化，而且可以提供双二次插值函数，这样可以保证所求解的数值解具有双二阶精度。通过Delaunay等参三角网格单元模型剖分获得的网格单元，每一个网格单元包括6个节点，其中 3个角节点，对应着三角形的三个顶点；另外3个是边节点，对应三角形三条边的中点。在网格剖分的过程中，处理单元20实时进行采集单元10获取的人工裂缝的空间形态扫描，包括裂缝的空间位置、裂缝孔径、裂缝面光滑程度以及裂缝走向。其中，裂缝的空间位置用于确认网格单元在空间中的位置，裂缝孔径以及裂缝表面光滑程度则用以确定网格单元的大小尺寸。裂缝孔径越小、裂缝面粗糙程度越大，则网格的尺寸要求也就更精细。而裂缝的走向则用以确认三角形内角的角度分布。

步骤S20：扫描所有网格单元，筛选出其中与裂缝开口相关的网格单元，并确定所述与裂缝开口相关的网格单元的属性参数。

具体的，处理单元20获取剖分后的所有网格单元组合形成的网格单元模型，然后进行网格单元模型筛选，过滤掉其中远离裂缝开口的网格单元，仅保留与裂缝开口相关的网格单元。然后进行这些与裂缝开口相关网格单元的属性参数获取。在步骤S10中，进行三角网格剖分时，为避免干扰，获得裂缝开口相关网格单元时，需要保证与裂缝开口相关网格单元有且只有一条边与裂缝开口重合，以便于后续在裂缝开口处加载位移载荷。处理单元20 将筛选出来的与裂缝开口相关的网格单元进行整合，然后逐一获取这些网格单元中与裂缝开口重合的那一条边，以获取这些边上所有的网格单元节点的编号和坐标。这些参数在获取进行位移载荷计算时，需要根据坐标位置进行向量确定，以保证位移载荷计算准确。

步骤S30：根据所述与裂缝开口相关的网格单元的属性参数和预设位移载荷计算公式逐一计算各与裂缝开口相关的网格单元对应裂缝开口处的位移载荷。

具体的，位移载荷的加载原则为，与裂缝形成原理相同，在裂缝开口面上，从裂缝的几何中心向两个侧裂缝壁的位置，位移载荷的值递减。为了表征这种递减关系，需要在载荷方程中定义一个空间尺度函数。该空间尺度函数值在裂缝中心处达到峰值1，而在两个侧裂缝壁的位置，则为0。在介于裂缝中心和两个侧裂缝壁之间的位置，尺度函数对称光滑递交。通过空间尺度函数建立，体现了位移载荷从裂缝中心向两个侧裂缝壁递减的特性。优选的，本发明根据上述规则建立的空间尺度函数关系式为：

其中，为落在开口处节点相对于裂缝中心点的位置向量；为表征裂缝开口的向量。其中，定义为由裂缝中心点指向该节点的向量，其具体表达式为：

其中，(x₀,y₀)为裂缝的中心店坐标；(x_node,y_node)为落在裂缝开口处的节点坐标。为表征裂缝开口的向量，由裂缝开口的起点志向裂缝开口终点的向量，其表达式为：

其中，为裂缝开口按照逆时针方向的起点；为裂缝开口按照逆时针方向的终点。除了需要空间尺度函数，Krauklis波的生成和传输需要用时间域来进行表征，时间域的时间序列的性质与特征则是各种数学滤波的理论基础，因此，对地震勘探资料作时间域内的研究，在现代地震勘探中具有重要意义。裂缝在开裂时会激发地震波，Krauklis波就是在这过程中激发的一种导波。所以需要定义一个时间与子波函数，优选的，本发明提出的时间域子波函数表达式为：

其中，A₀为振幅系数；τ为子波函数的时间尺度因子；t₀为子波函数的中心时刻；t为时间序列。而位移载荷方程便是空间尺度函数与时间域子波函数相乘获得，所以最终的位移载荷方程为：F_time＝F_space·F_time，即：

通过上述关系式，逐一进行各与裂缝开口相关网格单元中，开口线上的位移载荷计算。

步骤S40：根据位移载荷计算结果对所有网格单元进行位移载荷进行装配，并将装配完成的位移载荷加到常规有限元的位移项上，获得最终有限元方程。

具体的，完成每一个与裂缝开口相关网格单元位移载荷计算后，需要将这些计算获得的位移载荷装配到整个网格单元模型中，获得完整的人工裂缝的位移载荷。有限元节点装配方式就是将创建的有限元模型整合装配形成完整部件或装配体。所以本发明选用有限元节点装配方式进行位移载荷装配。其计算公式如下：

其中，F为最终的位移载荷矩阵；N_opening为落在裂缝开口处的节点个数； C_i为描述第i个落在裂缝开口处节点的编号位置向量；F_i为第i个落在裂缝开口处节点的位移载荷。在一种可能的实施方式中，位移载荷不局限于仅在与开口裂缝相关网格单元中进行，可以在所有网格单元中进行。获取完整的网格单元模型，其中，与裂缝开口相关的网格单元通过上述装配公式进行装配结果计算。而处理单元20识别到的与裂缝开口无关的网格单元，直接默认结果为0，将所有结果整合到一起，便获得了最后的位移载荷。完成位移载荷装配后，直接将最终的位移载荷加到有限元方程的位移项上，作为有限元方程的已知条件。

步骤S50：求解所述有限元方程，获得并输出模拟结果。

具体的，获得最终的位移载荷后，在有限元方程中，位移项变为已知条件。输出单元30将装配获得的位移载荷值直接放入到有限元方程中进行求解，获得求解结果，当前结果便为Kraukls波的模拟结果。

本发明实施方式还提供一种计算机可读储存介质，该计算机可读存储介质上储存有指令，其在计算机上运行时使得计算机执行上述的水力压裂激发 Kraukls波模拟方法。

本领域技术人员可以理解实现上述实施方式的方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序存储在一个存储介质中，包括若干指令用以使得单片机、芯片或处理器(processor)执行本发明各个实施方式所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM， Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

以上结合附图详细描述了本发明的可选实施方式，但是，本发明实施方式并不限于上述实施方式中的具体细节，在本发明实施方式的技术构思范围内，可以对本发明实施方式的技术方案进行多种简单变型，这些简单变型均属于本发明实施方式的保护范围。另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明实施方式对各种可能的组合方式不再另行说明。

此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明实施方式的思想，其同样应当视为本发明实施方式所公开的内容。