CN114016989B - Completion fracturing staged optimization method based on element logging interpretation data - Google Patents

Completion fracturing staged optimization method based on element logging interpretation data Download PDF

Info

Publication number
CN114016989B
CN114016989B CN202111294308.8A CN202111294308A CN114016989B CN 114016989 B CN114016989 B CN 114016989B CN 202111294308 A CN202111294308 A CN 202111294308A CN 114016989 B CN114016989 B CN 114016989B
Authority
CN
China
Prior art keywords
fracturing
parameters
max
parameter
index
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
Application number
CN202111294308.8A
Other languages
Chinese (zh)
Other versions
CN114016989A (en
Inventor
陈远林
聂林
林春容
田野
叶鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Moyan Technology Co ltd
Original Assignee
Chengdu Yiyetulan Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Yiyetulan Technology Co ltd filed Critical Chengdu Yiyetulan Technology Co ltd
Priority to CN202111294308.8A priority Critical patent/CN114016989B/en
Publication of CN114016989A publication Critical patent/CN114016989A/en
Application granted granted Critical
Publication of CN114016989B publication Critical patent/CN114016989B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25Methods for stimulating production
    • E21B43/26Methods for stimulating production by forming crevices or fractures
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Animal Husbandry (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • General Business, Economics & Management (AREA)
  • Agronomy & Crop Science (AREA)
  • Geophysics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a well completion fracturing staged optimization method based on element logging interpretation data, which comprises the following steps of: collecting interpretation data of element logging as a basic data set, wherein the basic data set comprises a plurality of parameters of shale gas well later fracturing; minimum assignment sigma using parameter i imin Sum maximum assignment sigma imax Calculating the relative fracturing index FR of each parameter i The method comprises the steps of carrying out a first treatment on the surface of the Calculating a comprehensive relative fracturing index FR of a plurality of parameters n : calculating the difference beta between the maximum value FRnmax 'of the input parameter FRn' and the input parameter FRn max : sequentially adding the difference beta max Dividing the steel into s grades, wherein the section length of each grade is Ls; sequentially searching a number Lsx which is between the minimum fracturing segment length Lmin and the maximum fracturing segment length Lmax and can be divided by the segment length Ls by using a remainder function: the fracturing index is segmented. The method is a well completion optimization technology based on element interpretation results and suitable for shale gas horizontal wells, can effectively improve the accuracy and effectiveness of horizontal well sectional design, and provides important support for oil testing geological design scheme and post-compaction evaluation.

Description

Completion fracturing staged optimization method based on element logging interpretation data
Technical Field
The invention relates to the technical field of shale gas development, in particular to a well completion fracturing staged optimization method based on element logging interpretation data.
Background
Over ten years, the exploration and development of domestic unconventional shale oil and gas resources are greatly advanced, a series of matched technologies for scale benefit development are formed, and a development mode of exploration and deployment integration, drilling and completion integration and production operation intensification is formed. Because the unconventional oil and gas reservoir has the characteristics of low holes and low permeability, the utilization of the horizontal well staged fracturing technology is an important means for realizing commercial development, wherein how to integrate multiple data sources according to logging or logging interpretation results is a key technical link in the implementation of overall staged optimization.
At present, after shale gas reservoirs are exploited in a large scale, in order to save cost, only part of well logging interpretation work is carried out on the same platform well, and how the well which is not subjected to logging interpretation further carries out horizontal well completion optimization design according to logging results is very important. In the conventional shale gas horizontal well completion design, the fracturing segmentation design is mainly based on logging interpretation results, so that the characteristics of heterogeneity and metamorphosis of horizontal segment reservoir parameters and completion parameters are difficult to consider more carefully, and the fracturing segmentation design is qualitative evaluation segmentation, so that the accuracy of later fracturing effect evaluation is affected.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a well completion fracturing staged optimization method with high precision based on element logging interpretation data.
In order to achieve the aim of the invention, the invention adopts the following technical scheme:
the well completion fracturing staged optimization method based on the element logging interpretation data comprises the following steps:
s1: collecting interpretation data of element logging as a basic data set, wherein the basic data set comprises a plurality of parameters of shale gas well later fracturing;
s2: minimum assignment sigma using parameter i imin Sum maximum assignment sigma imax Calculating the relative fracturing index FR of each parameter i
S3: calculating a relative fracturing index FR for each parameter in the base data set i Obtaining comprehensive relative fracturing index indexes FR of a plurality of parameters in the basic data set:
wherein n is the number of parameters in the basic data set;
s4: calculating a comprehensive relative fracturing index FR of a plurality of parameters by utilizing the comprehensive relative fracturing index FR n
S5: comprehensive relative fracturing index FR using weighting function n Smoothing to obtain input parameter FR n ';
S6: calculate input parameter FR n Maximum value FR of nmax ' and input parameter FR n ' difference beta max
β max =FR nmax '-FR n '
S7: sequentially adding the difference beta max Is divided into s grades, and the segment length of each grade is L s
S8: sequentially searching the minimum fracturing segment length L by using remainder function min To maximum fracture length L max Can be divided into sections of L s Number of integer divisors L sx
S9: judging whether the remainder is equal to 0, ifThe fracturing index is segmented into: (D) max -L ix )~(D max ),(D max -2L ix )~(D max -L ix ),···,(D max -sL ix )~(D max -(s-1)L ix ) The method comprises the steps of carrying out a first treatment on the surface of the Wherein D is max The bottom depth of the logging data;
s10: if it isReturning to step S8, the available segment length L is found again s Number of integer divisors L sx
Further, the plurality of shale gas well later fracturing parameters comprise positive parameters and negative parameters, wherein the positive parameters are parameters with positive effects on shale gas well later fracturing, and the larger the positive parameters are, the better the positive parameters are; the negative parameters are parameters with negative effects on the later fracturing of the shale gas well and serve as negative parameters; the better the negative going parameter.
Further, step S7 includes:
s71: sequentially adding the difference beta max The method comprises the steps of dividing the system into s grades, wherein the s grades sequentially comprise: first stageSecond stageS. level->
S72: calculating the length L of each grade segment s
Further, relative fracturing index FR i The calculation method comprises the following steps:
wherein W is i As the weight coefficient of parameter I, I i p2 is the percentile of 2% of parameter I i p98 is the percentile of 98% of parameter i.
The beneficial effects of the invention are as follows: the method is a well completion optimization technology based on element interpretation results and suitable for shale gas horizontal wells, can effectively improve the accuracy and effectiveness of horizontal well sectional design, and provides important support for oil testing geological design scheme and post-compaction evaluation. The invention can integrate data from multiple aspects such as drilling data, logging, ground logging, petrophysical and geomechanical interpretation, core, cuttings analysis, production data and the like, reduce the variation of reservoir heterogeneity in the section, and provide important support for oil testing geological design and post-pressure evaluation.
Drawings
FIG. 1 is a flow chart of a completion fracture staging optimization method based on elemental log interpretation data.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and all the inventions which make use of the inventive concept are protected by the spirit and scope of the present invention as defined and defined in the appended claims to those skilled in the art.
As shown in fig. 1, the well completion fracturing staged optimization method based on the element logging interpretation data of the scheme comprises the following steps:
s1: collecting interpretation data of element logging as a basic data set, wherein the basic data set comprises a plurality of parameters of shale gas well later fracturing; the shale gas well later fracturing parameters comprise positive parameters and negative parameters, wherein the positive parameters are parameters with positive effects on shale gas well later fracturing, and the larger the positive parameters are, the better the positive parameters are; such as porosity, TOC, total gas content, total hydrocarbon content, methane content, brittle mineral content, young's modulus, natural fracture density, etc. The negative parameters are parameters with negative effects on the later fracturing of the shale gas well and serve as negative parameters; the better the negative parameters are; for example: density, water saturation, formation fracture pressure, minimum level principal stress, poisson's ratio, etc.
S2: minimum assignment sigma using parameter i imin Sum maximum assignment sigma imax Calculating the relative fracturing index FR of each parameter i The method for calculating the relative fracturing index FRi comprises the following steps:
wherein W is i As the weight coefficient of parameter I, I i p2 is the percentile of 2% of parameter I i p98 is the percentile of 98% of parameter i.
S3: calculating a relative fracturing index FR for each parameter in the base data set i Obtaining comprehensive relative fracturing index indexes FR of a plurality of parameters in the basic data set:
wherein n is the number of parameters in the basic data set;
s4: calculating a comprehensive relative fracturing index FR of a plurality of parameters by utilizing the comprehensive relative fracturing index FR n
S5: comprehensive relative fracturing index FR using weighting function n Smoothing to obtain input parameter FR n ';
S6: calculate input parameter FR n Maximum value FR of nmax ' and input parameter FR n ' difference beta max
β max =FR nmax '-FR n '
S7: sequentially adding the difference beta max Is divided into s grades, and the segment length of each grade is L s The method comprises the steps of carrying out a first treatment on the surface of the Comprising the following steps:
s71: sequentially adding the difference beta max The method comprises the steps of dividing the system into s grades, wherein the s grades sequentially comprise: first stageSecond stageS. level->
S72: calculating the length L of each grade segment s
S8: sequentially searching the minimum fracturing segment length L by using remainder function min To maximum fracture length L max Can be divided into sections of L s Number of integer divisors L sx
S9: judging whether the remainder is equal to 0, ifThe fracturing index is segmented into: (D) max -L ix )~(D max ),(D max -2L ix )~(D max -L ix ),···,(D max -sL ix )~(D max -(s-1)L ix ) The method comprises the steps of carrying out a first treatment on the surface of the Wherein D is max The bottom depth of the logging data;
s10: if it isReturning to step S8, the available segment length L is found again s Number of integer divisors L sx
The method is a well completion optimization technology based on element interpretation results and suitable for shale gas horizontal wells, can effectively improve the accuracy and effectiveness of horizontal well sectional design, and provides important support for oil testing geological design scheme and post-compaction evaluation. The invention can integrate data from multiple aspects such as drilling data, logging, ground logging, petrophysical and geomechanical interpretation, core, cuttings analysis, production data and the like, reduce the variation of reservoir heterogeneity in the section, and provide important support for oil testing geological design and post-pressure evaluation.

Claims (4)

1. A well completion fracturing staged optimization method based on element logging interpretation data is characterized by comprising the following steps of:
s1: collecting interpretation data of element logging as a basic data set, wherein the basic data set comprises a plurality of parameters of shale gas well later fracturing;
s2: minimum assignment sigma using parameter i imin Sum maximum assignment sigma imax Calculating the relative fracturing index FR of each parameter i
S3: calculating a relative fracturing index FR for each parameter in the base data set i Obtaining comprehensive relative fracturing index indexes FR of a plurality of parameters in the basic data set:
wherein n is the number of parameters in the basic data set;
s4: calculating a comprehensive relative fracturing index FR of a plurality of parameters by utilizing the comprehensive relative fracturing index FR n
S5: comprehensive relative fracturing index FR using weighting function n Smoothing to obtain input parameter FR n ';
S6: calculate input parameter FR n Maximum value FR of nmax ' and input parameter FR n ' difference beta max
β max =FR nmax '-FR n '
S7: sequentially adding the difference beta max Is divided into s grades, and the segment length of each grade is L s
S8: sequentially searching the minimum fracturing segment length L by using remainder function min To maximum fracture length L max Can be divided into sections of L s Number of integer divisors L sx
S9: judging whether the remainder is equal to 0, ifThe fracturing index is segmented into: (D) max -L ix )~(D max ),(D max -2L ix )~(D max -L ix ),···,(D max -sL ix )~(D max -(s-1)L ix ) The method comprises the steps of carrying out a first treatment on the surface of the Wherein D is max The bottom depth of the logging data;
s10: if it isReturning to step S8, the available segment length L is found again s Number of integer divisors L sx
2. The staged optimization method for well completion and fracturing based on element logging interpretation data according to claim 1, wherein the plurality of shale gas well post-fracturing parameters comprise positive parameters and negative parameters, the positive parameters are parameters with positive effects on shale gas well post-fracturing, and the larger the positive parameters are, the better the positive parameters are; the negative parameters are parameters with negative effects on the later fracturing of the shale gas well and serve as negative parameters; the better the negative going parameter.
3. The method for staged optimization of completion fracturing based on elemental log interpretation data of claim 1, wherein step S7 comprises:
s71: sequentially adding the difference beta max The method comprises the steps of dividing the system into s grades, wherein the s grades sequentially comprise: first stageSecond stageS. level->
S72: calculating the length L of each grade segment s
4. The method of staged optimization of completion fracturing based on elemental log interpretation data of claim 1, wherein the relative fracturing index FR i The calculation method comprises the following steps:
wherein W is i As the weight coefficient of parameter I, I i p2 is the percentile of 2% of parameter I i p98 is the percentile of 98% of parameter i.
CN202111294308.8A 2021-11-03 2021-11-03 Completion fracturing staged optimization method based on element logging interpretation data Active CN114016989B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111294308.8A CN114016989B (en) 2021-11-03 2021-11-03 Completion fracturing staged optimization method based on element logging interpretation data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111294308.8A CN114016989B (en) 2021-11-03 2021-11-03 Completion fracturing staged optimization method based on element logging interpretation data

Publications (2)

Publication Number Publication Date
CN114016989A CN114016989A (en) 2022-02-08
CN114016989B true CN114016989B (en) 2024-02-20

Family

ID=80060297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111294308.8A Active CN114016989B (en) 2021-11-03 2021-11-03 Completion fracturing staged optimization method based on element logging interpretation data

Country Status (1)

Country Link
CN (1) CN114016989B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110163533A (en) * 2019-06-03 2019-08-23 西南石油大学 A kind of comprehensive compressibility Slope map of pixels of shale air cleft net

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10677036B2 (en) * 2018-03-29 2020-06-09 Baker Hughes, A Ge Company, Llc Integrated data driven platform for completion optimization and reservoir characterization

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110163533A (en) * 2019-06-03 2019-08-23 西南石油大学 A kind of comprehensive compressibility Slope map of pixels of shale air cleft net

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
钟光海 ; 陈丽清 ; 廖茂杰 ; 王广耀 ; 杨杨 ; 高翔 ; .页岩气储层品质测井综合评价.天然气工业.(02),全文. *

Also Published As

Publication number Publication date
CN114016989A (en) 2022-02-08

Similar Documents

Publication Publication Date Title
CN112502701B (en) Classification evaluation method for comprehensive geology-engineering of low-permeability reservoir
CN105426650A (en) Method for optimizing fracturing reconstruction well section of shale gas well by using logging information
CN105257252A (en) Method for optimizing shale gas horizontal well clustering perforation well section by using logging information
CN113312785B (en) Deep shale gas reservoir fracturing optimization method based on compressibility evaluation
CN103616724B (en) Unconventional compact and shale core lithology fine and accurate description method
CN106570262A (en) Reservoir configuration structure description method
CN112196513A (en) Longmaxi group shale gas well productivity prediction method based on horizontal well trajectory evaluation
CN116122801A (en) Shale oil horizontal well volume fracturing compressibility comprehensive evaluation method
CN111487691B (en) Method for quantitatively identifying lithology and sand body structure of tight sandstone reservoir
CN108763751A (en) The bearing calibration of total hydrocarbon in gas logging data in a kind of petroleum geology well logging
CN114016989B (en) Completion fracturing staged optimization method based on element logging interpretation data
CN112832738B (en) Determination method for strength of hydrocarbon accumulated in clastic rock and identification and evaluation method for dessert layer
CN111749688B (en) Method for predicting development position and direction of dominant seepage channel
CN111764895B (en) Logging evaluation method suitable for shale gas reservoir geological model
CN108241785A (en) A kind of fine characterizing method of Heterogeneous reservoir saturation field
CN110410068B (en) Logging method for determining stratum fracture pressure gradient
CN115434698B (en) Stratum calcareous cement content calculation method based on photoelectric absorption section index
CN116146176A (en) Geological-engineering factor-based quantitative classification method for deep coalbed methane reservoir logging
CN116522583A (en) Horizontal well stratum pressure prediction method based on linear regression method
CN115861551A (en) Digital well construction method for in-situ leaching uranium mining
CN114109374B (en) Shale gas reservoir target window position determination method and device
CN115099014A (en) Natural gas well geological exploration reserves calculation method based on logging while drilling
CN114488298A (en) Method for predicting microbial rock dessert reservoir
CN113514879A (en) 'ear' layer identification method
CN106909717B (en) Method for determining main oil production layer of lake-phase low-abundance and multi-production-layer compact oil field

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
TR01 Transfer of patent right

Effective date of registration: 20240527

Address after: No. 18 Yaobazi Street, Hanpo Township, Nanbu County, Nanchong City, Sichuan Province, 637300

Patentee after: Tian Ye

Country or region after: China

Address before: No. 2420, floor 24, building 5, No. 46, erxianqiao East Road, Chenghua District, Chengdu, Sichuan 610000

Patentee before: Chengdu yiyetulan Technology Co.,Ltd.

Country or region before: China

TR01 Transfer of patent right
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240618

Address after: Room 1230, 12th Floor, Unit 1, Building 9, No. 1700, Tianfu Avenue North Section, Chengdu High tech Zone, China (Sichuan) Pilot Free Trade Zone, Chengdu, Sichuan Province, 610000

Patentee after: Sichuan Moyan Technology Co.,Ltd.

Country or region after: China

Address before: No. 18 Yaobazi Street, Hanpo Township, Nanbu County, Nanchong City, Sichuan Province, 637300

Patentee before: Tian Ye

Country or region before: China

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240823

Address after: No. 2, Building 4, No. 76 Shuyuan Road, Wenlin Town, Renshou County, Meishan City, Sichuan Province 620000

Patentee after: Peng Mengqun

Country or region after: China

Patentee after: Zhang Guiliang

Address before: Room 1230, 12th Floor, Unit 1, Building 9, No. 1700, Tianfu Avenue North Section, Chengdu High tech Zone, China (Sichuan) Pilot Free Trade Zone, Chengdu, Sichuan Province, 610000

Patentee before: Sichuan Moyan Technology Co.,Ltd.

Country or region before: China

TR01 Transfer of patent right

Effective date of registration: 20240918

Address after: Room 1230, 12th Floor, Unit 1, Building 9, No. 1700, Tianfu Avenue North Section, Chengdu High tech Zone, China (Sichuan) Pilot Free Trade Zone, Chengdu, Sichuan Province, 610000

Patentee after: Sichuan Moyan Technology Co.,Ltd.

Country or region after: China

Address before: No. 2, Building 4, No. 76 Shuyuan Road, Wenlin Town, Renshou County, Meishan City, Sichuan Province 620000

Patentee before: Peng Mengqun

Country or region before: China

Patentee before: Zhang Guiliang