[发明专利]地下开采引起的地下水变化分析方法和系统

权利要求书

1.一种地下开采引起的地下水变化分析方法，其特征在于，包括：

获取待分析区域的地下水位信息，所述待分析区域的时移地震数据和所述待分析区域的地下含水层及隔水层的分布信息；所述地下水位信息为所述待分析区域在地下开采的不同时期的地下水位变化信息；

基于所述时移地震数据，对所述待分析区域进行断裂属性的叠前反演计算，得到所述待分析区域的断裂发育特征信息和所述待分析区域的断裂含水信息；

基于所述地下水位信息、所述地下含水层及隔水层的分布信息、所述断裂发育特征信息和所述断裂含水信息，分析得到所述待分析区域的地下水变化信息；所述地下水变化信息包括：地下水位变化原因和地下水位随地下开采进行的周期变化信息。

2.根据权利要求1所述的方法，其特征在于，获取待分析区域的地下水位信息，包括：

获取所述待分析区域内多个探测井中的多个初始地下水位信息；

对所述多个初始地下水位信息进行插值计算，得到所述待分析区域的地下水位信息。

3.根据权利要求1所述的方法，其特征在于，获取所述待分析区域的地下含水层及隔水层的分布信息，包括：

获取所述待分析区域的钻井及测井资料信息；

基于所述钻井及测井资料信息，采用地质建模的方法，得到所述待分析区域的地下含水层及隔水层的分布信息。

4.根据权利要求1所述的方法，其特征在于，基于所述时移地震数据，对所述待分析区域进行断裂属性的叠前反演计算，得到所述待分析区域的断裂发育特征信息和所述待分析区域的断裂含水信息，包括：

通过算式组：

R_p(θ)＝ωS_TT(θ)+ωS_NN(θ)，

μ＝ρV_S²，

对所述待分析区域进行断裂属性的叠前反演计算，得到地层参数信息；其中，P为关于所述时移地震数据的射线参数，θ为关于所述时移地震数据的入射角，ω为关于所述时移地震数据的角频率，R_p为关于所述时移地震数据的反射系数；所述地层参数信息包括：法向断裂参数S_T，切向断裂参数S_N，纵波速度V_P，横波速度V_S，质量密度ρ；

基于所述地层参数信息，得到所述待分析区域的断裂发育特征信息和所述待分析区域的断裂含水信息。

5.根据权利要求4所述的方法，其特征在于，基于所述时移地震数据，对所述待分析区域进行断裂属性的叠前反演计算，还包括：

基于所述时移地震数据，获取断裂属性反演的先验信息；

基于所述先验信息对断裂属性反演过程进行约束，其中，约束算式为：其中，Q为反演目标函数，S^obs为所述时移地震数据，m为所述地层参数信息，m₁为所述先验信息，S为基于所述地层参数信息进行正演模拟之后得到的正演数据，λ为正则化参数。

6.根据权利要求5所述的方法，其特征在于，基于所述时移地震数据，获取断裂属性反演的先验信息，包括：

基于所述时移地震数据，获取预设纵波速度、预设横波速度和预设质量密度；

获得等效断裂孔隙度模型，然后通过如下算式计算得到预设法向断裂参数S_N'和预设切向断裂参数S_T'：

其中，预设法向断裂参数和预设切向断裂参数的单位为m/Gpa,φ_f为关于所述时移地震数据的孔隙度；

将所述预设纵波速度、预设横波速度、预设质量密度、预设法向断裂参数和预设切向断裂参数确定为断裂属性反演的先验信息。