《煤层气数值模拟技术进展03for本科讲座详解说课讲解.ppt》由会员分享，可在线阅读，更多相关《煤层气数值模拟技术进展03for本科讲座详解说课讲解.ppt（37页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、煤层气数值模拟技术进展200703for本科讲座详解何谓数模何谓数模何谓数模何谓数模 煤层气数值模拟技术，是一项利用现代数煤层气数值模拟技术，是一项利用现代数煤层气数值模拟技术，是一项利用现代数煤层气数值模拟技术，是一项利用现代数值方法，采用系列偏微分方程组来描述煤层气值方法，采用系列偏微分方程组来描述煤层气值方法，采用系列偏微分方程组来描述煤层气值方法，采用系列偏微分方程组来描述煤层气及孔隙水在煤储层中的渗流过程，再通过离散及孔隙水在煤储层中的渗流过程，再通过离散及孔隙水在煤储层中的渗流过程，再通过离散及孔隙水在煤储层中的渗流过程，再通过离散化方法把连续函数转变成离散函数，进一步求化方法把连
2、续函数转变成离散函数，进一步求化方法把连续函数转变成离散函数，进一步求化方法把连续函数转变成离散函数，进一步求解偏微分方程组，从而模拟煤层气的产出过程解偏微分方程组，从而模拟煤层气的产出过程解偏微分方程组，从而模拟煤层气的产出过程解偏微分方程组，从而模拟煤层气的产出过程及产出数量。及产出数量。及产出数量。及产出数量。12/14/20222中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心一种有效工具的用途概述一种有效工具的用途概述 储层模拟提供了一种说明煤层气解吸、扩储层模拟提供了一种说明煤层气解吸、扩散和渗
3、流复杂机理的统一可靠的方法。储层散和渗流复杂机理的统一可靠的方法。储层模拟器也提供了将现场数据和实验室数据综模拟器也提供了将现场数据和实验室数据综合成一种简单的地质模型和储层模型的手段，合成一种简单的地质模型和储层模型的手段，以便评价勘探草案和长期开发策略。以便评价勘探草案和长期开发策略。数摸的用途数摸的用途数摸的用途数摸的用途优优 点点1)1)可以重复进行，能进行所谓的可以重复进行，能进行所谓的“多次开发多次开发”2)2)可以模拟各种非均质情况及复杂流体流动可以模拟各种非均质情况及复杂流体流动3)3)可以在短时间内进行反复试验，成本较低可以在短时间内进行反复试验，成本较低缺缺 点点1)1)模
4、拟精度依赖于对储层描述的精度和生产动态模拟精度依赖于对储层描述的精度和生产动态2)2)模型本身有一定的假设条件，有一定的误差模型本身有一定的假设条件，有一定的误差数模的优缺点数模的优缺点数模的优缺点数模的优缺点12/14/20224中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心数模的实现过程数模的实现过程数模的实现过程数模的实现过程建立地质模型建立地质模型建立数值模型建立数值模型建立计算机模型（软件）建立计算机模型（软件）建立数学模型建立数学模型模拟计算模拟计算12/14/20225中国石油大学（北京）煤
5、层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心 一维、二维、三维一维、二维、三维 单相、两相、单相、两相、三相三相 单组分、两组分、单组分、两组分、N组分组分 双重介质、三重介质双重介质、三重介质 直井、直井、水平井、水平井、ECBM按空间维数按空间维数按流体相数按流体相数按流体组分按流体组分按岩石类型按岩石类型地质模型地质模型地质模型地质模型按模型功能按模型功能12/14/20226中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心用于储层
6、模拟的典型网格几何形状用于储层模拟的典型网格几何形状 Tank 罐1D 一维网格 1D Radial 一维径向网格 2D Cross-Sectional 二维横截面网格 2D Areal 二维平面网格 2D Radial Cross-Sectional 二维径向横截面网格 3D 三维网格 网格几何形状网格几何形状 1 1维，维，2 2维维 或或 3 3维维利用笛卡尔利用笛卡尔(x-y-z)坐标坐标网格坐标系网格坐标系 网格坐标系网格坐标系 笛卡尔坐标、极坐标、不规则坐标、笛卡尔坐标、极坐标、不规则坐标、voronoivoronoi坐标等坐标等径向网格坐标系径向网格坐标系(r-z)(r-z)坐标
7、坐标利用利用r-r-q-q-z z 坐标坐标网格坐标系平面投影网格坐标系平面投影 网格坐标系网格坐标系 笛卡尔坐标、极坐标、不规则坐标、笛卡尔坐标、极坐标、不规则坐标、voronoivoronoi坐标等坐标等储层孔渗模型储层孔渗模型储层孔渗模型储层孔渗模型8.8.2 2 地质模型与数学模型地质模型与数学模型地质模型与数学模型地质模型与数学模型12/14/202210中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心8.8.2 2 地质模型与数学模型地质模型与数学模型地质模型与数学模型地质模型与数学模型储层孔渗
8、模型储层孔渗模型储层孔渗模型储层孔渗模型12/14/202211中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心8.8.2 2 地质模型与数学模型地质模型与数学模型地质模型与数学模型地质模型与数学模型煤层气产出模型煤层气产出模型煤层气产出模型煤层气产出模型12/14/202212中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心建立一套描述储层中流体渗流建立一套描述储层中流体渗流的偏微分方程组及其定解条件的偏微分方程组及其
9、定解条件 (初始条件、边界条件初始条件、边界条件)。守恒关系式守恒关系式运动方程运动方程状态方程状态方程辅助方程辅助方程物质平衡关系物质平衡关系能量平衡关系能量平衡关系解吸解吸Langmuir方程方程扩散扩散Fick定律定律渗流渗流Darcy定律定律流体状态方程流体状态方程岩石状态方程岩石状态方程流动辅助方程流动辅助方程参数辅助方程参数辅助方程化学辅助方程化学辅助方程物理辅助方程物理辅助方程质量守恒方程质量守恒方程(组组)能量守恒方程能量守恒方程偏微分方程偏微分方程(组组)数学模型数学模型数学模型数学模型12/14/202213中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中
10、心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心数学模型数学模型数学模型数学模型解吸模型解吸模型Langmuir方程方程式中：式中：C(p)C(p)吸附量，吸附量，ftft3 3/t/t；V VL L 兰兰氏体氏体积积，ftft3 3/t/t；P P 地地层压层压力（力（psipsi）；）；P PL L 兰兰氏氏压压力（力（psipsi）。12/14/202214中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心式中：式中：q qm m 为为煤基煤基质质中甲中甲烷扩烷扩散量，散量，m3/
11、daym3/day；D D 为扩为扩散系数，散系数，m2/daym2/day；为为形状因子，形状因子，m-2m-2；g g 为为甲甲烷烷的密度，的密度，t/m3t/m3；Vm Vm 为为煤基煤基质块质块的体的体积积，m3;m3;C(t)C(t)为为煤基煤基质质中甲中甲烷烷的平均的平均浓浓度，度，m3/t;m3/t;C(C(P P)为为基基质质-割理割理边边界上的平衡甲界上的平衡甲烷浓烷浓度，度，m m3 3/t/t。q q数学模型数学模型数学模型数学模型扩散模型扩散模型Fick定律定律12/14/202215中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京
12、）煤层气研究中心中国石油大学（北京）煤层气研究中心式中：式中：V Vl l为为l l相的渗流速度，相的渗流速度，m/sm/s；l l 为为l l相的粘滞系数，相的粘滞系数，MpaMpas s；P Pl l为为l l相的相的压压差，差，MPaMPa；L L 为为渗流途径的渗流途径的长长度，度，m m；K Kl l为为l l相的有效渗透率，相的有效渗透率，10-3m2；K K 为为多孔介多孔介质质的的绝对绝对渗透率，渗透率，10-3m2；K Kr rl l为为l l相的相相的相对对渗透率，渗透率，10-3m2。Kl=K Krl数学模型数学模型数学模型数学模型渗流模型渗流模型Darcy定律定律12/
13、14/202216中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心Tau()=1/(D*)式中：式中：吸附时间（天）吸附时间（天）基质单元形状因子基质单元形状因子 D D 扩散系数扩散系数吸附时间（吸附时间（）的确定）的确定“63的甲烷分子从微孔单元中央的甲烷分子从微孔单元中央运动到割理中所需的时间运动到割理中所需的时间”煤层气模拟煤层气模拟所需数据所需数据 储层描述数据 割理绝对渗透率割理渗透率方向 垂直渗透率 孔隙度 初始气含量 等温吸附曲线 解吸压力 吸附时间 扩散系数 割理间距 孔隙体积压缩性 基质
14、收缩压缩性 储层几何形状 构造高程（倾向）深度 净厚度 层理（层）灰分含量 井的排气面积 初始储层压力初始水饱和度 气水相对渗透率 气水毛细压力 含水层岩石性质流体PVT数据 气体地层体积系数 气体粘度 气体比重 气体组成 水地层体积系数 水粘度 水储罐密度 气体在水中的溶解度再现数据 最小时间步长 最大时间步长 时间步长增量 水产量与时间 气产量（注入量）与时间 井底（井口）压力与时间 井产能指数 表皮因子 随时间步长变化最大饱和度 随时间步长变化最大压力 有限差分求解公差 最大容许的水产量 最大容许的气产量 最小容许的井底压力 井筒半径 压裂裂隙长度12/14/202218中国石油大学（北
15、京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心数 据 项主 要 来 源 渗透率 试井 初始压力 试井 初始水饱和度 试井 气体解吸压力 试井 孔隙压缩性 试井 解吸等温线 岩心测试 吸附气含量 岩心测试 解吸时间 岩心测试 相对渗透率 生产数据和岩心测试 孔隙度 岩心测试 净产层厚度 测井和岩心测试 温度 测井 气体PVT特性 气体分析 水PVT特性 水分析 完井效果 试井 井抽排面积（间距）地质描述数数 据据来来 源源 数模技术的发展数模技术的发展 储层模型已由方糖模型发展到全三维模型（储层模型已由方糖模型发展到全三维
16、模型（储层模型已由方糖模型发展到全三维模型（储层模型已由方糖模型发展到全三维模型（Fully 3DFully 3DFully 3DFully 3D）储层孔隙模型由双重孔隙模型（裂隙系统和吸附气体）发展储层孔隙模型由双重孔隙模型（裂隙系统和吸附气体）发展储层孔隙模型由双重孔隙模型（裂隙系统和吸附气体）发展储层孔隙模型由双重孔隙模型（裂隙系统和吸附气体）发展为三重孔隙模型（基质孔隙与割理孔隙及吸附气）为三重孔隙模型（基质孔隙与割理孔隙及吸附气）为三重孔隙模型（基质孔隙与割理孔隙及吸附气）为三重孔隙模型（基质孔隙与割理孔隙及吸附气）为进行为进行为进行为进行ECBMECBMECBMECBM评价，将三重
17、孔隙度模型转换成双孔隙度模型评价，将三重孔隙度模型转换成双孔隙度模型评价，将三重孔隙度模型转换成双孔隙度模型评价，将三重孔隙度模型转换成双孔隙度模型储层孔隙模型也由一成不变的孔隙模型加入了基质收缩与孔储层孔隙模型也由一成不变的孔隙模型加入了基质收缩与孔储层孔隙模型也由一成不变的孔隙模型加入了基质收缩与孔储层孔隙模型也由一成不变的孔隙模型加入了基质收缩与孔隙膨胀模型（隙膨胀模型（隙膨胀模型（隙膨胀模型（matrix swellingmatrix swellingmatrix swellingmatrix swelling），目前已发展到所谓的微分），目前已发展到所谓的微分），目前已发展到所谓的微
18、分），目前已发展到所谓的微分膨胀模型（膨胀模型（膨胀模型（膨胀模型（differential swellingdifferential swellingdifferential swellingdifferential swelling）。）。）。）。12/14/202220中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心由于孔隙压缩、收缩和膨胀，渗由于孔隙压缩、收缩和膨胀，渗透率受到孔隙度变化的强烈影响透率受到孔隙度变化的强烈影响k=ki(/i)n 式中：式中：n n渗透率指数，通常为渗透率指数，通常为3
19、 3。压缩和基质收缩对煤的压缩和基质收缩对煤的渗透率的影响渗透率的影响 式中：式中：C Cp p孔隙压缩系数孔隙压缩系数 C Cm m基质收缩压缩系数。基质收缩压缩系数。=i i cp(Pi P)+cm(1-i)dPi(Ci-C)dCi数模技术的发展数模技术的发展 煤层气储层的渗透率模型也由单一渗透率模型（裂隙渗透率）煤层气储层的渗透率模型也由单一渗透率模型（裂隙渗透率）煤层气储层的渗透率模型也由单一渗透率模型（裂隙渗透率）煤层气储层的渗透率模型也由单一渗透率模型（裂隙渗透率）发展成双重渗透率（裂隙渗透率和基质孔隙渗透率）；渗透发展成双重渗透率（裂隙渗透率和基质孔隙渗透率）；渗透发展成双重渗透
20、率（裂隙渗透率和基质孔隙渗透率）；渗透发展成双重渗透率（裂隙渗透率和基质孔隙渗透率）；渗透率模型还加进了应力敏感模型。率模型还加进了应力敏感模型。率模型还加进了应力敏感模型。率模型还加进了应力敏感模型。煤层气解吸模型也已由单组分（煤层气解吸模型也已由单组分（煤层气解吸模型也已由单组分（煤层气解吸模型也已由单组分（CH4CH4CH4CH4）的）的）的）的LangmuirLangmuirLangmuirLangmuir方程发展成方程发展成方程发展成方程发展成多组分（多组分（多组分（多组分（CH4CH4CH4CH4、CO2CO2CO2CO2、N2N2N2N2）扩展的）扩展的）扩展的）扩展的Langm
21、uirLangmuirLangmuirLangmuir方程。方程。方程。方程。为满足为满足为满足为满足ECBMECBMECBMECBM技术研发的需要，技术研发的需要，技术研发的需要，技术研发的需要，COMET 3COMET 3COMET 3COMET 3（研发者（研发者（研发者（研发者Advanced Advanced Advanced Advanced Resources InternationalResources InternationalResources InternationalResources International）、）、）、）、GEMGEMGEMGEM（研发者（研发者（
22、研发者（研发者Computer Computer Computer Computer Modelling Group Ltd.Modelling Group Ltd.Modelling Group Ltd.Modelling Group Ltd.）、）、）、）、ECLIPSEECLIPSEECLIPSEECLIPSE（研发者（研发者（研发者（研发者SchlumbergerSchlumbergerSchlumbergerSchlumberger）、）、）、）、SIMED IISIMED IISIMED IISIMED II（研发者（研发者（研发者（研发者CSIROCSIROCSIROCSIRO）
23、等煤层气数值模拟软件陆续加入了）等煤层气数值模拟软件陆续加入了）等煤层气数值模拟软件陆续加入了）等煤层气数值模拟软件陆续加入了ECBMECBMECBMECBM模拟功能。模拟功能。模拟功能。模拟功能。12/14/202223中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心Enhanced Coalbed Methane(ECBM)RecoveryGreen House Gas(GHG)SequestrationECBM MechanismsCO2Deep CoalbedCH4CH4CH4CH4 to Sale
24、sN2 CoalFlue GasCO2N2InjectionGreen Power PlantSeparation12/14/202224中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心模拟网格精确化模拟网格精确化模拟网格精确化模拟网格精确化应用软件一体化应用软件一体化应用软件一体化应用软件一体化 前后处理可视化前后处理可视化前后处理可视化前后处理可视化 数值计算并行化数值计算并行化数值计算并行化数值计算并行化 软件技术网络化软件技术网络化软件技术网络化软件技术网络化模拟技术工程化模拟技术工程化模拟技术工程
25、化模拟技术工程化数模技术的发展数模技术的发展 12/14/202225中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心数模技术的发展数模技术的发展 12/14/202226中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心计算机模型计算机模型计算机模型计算机模型COMET 3（研发者（研发者Advanced Resources International）GEM（研发者（研发者Computer Modelling Grou
26、p Ltd.）ECLIPSE（研发者（研发者Schlumberger）SIMED II（研发者（研发者CSIRO）等煤层气数值模拟软件等煤层气数值模拟软件 12/14/202227中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心 通过离散化，将连通过离散化，将连续的偏微分方程组转换续的偏微分方程组转换成离散的有限差分方程成离散的有限差分方程组，再用多种方法将非组，再用多种方法将非线性系数线性化，成为线性系数线性化，成为线性代数方程组，然后线性代数方程组，然后求解线性代数方程组。求解线性代数方程组。偏微分方程
27、组偏微分方程组线性代数方程组线性代数方程组得到压力、饱和度等得到压力、饱和度等有限差分方程组有限差分方程组离散化离散化线性化线性化解方程组解方程组 求解技术求解技术 8.38.3方程求解技术方程求解技术方程求解技术方程求解技术12/14/202228中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心 离散化的概念离散化的概念 对储层数值模拟来说，它的数学模型是一组对储层数值模拟来说，它的数学模型是一组偏微分方程，其自变量是空间和时间。偏微分方程，其自变量是空间和时间。离散空间离散空间即把储层这个连续空间变量即
28、把储层这个连续空间变量 离散离散成若干个小单元。成若干个小单元。离散时间离散时间即把在所研究的时间范围内离散即把在所研究的时间范围内离散成一定数量的时间段。成一定数量的时间段。8.38.3方程求解技术方程求解技术方程求解技术方程求解技术12/14/202229中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心有限差分方程组的线性化方法有限差分方程组的线性化方法 IMPES方法方法(Implicit Pressure Explicit Saturation)半隐式方法半隐式方法(Semi-implicit me
29、thod)全隐式方法全隐式方法(Fully Implicit method)SEQ方法方法(Sequencial method)自适应隐式方法自适应隐式方法(Adaptive Implicit method)8.38.3方程求解技术方程求解技术方程求解技术方程求解技术12/14/202230中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心线性方程组的求解方法线性方程组的求解方法 直接解法直接解法 高斯消去法、高斯消去法、LU分解法分解法 迭代解法迭代解法 线松弛法线松弛法(LSOR)、面松弛法、面松弛法(P
30、SOR)、预处理共轭梯度法、预处理共轭梯度法 直接解法占用内存多，但计算速度快；迭代解法占直接解法占用内存多，但计算速度快；迭代解法占用内存少，但由于迭代次数多，而降低计算速度。用内存少，但由于迭代次数多，而降低计算速度。预处理共轭梯度法预处理共轭梯度法 在在80年代兴起，该方法适用于解大型稀疏矩阵。预处理是将年代兴起，该方法适用于解大型稀疏矩阵。预处理是将稀疏矩阵不完全稀疏矩阵不完全LU分解成近似阵，然后用正交极小化使迭代过程分解成近似阵，然后用正交极小化使迭代过程沿着最快的方向收敛。沿着最快的方向收敛。8.38.3方程求解技术方程求解技术方程求解技术方程求解技术12/14/202231中国
31、石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心双重孔隙模型与三重孔隙模型对比双重孔隙模型与三重孔隙模型对比 气产量历史拟合气产量历史拟合三重孔隙度气产量双孔隙度气产量三重孔隙度历史气产量数据 12/14/202232中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心双重孔隙模型与三重孔隙模型对比双重孔隙模型与三重孔隙模型对比 水产量历史拟合水产量历史拟合三重孔隙度水产量双孔隙度水产量三重孔隙度历史水产量数据 12/14/20
32、2233中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心双孔隙度模型与三重孔隙度模型对比双孔隙度模型与三重孔隙度模型对比 预测气产量预测气产量三重孔隙度气产量双孔隙度气产量三重孔隙度历史气产量数据 双孔隙度模型与三重孔隙度模型对比双孔隙度模型与三重孔隙度模型对比 预测水产量预测水产量三重孔隙度水产量双孔隙度水产量三重孔隙度历史水产量数据 张遂安张遂安张遂安张遂安张遂安张遂安 教授教授教授教授教授教授12/14/202236中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心中国石油大学（北京）煤层气研究中心此课件下载可自行编辑修改，仅供参考！此课件下载可自行编辑修改，仅供参考！感谢您的支持，我们努力做得更好！谢谢感谢您的支持，我们努力做得更好！谢谢