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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005457 optimization Methods 0.000 title claims abstract description 17
- 230000000694 effects Effects 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 3
- 230000008092 positive effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 12
- 238000011156 evaluation Methods 0.000 abstract description 7
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000005056 compaction Methods 0.000 abstract description 3
- 229910000831 Steel Inorganic materials 0.000 abstract 1
- 239000010959 steel Substances 0.000 abstract 1
- 238000011161 development Methods 0.000 description 5
- 230000018109 developmental process Effects 0.000 description 5
- 238000005553 drilling Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000029052 metamorphosis Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- 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
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.
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)
| 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)
| 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 |
-
2021
- 2021-11-03 CN CN202111294308.8A patent/CN114016989B/en active Active
Patent Citations (1)
| 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)
| 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 | |
| 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 | |
| CN112443322A (en) | Hydrocarbon source rock logging evaluation method based on equivalent saturation | |
| 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 | |
| CN110410068B (en) | Logging method for determining stratum fracture pressure gradient | |
| 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 |
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 |