CN117035508B - 基于多因素的致密气藏河道砂体连通性评价方法 - Google Patents
基于多因素的致密气藏河道砂体连通性评价方法 Download PDFInfo
- Publication number
- CN117035508B CN117035508B CN202310958006.9A CN202310958006A CN117035508B CN 117035508 B CN117035508 B CN 117035508B CN 202310958006 A CN202310958006 A CN 202310958006A CN 117035508 B CN117035508 B CN 117035508B
- Authority
- CN
- China
- Prior art keywords
- river channel
- deposition
- river
- sand
- type
- 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.)
- Active
Links
- 239000004576 sand Substances 0.000 title claims abstract description 184
- 238000011156 evaluation Methods 0.000 title claims abstract description 35
- 238000011161 development Methods 0.000 claims abstract description 42
- 238000004458 analytical method Methods 0.000 claims abstract description 32
- 238000011160 research Methods 0.000 claims abstract description 32
- 239000013049 sediment Substances 0.000 claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000004044 response Effects 0.000 claims abstract description 18
- 230000007246 mechanism Effects 0.000 claims abstract description 12
- 238000000151 deposition Methods 0.000 claims description 114
- 230000008021 deposition Effects 0.000 claims description 110
- 239000010410 layer Substances 0.000 claims description 53
- 238000004891 communication Methods 0.000 claims description 39
- 230000000704 physical effect Effects 0.000 claims description 23
- 239000011435 rock Substances 0.000 claims description 16
- 239000002699 waste material Substances 0.000 claims description 16
- 230000014759 maintenance of location Effects 0.000 claims description 15
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 15
- 230000033764 rhythmic process Effects 0.000 claims description 12
- 239000011229 interlayer Substances 0.000 claims description 10
- 208000035126 Facies Diseases 0.000 claims description 7
- 238000005520 cutting process Methods 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 4
- 230000002349 favourable effect Effects 0.000 claims description 3
- 239000005416 organic matter Substances 0.000 claims description 3
- 238000009991 scouring Methods 0.000 claims description 3
- 238000004898 kneading Methods 0.000 claims 1
- 230000004927 fusion Effects 0.000 abstract description 5
- 239000007789 gas Substances 0.000 description 30
- 230000003068 static effect Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 238000004519 manufacturing process Methods 0.000 description 9
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000011158 quantitative evaluation Methods 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 239000003345 natural gas Substances 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000004568 cement Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 241001415846 Procellariidae Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011449 brick Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- -1 pressure test Substances 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000004936 stimulating effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Human Resources & Organizations (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Strategic Management (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Economics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Primary Health Care (AREA)
- Operations Research (AREA)
- Marine Sciences & Fisheries (AREA)
- Game Theory and Decision Science (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Quality & Reliability (AREA)
- Revetment (AREA)
Abstract
本发明提供基于多因素的致密气藏河道砂体连通性评价方法,包括(1)沉积微相研究、(2)河道充填样式、(3)河道成因类型及特征、(4)砂体连通类型及特征、(5)砂体连通性识别。本发明通过对井上数据的分析,高效快速完成单井点的河道砂体连通性评价;对实现高效开采单井点致密砂岩气提供理论基础,采用地震约束响应机制,分析地震反演图属性,完成整条河道砂体连通性评价;将地震响应机制应用在油气勘探开发领域,实现了多学科融合交叉发展,同时快速完成河道砂体连通性评价,为高效开发致密砂岩气提供完整的知识理论体系,为早日实现川渝气大庆奠定稳定的理论根基。
Description
技术领域
本发明属于地质分析技术领域,尤其涉及基于多因素的致密气藏河道连通性评价方法。
背景技术
行业发展中补贴的延续将对国内非常规天然气勘探开发起到很大的刺激作用。由于特殊的地质结构与地理条件,相比较于其他的平原国家,中国的致密气开采更具有发展潜力,尤其是在当前国家经济腾飞的新形势下,致密气的大力勘探与开发对中国工业发展以及社会发展具有重要的推动作用,同时,要深入推进能源革命,加大油气资源勘探开发和增储上产力度。中国致密砂岩气资源量约为112万亿立方米(部分与常规气在资源量上存在着交叉),广泛分布于鄂尔多斯、四川、松辽、渤海湾、柴达木、塔里木及准噶尔等10余个盆地,其中鄂尔多斯和四川盆地最为丰富。据2022年天然气年产量报告分析,四川年产气量452.4亿立方米,年产气量排名全国第二,截止2022年9月30日,川中致密气日产量达到650万方,年累产气达到9亿方,与去年同期相比增长401%。致密气已成为我国天然气增储上产的重要领域,在天然气工业发展中占有非常重要的地位。我国致密气完全动用难度很大,主要原因是关键技术尚需进一步攻关,对于难以动用的储量,高投入产生的高成本与气价不匹配,完全没有经济效益。因此,关键技术的突破能快速实现致密气的高效开采的发展目标,为实现我国增储上产添砖加瓦。
砂体连通性是指各成因单元砂体在横向上和纵向上相互接触方式及连通程度,它决定了储层的连续性,是导致储层平面非均质性的主要因素。根据所用资料的不同,砂体连通性的研究方法包括动态分析法和静态分析法。动态分析法主要利用动态资料,比如示踪剂、压力测试及水淹状况分析等,可以直观地通过动态资料分析砂体的连通状况,其连通性评价结果有较高的准确性和可信性;静态分析法主要利用静态资料,包括岩芯、测井及地震资料等,在地下砂体展布刻画的基础上,通过综合分析砂体的沉积类型、发育程度及规模、相变接触关系来判断各类砂体间的连通状况。
首先要谈到的是静态的分析法,前人利用“砂体分析”法,从河流的类型、规模、河道的充填程度、砂体接触关系及相变特征入手,分析研究了河流沉积中河道砂与河间砂、不同河道砂体间的连通关系,认为河道砂与被切割的河间砂的连通程度取决于其间的废弃河道泥质充填程度,而废弃河道的发育程度又受河型和具体的沉积部位控制;不同河道砂体间的连通关系则不仅与河型及沉积部位有关,还与砂体发育的期次和相互接触关系有关。同一河道砂体的连通性主要受砂体内部的不连续薄夹层及物性分布规律影响;不同河道砂体之间的连通性与两河道的接触位置有关,由于好的点坝砂发育在河道的凸岸,而废弃河道的发育程度又受河型和具体的沉积部位控制;不同河道砂体间的连通关系则不仅与河型及沉积部位有关,还与砂体发育的期次和相互接触关系有关。同一河道砂体的连通性主要受砂体内部的不连续薄夹层及物性分布规律影响;不同河道砂体之间的连通性与两河道的接触位置有关,由于好的点坝砂发育在河道的凸岸,而废弃河道的泥质充填发育在河道的凹岸,所以两河道以凸岸相接,则连通性好,两河道以凹岸相接,则连通性一般较差,两河道相接位置是“一凸一凹”,则连通性取决于中间泥质充填的保留程度;至于河间砂与河道砂之间的连通状况,不仅取决于河间砂的发育位置,还与他们之间的层位关系有关,处于河道凸岸一侧的河间砂与主河道的连通性优于处于凹岸侧的,河间砂与河道砂层位相互对应时的连通性优于层位不对应的。而动态分析法根据所用动态资料的不同,分析的难易程度及可信程度也不同。一些学者等将油藏近似为一个动态平衡系统,应用系统分析方法,建立了油藏井间动态连通性反演模型。根据该模型可利用注采资料建立各注采井之间连通系数的多元一次方程组,然后通过拟牛顿迭代解方程,得到井与井之间的连通系数,该连通系数可为砂体连通性评价提供理论依据。关于砂体连通性影响因素研究,学者Allen基于冲积环境中规则河道的研究,认为砂体连通程度与砂岩密度有较高的相关性,当砂岩密度大于50%时,砂体的连通程度会迅速提高。Hovadik and Larue指出砂地比对河道砂体连通性的影响作用最大,他们通过对河流相储层砂体连通规模的三维建模研究发现,砂体连通程度与砂地比的关系曲线呈“S”形变化,当砂地比大于35%时,砂体连通程度随着砂地比的升高出现瀑布式陡增,而其他的一些因素,例如宽厚壁、河漫滩的沉积率、连续的泥岩层及河道之间的角度等,只是使该“S”型曲线左右偏移。Pranter和Sommer以Piceance盆地晚白垩世Williams Fork地层为例,通过野外露头和地质模拟的方法,分析了河道砂体静态连通性的控制因素,认为砂体宽度、砂地比和井距堆砂体的连通程度影响较大,当井距较大时(大于800m),砂地比小于20%时,砂体连通性较差,当砂地比增加时,连通程度逐渐增加,当砂地比增加到30%后。砂体连通性趋于平稳;当井距较小时(小于200m),砂体间的连通程度主要由砂体发育的宽度决定,此时如果河道砂体发育的宽度较大,井距也足够小,砂体间的连通程度将非常高。
目前来看,国内外关于井间动态连通性反演模型的研究评价方法较为单一,没有考虑地层砂体连通条件对井间动态连通的影响,缺少结合地层砂体叠置模式的井间动静态一体化对井间连通进行识别的方法。
现有技术一的方案
基于动静态信息融合的致密气藏河道砂体连通性验证方法,首先确定河道砂体的砂体叠置样式,然后确定河道砂体的垂直构型组合类型,通过确定河道砂体的不同砂体构型的砂泥组合类型来确定了河道砂体的构型单元组合类型,然后确定河道砂体岩相,之后是结合地震属性分析再根据前面的研究基础来根据河道砂体的参数确定河道砂体的砂体连通性。最后通过结合生产井的动态数据验证以及通过静态信息融合开展河道砂体连通性验证,从而高效且准确地评价储层非均质性强弱,为我国致密气藏生产开发提供充分的理论依据。
现有技术一的缺点
该技术仅对河道砂体连通性进行验证,但并未对河道砂体连通性进行定性或定量的评价。
现有技术二的方案
基于沉积环境和评价尺度的井间砂体静态连通性评价方法,申请号为CN202111471432.7,其包括:单一砂体输导层进行识别和划分;将地层砂地比、沉积环境、评价尺度、砂体几何学参数作为输导层评价主控因素,以其中某一评价主控因素为对象,采用因素控制阀,分别明确其他三个评价主控因素对疏导层静态连通性的控制作用;建立单一砂体疏导层三维静态连通模型,采用Petrel地质建模软件实现建模,并实现定量化输出。
该技术主要是根据砂体几何学参数、地层砂地比等作为主控因素,研究井间的河道砂体连通性。
发明内容
本发明的目的在于解决上述现有技术存在的缺陷,提供基于多因素的致密气藏河道砂体连通性评价方法,综合考虑测井相标志以及岩芯的岩石学特征确定沉积微相,然后根据沉积微相在垂向上的序列分析河道充填样式,判断河道成因类型,最后根据地震响应约束机制等形成一种基于“一种基于地质成因+钻井约束+地震响应”的致密气藏河道砂体连通性评价新技术,为国家增储上产发展提供理论支撑。
本发明通过采用地震响应机制法和井间砂体对比法,通过控制变量法对比研究了宽厚比,砂泥比、砂层厚度等参数对砂体连通性的影响,应用Fisher线性判定函数,确定了砂体连通性判断的连通系数。在沉积、成岩的研究基础上,按照“储层非均质性-砂体构型-河道连通性”提出了致密气藏河道砂体连通性评价技术。在一定程度上解决地层砂体连通条件对井间连通的影响判断机制,同时也能够快速应用砂体连通系数对致密气藏河道砂体连通性进行评价,同时对多个单井进行分析,在短时间内可以对整个河道的砂体连通性进行分析,达到短时高效的研究目的。
本发明采用如下技术方案:
基于多因素的致密气藏河道砂体连通性评价方法,包括:
(1)沉积微相研究
结合单井的测井曲线,分析其测井响应特征以及测井相特征,按照测井相特征中的自然伽马电位、密度、孔隙度以及渗透率等参数将沉积微相分为边滩、河床滞留沉积、废弃河道、天然堤+决口扇和泛滥平原五种类型,基于GR曲线特征分析,边滩的GR曲线整体高幅、呈箱式;河床滞留沉积的GR曲线整体中高幅、呈微齿状;决口扇的GR曲线整体呈低中幅、指状。泛滥平原的GR曲线整体低幅、呈锯齿状且平直。其次基于岩芯,结合其岩芯的颜色特征、沉积构造、粒度对单井点的岩性进行判断,其中颜色是反映沉积环境的最直观特征。位于研究1区的岩性主要以灰色-灰白色中-细砂岩、灰色-灰绿色粉砂岩、红褐色-灰绿色泥岩为主,可见暗色有机质沉积。通过对研究区的岩石的岩性判断,然后结合研究1区的岩石粒度大小以及沉积层序来初步确定沉积微相。粒度大小所绘制出的曲线能够在一定程度上反演岩石组合类型,而沉积层序有利于识别不同的沉积旋回和沉积构造,最后结合沉积构造来开展沉积微相的研究,水动力条件的强弱、沉积速度、水流作用方式都是决定沉积构造类型的重要因素。不同介质条件下,上述三种特征不相同,所形成的沉积构造截然不同。分流河道、水下分流河道沉积多发育小型槽状交错层理、小型楔状交错层理、平行层理、沙纹层理、滑塌构造和冲刷构造;河口坝沉积多发育小型槽状交错层理,波状层理;分流间湾中见揉皱构造、水平层理。这些层理规模相对较小,水动力相对较强,多出现在浅水三角洲沉积内。而研究1区中出现浅水三角洲沉积体系,包括浅水三角洲平原、内前缘沉积,这种沉积体系中主要发育块状层理、大型的交错层理、小型交错层理、平行层理、水平层理等沉积构造,可通过其所形成的沉积构造对沉积微相进行深度研究。
(2)河道充填样式
基于沉积微相的识别,分析沉积微相在垂向上的组合序列,可以把河道充填样式分为七大类,分别为废弃河道+边滩、边滩+河床滞留、边滩+边滩+边滩、水下分流河道+水下分流河道、水下分流间湾+水下天然堤+水下分流间湾+水下分流河道+水下分流间湾、水下分流河道+河口坝、泛滥平原+天然堤+决口扇+边滩+泛滥平原。
(3)河道成因类型及特征
基于河道充填样式,分析河道的砂体物性特征和沉积特征将河道成因类型分为五类,包括多期叠加型、深切型、废弃型、无期叠加型和进叠加型。多期叠加型河道沉积类型为边滩沉积加边滩沉积或水下分流河道沉积加水下分流间湾沉积以及水下天然堤沉积;深切型河道沉积类型为边滩沉积加河床滞留沉积;废弃型河道沉积类型为废弃河道沉积加边滩沉积;无期叠加型河道沉积类型为边滩沉积、决口扇沉积、天然堤沉积加泛滥平原沉积或水下分流河道沉积、水下分流间湾沉积加水下天然堤沉积,进叠加型河道沉积类型为水下分流河道沉积加河口坝沉积。从对河道物性特征展开的研究中可以得到,发育边滩+边滩构型单元组合的砂体物性最好;边滩+河床滞留沉积、边滩+废弃河道、水下分流河道+水下分流河道、水下分流河道+河口坝次之;边滩+决口扇+天然堤、水下分流河道+水下天然堤最差。在沉积特征研究方面,边滩、水下分流河道发育块状层理、大型交错层理岩相为主,分别发育均匀韵律和正韵律结构,物性整体较好;河口坝发育反韵律物性次之;天然堤、决口扇等薄层砂体复合韵律发育,物性整体较差。
(4)砂体连通类型及特征
通过河道成因类型的分类,将河道砂体连通性分为通畅、局部通畅以及不通畅三类,通畅的河道成因类型包括多期叠加型和深切型,二者的砂泥比大于0.94,宽深比大于40,砂体厚度大于19m,隔夹层厚度小于28.32m,地震振幅均方根值(RMS)大于8500,纵横波速比值小于1.65,根据参数分析为通畅性好;局部通畅的河道包括进积叠加型和无叠加单期次型,二者的砂泥比处于0.59-0.94之间,宽深比33-40,砂体厚度16-19m,隔夹层厚度28.32-52m,地震振幅均方根值(RMS)为8000-8500,纵横波速比值大于1.65,根据参数分析为局部通畅;不通畅的河道成因类型为废弃型,其砂泥比小于0.59,宽深比小于33,砂体厚度小于16m,隔夹层厚度小于52m。地震振幅均方根值(RMS)处于7000-8500,纵横波速比值处于1.6左右,根据参数分析为通畅性差。
(5)砂体连通性识别
应用岩石学特征、测井相特征等对河道沉积微相类型判断,通过沉积微相在垂向上的组合序列确定河道充填样式,基于七大类河道充填样式将河道成因类型分为五类,基于单井河道类型识别和砂泥比、宽深比等参数应用Fisher函数确定多个参数对河道砂体连通性的影响,通过研究发现、砂泥比的比值越大,河道砂体连通性越好;宽厚比的比值越大,河道砂体越通畅;夹层厚度越小,砂体越通畅。根据多个参数的影响,采用Fisher线性判别函数确定河道的连通系数,完成单井点河道砂体连通程度评价,再结合地震振幅响应约束,识别地震振幅均方根图(RMS)和地震响应约束机制完成各条河道砂体整体连通性识别和评价。
本发明的有益效果:
本发明与基于动静态信息融合的致密气藏河道砂体连通性验证方法相比,提出了河道砂体的连通性的评价方法,能够实现河道砂体连通性的半定量评价。
本发明与基于沉积环境和评价尺度的井间砂体静态连通性评价方法相比,本发明侧重于根据砂泥比、宽深比、砂体厚度等参数,半定量研究整条河道砂体连通性。
前人主要基于岩石学特征中的岩芯颜色和沉积构造类型方面来对沉积微相进行研究,在前人研究的基础之上,本发明还结合了测井相特征,更为动态和精确的对沉积微相进行评价。同时,前人主要根据河道成因类型、河道砂体几何参数等对河道连通性完成定性评价,本发明基于各类参数,采用Fisher线性判别函数计算河道连通系数,结合地震响应约束机制,分析地震均方根图,首次完成整条河道的半定量评价。
本发明基于对河道沉积微相、河道充填样式以及河道成因类型的研究,首次将测井曲线与致密砂岩气结合起来,评价其河道成因类型,确定河道砂体连通性,计算连通系数,首次完成了对河道砂体连通性的半定量评价。首次结合地震反演图的特征对河道砂体连通性评价,同时首次提出用地震响应机制对河道砂体连通性进行评价,完成了从单井到井间最后到整条河道的砂体连通性评价,形成一种基于“地质成因+钻井约束+地震响应”的致密气藏河道砂体连通性评价方法,为我国致密气藏生产开发提供充分的理论依据。
本发明通过对井上数据的分析,高效快速完成单井点的河道砂体连通性评价;对实现高效开采单井点致密砂岩气提供理论基础,采用地震约束响应机制,分析地震反演图属性,完成整条河道砂体连通性评价;将地震响应机制应用在油气勘探开发领域,实现了多学科融合交叉发展,同时快速完成河道砂体连通性评价,为高效开发致密砂岩气提供完整的知识理论体系,为早日实现川渝气大庆奠定稳定的理论根基。
附图说明
图1为三角洲平原沉积微相测井相识别图版;
图2为微相纵向发育序列与对应河道成因类型;
图3为粒度分布参数统计图Ⅰ;
图4为粒度分布参数统计图Ⅱ;
图5(a)、图5(b)为通过沉积微相序列判断河道充填样式图;
图6为通过河道充填样式判断河道成因类型图;
图7为边滩+边滩物性特征分析图;
图8为边滩+河床滞留物性特征图;
图9为废弃河道+边滩物性特征;
图10为水下分流河道+水下分流河道物性分析图;
图11为水下分流河道+河口坝河道物性分析图;
图12为边滩+决口扇+天然堤河道物性分析图;
图13为水下分流河道+天然堤河道物性分析图;
图14为砂泥比对地震均方根值的影响;
图15为宽深比对地震均方根值的影响;
图16为隔夹层厚度对地震均方根值的影响;
图17为地震均方根值对纵横波速比的影响;
图18为砂泥比对纵横波速比的影响;
图19为宽深比对纵横波速比的影响;
图20为河道连通性半定量评价标准;
图21为Fisher判别函数确定的分类系数;
图22为河道砂体连通性识别流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明通过岩石学特征、沉积相标志和测井相特征判断沉积微相,然后根据纵向上的沉积微相序列来判断河道充填样式,通过不同的河道充填样式来确定河道成因类型,依据河道的砂泥比、宽厚比、砂体厚度等参数,应用Fisher线性判别函数分析多个参数对河道砂体连通性的影响,确定砂体连通性系数,最后结合地震反演图的属性如RMS、纵横波速比等完成整条河流的砂体连通性评价。
如图22所示,本发明的基于多因素的致密气藏河道砂体连通性评价方法,包括:
(1)沉积微相研究
结合单井的测井曲线,分析其测井响应特征以及测井相特征,按照测井相特征中的自然伽马电位、密度、孔隙度以及渗透率等参数将沉积微相分为边滩、河床滞留沉积、废弃河道、天然堤+决口扇和泛滥平原五种类型,基于GR曲线特征分析,边滩的GR曲线整体高幅、呈箱式;河床滞留沉积的GR曲线整体中高幅,呈微齿化锯齿状;决口扇的GR曲线整体呈低中幅、指状。泛滥平原的GR曲线整体低幅、呈锯齿状且平直。(如图1所示)其次基于岩芯,结合其岩芯的颜色特征、沉积构造、粒度对单井点的岩性进行评价,其中颜色是反映沉积环境的最直观特征。位于研究1区的岩性主要以灰色-灰白色种-细砂岩、灰色-灰绿色粉砂岩、红褐色-灰绿色泥岩为主,可见暗色有机质沉积。通过对研究区的岩石的岩性判断,然后结合研究区1区的岩石粒度大小(如图2、图3所示)以及沉积层序来初步确定沉积微相。粒度大小所绘制出的曲线能够在一定程度上反演岩石组合类型,而沉积层序有利于识别不同的沉积旋回和沉积构造,最后结合沉积构造来开展沉积微相的研究,水动力条件的强弱、沉积速度、水流作用方式都是决定沉积构造类型的重要因素(如图4所示)。不同介质条件下,上述三种特征不相同,所形成的沉积构造截然不同。分流河道、水下分流河道沉积多发育小型槽状交错层理、小型楔状交错层理、平行层理、沙文层理、滑塌构造和冲刷构造;河口坝沉积多发育小型槽状交错层理、波纹层理;分流间湾中见揉皱构造、水平层理。这些层理规模相对较小,水动力相对较强,多出现在浅水三角洲沉积内。而研究1区中出现浅水三角洲沉积体系,包括浅水三角洲平原、内前缘沉积,这种沉积体系中主要发育块状层理、大型的交错层理、小型交错层理、平行层理、水平层理等沉积构造,可通过其所形成的沉积构造对沉积微相进行深度研究。
(2)河道充填样式
基于沉积微相的识别,分析沉积微相在垂向上的组合序列,可以把河道充填样式分为七大类(如图5(a)、图5(b)所示),分别为废弃河道+边滩、边滩+河床滞留、边滩+边滩+边滩、水下分流河道+水下分流河道、水下分流间湾+水下天然堤+水下分流间湾+水下分流河道+水下分流间湾、水下分流河道+河口坝、泛滥平原+天然堤+决口扇+边滩+泛滥平原。
(3)河道成因类型及特征
基于河道充填样式,分析河道的孔隙度、渗透率、饱和度等砂体物性特征和沉积特征将河道成因类型分为五类,包括多期叠加型、深切型、废弃型、无期叠加型和进叠加型(如图6所示)。多期叠加型河道沉积类型为边滩沉积加边滩沉积或水下分流河道沉积加水下分流间湾沉积以及水下天然堤沉积;深切型河道沉积类型为边滩沉积加河床滞留沉积;废弃型河道沉积类型为废弃河道沉积加边滩沉积;无期叠加型河道沉积类型为边滩沉积、决口扇沉积、天然堤沉积加泛滥平原沉积或水下分流河道沉积、水下分流间湾沉积加水下天然堤沉积,进叠加型河道沉积类型为水下分流河道沉积加河口坝沉积。从对河道物性特征展开的研究中可以得到,发育边滩+边滩构型单元组合的砂体物性最好(如图7所示);边滩+河床滞留沉积(如图8所示)、边滩+废弃河道(如图9所示)、水下分流河道+水下分流河道(如图10所示)、水下分流河道+河口坝次之(如图11);边滩+决口扇+天然堤(如图12所示)、水下分流河道+水下天然堤(如图13)最差。在沉积特征研究方面,边滩、水下分流河道发育块状层理、大型交错层理岩相为主,分别发育均匀韵律和正韵律结构,物性整体较好;河口坝发育反韵律物性次之;天然堤、决口扇等薄层砂体复合韵律发育,物性整体较差。
(4)砂体连通类型及特征
通过河道成因类型的分类,将河道砂体连通性分为通畅、局部通畅以及不通畅三类,通畅的河道成因类型包括多期叠加型和深切型,二者的砂泥比大于0.94,宽深比大于40,砂体厚度大于19m,隔夹层厚度小于28.32m,地震振幅均方根值(RMS)大于8500,纵横波速比值小于1.65,根据参数分析为通畅性好(图14、图15、图16、图17、图18、图19);局部通畅的河道包括进积叠加型和无叠加单期次型,二者的砂泥比处于0.59-0.94之间,宽深比33-40,砂体厚度16-19m,隔夹层厚度28.32-52m,地震振幅均方根值(RMS)为8000-8500,纵横波速比值大于1.65,根据参数分析为局部通畅(图14、图15、图16、图17、图18、图19);不通畅的河道成因类型为废弃型,其砂泥比小于0.59,宽深比小于33,砂体厚度小于16m,隔夹层厚度小于52m。地震振幅均方根值(RMS)处于7000-8500,纵横波速比值处于1.6左右,根据参数分析为通畅性差(图14、图15、图16、图17、图18、图19)。
(5)砂体连通类型及特征
通过河道成因类型的分类,将河道砂体连通性分为通畅、局部通畅以及不通畅三类,通畅的河道成因类型包括多期叠加型和深切型,二者的砂泥比大于0.94,宽深比大于40,砂体厚度大于19m,隔夹层厚度小于28.32m,地震振幅均方根值(RMS)大于8500,纵横波速比值小于1.65,根据参数分析为通畅性好(图14、图15、图16、图17、图18、图19);局部通畅的河道包括进积叠加型和无叠加单期次型,二者的砂泥比处于0.59-0.94之间,宽深比33-40,砂体厚度16-19m,隔夹层厚度28.32-52m,地震振幅均方根值(RMS)为8000-8500,纵横波速比值大于1.65,根据参数分析为局部通畅(图14、图15、图16、图17、图18、图19);不通畅的河道成因类型为废弃型,其砂泥比小于0.59,宽深比小于33,砂体厚度小于16m,隔夹层厚度小于52m。地震振幅均方根值(RMS)处于7000-8500,纵横波速比值处于1.6左右,根据参数分析为通畅性差(图14、图15、图16、图17、图18、图19)。
通过自然伽马曲线对河道沉积微相类型评价,结合岩芯确定河道充填样式,基于单井河道类型识别和砂泥比、宽深比等参数应用Fisher函数确定多个参数对河道砂体连通性的影响,通过研究发现,砂泥比的比值越大,河道砂体连通性越好;宽厚比的比值越大,河道砂体越通畅;夹层厚度越小,砂体越通畅。根据多个参数的影响,采用Fisher线性判别函数,如图20、图21所示,确定河道的连通系数,完成单井点河道砂体连通程度评价,再结合地震振幅响应约束,识别地震振幅均方根图(RMS)和地震响应约束机制完成各条河道砂体整体连通性识别和评价。
本发明首先根据岩石学特征、测井相特征来确定沉积微相并刻画岩性曲线,然后根据纵向上的沉积微相序列确定河道充填样式,根据七类河道充填样式确定五类河道成因类型,基于对河道成因类型的研究,分析砂泥比、宽深比、砂体厚度等影响因素,应用Fisher线性判别函数综合计算砂体连通系数,结合地震响应约束机制和地震均方根图,完成对河道砂体连通性的评价,形成一种基于“地质成因+钻进约束+地震响应”的致密气藏河道砂体连通性评价方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (6)
1.基于多因素的致密气藏河道砂体连通性评价方法,其特征在于,包括:
(1)沉积微相研究
结合单井的测井曲线,分析其测井响应特征以及测井相特征,按照测井相特征中的自然伽马电位、密度、孔隙度以及渗透率将沉积微相分为边滩、河床滞留沉积、废弃河道、天然堤+决口扇和泛滥平原五种类型,基于GR曲线特征分析,边滩的GR曲线整体高幅、呈箱式;河床滞留沉积的GR曲线整体中高幅、呈微齿状;决口扇的GR曲线整体呈低中幅、指状;泛滥平原的GR曲线整体低幅、呈锯齿状且平直;
其次基于岩芯,结合其岩芯的颜色特征、沉积构造、粒度对单井点的岩性进行判断,粒度大小所绘制出的曲线能够反演岩石组合类型,而沉积层序有利于识别不同的沉积旋回和沉积构造;
最后结合沉积构造来开展沉积微相的研究,水动力条件的强弱、沉积速度、水流作用方式都是决定沉积构造类型的重要因素,不同介质条件下,上述三种特征不相同,所形成的沉积构造截然不同,分流河道、水下分流河道沉积多发育小型槽状交错层理、小型楔状交错层理、平行层理、沙纹层理、滑塌构造和冲刷构造;河口坝沉积多发育小型槽状交错层理,波状层理;分流间湾中见揉皱构造、水平层理;
(2)河道充填样式
基于沉积微相的识别,分析沉积微相在垂向上的组合序列,把河道充填样式分为七大类,分别为废弃河道+边滩、边滩+河床滞留、边滩+边滩+边滩、水下分流河道+水下分流河道、水下分流间湾+水下天然堤+水下分流间湾+水下分流河道+水下分流间湾、水下分流河道+河口坝、泛滥平原+天然堤+决口扇+边滩+泛滥平原;
(3)河道成因类型及特征
基于河道充填样式,分析河道的砂体物性特征和沉积特征将河道成因类型分为五类,包括多期叠加型、深切型、废弃型、无期叠加型和进叠加型;
多期叠加型河道沉积类型为边滩沉积加边滩沉积或水下分流河道沉积加水下分流间湾沉积以及水下天然堤沉积;深切型河道沉积类型为边滩沉积加河床滞留沉积;废弃型河道沉积类型为废弃河道沉积加边滩沉积;无期叠加型河道沉积类型为边滩沉积、决口扇沉积、天然堤沉积加泛滥平原沉积或水下分流河道沉积、水下分流间湾沉积加水下天然堤沉积;进叠加型河道沉积类型为水下分流河道沉积加河口坝沉积;
(4)砂体连通类型及特征
通过河道成因类型的分类,将河道砂体连通性分为通畅、局部通畅以及不通畅三类,通畅的河道成因类型包括多期叠加型和深切型,二者的砂泥比大于0.94,宽深比大于40,砂体厚度大于19m,隔夹层厚度小于28.32m,地震振幅均方根值大于8500,纵横波速比值小于1.65,根据参数分析为通畅性好;
局部通畅的河道包括进积叠加型和无叠加单期次型,二者的砂泥比处于0.59-0.94之间,宽深比33-40,砂体厚度16-19m,隔夹层厚度28.32-52m,地震振幅均方根值为8000-8500,纵横波速比值大于1.65,根据参数分析为局部通畅;不通畅的河道成因类型为废弃型,其砂泥比小于0.59,宽深比小于33,砂体厚度小于16m,隔夹层厚度小于52m,地震振幅均方根值处于7000-8500,纵横波速比值处于1.6,根据参数分析为通畅性差;
(5)砂体连通性识别
应用岩石学特征、测井相特征对河道沉积微相类型判断,通过沉积微相在垂向上的组合序列确定河道充填样式,基于七大类河道充填样式将河道成因类型分为五类,基于单井河道类型识别和砂泥比、宽深比应用Fisher函数确定多个参数对河道砂体连通性的影响,通过研究发现、砂泥比的比值越大,河道砂体连通性越好;
宽厚比的比值越大,河道砂体越通畅;夹层厚度越小,砂体越通畅;根据多个参数的影响,采用Fisher线性判别函数确定河道的连通系数,完成单井点河道砂体连通程度评价,再结合地震振幅响应约束,识别地震振幅均方根图和地震响应约束机制完成各条河道砂体整体连通性识别和评价。
2.根据权利要求1所述的基于多因素的致密气藏河道砂体连通性评价方法,其特征在于,(1)沉积微相研究中,岩芯颜色是反映沉积环境的最直观特征。
3.根据权利要求1所述的基于多因素的致密气藏河道砂体连通性评价方法,其特征在于,(1)沉积微相研究中,当位于研究1区的岩性为灰色-灰白色中-细砂岩、灰色-灰绿色粉砂岩、红褐色-灰绿色泥岩时,见暗色有机质沉积,通过对研究区的岩石的岩性判断,然后结合研究1区的岩石粒度大小以及沉积层序来初步确定沉积微相。
4.根据权利要求1所述的基于多因素的致密气藏河道砂体连通性评价方法,其特征在于,(1)沉积微相研究中,当研究1区中出现浅水三角洲沉积体系,包括浅水三角洲平原、内前缘沉积,这种沉积体系中发育块状层理、大型的交错层理、小型交错层理、平行层理、水平层理的沉积构造,通过其所形成的沉积构造对沉积微相进行深度研究。
5.根据权利要求1所述的基于多因素的致密气藏河道砂体连通性评价方法,其特征在于,(3)河道成因类型及特征中,从对河道物性特征展开的研究中得到,发育边滩+边滩构型单元组合的砂体物性最好;边滩+河床滞留沉积、边滩+废弃河道、水下分流河道+水下分流河道、水下分流河道+河口坝次之;边滩+决口扇+天然堤、水下分流河道+水下天然堤最差。
6.根据权利要求1所述的基于多因素的致密气藏河道砂体连通性评价方法,其特征在于,(3)河道成因类型及特征中,在沉积特征研究方面,边滩、水下分流河道发育块状层理、大型交错层理岩相为主,分别发育均匀韵律和正韵律结构,物性整体较好;河口坝发育反韵律物性次之;天然堤、决口扇的薄层砂体复合韵律发育,物性整体较差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310958006.9A CN117035508B (zh) | 2023-08-01 | 2023-08-01 | 基于多因素的致密气藏河道砂体连通性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310958006.9A CN117035508B (zh) | 2023-08-01 | 2023-08-01 | 基于多因素的致密气藏河道砂体连通性评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117035508A CN117035508A (zh) | 2023-11-10 |
CN117035508B true CN117035508B (zh) | 2024-06-18 |
Family
ID=88639620
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310958006.9A Active CN117035508B (zh) | 2023-08-01 | 2023-08-01 | 基于多因素的致密气藏河道砂体连通性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117035508B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109184677A (zh) * | 2018-09-26 | 2019-01-11 | 西南石油大学 | 用于非均质互层砂体的储层评价方法 |
CN111784160A (zh) * | 2020-06-30 | 2020-10-16 | 中国水利水电科学研究院 | 一种河流水文情势变化的评估方法及系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112084660B (zh) * | 2020-09-10 | 2022-05-31 | 西南石油大学 | 基于岩电解释模型对深层/超深层碳酸盐岩沉积微相精细划分的方法 |
CN113159994B (zh) * | 2021-04-25 | 2024-01-23 | 黄河勘测规划设计研究院有限公司 | 多泥沙河流的船闸淤积处理方法及装置 |
CN115097525A (zh) * | 2022-07-11 | 2022-09-23 | 成都理工大学 | 一种基于多信息融合的致密气藏砂体构型分析方法 |
CN114994789A (zh) * | 2022-07-11 | 2022-09-02 | 成都理工大学 | 基于动静态信息融合的致密气藏河道砂体连通性验证方法 |
CN114966886B (zh) * | 2022-07-11 | 2023-05-05 | 成都理工大学 | 一种基于构型的致密气藏砂体连通性等级判定方法 |
CN116341925B (zh) * | 2023-02-23 | 2024-06-18 | 成都理工大学 | 致密砂岩优质储层连续性评价系统、方法、设备及终端 |
-
2023
- 2023-08-01 CN CN202310958006.9A patent/CN117035508B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109184677A (zh) * | 2018-09-26 | 2019-01-11 | 西南石油大学 | 用于非均质互层砂体的储层评价方法 |
CN111784160A (zh) * | 2020-06-30 | 2020-10-16 | 中国水利水电科学研究院 | 一种河流水文情势变化的评估方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN117035508A (zh) | 2023-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106569287B (zh) | 一种基于复合砂体构型模式的扇三角洲前缘储层预测方法 | |
Gani et al. | Basic building blocks and process variability of a Cretaceous delta: internal facies architecture reveals a more dynamic interaction of river, wave, and tidal processes than is indicated by external shape | |
Leckie et al. | Storm-and tide-dominated shorelines in Cretaceous Moosebar-Lower Gates interval—outcrop equivalents of deep basin gas trap in western Canada | |
CN105373648B (zh) | 一种砂泥岩互层型沉积体储层构型的建模方法 | |
CN110838175B (zh) | 注气开发油藏的地质模型建立方法 | |
CN111706317B (zh) | 一种确定加密调整区低渗储层剩余油分布状况的方法 | |
CN108828666B (zh) | 一种正断层附近地层埋藏史的恢复方法 | |
CN103901478A (zh) | 一种井震信息联合确定储层沉积特征和分布的方法 | |
CN114966886B (zh) | 一种基于构型的致密气藏砂体连通性等级判定方法 | |
Handford et al. | Regional facies relationships and sequence stratigraphy of a super-giant reservoir (Arab-D Member), Saudi Arabia | |
CN107808068A (zh) | 一种用于致密砂岩气藏高产富集评价的方法 | |
CN115097525A (zh) | 一种基于多信息融合的致密气藏砂体构型分析方法 | |
CN114994789A (zh) | 基于动静态信息融合的致密气藏河道砂体连通性验证方法 | |
CN117035508B (zh) | 基于多因素的致密气藏河道砂体连通性评价方法 | |
Liu et al. | The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs | |
Lindsay et al. | Depositional facies, diagenesis, and reservoir character of Mississippian cyclic carbonates in the Mission Canyon Formation, Little Knife field, Williston basin, North Dakota | |
CN115577213B (zh) | 一种水侵方向预测方法 | |
CN111624678A (zh) | 一种利用同期河道发育高程来划分高精度层序地层的方法 | |
CN113433589A (zh) | 一种基于数理统计的风化壳储层底界面识别方法 | |
CN113775337A (zh) | 单砂体地质模型的建立方法、装置和计算机可读存储介质 | |
Watts et al. | Application of reservoir geology to enhanced oil recovery from Upper Devonian Nisku reefs, Alberta, Canada | |
Shi et al. | Case Stories of Horizontal Wells in Offshore Fluvial Oil Reservoir | |
Qian et al. | Numerical simulation study of water injection development in an extra-low-permeability fractured reservoir, Xiaoguai oilfield | |
Best et al. | Contribution of outcrop data to improve understanding of field performance: rock exposures at Eight Foot Rapids tied to the Aneth field | |
CN118228614B (zh) | 一种基于多代理模型的裂缝性气藏智能化提采方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |