大多数初始地应力场的反演均将对象岩体简化为各向同性弹性体，但是对于具有层状构造、倾角多变、力学性质方向差异较大的含煤地层，这样简化就会导致反演结果误差很大，同时这种多层交互岩体在计算建模、参数取值等过程中也存在难度。为了提高具有复杂地质条件的煤系地层初始地应力场反演精度，提出了基于FLAC3D横观各向同性弹性模型和各向同性弹性模型相结合的地应力反演方法。同时为解决多层交互、倾角多变等煤系地层的计算建模、参数选取等问题，根据岩层交互规律、岩性特征，通过对复杂地层进行地质概化，建立与实际岩层基本一致的计算模型，进行数值试验，得到了用于数值计算的多层岩体等效力学参数。根据数值试验结果，建立了局部区域的横观各向同性弹性模型和各向同性弹性模型相结合的计算模型，利用偏最小二乘回归方法拟合地应力场并将地应力场精确加载至计算模型并实现力的平衡，最终完成初始地应力场的加载。结果表明，采用建立数值模型来计算各层岩体的等效力学参数的方法优化了建模过程，提高了数值计算效率。对比两种方法反演结果，采用横观各向同性模型进行计算更能体现具有层状构造的煤系地层力学性质的方向差异性，能有效提高局部区域的地应力拟合精度，其地应力反演结果更符合实测值。

In most inversion on initial ground stress field,the object rock mass is simplified to isotropic elastomer,but for coal-bearing strata with layered structure,varied dip angle and different directions of mechanical properties,such simplification will lead to a large inversion error. Also,it is difficult to model and select parameters for such multilayer interactive rock mass. In order to improve the inversion accuracy of initial ground stress field in coal measure strata with complicated geological conditions,a method of ground stress inversion based on FLAC3D transverse isotropic elas- tic model and isotropic elastic model was proposed. At the same time,to solve the problems of computational modeling and parameter selection in the coal strata with multiple layer interactions and changeable angles,according to the laws of rock interbed,lithologic characteristics,through the study of the geological sketch of complex formation,a calculation model is established,which is basically consistent with the actual rock situation. A numerical experiment is conducted and the multi-layer equivalent rock mass mechanics parameters for the numerical calculation is obtained. According to the numerical test results,the transverse isotropic elastic model and the isotropic elastic model of the local region are established,and the partial least square regression method is used to fit the ground stress field,and the ground stress field is precisely loaded into the calculation model and the force balance is realized,and finally the initial ground stress field is loaded. The results show that the method of constructing numerical model to calculate the equivalent mechani- cal parameters of each layer of rock mass optimizes the modeling process and improves the efficiency of numerical cal- culation. By comparing the two inversion results,the horizontal isotropic model can better reflect the directional differ- ence of the mechanical properties of coal measures strata with layered structure,and can effectively improve the ground stress fitting accuracy in local areas,and the ground stress inversion results are more in line with the measured values.

1 基本原理

　　1.1 多元回归原理

　　1.2 FLAC3D横观各向同性模型

2 应力场加载至模型

3 潘一矿东区地应力反演分析

　　3.1 工程概况

　　3.2 地质模型概化

　　3.3 材料等效参数计算

　　3.4 地应力反演

4 结论