[发明专利]一种模拟CO2

权利要求书

1.一种模拟CO2注入过程的被动地震监测实验分析方法，其特征在于，包括：

根据预设的多个接收模块对模拟试件进行连续声发射数据采集，获得模拟试件表面多个连续地震记录；

对所述连续声发射数据进行预处理，获得预处理后的预处理数据；

将所述预处理数据分割为预设时间窗长度的片段数据；

对任意两个接收模块之间的所述片段数据进行互相关计算，获得任意接收模块对之间的所述预处理数据对应的互相关函数；

将模拟试件表面范围划分为若干个微小单元，并假设每个所述微小单元为一个假设源，计算每个假设源与任意两个接收模块之间校正过的互相关函数；

将所有接收模块对之间校正过的互相关函数进行叠加，取零时间附近窗口内的最大值作为所述微小单元的的偏移幅值；

计算所有所述微小单位内偏移幅值，获得模拟试件内噪声源分布情况；

计算所有所述预设时间窗长度片段数据的所有微小单元的偏移幅值，获得模拟试件内噪声源时空分布情况。

2.根据权利要求1所述的模拟CO2注入过程的被动地震监测实验分析方法，所述对所述计算每个假设源与任意两个接收模块之间校正过的互相关函数的步骤，包括：

将所述假设源与接收模块之间之间的距离除以区域背景速度v，获得所述假设源信号的预计到达时间，并作为每个接收模块对之间面/体波的理论滞后时间；

根据时移从所述接收模块对之间的互相关函数中除去理论滞后时间，获得所述校正过的互相关函数。

3.根据权利要求1所述的模拟CO2注入过程的被动地震监测实验分析方法，所述将所有接收模块对之间校正过的互相关函数进行叠加，取零时间附近窗口内的最大值作为所述微小单元的的偏移幅值的步骤，包括：

将所述将所有接收模块对之间校正过的互相关函数进行叠加；

取零滞后附近窗口内的最大值作为所述微小单元的偏移幅值：

其中N为台站数，S为第i和j个台站归一化到0-1范围的互相关函数包络幅值，d表示所述假设源单元位置(x，y)与接收模块位置之间的距离。

4.根据权利要求1所述的模拟CO2注入过程的被动地震监测实验分析方法，其特征在于，所述对所述噪声数据进行预处理，获得预处理数据的步骤，包括：去均值、去线性趋势、归一化、谱白化和锥形化处理。

5.根据权利要求1所述的模拟CO2注入过程的被动地震监测实验分析方法，其特征在于，包括：所述预设时长为1秒，包括0.5秒的窗重叠。

6.根据权利要求1所述的模拟CO2注入过程的被动地震监测实验分析方法，其特征在于，所述将所述预处理数据分割为预设时间窗长度的片段数据的步骤还包括：

将所述预处理数据分割为多个预设时长的连续数据；

将所述连续数据分割为若干个预设时间窗长度的所述片段数据。

7.根据权利要求6所述的模拟CO2注入过程的被动地震监测实验分析方法，其特征在于，所述预设时间窗长度为1个小时。

8.根据权利要求1所述的模拟CO2注入过程的被动地震监测实验分析方法，其特征在于，所述模拟试件用于模拟岩心，包括第一砂岩层、第二砂岩层和片岩层，所述片岩层设置于所述第一砂岩层与所述第二砂岩层之间。

9.根据权利要求1所述的模拟CO2注入过程的被动地震监测实验分析方法，其特征在于，所述接收模块包括声发射采集模块和流体压力采集模块，所述声发射采集模块，用于采集所述连续声发射数据；所述流体压力采集模块用于采集注入CO₂流量和流体压力。