1.一种煤层瓦斯压力和渗透率同步反演方法，其特征在于，在井下施工钻孔并封孔，测试孔口瓦斯流量，构建双重孔隙/双重渗透率气固耦合数学模型，并通过试验确定双重孔隙/双重渗透率气固耦合数学模型的参数，生成地质模拟几何模型结合地质模拟几何模型的边界条件计算得出对应钻孔瓦斯流量的模拟值，利用钻孔瓦斯流量的实测值和模拟值构建目标函数，通过代理优化算法反演得出煤层原始瓦斯压力和渗透率；

具体步骤如下：

a.在目标煤层施工测试钻孔，钻孔施工完成后对钻孔进行密封，通过孔口流量计测试不同时间的钻孔瓦斯流量值；

b.构建双重孔隙/双重渗透率气固耦合数学模型，并基于实验室和现场试验，确定双重孔隙/双重渗透率气固耦合数学模型的所有参数；

c.根据钻孔布置和钻孔周围煤岩层赋存情况，建立目标煤层相应的地质模拟几何模型并确定地质模拟几何模型的边界条件；

d.将目标煤层原始瓦斯压力和渗透率作为输入变量，并根据上述双重孔隙/双重渗透率气固耦合数学模型、双重孔隙/双重渗透率气固耦合数学模型参数、地质模拟几何模型及地质模拟几何模型边界条件，调用COMSOL Multiphysics with MATLAB软件中的固体力学和达西定律模块仿真计算出对应钻孔瓦斯流量的模拟值；

e.根据钻孔瓦斯流量实测值和模拟值构建目标函数，给定优化参数约束条件和优化停止条件，调用COMSOL Multiphysics with MATLAB软件中的代理优化求解器反演得出煤层原始瓦斯压力和渗透率；

f.根据煤层原始瓦斯压力以及现场采集到的信息：煤层孔隙率、煤层密度、煤层温度和煤体吸附常数，进一步计算出煤层瓦斯含量，根据煤层原始瓦斯压力、煤层瓦斯含量预测出地质模拟几何模型范围内煤层的突出危险性。

2.根据权利要求1所述煤层瓦斯压力和渗透率同步反演方法，其特征在于，对目标煤层施工的顺层钻孔或穿层钻孔，需要保证施工时钻孔50m范围内的煤层无地质构造、含水层、溶洞，且不受采动、瓦斯抽放及其它人为卸压的影响，顺层钻孔长度大于50m，且钻孔密封后钻孔前端保留1‑2m的气室。

3.根据权利要求1所述煤层瓦斯压力和渗透率同步反演方法，其特征在于，所述的双重孔隙/双重渗透率气固耦合数学模型为：将煤视为包含裂隙和基质的双重孔隙介质，煤体的变形为弹性小变形，且服从广义胡克定律，瓦斯为理想气体且以游离态和吸附态同时赋存于裂隙和基质中，钻孔施工后，裂隙中的大部分瓦斯和基质中的少部分瓦斯均以渗流的形式进入钻孔，其余瓦斯以扩散的形式在裂隙和基质间流动，其数学表达式如下：式中，E为煤体弹性模量，MPa；v为煤体泊松比；ui,jj和uj,ji为煤体位移的张量形式，其中第一个下标分别表示i和j方向上的位移分量ui和uj，第二个下标分别表示ui和uj在j方向上的偏导数ui,j和uj,j，第三个下标分别表示ui,j和uj,j在j和i方向上的偏导数；αm和αf分别为基质和裂隙的比奥系数，αm＝K/Km‑K/Ks，αf＝1‑K/Km；K、Km和Ks分别为煤体、煤基质和煤骨架的体积模量，MPa；pm,i和pf,i分别为基质瓦斯压力pm和裂隙瓦斯压力pf在i方向上的偏导数，MPa；eL为朗格缪尔应变常数；PL为朗格缪尔压力常数，MPa；fi为i方向上的体积力分量，MPa；

3

φm和φf分别为基质和裂隙的孔隙率；ρga为标准状态下的瓦斯密度，kg/m；ρc为煤体密度，

3 3

kg/m；VL为朗格缪尔体积常数，m /kg；R为理想气体常数，J/(mol·K)；T为煤体温度，K；Mg为瓦斯摩尔质量，kg/mol；为哈密顿算子；μ为瓦斯动力粘度，Pa·s；km和kf分别为基质和裂

2 2

隙的渗透率，mD；Df0、Da0为初始扩散系数，m/s；Dr为残余扩散系数，m/s；λf和λa为衰减系数，‑1s ；Lm为基质间距，m；φm0和φf0分别为基质和裂隙的初始孔隙率； 为平均应力，MPa； 为初始平均应力，MPa；p0为煤层原始瓦斯压力，MPa；f1和f2为吸附系数；ξ为无量纲系数，其表示基质初始渗透率与煤层原始渗透率的比值；k0为煤层原始渗透率，mD。

4.根据权利要求1所述煤层瓦斯压力和渗透率同步反演方法，其特征在于，将钻孔瓦斯流量的数值计算定义为封装函数，在该函数中，煤层瓦斯压力和渗透率为输入变量，钻孔瓦斯流量为输出变量，输出变量定义的时间点与现场观测数据的时间点保持一致；

地质模拟几何模型为二维地质几何模型或者三维地质几何模型，二维地质几何模型的计算量明显小于三维地质几何模型，在计算结果相近的情况下选择二维地质几何模型则有效提高计算效率；钻孔瓦斯流量的计算公式如下：cal

当构建二维地质几何模型时，钻孔瓦斯流量模拟值q 利用钻孔边界达西速度的线积分计算：cal

当构建三维地质几何模型时，钻孔瓦斯流量模拟值q 利用钻孔边界达西速度的面积分计算：式中，L为气室长度，m；Ω为钻孔边界。

5.根据权利要求1所述煤层瓦斯压力和渗透率同步反演方法，其特征在于，目标函数通过钻孔瓦斯流量实测值和模拟值的平均绝对百分比误差MAPE表示：3

式中，N为钻孔瓦斯流量测量点的个数； 为第i时刻的钻孔瓦斯流量实测值，m/s。

6.根据权利要求1所述煤层瓦斯压力和渗透率同步反演方法，其特征在于，所述的优化参数约束条件为煤层初始瓦斯压力和渗透率的取值范围，其根据测试地点所在煤层或者邻近煤层的历史观测数据确定。

7.根据权利要求1所述煤层瓦斯压力和渗透率同步反演方法，其特征在于，优化停止条件为设定的最大函数计算数，当函数计算数满足停止条件时，优化停止，输出最优解，否则，改变参数值继续优化。

8.根据权利要求1所述煤层瓦斯压力和渗透率同步反演方法，其特征在于，煤层瓦斯含量m的计算公式如下：