COMSOL能源开采仿真：瓦斯抽采与注气的深度探索

COMSOL能源开采仿真基质中瓦斯扩散、裂隙中瓦斯渗流分析不同工况条件下渗透率演化、有效抽采半径、抽采产量。 使用模块PDE基质瓦斯扩散达西定律/PDE裂隙瓦斯渗流固体力学/PDE煤体变形控制方程。 工况条件双重孔隙模型、仅考虑裂隙渗流模型。 采用时变扩散模型、双孔扩散模型分析煤粒瓦斯解吸扩散特性。 耦合不同扩散模型与裂隙中瓦斯渗流分析瓦斯渗流特性。 煤粒中瓦斯解吸扩散采用PDE模块采用吸附与解吸模型分析不同扩散模型下的吸附平衡时间。 瓦斯抽采过程巷道工作面中空气流入煤层中钻孔裂隙圈周围空气流入钻孔中降低瓦斯抽采浓度。 物理场方程选用达西定律与固体力学方程瓦斯的压力梯度影响空气流动混合气体抽采。 注气井中注入CO2/N2与甲烷发生竞争吸附作用促进甲烷解吸起到增产增流的目的。 注入井中可注入CO2/N2或者同时注入CO2、N2分析不同注入压力、注入速率下甲烷的抽采产量、流动速率。在能源开采领域瓦斯的有效抽采一直是关键课题。今天咱们就来唠唠利用COMSOL进行瓦斯抽采相关的仿真这里面涉及到基质中瓦斯扩散、裂隙中瓦斯渗流还有不同工况条件下各种参数的分析可有意思了。一、涉及的模块与模型咱们使用的模块有PDE基质瓦斯扩散达西定律/PDE裂隙瓦斯渗流固体力学/PDE煤体变形控制方程。就像搭建乐高一样每个模块负责不同的功能最后拼出整个瓦斯抽采仿真的“大厦”。一扩散模型时变扩散模型与双孔扩散模型采用时变扩散模型、双孔扩散模型能帮助我们分析煤粒瓦斯解吸扩散特性。比如说在COMSOL中设置时变扩散模型代码大概像这样model createpde(); geometryFromEdges(model, gdm); % gdm是自定义的几何函数 specifyCoefficients(model, m, 0, d, 1, c, 1, a, 0, f, 0); applyBoundaryCondition(model, Dirichlet, Edge, 1:model.Geometry.NumEdges, u, 0); generateMesh(model, Hmax, 0.1); tlist 0:0.01:1; u0 zeros(model.NumNodes, 1); result solvepde(model, tlist, u0);这段代码创建了一个PDE模型设定了系数施加了边界条件划分了网格最后求解时变问题。在瓦斯扩散仿真里这里的u就可以代表瓦斯浓度随时间的变化通过这个模型咱们能直观看到瓦斯在基质中如何随着时间扩散。COMSOL能源开采仿真基质中瓦斯扩散、裂隙中瓦斯渗流分析不同工况条件下渗透率演化、有效抽采半径、抽采产量。 使用模块PDE基质瓦斯扩散达西定律/PDE裂隙瓦斯渗流固体力学/PDE煤体变形控制方程。 工况条件双重孔隙模型、仅考虑裂隙渗流模型。 采用时变扩散模型、双孔扩散模型分析煤粒瓦斯解吸扩散特性。 耦合不同扩散模型与裂隙中瓦斯渗流分析瓦斯渗流特性。 煤粒中瓦斯解吸扩散采用PDE模块采用吸附与解吸模型分析不同扩散模型下的吸附平衡时间。 瓦斯抽采过程巷道工作面中空气流入煤层中钻孔裂隙圈周围空气流入钻孔中降低瓦斯抽采浓度。 物理场方程选用达西定律与固体力学方程瓦斯的压力梯度影响空气流动混合气体抽采。 注气井中注入CO2/N2与甲烷发生竞争吸附作用促进甲烷解吸起到增产增流的目的。 注入井中可注入CO2/N2或者同时注入CO2、N2分析不同注入压力、注入速率下甲烷的抽采产量、流动速率。双孔扩散模型稍微复杂点它考虑了煤体中不同尺度孔隙的扩散差异。在代码实现上可能需要额外定义不同孔隙尺度的参数像孔隙率、扩散系数等然后分别设置它们在模型中的作用。吸附与解吸模型煤粒中瓦斯解吸扩散采用PDE模块并且采用吸附与解吸模型来分析不同扩散模型下的吸附平衡时间。这时候在PDE系数设置里就需要考虑吸附和解吸相关的参数比如吸附速率、解吸速率等。代码示例model createpde(); geometryFromEdges(model, gdm); % 这里假设吸附速率为ka解吸速率为kd specifyCoefficients(model, m, 0, d, 1, c, 1, a, -ka kd, f, 0); applyBoundaryCondition(model, Dirichlet, Edge, 1:model.Geometry.NumEdges, u, 0); generateMesh(model, Hmax, 0.1); tlist 0:0.01:1; u0 zeros(model.NumNodes, 1); result solvepde(model, tlist, u0);这里的a系数就结合了吸附和解吸的影响通过调整ka和kd的值咱们就能观察不同吸附解吸条件下瓦斯浓度随时间的变化进而分析吸附平衡时间。二、工况条件分析一双重孔隙模型与仅考虑裂隙渗流模型双重孔隙模型双重孔隙模型考虑了基质孔隙和裂隙孔隙的双重作用。在COMSOL仿真里就需要分别对基质和裂隙设置不同的参数比如基质的渗透率低裂隙的渗透率高。代码里可能会有类似这样的设置% 假设k1为基质渗透率k2为裂隙渗透率 k1 1e-15; k2 1e-12;通过这样的设置就能模拟瓦斯在双重孔隙结构中的流动分析渗透率演化等特性。仅考虑裂隙渗流模型这个模型相对简单只关注裂隙中的瓦斯渗流。代码实现上可能就不需要过多考虑基质相关的参数设置主要精力放在裂隙的参数像裂隙宽度、裂隙方向等对瓦斯渗流的影响。三、瓦斯渗流特性分析耦合不同扩散模型与裂隙中瓦斯渗流这可是个关键步骤。在实际瓦斯抽采过程中瓦斯在基质中扩散出来后要通过裂隙流向抽采钻孔。咱们来看看代码实现思路% 假设已经分别设置好基质扩散模型和裂隙渗流模型 % 耦合部分 modelCoupled createMultiphysics(); modelDiffusion createpde(); % 基质扩散模型 modelSeepage createpde(); % 裂隙渗流模型 addPhysics(modelCoupled, modelDiffusion, Diffusion); addPhysics(modelCoupled, modelSeepage, DarcyFlow);这里创建了多物理场模型modelCoupled把基质扩散和裂隙渗流两个物理场添加进去实现耦合。通过这样的耦合模型咱们可以分析瓦斯渗流特性比如有效抽采半径、抽采产量等。四、瓦斯抽采中的空气流入问题瓦斯抽采过程巷道工作面中空气流入煤层中钻孔裂隙圈周围空气流入钻孔中这会降低瓦斯抽采浓度。物理场方程选用达西定律与固体力学方程瓦斯的压力梯度影响空气流动混合气体抽采。在COMSOL里就需要设置多相流相关的参数代码示例model createpde(); addPhysics(model, DarcyFlow); addPhysics(model, SolidMechanics); % 设置多相流参数比如空气和瓦斯的密度、粘度等 rhoAir 1.2; rhoGas 0.7; muAir 1.8e-5; muGas 1.1e-5; setCoefficients(model, darcy, {rho, rhoAir rhoGas,mu, muAir muGas});这样设置后就能模拟空气和瓦斯混合流动的情况分析空气流入对瓦斯抽采浓度的影响。五、注气增产分析注气井中注入CO2/N2与甲烷发生竞争吸附作用促进甲烷解吸起到增产增流的目的。注入井中可注入CO2/N2或者同时注入CO2、N2分析不同注入压力、注入速率下甲烷的抽采产量、流动速率。在COMSOL里实现这个过程代码大概如下model createpde(); addPhysics(model, ReactionEngineering); % 设置CO2、N2和甲烷的反应参数比如吸附反应速率 k1 0.1; % CO2吸附反应速率 k2 0.05; % N2吸附反应速率 k3 0.2; % 甲烷解吸反应速率 reactionRates [-k1 * cCO2 k3 * cCH4, -k2 * cN2 k3 * cCH4, k1 * cCO2 k2 * cN2 - 2 * k3 * cCH4]; setReaction(model, Rates, reactionRates);这里通过ReactionEngineering模块设置了注气和甲烷解吸的反应参数通过调整这些参数就能分析不同注入情况下甲烷的抽采产量和流动速率。总的来说利用COMSOL进行能源开采瓦斯抽采仿真通过不同模块、模型和工况条件的设置结合代码实现能帮助我们深入了解瓦斯抽采过程中的各种特性为实际开采提供有力的理论支持。