一种用于复杂裂缝系统波场正演模拟的跨尺度数值方法技术方案

本发明专利技术涉及一种复杂裂缝系统波场正演模拟的跨尺度数值方法，包括以下步骤：1)生成并获取RVE晶格节点输入参数；2)获取RVE弹性模量并量化RVE尺度特征参数；3)模拟复杂裂缝系统波场；4)对照超声实验探索尺度效应关联机制。本发明专利技术可以实现一种能够计及更小样本尺度微观裂缝诱发波场异常特征的跨尺度数值方法，对阐明不同样本尺度裂缝系统波场特征的关联机制与叠合效应具有重要的理论意义，对推动非常规油气藏的勘探与开发具有广泛的实际应用价值。值。值。

【技术实现步骤摘要】
一种用于复杂裂缝系统波场正演模拟的跨尺度数值方法


[0001]专利技术涉属于含油气裂缝储层地震波传播及各向异性分析
，具体的说，涉及一种复杂裂缝系统波场正演模拟的跨尺度数值方法。

技术介绍

[0002]裂缝型储层中多构造期次发育天然裂缝与后期诱发新生裂缝等构成多成因、多尺度复杂裂缝系统。以页岩气储层为代表的典型复杂裂缝系统，主要包含纳米
‑
微米级粒间孔隙、微米级
‑
毫米级天然裂缝以及更大尺度的水力压裂缝，尺度横跨纳米到千米等多个量级。然而，实际测量只能针对各自所在的样本尺度，如何综合不同样本尺度的测量结果以得到更全面的复杂裂缝系统波场信息，是业界亟待解决的关键基础科学问题。因此，如果能够从孔隙尺度展开研究，建立跨微观到宏观多个样本尺度的桥梁，提出一种跨尺度研究方法，则具有重要的理论意义与广泛的实际应用价值。
[0003]但相关传统研究存在如下局限性：(1)裂缝仿真多基于VTI(垂直横向各向同性)、HTI(水平横向各向同性)与TTI(倾斜横向各向同性)等效裂缝介质或理想规则形态裂隙结构，过度简化了实际裂缝系统的复杂非均质性与多尺度性；(2)忽略了不同样本尺度裂缝系统波场特征的关联机制与叠合效应，缺乏一种理想的跨尺度数值方法。

技术实现思路

[0004]针对上述问题，本专利技术的目的是提供一种跨尺度数值方法，深入探索不同样本尺度裂缝系统波场特征的关联机制与叠合效应，发展相关前沿理论并为实际工程应用提供必要理论基础与技术支持。
[0005]为实现上述目的，本专利技术采取以下技术方案：一种复杂裂缝系统波场正演模拟的跨尺度数值方法，包括以下步骤：
[0006]1)生成并获取代表体积单元(RVE)晶格节点输入参数：选取毫米级的圆柱页岩样本开展微米CT原位单轴压缩
‑
拉伸实验，获取毫米级的圆柱页岩样本的内部微结构快照与力学特征响应曲线，将毫米级的圆柱页岩样本离散为微米级RVE，根据内部微结构快照提取微结构信息，建立对应的修正晶格弹簧(LSM)
‑
离散裂缝网络(DFN)数值模型，模拟单轴压缩
‑
拉伸实验，基于最小二乘法构造数值模拟数据与微米CT原位单轴压缩
‑
拉伸实验数据的目标函数从而获取RVE晶格节点输入参数；
[0007]2)获取RVE弹性模量并量化RVE尺度特征参数：建立每个微米级代表体积单元尺度(RVE)的修正晶格弹簧
‑
离散裂缝网络模型，开展三轴压缩和双轴压缩数值模拟分步识别量化RVE弹性模量与尺度特征参数；
[0008]3)模拟复杂裂缝系统波场：选取厘米级的圆柱页岩试样作为开展超声实验的研究对象，提取裂缝结构生成DFN，据此建立该岩样的修正晶格弹簧
‑
离散裂缝网络模型(LSM
‑
DFN模型)，反解出RVE等效晶格单胞中晶格节点，即步骤1中的跨尺度模型中晶格节点的输入参数，以雷克子波作为震源函数模拟超声信号，通过求解晶格节点的局部广义力与广义
位移，进而得到整个研究区域的波场信息，模拟复杂裂缝系统波场；
[0009]4)对照超声实验探索尺度效应关联机制：取厘米级的圆柱页岩试样开展不同频率的超声实验，对照数值模拟数据，综合波场记录，基于机械工程学科中的信号处理分析方法深入探索裂缝诱发波场异常特征的尺度效应与关联机制。
[0010]其中，所述步骤1)的具体过程如下：
[0011]11)选取毫米级的圆柱页岩样本开展微米CT原位单轴压缩
‑
拉伸实验，获取毫米级样本的内部微结构快照与力学特征响应曲线；
[0012]12)将毫米级的圆柱页岩样本离散为微米级代表体积单元，根据内部微结构快照提取微结构信息，建立微米级代表体积单元尺度的修正晶格弹簧
‑
离散裂缝网络(LSM
‑
DFN)数值模型；
[0013]13)建立对应毫米级的圆柱页岩的修正晶格弹簧
‑
离散裂缝网络(LSM
‑
DFN)数值模型，模拟单轴压缩
‑
拉伸实验，基于最小二乘法构造数值模拟数据与微米CT原位单轴压缩
‑
拉伸实验数据的目标函数如下式所示，
[0014][0015]式中，x是待校准晶格节点的输入参数集，N是用于校准的已有实验数据数量，d
i
是第i个数据点，w
i
是对应第i个数据点的权重因子，s
i
(x)是对应第i个数据点的数据集x作为输入参数得到的模拟结果，r
i
是对应第i个数据点的模拟结果与实验数据的残差；
[0016]14)采用NL2SOL算法最小化目标函数，经过多次迭代校准得到最优输入参数数据集x。
[0017]其中，所述步骤2)的具体过程如下：
[0018]21)建立每个RVE的修正LSM
‑
DFN模型；
[0019]22)开展三轴压缩数值模拟，通过以下公式得到每个RVE的λ、μ、等传统弹性参数，
[0020][0021][0022]式中，σ是应力张量，是应变张量的对称部分，I是单位矩阵，是临界体积应变，是RVE加载过程中晶格节点体积最大占比，是初始RVE的晶格节点体积占比，ε
V
是体积应变；
[0023]23)开展双轴压缩数值模拟，通过以下公式识别材料尺度特征参数，
[0024][0025]式中，下标i表示第i个数据点，N
σ
表示加载数据数目，σ0是参考应力值，用于无量纲化目标函数，σ
11i
(s)和分别为第i个数据点的轴向应力的模拟值和实验数据值；N
H
表示
出现剪切带数据值数目，H0是剪切带厚度，H
i
(s)和分别为第i个数据点的剪切带厚度的模拟值和实验数据值。
[0026]其中，所述步骤3)的具体过程如下：
[0027]31)选取厘米级的圆柱页岩试样作为开展超声实验的研究对象，基于裂缝结构生成DFN，据此建立该岩样的修正LSM
‑
DFN模型，其中LSM的晶键长度与RVE尺寸需满足一定比例关系以保证晶格单胞与RVE体积相等；
[0028]32)基于能量守恒原理与包含材料尺度特征参数的广义连续介质力学偶应力理论，建立修正LSM晶格节点的输入参数与对应等效晶格单胞材料尺度特征参数间的量化关系，具体公式如下所示，
[0029]U
cell
＝U
continuum
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0030][0031][0032][0033]m＝2l2μχ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0034][0035][0036]式中，U
cell
、U
contin本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种复杂裂缝系统波场正演模拟的跨尺度数值方法，其特征在于，包括以下步骤：1)生成并获取RVE晶格节点输入参数：选取毫米级的圆柱页岩样本，获取毫米级的圆柱页岩样本的内部微结构快照与力学特征响应曲线，将毫米级的圆柱页岩样本离散为微米级代表体积单元(RVE)，根据内部微结构快照提取微结构信息，建立毫米级的修正晶格弹簧
‑
离散裂缝网络数值模型，从而获取RVE晶格节点输入参数；2)获取RVE弹性模量并量化RVE尺度特征参数：基于步骤1)中的微米级代表体积单元，建立每个微米级代表体积单元尺度的修正晶格弹簧
‑
离散裂缝网络模型，开展三轴压缩和双轴压缩数值模拟分步识别量化RVE弹性模量与尺度特征参数；3)模拟复杂裂缝系统波场：选取厘米级的圆柱页岩试样作为开展超声实验的研究对象，据此建立厘米级的修正晶格弹簧
‑
离散裂缝网络模型，反解出RVE等效晶格单胞中晶格节点，所述RVE等效晶格单胞中晶格节点与的跨尺度模型中晶格节点的输入参数相一致，以雷克子波作为震源函数模拟超声信号，通过求解晶格节点的局部广义力与广义位移，进而得到整个研究区域的波场信息，模拟复杂裂缝系统波场；4)对照超声实验探索尺度效应关联机制：取厘米级的圆柱页岩试样开展不同频率的超声实验，对照数值模拟数据，综合波场记录，基于机械工程学科中的信号处理分析方法深入探索裂缝诱发波场异常特征的尺度效应与关联机制。2.根据权利要求1所述的复杂裂缝系统波场正演模拟的跨尺度数值方法，其特征在于，所述步骤1)的具体过程如下：11)选取毫米级的圆柱页岩样本开展微米CT原位单轴压缩
‑
拉伸实验，获取毫米级样本的内部微结构快照与力学特征响应曲线；12)将毫米级的圆柱页岩样本离散为微米级代表体积单元，根据内部微结构快照提取微结构信息，建立微米级代表体积单元尺度的修正晶格弹簧
‑
离散裂缝数值模型；13)建立对应毫米级的圆柱页岩的修正晶格弹簧
‑
离散裂缝网络数值模型，模拟单轴压缩
‑
拉伸实验，基于最小二乘法构造数值模拟数据与微米CT原位单轴压缩
‑
拉伸实验数据的目标函数如下式所示，式中，x是待校准晶格节点的输入参数集，N是用于校准的已有实验数据数量，d
i
是第i个数据点，w
i
是对应第i个数据点的权重因子，s
i
(x)是对应第i个数据点的数据集x作为输入参数得到的模拟结果，r
i
是对应第i个数据点的模拟结果与实验数据的残差；14)采用NL2SOL算法最小化目标函数，经过多次迭代校准得到最优输入参数数据集x。3.根据权利要求1所述的复杂裂缝系统波场正演模拟的跨尺度数值方法，其特征在于，所述步骤2)的具体过程如下：21)建立每个RVE的修正LSM
‑
DFN模型；22)开展三轴压缩数值模拟，通过以下公式得到每个RVE的λ、μ、等传统弹性参数，
式中，σ是应力张量，是应变张量的对称部分，I是单位矩阵，是临界体积应变，是RVE加载过程中晶格节点体积最大占比，是初始R...

【专利技术属性】
技术研发人员：刘宁，李吴霜，
申请(专利权)人：北京化工大学，
类型：发明
国别省市：