CN106251232A - 确定页岩含气量的方法和装置 - Google Patents
确定页岩含气量的方法和装置 Download PDFInfo
- Publication number
- CN106251232A CN106251232A CN201610663994.4A CN201610663994A CN106251232A CN 106251232 A CN106251232 A CN 106251232A CN 201610663994 A CN201610663994 A CN 201610663994A CN 106251232 A CN106251232 A CN 106251232A
- Authority
- CN
- China
- Prior art keywords
- unit
- shale
- gas content
- represent
- centerdot
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 30
- 230000015572 biosynthetic process Effects 0.000 claims description 20
- 239000011148 porous material Substances 0.000 claims description 18
- 238000010521 absorption reaction Methods 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 7
- 238000005859 coupling reaction Methods 0.000 claims description 7
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 238000010219 correlation analysis Methods 0.000 claims description 6
- 230000002596 correlated effect Effects 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 7
- 239000007789 gas Substances 0.000 description 110
- 238000001179 sorption measurement Methods 0.000 description 26
- 238000000611 regression analysis Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000000556 factor analysis Methods 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000004459 forage Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 239000006104 solid solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种确定页岩含气量的方法和装置,其中,该方法包括:确定当前工区页岩吸附气含量的第一主控因素;获取页岩游离气含量的第二主控因素;基于所述第一主控因素,建立所述当前工区的页岩吸附气含量预测公式;基于所述第二主控因素,建立所述当前工区的页岩游离气含量预测公式;根据所述页岩吸附气含量预测公式和所述页岩游离气含量预测公式,确定所述当前工区的含气量。本发明实施例解决了现有技术中仅定性地对页岩含气量进行预测,从而导致的预测出的含气量准确性不高的技术问题,达到了有效提高页岩含气量预测结果准确度的技术效果。
Description
技术领域
本发明涉及地质勘探技术领域,特别涉及一种确定页岩含气量的方法和装置。
背景技术
页岩气作为一种新型非常规天然气能源其勘探开发越来越受重视,页岩气的赋存形式包括:游离态、吸附态、溶解态和固溶态,且以游离态和吸附态为主。因此,预测游离气和吸附气含量对评价页岩气的资源潜力具有重要的意义。
作为与常规天然气存在重大差异的吸附态页岩气,现在预测页岩吸附气含量的模型主要有:Langmuir吸附模型、DR吸附模型、DA吸附模型、BET吸附模型等,然而,这些模型仅是定性地对页岩含气量进行预测,从而导致预测出的含气量准确性不高,难以对实际生产进行指导生产。
针对上述问题,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种确定页岩含气量的方法,以达到定性实现页岩含气量预测的技术效果,该方法包括:
确定当前工区页岩吸附气含量的第一主控因素;
确定页岩游离气含量的第二主控因素;
基于所述第一主控因素,建立所述当前工区的页岩吸附气含量预测公式;
基于所述第二主控因素,建立所述当前工区的页岩游离气含量预测公式;
根据所述页岩吸附气含量预测公式和所述页岩游离气含量预测公式,确定所述当前工区的含气量。
在一个实施方式中,确定当前工区页岩吸附气含量的第一主控因素,包括:
通过单因素相关性分析和多因素耦合的方法,确定出多个第一主控因素,其中,所述第一主控因素包括以下至少之一:页岩储层的地化参数、无机矿物参数、和孔隙结构参数。
在一个实施方式中,所述页岩吸附气含量预测公式为:
其中,Qa表示吸附气含量,单位为m3/t,P表示地层压力,单位为MPa,T表示地层温度,单位为℃,no表示兰氏体积,单位为m3/t,KT表示吸附常数,是温度T的函数,Xi表示第一主控因素,Ai表示Xi的待定系数。
在一个实施方式中,所述第二主控因素包括:孔隙度、地层压力和地层温度。
在一个实施方式中,所述页岩游离气含量预测公式为:
其中,Qf表示实际地质条件下的游离气含量,单位为cm3/g,Vp表示总孔体积,单位为ml/g,表示甲烷分子直径,单位为nm,S表示比表面积,单位为m2/g,Sw表示孔隙的含水饱和度,单位为%,So表示含油饱和度,单位为%,单位为cm3/g,T表示地层温度,单位为℃,P表示地层压力,单位为MPa,Tsc表示地面的标准温度,单位为293K,Psc表示地面的标准压力,单位为0.101MPa,Z表示原始气体压缩因子。
本发明实施例还提供了一种确定页岩含气量的装置,以达到定性实现页岩含气量预测的技术效果,该装置包括:
第一确定模块,用于确定当前工区页岩吸附气含量的第一主控因素;
获取模块,用于获取页岩游离气含量的第二主控因素;
第一建立模块,用于基于所述第一主控因素,建立所述当前工区的页岩吸附气含量预测公式;
第二建立模块,用于基于所述第二主控因素,建立所述当前工区的页岩游离气含量预测公式;
第二确定模块,用于根据所述页岩吸附气含量预测公式和所述页岩游离气含量预测公式,确定所述当前工区的含气量。
在一个实施方式中,所述第一确定模块具体用于通过单因素相关性分析和多因素耦合的方法,确定出多个第一主控因素,其中,所述第一主控因素包括以下至少之一:页岩储层的地化参数、无机矿物参数、和孔隙结构参数。
在一个实施方式中,所述页岩吸附气含量预测公式为:
其中,Qa表示吸附气含量,单位为m3/t,P表示地层压力,单位为MPa,T表示地层温度,单位为℃,no表示兰氏体积,单位为m3/t,KT表示吸附常数,是温度T的函数,Xi表示第一主控因素,Ai表示Xi的待定系数。
在一个实施方式中,所述第二主控因素包括:孔隙度、地层压力和地层温度。
在一个实施方式中,所述页岩游离气含量预测公式为:
其中,Qf表示实际地质条件下的游离气含量,单位为cm3/g,Vp表示总孔体积,单位为ml/g,表示甲烷分子直径,单位为nm,S表示比表面积,单位为m2/g,Sw表示孔隙的含水饱和度,单位为%,So表示含油饱和度,单位为%,单位为cm3/g,T表示地层温度,单位为℃,P表示地层压力,单位为MPa,Tsc表示地面的标准温度,单位为293K,Psc表示地面的标准压力,单位为0.101MPa,Z表示原始气体压缩因子。
在本发明实施例中,建立了当前工区的页岩吸附气含量预测公式和页岩游离气含量预测公式,结合页岩吸附气含量预测公式和页岩游离气含量预测公式来确定当前工区的含气量,通过上述方式解决了现有技术中仅定性地对页岩含气量进行预测,从而导致的预测出的含气量准确性不高的技术问题,达到了有效提高页岩含气量预测结果准确度的技术效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是根据本发明实施例的确定页岩含气量的方法流程图;
图2是根据本发明实施例的确定页岩含气量的装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
发明人考虑到对于页岩储层而言,不仅具有吸附气还具有游离气,因此,如果希望得到较为准确的含气量,不仅需要计算吸附气含量,还需要计算游离气含量,对此,在本例中,提供了一种确定页岩含气量的方法,如图1所示,该方法包括:
步骤101:确定当前工区页岩吸附气含量的第一主控因素;
在确定第一主控因素的过程中,可以先将实测过剩吸附量校正为绝对吸附量,具体地,可以根据选出的具有代表性的目的层岩心,进行有机地化分析、储层特征分析、等温吸附实验等。其中,等温吸附实验测得的吸附量是在吸附平衡后,根据气体状态方程计算气体物质的量的差值,该过程认为体积不变,固体表面吸附层中的气体分子中按气相密度分布于吸附相空间的气体分子与气/固分子间的作用力无关,因此,等温吸附实验得到的页岩吸附量对应于吸附相中超过气相密度的过剩量,而非页岩的真实吸附量即绝对吸附量,因此,在进行模型建立之前,可以将实测的过剩吸附量校正为绝对吸附量。具体地,可以通过以下方式进行校正:
实测过剩吸附(nexcess,单位为mmol/g)与绝对吸附(nabs,单位为mmol/g)的关系表达式为:
其中,nexcess表示实测过剩吸附量,单位为mmol/g,nabs表示绝对吸附量,单位为mmol/g,ρg表示甲烷游离相密度,单位为mol/L,可通过气体状态方程求取,ρads表示甲烷吸附相密度,单位为mol/L。
进一步的,假设甲烷在页岩中的吸附为单层吸附,则甲烷的绝对吸附量可以表示为:
其中,no表示兰氏体积,单位为m3/t,P表示地层压力,单位为MPa,T表示地层温度,单位为℃,KT表示吸附常数,是温度T的函数。
将nabs的公式代入nexcess的公式,可以得到:
在该公式中,未知参数包括:no、KT及ρads,通过将已知实测数据nexcess、T及P代入该公式进行拟合,可以获得上述三个未知参数。
将拟合得到的上述三个未知参数代入nexcess公式,就可以获得实际页岩样品的绝对吸附量,将各个样品的绝对吸附量与页岩储层的内部地质参数进行单因素相关性分析,其中,这些地质参数可以包括:页岩储层的地化参数、无机矿物参数、孔隙结构参数等。在单因素相关分析的基础上,可以建立如下的主控因素的优选原则:
1)选取比较常规而且容易获得的参数;
2)选取数据量多的参数,更能反应准确的规律;
3)彼此相关的参数选择其中与含气量关系更加明显的。
在确定页岩吸附含量的主控因素的过程中,可以通过单因素分析多因素耦合的方法确定,单因素分析就是逐个因素确定该因素对吸附气含量的影响,多因素耦合主要是确定出较为关键的因素。
进一步的,之所以需要进行单因素分析多因素耦合的方式确定主控因素,是因为不同的工区影响吸附气含量的主控因素是不同的,因此在确定主控因素和进行后续的公式确定的过程中,所采用的数据都是该工区的数据。
具体地,对于吸附气含量预测公式需要确定主控内因,即,确定主控因素,具体地,可以按照以下规则建立主控因素:
1)选取比较常规而且容易获得的参数;
2)选取数据量多的参数,更能反应准确的规律;
3)彼此相关的参数选择其中与含气量关系更加明显的参数。
即,无论对于哪个研究区,温度和压力都是必须的第一主控因素,而其它的因素,例如:页岩储层的地化参数、无机矿物参数、孔隙结构参数等是否作为第一主控因素都是需要根据实际的研究区的情况确定的。
步骤102:获取页岩游离气含量的第二主控因素;
页岩游离态天然气的赋存与常规油气赋存机理相似,同样服从气体方程,取决于孔隙结构和孔隙体积、地层温度以及气体压力和压缩系数。
尽管页岩储层具有低孔低渗的特征,在扫描电镜下依然可以观察到不同大小的孔隙、喉道、微裂缝等组成的复杂多孔系统,并具有网状连通的特征。在页岩储层中的这些孔隙是游离态页岩气主要的赋存空间,因此,孔隙度的大小直接控制着游离态页岩气的含量。一般来说,当页岩中孔隙直径较大时,其内部主要以游离态赋存,孔隙度越大,所含游离气含量就越高。例如:当孔隙度从0.5%提升到4.2%时,页岩中游离气含量也会从开始的5%上升到50%。因此,影响页岩游离气含量的其中一个主要因素是孔隙度。
进一步的,游离气的赋存也与地层的温压条件密切相关。一般情况下,在一定温度内,压力的增加会使游离气含量平稳增大,而压力往往随着埋深的增加而增大,因此页岩埋深较大时,其中游离气量也逐步增加。1150m的页岩储层中游离气含量与吸附气含量基本相当,随着埋深逐渐增大,到达2800m左右时,页岩中游离气含量明显增大,至吸附气含量的2倍以上。因此,影响页岩游离气含量的另外的主控因素是地层压力和地层温度。
由此,可以确定页岩游离气含量的第二主控因素就可以包括但不限于:孔隙度、地层温度和地层压力。
步骤103:基于所述第一主控因素,建立所述当前工区的页岩吸附气含量预测公式;
在具体实现的时候,可以通过数学地质中的多元回归分析法建立最终的页岩含气量预测公式。
其中,所谓的多元回归分析法,就是当某一个特定地质变量y与其他多个地质变量xi(i=1,2,…,m)间存在一定的依赖性,而它们之间的数量关系却是不确定的情况下,需要通过这些相关变量间的观测值,建立出一个可以表达其相关关系的数学表达式,这一统计分析方法即为多元回归分析。
在地质勘探研究过程中,多元回归分析总体上可以解决三方面难题:
1)确定特定地质变量y与其他变量Xi(i=1,2,…,p,其中,p≤m)是否存在相关关系,若存在,寻求彼此间相关关系的数学表达式;
2)根据相关变量Xi(i=1,2,…,p,其中,p≤m)的观测值,利用确定的数学表达式可以预测该特定变量的估计值,同时也可以确定出其结果的精确度;
3)通过多元回归分析也可以优选出对该特定变量y作用大的相关地质变量,剔除对其影响无足轻重的地质变量,从而简化地质研究。
按回归模型类型可将多元回归分析划分为:线性回归分析和非线性回归分析,在本例中,采用了非线性回归分析,建立了相关变量的回归方程为:
y=A1+A2*f(X1)+A3*f(X2)+…
其中,y表示因变量,f(X1)、f(X2)、…为选出的自变量参数,A1、A2和A3、…为回归系数。
具体地,结合实际的地层温度、压力等主控外因进行分析后,发现吸附气含量与确定的各第一主控因素呈线性关系,与温度及压力的关系满足Langmuir吸附模型,因此,可以建立如下的页岩吸附气含量预测公式:
其中,Qa表示吸附气含量,单位为m3/t,P表示地层压力,单位为MPa,T表示地层温度,单位为℃,no表示兰氏体积,单位为m3/t,KT表示吸附常数,是温度T的函数,Xi表示第一主控因素,Ai表示Xi的待定系数。
步骤104:基于所述第二主控因素,建立所述当前工区的页岩游离气含量预测公式;
游离气含量的预测公式的建立主要是以页岩中孔隙发育的特征研究为基础,假设吸附气含量为实际地质条件下的最大吸附量,即,整个孔隙抽象为饱和吸附状态下布满单层甲烷分子的孔隙空间,其剩余的孔隙体积即为游离气体积,其中,吸附气的体积通过甲烷分子直径与吸附面积来表征,从而建立地表条件下游离气的单因素理论预测模型:
Vp=EXP(a·h+b)
其中,Vs表示地表条件下的游离气含量,单位为cm3/g,Vp表示为总孔体积,单位为ml/g,表示甲烷分子的直径,单位为nm,S表示比表面积,单位为m2/g,Sw表示孔隙含水饱和度,单位为%,So表示含油饱和度,单位为%,h表示深度,单位为m,a、b表示系数,EXP表示指数函数,EXP{F(x)}EXP{F(X)}表示e的F(x)次方。
由于实际地质条件与地表条件温压之间是存在差异的,地表条件下的游离气含量预测模型与实际地质条件下的差异也很大,因此,需要在游离气的理论预测模型建立之后,恢复实际地层埋藏条件下的温度和压力,为此,可以建立如下相应的游离气含量的预测公式:
其中,Qf表示实际地质条件下的游离气含量,单位为cm3/g,Vs表示地表条件下的游离气含量,单位为cm3/g,T表示地层温度,单位为℃,P表示地层压力,单位为MPa,Tsc地面的标准温度,单位为293K,Psc表示地面的标准压力,单位为0.101MPa,Z表示原始气体压缩因子。
将上述Vs的公式代入该Qf的公式即可得到最终的游离气含量的预测公式:
其中,Qf表示实际地质条件下的游离气含量,单位为cm3/g,Vp表示总孔体积,单位为ml/g,表示甲烷分子直径,单位为nm,S表示比表面积,单位为m2/g,Sw表示孔隙的含水饱和度,单位为%,So表示含油饱和度,单位为%,单位为cm3/g,T表示地层温度,单位为℃,P表示地层压力,单位为MPa,Tsc表示地面的标准温度,单位为293K,Psc表示地面的标准压力,单位为0.101MPa,Z表示原始气体压缩因子。
步骤105:根据所述页岩吸附气含量预测公式和所述页岩游离气含量预测公式,确定所述当前工区的含气量。
结合上述的页岩吸附气含量预测公式和页岩游离气含量预测公式,可以得到最终确定页岩含气量的公式:
其中,Q表示页岩含气量,单位为m3/t,P表示地层压力,单位为MPa,T表示地层温度,单位为℃,no表示兰氏体积,单位为m3/t,KT表示吸附常数,是温度T的函数,Xi表示第一主控因素,Ai表示Xi的待定系数,Vp表示总孔体积,单位为ml/g,表示甲烷分子直径,单位为nm,S表示比表面积,单位为m2/g,Sw表示孔隙的含水饱和度,单位为%,So表示含油饱和度,单位为%,单位为cm3/g,T表示地层温度,单位为℃,P表示地层压力,单位为MPa,Tsc表示地面的标准温度,单位为293K,Psc表示地面的标准压力,单位为0.101MPa,Z表示原始气体压缩因子。
在上例中,不仅计算了页岩吸附气含量,也计算了页岩游离气含量,将两者结合作为最终的页岩总的含气量,从而使得得到的页岩总的含气量的结果更为准确。
基于同一发明构思,本发明实施例中还提供了一种确定页岩含气量的装置,如下面的实施例所述。由于确定页岩含气量的装置解决问题的原理与确定页岩含气量的方法相似,因此确定页岩含气量的装置的实施可以参见确定页岩含气量的方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图2是本发明实施例的确定页岩含气量的装置的一种结构框图,如图2所示可以包括:第一确定模块201、获取模块202、第一建立模块203、第二建立模块204和第二确定模块205,下面对该结构进行说明。
第一确定模块201,用于确定当前工区页岩吸附气含量的第一主控因素;
获取模块202,用于获取页岩游离气含量的第二主控因素;
第一建立模块203,用于基于所述第一主控因素,建立所述当前工区的页岩吸附气含量预测公式;
第二建立模块204,用于基于所述第二主控因素,建立所述当前工区的页岩游离气含量预测公式;
第二确定模块205,用于根据所述页岩吸附气含量预测公式和所述页岩游离气含量预测公式,确定所述当前工区的含气量。
在一个实施方式中,第一确定模块201具体可以用于通过单因素相关性分析和多因素耦合的方法,确定出多个第一主控因素,其中,所述第一主控因素包括以下至少之一:页岩储层的地化参数、无机矿物参数、和孔隙结构参数。
在一个实施方式中,页岩吸附气含量预测公式可以表示为:
其中,Qa表示吸附气含量,单位为m3/t,P表示地层压力,单位为MPa,T表示地层温度,单位为℃,no表示兰氏体积,单位为m3/t,KT表示吸附常数,是温度T的函数,Xi表示第一主控因素,Ai表示Xi的待定系数。
在一个实施方式中,所述第二主控因素可以包括但不限于以下至少之一:孔隙度、地层压力和地层温度。
在一个实施方式中,所述页岩游离气含量预测公式可以表示为:
其中,Qf表示实际地质条件下的游离气含量,单位为cm3/g,Vp表示总孔体积,单位为ml/g,表示甲烷分子直径,单位为nm,S表示比表面积,单位为m2/g,Sw表示孔隙的含水饱和度,单位为%,So表示含油饱和度,单位为%,单位为cm3/g,T表示地层温度,单位为℃,P表示地层压力,单位为MPa,Tsc表示地面的标准温度,单位为293K,Psc表示地面的标准压力,单位为0.101MPa,Z表示原始气体压缩因子。
从以上的描述中,可以看出,本发明实施例实现了如下技术效果:建立了当前工区的页岩吸附气含量预测公式和页岩游离气含量预测公式,结合页岩吸附气含量预测公式和页岩游离气含量预测公式来确定当前工区的含气量,通过上述方式解决了现有技术中仅定性地对页岩含气量进行预测,从而导致的预测出的含气量准确性不高的技术问题,达到了有效提高页岩含气量预测结果准确度的技术效果。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种确定页岩含气量的方法,其特征在于,包括:
确定当前工区页岩吸附气含量的第一主控因素;
获取页岩游离气含量的第二主控因素;
基于所述第一主控因素,建立所述当前工区的页岩吸附气含量预测公式;
基于所述第二主控因素,建立所述当前工区的页岩游离气含量预测公式;
根据所述页岩吸附气含量预测公式和所述页岩游离气含量预测公式,确定所述当前工区的含气量。
2.如权利要求1所述的方法,其特征在于,确定当前工区页岩吸附气含量的第一主控因素,包括:
通过单因素相关性分析和多因素耦合的方法,确定出多个第一主控因素,其中,所述第一主控因素包括以下至少之一:页岩储层的地化参数、无机矿物参数、和孔隙结构参数。
3.如权利要求1或2所述的方法,其特征在于,所述页岩吸附气含量预测公式为:
其中,Qa表示吸附气含量,单位为m3/t,P表示地层压力,单位为MPa,T表示地层温度,单位为℃,no表示兰氏体积,单位为m3/t,KT表示吸附常数,是温度T的函数,Xi表示第一主控因素,Ai表示Xi的待定系数。
4.如权利要求1所述的方法,其特征在于,所述第二主控因素包括:孔隙度、地层压力和地层温度。
5.如权利要求1或4所述的方法,其特征在于,所述页岩游离气含量预测公式为:
其中,Qf表示实际地质条件下的游离气含量,单位为cm3/g,Vp表示总孔体积,单位为ml/g,表示甲烷分子直径,单位为nm,S表示比表面积,单位为m2/g,Sw表示孔隙的含水饱和度,单位为%,So表示含油饱和度,单位为%,单位为cm3/g,T表示地层温度,单位为℃,P表示地层压力,单位为MPa,Tsc表示地面的标准温度,单位为293K,Psc表示地面的标准压力,单位为0.101MPa,Z表示原始气体压缩因子。
6.一种确定页岩含气量的装置,其特征在于,包括:
第一确定模块,用于确定当前工区页岩吸附气含量的第一主控因素;
获取模块,用于获取页岩游离气含量的第二主控因素;
第一建立模块,用于基于所述第一主控因素,建立所述当前工区的页岩吸附气含量预测公式;
第二建立模块,用于基于所述第二主控因素,建立所述当前工区的页岩游离气含量预测公式;
第二确定模块,用于根据所述页岩吸附气含量预测公式和所述页岩游离气含量预测公式,确定所述当前工区的含气量。
7.如权利要求6所述的装置,其特征在于,所述第一确定模块具体用于通过单因素相关性分析和多因素耦合的方法,确定出多个第一主控因素,其中,所述第一主控因素包括以下至少之一:页岩储层的地化参数、无机矿物参数、和孔隙结构参数。
8.如权利要求6或7所述的装置,其特征在于,所述页岩吸附气含量预测公式为:
其中,Qa表示吸附气含量,单位为m3/t,P表示地层压力,单位为MPa,T表示地层温度,单位为℃,no表示兰氏体积,单位为m3/t,KT表示吸附常数,是温度T的函数,Xi表示第一主控因素,Ai表示Xi的待定系数。
9.如权利要求6所述的装置,其特征在于,所述第二主控因素包括:孔隙度、地层压力和地层温度。
10.如权利要求6或9所述的装置,其特征在于,所述页岩游离气含量预测公式为:
其中,Qf表示实际地质条件下的游离气含量,单位为cm3/g,Vp表示总孔体积,单位为ml/g,表示甲烷分子直径,单位为nm,S表示比表面积,单位为m2/g,Sw表示孔隙的含水饱和度,单位为%,So表示含油饱和度,单位为%,单位为cm3/g,T表示地层温度,单位为℃,P表示地层压力,单位为MPa,Tsc表示地面的标准温度,单位为293K,Psc表示地面的标准压力,单位为0.101MPa,Z表示原始气体压缩因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610663994.4A CN106251232B (zh) | 2016-08-12 | 2016-08-12 | 确定页岩含气量的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610663994.4A CN106251232B (zh) | 2016-08-12 | 2016-08-12 | 确定页岩含气量的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106251232A true CN106251232A (zh) | 2016-12-21 |
CN106251232B CN106251232B (zh) | 2019-09-06 |
Family
ID=57591582
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610663994.4A Active CN106251232B (zh) | 2016-08-12 | 2016-08-12 | 确定页岩含气量的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106251232B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107271322A (zh) * | 2017-07-27 | 2017-10-20 | 中国石油大学(华东) | 一种定量评价泥页岩页岩气资源量及表征吸附气和游离气转化规律的方法 |
CN107818379A (zh) * | 2017-09-28 | 2018-03-20 | 中国石油天然气股份有限公司 | 吸附气量确定方法和装置 |
CN108150161A (zh) * | 2017-12-27 | 2018-06-12 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 页岩含气性评价方法装置 |
CN108593493A (zh) * | 2018-01-10 | 2018-09-28 | 中国石油天然气股份有限公司 | 一种确定致密储层含气量的方法和装置 |
CN108717202A (zh) * | 2018-05-31 | 2018-10-30 | 中国地质大学(武汉) | 一种考虑流体温度的页岩气异常地层压力预测方法 |
CN108846218A (zh) * | 2018-06-22 | 2018-11-20 | 中国石油大学(北京) | 预测泥页岩储层不同类型孔隙的内比表面积的方法及系统 |
CN110018293A (zh) * | 2019-05-05 | 2019-07-16 | 西安石油大学 | 一种考虑水敏性的页岩黏土多相耦合甲烷含气量计算方法 |
CN110579797A (zh) * | 2019-09-02 | 2019-12-17 | 长江师范学院 | 一种泥页岩储层含气量地球物理定量预测方法 |
CN110672813A (zh) * | 2018-07-02 | 2020-01-10 | 中国石油化工股份有限公司 | 一种页岩含气量计算方法 |
CN111175471A (zh) * | 2018-11-12 | 2020-05-19 | 中石化石油工程技术服务有限公司 | 一种岩屑录井方法 |
CN111188608A (zh) * | 2020-03-10 | 2020-05-22 | 中国石油天然气股份有限公司 | 页岩气井压裂改造方案的确定方法、装置及存储介质 |
CN111307652A (zh) * | 2018-12-11 | 2020-06-19 | 中国石油天然气股份有限公司 | 吸附气量的确定方法及装置 |
CN112304843A (zh) * | 2020-10-12 | 2021-02-02 | 四川省科源工程技术测试中心 | 一种泥页岩中页岩气吸附量定量表征方法 |
CN112613174A (zh) * | 2020-12-18 | 2021-04-06 | 西南石油大学 | 考虑多重吸附机制的页岩吸附甲烷能力评估方法 |
CN116773396A (zh) * | 2023-06-07 | 2023-09-19 | 中国地质调查局油气资源调查中心 | 地层温压条件下页岩总含气量获取方法和系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707333A (zh) * | 2012-06-06 | 2012-10-03 | 陕西延长石油(集团)有限责任公司研究院 | 页岩气资源量/储量的测量方法 |
CN104573344A (zh) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | 一种通过测井数据获取页岩储层含气量的方法 |
CN105184106A (zh) * | 2015-10-19 | 2015-12-23 | 中国石油大学(北京) | 一种确定泥页岩储层中吸附气含量的方法及装置 |
CN105842416A (zh) * | 2015-01-16 | 2016-08-10 | 中石化西南石油工程有限公司 | 一种基于兰格缪尔模型的页岩吸附气含量动态计算方法 |
-
2016
- 2016-08-12 CN CN201610663994.4A patent/CN106251232B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707333A (zh) * | 2012-06-06 | 2012-10-03 | 陕西延长石油(集团)有限责任公司研究院 | 页岩气资源量/储量的测量方法 |
CN104573344A (zh) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | 一种通过测井数据获取页岩储层含气量的方法 |
CN105842416A (zh) * | 2015-01-16 | 2016-08-10 | 中石化西南石油工程有限公司 | 一种基于兰格缪尔模型的页岩吸附气含量动态计算方法 |
CN105184106A (zh) * | 2015-10-19 | 2015-12-23 | 中国石油大学(北京) | 一种确定泥页岩储层中吸附气含量的方法及装置 |
Non-Patent Citations (1)
Title |
---|
PENG LI 等: "Prediction model for gas adsorption capacity of the Lower Ganchaigou Formation in the Qaidam Basin, China", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107271322A (zh) * | 2017-07-27 | 2017-10-20 | 中国石油大学(华东) | 一种定量评价泥页岩页岩气资源量及表征吸附气和游离气转化规律的方法 |
CN107818379A (zh) * | 2017-09-28 | 2018-03-20 | 中国石油天然气股份有限公司 | 吸附气量确定方法和装置 |
CN108150161A (zh) * | 2017-12-27 | 2018-06-12 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 页岩含气性评价方法装置 |
CN108593493A (zh) * | 2018-01-10 | 2018-09-28 | 中国石油天然气股份有限公司 | 一种确定致密储层含气量的方法和装置 |
CN108717202A (zh) * | 2018-05-31 | 2018-10-30 | 中国地质大学(武汉) | 一种考虑流体温度的页岩气异常地层压力预测方法 |
CN108846218B (zh) * | 2018-06-22 | 2020-07-28 | 中国石油大学(北京) | 预测泥页岩储层不同类型孔隙的内比表面积的方法及系统 |
CN108846218A (zh) * | 2018-06-22 | 2018-11-20 | 中国石油大学(北京) | 预测泥页岩储层不同类型孔隙的内比表面积的方法及系统 |
CN110672813B (zh) * | 2018-07-02 | 2021-09-14 | 中国石油化工股份有限公司 | 一种页岩含气量计算方法 |
CN110672813A (zh) * | 2018-07-02 | 2020-01-10 | 中国石油化工股份有限公司 | 一种页岩含气量计算方法 |
CN111175471A (zh) * | 2018-11-12 | 2020-05-19 | 中石化石油工程技术服务有限公司 | 一种岩屑录井方法 |
CN111175471B (zh) * | 2018-11-12 | 2022-06-24 | 中石化石油工程技术服务有限公司 | 一种岩屑录井方法 |
CN111307652A (zh) * | 2018-12-11 | 2020-06-19 | 中国石油天然气股份有限公司 | 吸附气量的确定方法及装置 |
CN111307652B (zh) * | 2018-12-11 | 2022-07-05 | 中国石油天然气股份有限公司 | 吸附气量的确定方法及装置 |
CN110018293A (zh) * | 2019-05-05 | 2019-07-16 | 西安石油大学 | 一种考虑水敏性的页岩黏土多相耦合甲烷含气量计算方法 |
CN110018293B (zh) * | 2019-05-05 | 2021-09-21 | 西安石油大学 | 一种考虑水敏性的页岩黏土多相耦合甲烷含气量计算方法 |
CN110579797A (zh) * | 2019-09-02 | 2019-12-17 | 长江师范学院 | 一种泥页岩储层含气量地球物理定量预测方法 |
CN111188608A (zh) * | 2020-03-10 | 2020-05-22 | 中国石油天然气股份有限公司 | 页岩气井压裂改造方案的确定方法、装置及存储介质 |
CN112304843A (zh) * | 2020-10-12 | 2021-02-02 | 四川省科源工程技术测试中心 | 一种泥页岩中页岩气吸附量定量表征方法 |
CN112304843B (zh) * | 2020-10-12 | 2022-09-30 | 四川省科源工程技术测试中心 | 一种泥页岩中页岩气吸附量定量表征方法 |
CN112613174A (zh) * | 2020-12-18 | 2021-04-06 | 西南石油大学 | 考虑多重吸附机制的页岩吸附甲烷能力评估方法 |
CN116773396A (zh) * | 2023-06-07 | 2023-09-19 | 中国地质调查局油气资源调查中心 | 地层温压条件下页岩总含气量获取方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106251232B (zh) | 2019-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106251232A (zh) | 确定页岩含气量的方法和装置 | |
CN105927218B (zh) | 一种陆相泥页岩储层含气量预测方法及装置 | |
Wang et al. | Quantitative characterization of reservoir space in the Lower Silurian Longmaxi Shale, southern Sichuan, China | |
US10466386B2 (en) | Methods for determining a saturation-height function in oil and gas reservoirs | |
CN105134195A (zh) | 一种基于测井资料的页岩气储层品质评价方法 | |
Yu et al. | Porosity estimation in kerogen-bearing shale gas reservoirs | |
CN101887132A (zh) | 一种井震结合定量预测砂岩储层流体饱和度的方法 | |
CN105487122B (zh) | 用于确定储层流体识别因子的敏感性的方法及系统 | |
CN104819923A (zh) | 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 | |
US20050216223A1 (en) | Method of evaluating the capillary pressure curve of an underground deposit rocks based on rock cuttings measurements | |
US20130229892A1 (en) | Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks | |
CN105093331B (zh) | 获取岩石基质体积模量的方法 | |
CN107102377B (zh) | 定量预测致密砂岩油气有利勘探区的方法 | |
CN107515957A (zh) | 泥页岩层序地层划分方法 | |
CN103630939A (zh) | 一种气层识别评价方法 | |
Wood | On the small strain stiffness of some scandinavian clays and impact on deep excavation | |
US8384389B2 (en) | Method of monitoring behavior of carbon dioxide in porous aquifer by using marine controlled-sources electromagnetic survey | |
CN103412335B (zh) | 一种利用地震物相体预测储层的方法 | |
CN106443800A (zh) | 一种中基性火山岩的核磁孔隙度校正方法 | |
Torrese et al. | 3D ERT imaging of the fractured-karst aquifer underlying the experimental site of Poitiers (France): comparing Wenner-Schlumberger, Pole-Dipole and hybrid arrays | |
Zhang et al. | Probing water partitioning in unsaturated weathered rock using nuclear magnetic resonance | |
CN103197348A (zh) | 利用各层内部样品进行加权编制测井交会图的方法 | |
CN101798922A (zh) | 天然气充盈度判断方法及其在低丰度气藏中的应用 | |
CN110954942A (zh) | 粘土含量识别的图版制作方法及系统 | |
CN104749618A (zh) | 泥页岩缓倾角裂缝叠后概率定量表征方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |