《Gas Science and Engineering》：Spatio-temporal modeling for assessing geoenergy resources: a workflow applied to gas in place variation in coal beds
编辑推荐：
研究人员提出了一种时空地理统计框架，用于估算煤层气储层中随时间变化的原位气量（gas-in-place, GIP）。该方法将压力与裂隙水饱和度（water saturation, Swf）的时空克里金（spatio-temporal kriging）与煤体属性
研究人员提出了一种时空地理统计框架，用于估算煤层气储层中随时间变化的原位气量（gas-in-place, GIP）。该方法将压力与裂隙水饱和度（water saturation, Swf）的时空克里金（spatio-temporal kriging）与煤体属性（厚度、密度、Langmuir体积VL与压力PL）的空间模型相融合，基于1996–2018年11个时点的92口排采井数据（1012组观测），量化了数十年间GIP演化，包括无直接测量的年份。GIP按年份计算：GIPyear=Ah(43560?f(1?Swf)/Bg(year)+1.359GCyearρ(1?fa?fm))，其中吸附气量由Langmuir等温式GCyear=VLPyear/(PL+Pyear)给出。研究人员采用R语言gstat与spacetime包拟合时空变异函数，压力与饱和度均优选sum-metric模型（MSE最小），各向异性因子κ分别为145 ft/day与117 ft/day。生成的时空GIP估计值与累计产气量及矿山甲烷排放报告进行了一致性验证。研究表明，该框架可作为完整油藏数值模拟的可行替代方案，适用于天然气管网与地质能源资源评价，在采样稀疏、计算受限场景下具有实用价值。
研究背景与立项依据
常规与非常规储层从发现到最终采收期间，压力、流体饱和度及孔隙度渗透率等属性随采出持续演化，控制着剩余气量（gas-in-place, GIP）与流动能力。实际油藏因不可达性与测量稀疏（仅井位有数据），工程师只能依赖地质与流动模型外推未采样时空。物理基数值油藏模拟（numerical reservoir simulator）虽经地理统计实现（geostatistical realization）辅助历史拟合（history matching）后预测能力增强，但全尺度场模拟在数据需求、算力与时间上代价高昂。已有研究用多点或两点地理统计生成时序GIP分布，但未在统一框架内显式建模时空联合协方差。水文等领域广泛使用的时空地理统计（spatio-temporal geostatistics）以变异函数/协方差结构同时刻画空间与时间依赖性，可在稀疏不规则采样下插值外推，而在煤储层排采后剩余GIP评价中应用有限。为此，研究人员在《Gas Science and Engineering》发表本文，以美国阿拉巴马州Black Warrior盆地Brookwood–Oak Grove煤层气田Mary Lee煤组（ML/BC连同New Castle与Jagger煤）为案例，构建统一时空地理统计工作流估算压力、水饱和度直至GIP，并与产气量和矿排甲烷校验。
主要关键技术方法
研究人员以1996、1998、2001、2003、2006、2008、2011、2012、2013、2015、2018共11期92口排采井历史拟合所得储层压力与裂隙水饱和度（1012观测）为时空变量，煤厚、密度、Langmuir参数（VL、PL）、裂隙孔隙度?f作空间不变或井点属性；在R中以gstat/spacetime建立东向–北向–年份三维时空域，比较separable、product-sum、metric、sum-metric、simple sum-metric五类时空变异函数，按MSE最小选sum-metric；以ordinary spatio-temporal kriging在336743网格节点（裁剪至Fig.1多边形内2997空间节点、23时间步）估算压力与Swf，代入GIP方程得各煤分层逐年GIP；最后将时空GIP体积变化与区块累计产气及长壁开采期（2010–2015）矿端上报甲烷排放对照验证。
研究结果
- 2.
Introduction
研究人员指出储层动态属性时空耦合演化对剩余气评价至关重要，物理模拟成本高，统一时空地理统计框架可弥补两点/多点统计未显式建模联合协方差的缺口，并以水文、流行病学等跨域应用佐证其稀疏数据适应性。
2.1. Description of the study area
研究区位于Black Warrior盆地Pennsylvanian Pottsville组Mary Lee煤组，ML/BC合采煤层厚4–11 ft（均6.6 ft）、深1200–1950 ft，VL≈791 scf/ton、PL≈455 psia；上覆New Castle（≈65 ft上）与下伏Jagger（≈41 ft下）煤分别贡献采动顶板/底板气，故排采井全煤组射孔，构成时空数据空间载体。
2.2. Equations for spatio-temporal estimation of GIP
给出GIPyear体积方程与Langmuir吸附气含量GCyear式，明确?f生产期假定恒定（因贫时间向孔隙度数据，且自由气占比小）、VL/PL取煤级均值（Table 1），GIP由压力Pyear与Swf两时空变量主导。
2.3. Data source and analysis
2.3.1. Spatio-temporal data analysis and modelling for pressure and water saturation
经验时空均值显示压力具西高东低（东向ρ=?0.561）、饱和度无强坐标趋势；时间均值表2显示1996→2018年压力167.3→106.0 psia、饱和度61.7%→50.2%单调降，压力空间标准差同步收窄。实验时空变异函数空间滞距1500 ft、时间滞距0.01–10（≈8–7240天），拟合结果表明sum-metric模型对压力（MSE=107869.7）与饱和度（MSE=902.9）均最优；Table 4参数中空间球状项范围约10109 ft（压力）/245 ft（饱和度），时间项范围10385 ft当量/10075 ft当量，联合项存在，κ=145/117 ft/day。ordinary kriging在裁剪网格输出Pyear、Swf及标准误，再逐节点算GCyear与GIPyear。研究人员据此获得ML/BC、New Castle、Jagger各煤1996–2018逐年GIP场，无实测年份由时间向插值填补；区块总量降速与92井累计产气一致，2010–2015长壁面内剩余GIP削减量与矿端通风甲烷量级吻合，证明时空地理统计估计在物理模拟缺席下仍可约束资源量与排放源。
讨论与结论翻译
研究人员讨论认为，所提统一时空地理统计框架并非替代物理基多井油藏模拟，而是在完整模拟不可行、算力受限或数据稀疏时提供实用互补：它把压力与裂隙水饱和度作为时空随机场，以sum-metric变异函数保留空间–时间非可分离依赖，将Langmuir吸附项与厚度/密度空间场耦合，输出带克里金标准误的GIP演化面。案例证实该方法在Mary Lee煤组可重现数十年GIP衰减、补齐无测年、并与产气/排放自洽。结论为：时空克里金（spatio-temporal kriging）融合煤属性空间模型，构成一套可迁移至其他气藏、井下甲烷治理及CO2封存监测的GIP评价工作流，当大模型不可行时，能以较低数据–计算成本支撑地质能源资源决策。
