CN116927740A - Three-dimensional development integral fracturing sequence optimization design method - Google Patents
Three-dimensional development integral fracturing sequence optimization design method Download PDFInfo
- Publication number
- CN116927740A CN116927740A CN202310632853.6A CN202310632853A CN116927740A CN 116927740 A CN116927740 A CN 116927740A CN 202310632853 A CN202310632853 A CN 202310632853A CN 116927740 A CN116927740 A CN 116927740A
- Authority
- CN
- China
- Prior art keywords
- fracturing
- well
- model
- sequence
- scheme
- 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 46
- 238000013461 design Methods 0.000 title claims abstract description 40
- 238000011161 development Methods 0.000 title claims abstract description 27
- 238000005457 optimization Methods 0.000 title claims abstract description 15
- 238000004088 simulation Methods 0.000 claims abstract description 23
- 239000010410 layer Substances 0.000 claims description 91
- 238000010276 construction Methods 0.000 claims description 12
- 239000011229 interlayer Substances 0.000 claims description 11
- 239000011435 rock Substances 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 10
- 230000035699 permeability Effects 0.000 claims description 4
- 239000011148 porous material Substances 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 6
- 238000012795 verification Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 11
- 238000011156 evaluation Methods 0.000 description 7
- 238000004519 manufacturing process Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 5
- 239000004744 fabric Substances 0.000 description 5
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004826 seaming Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000001360 synchronised effect Effects 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
-
- 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
-
- 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)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Business, Economics & Management (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)
- Geophysics (AREA)
- Agronomy & Crop Science (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to the technical field of oil and gas field development, in particular to a three-dimensional development integral fracturing sequence optimization design method; the method comprises the following steps: s1, constructing a geological engineering integrated model; s2, designing a fracturing sequence scheme, wherein the fracturing sequence scheme comprises same-layer sequential fracturing, staggered zipper fracturing and edge-center zipper fracturing; s3, comprehensively evaluating the designed fracturing sequence scheme, and evaluating the aspects of fracture network simulation, pressure field simulation, SRV prediction and yield prediction by adopting a numerical simulation method; firstly, a geological engineering integrated model is constructed, then, a plurality of fracturing sequence schemes are designed according to specific actual needs, finally, a numerical simulation method is adopted, the optimal scheme is selected according to the finally formed seam net complexity, SRV and 10 years of output, and the optimal scheme obtained through verification of the simulation method can achieve maximum improvement of fracturing effect.
Description
Technical Field
The invention relates to the technical field of oil and gas field development, in particular to a three-dimensional development integral fracturing sequence optimization design method.
Background
The unconventional reservoir has the geological characteristics of compactness and hypotonic, has the characteristics of long longitudinal span of an oil layer and overlapping development of multiple sets of oil layers, is mainly developed by adopting a novel mode of three-dimensional, horizontal well, large platform, group-type and factory-type integral fracturing construction in order to improve the longitudinal utilization degree of the reservoir, but in the integral fracturing construction of the large platform, a high-stress area caused by a broken section is monitored to generate a stress barrier and unbalanced expansion of cracks generated by an adjacent section, which causes insufficient integral transformation of the reservoir, weak production stabilizing capability after compression and low three-dimensional utilization degree, so that the integral fracturing construction sequence becomes critical in the integral fracturing construction of the large platform.
The prior art CN114737943A discloses a method, a device and equipment for transforming an unconventional reservoir three-dimensional well pattern, wherein the transforming method comprises the following steps: (1) According to the three-dimensional geological model of the target block, a three-dimensional ground stress model and a rock mechanics model are established; and performing three-dimensional reconstruction feasibility evaluation on the target block; (2) If the three-dimensional transformation feasibility evaluation result is feasible, carrying out cloth well cloth seam design according to the three-dimensional ground stress model and the rock mechanical model, wherein the cloth well cloth seam design comprises the steps of using trisection staggered seams between clusters and using W-shaped three-dimensional staggered cloth wells; (3) According to the well pattern fracturing modification process, the well pattern fracturing modification process is executed, and comprises the following steps: determining a single well design parameter, the single well parameter comprising at least one of: the horizontal well section is long, the well distance of the horizontal well, the number of perforation holes and the perforation aperture; and determining production parameters including at least one of: limit displacement, fracturing fluid viscosity, and fracturing scale; and executing a well pattern fracturing modification process according to the design parameters of the single well, the production parameters and the well distribution and seam distribution design, but the modification method does not mention the fracturing sequence of the horizontal well and the optimization method of the fracturing sequence.
Therefore, it is needed to provide a three-dimensional development integral fracturing sequence optimization design method, which guides and optimizes the fracturing sequence during multi-well synchronous fracturing and maximally improves the fracturing effect compared with the prior art.
Disclosure of Invention
In order to solve the problems, the invention provides a three-dimensional development integral fracturing sequence optimization design method.
In order to achieve the above purpose, the technical scheme of the invention is as follows:
a three-dimensional development integral fracturing sequence optimization design method comprises the following steps:
s1, constructing a geological engineering integrated model;
s2, designing a fracturing sequence scheme, wherein the fracturing sequence scheme comprises same-layer sequential fracturing, staggered zipper fracturing and edge-center zipper fracturing;
s3, comprehensively evaluating the designed fracturing sequence scheme, and evaluating the aspects of fracture network simulation, platform reservoir transformation volume prediction and yield prediction by adopting a numerical simulation method.
Further, S3 specifically includes the following steps:
s301, calculating single well crack expansion balance and platform crack expansion balance, and eliminating a fracturing sequence scheme with uneven crack expansion;
s302, calculating the transformation volume of the reservoir of the platform by using a Slab volumetric method for the rest fracturing sequence scheme eliminated in the step S301, and eliminating part of the fracturing sequence scheme according to the transformation volume of the reservoir of the platform;
s303, predicting the yield of the rest fracturing sequence schemes eliminated in the step S302, calculating the accumulated oil yield of the platform, and screening out the fracturing sequence scheme with the largest accumulated oil yield of the platform as an optimal scheme.
Further, the single well crack propagation balance is determined by a single well crack propagation balance factor Z N Determining Z N Calculated by the following formula:
in the above, Z N A single well crack propagation equalization factor; n is the number of three-dimensional development wells; n is the number of cracks; x is x n Represents the length of the nth slit;the average length of n cracks is shown.
Further, the platform crack propagation equilibrium is determined by a platform crack propagation equilibrium factor Z, which is calculated by the following formula:
in the above, Z 1 、Z 2 、Z 3 …Z N The single well crack propagation balance factors of the first well, the second well and the third well … Nth well are respectively represented, and N is the number of three-dimensional development wells.
Furthermore, the well arrangement mode in the step S2 adopts a W-shaped three-dimensional staggered well arrangement mode.
Furthermore, in the same-layer sequential fracturing, the fracturing is carried out from one end of the area, the single-section well-by-well fracturing or multi-section well-by-well fracturing mode is adopted for the horizontal wells of the same layer, and the single-section interlayer fracturing or multi-section interlayer fracturing mode is adopted for the horizontal wells of different layers.
Furthermore, in the staggered zipper fracturing, the fracturing is performed from one end of the area, and a single-section well-by-well alternating fracturing mode or a multi-section well-by-well alternating fracturing mode is adopted.
Furthermore, in the edge-center zipper fracturing, the fracturing is performed simultaneously from two ends of the area, and a single-section well-by-well alternating fracturing mode or a multi-section well-by-well alternating fracturing mode is adopted.
Further, S1 specifically includes the following steps:
s101, collecting original data of a target block;
s102, constructing a three-dimensional geological model of the target block by adopting a modeling method combining certainty and randomness, wherein the three-dimensional geological model comprises a construction model, a permeability model, a porosity model and an oil saturation model.
S103, constructing a three-dimensional ground stress and rock mechanics model of the target block on the basis of the geological model; the ground stress model specifically comprises a vertical stress model, a pore pressure model, a horizontal minimum principal stress model and a horizontal maximum principal stress model; the rock mechanics model specifically includes a young's modulus model, a poisson's ratio model, a shear modulus model, and a bulk modulus model.
Compared with the prior art, the invention has the beneficial effects that:
(1) According to the invention, a geological engineering integrated model is constructed, multiple fracturing sequence schemes are designed according to specific actual needs, and finally, a numerical simulation method is adopted, the optimal scheme is selected based on the finally formed seaming net complexity, SRV and 10 years yield, and the optimal scheme obtained through verification of the simulation method can realize the maximum promotion of the fracturing effect.
(2) According to the invention, the designed multiple fracturing sequence schemes are simulated one by one, so that the calculation results of factors influencing the crack expansion balance and factors influencing the SRV are more accurate, the actual fracturing sequence schemes can be more fit, and the comprehensive evaluation is more objective.
Drawings
FIG. 1 is a flow chart of the present invention.
Fig. 2 is a flowchart of the step S1 of the present invention.
FIG. 3 is a schematic diagram of an integrated model of the geological engineering of the present invention.
FIG. 4 is a schematic diagram showing a well pattern according to the present invention.
FIG. 5 is a graph showing horizontal well distribution in the case of two layers in the same-layer sequential fracturing of the present invention.
FIG. 6 is a graph of horizontal well distribution for a three-layer single stage condition in the same-layer sequential fracturing of the present invention.
FIG. 7 is a graph showing the horizontal profile of a three-layer multi-segment in a same-layer sequential fracture according to the present invention.
FIG. 8 is a graph showing the horizontal well profile for two layers of the staggered zipper fracturing of the present invention.
FIG. 9 is a graph showing the horizontal well profile for a three-layer single stage case in a staggered zipper fracture of the present invention.
FIG. 10 is a graph showing the horizontal well profile for a three-layer multi-segment case in a staggered zipper fracking in accordance with the present invention.
FIG. 11 is a graph showing the profile of a horizontal well for the two-layer case in a frac edge-center zipper of the present invention.
FIG. 12 is a graph of horizontal well profile for a three-layer single stage case in an edge-center zipper fracture of the present invention.
FIG. 13 is a graph showing the horizontal well profile for a three-layer multi-segment case in an edge-center zipper fracture of the present invention.
Fig. 14 is a schematic diagram of a design of the present invention in the case of a two-layer, 3-well horizontal well.
Fig. 15 is a schematic diagram of a design of the present invention in the case of a two-layer, 4-well horizontal well.
FIG. 16 is a schematic diagram of the horizontal well distribution in a double-layer well layout according to the present invention.
Fig. 17 is a schematic diagram of the simulation results of the seam network complexity of the designs 1, 2, and 3 of the block a of the present invention.
Fig. 18 is a schematic diagram of the simulation results of the seam complexity of the designs 4 and 5 of the block a of the present invention.
Fig. 19 is a schematic diagram of the simulation results of the seam network complexity of the designs 6, 7, and 8 of the block a of the present invention.
Fig. 20 is a schematic diagram of the simulation results of the seam complexity of the designs 9 and 10 of the block a according to the present invention.
FIG. 21 is a bar graph of the measured data values after simulation in all designs of block A of the present invention.
Fig. 22 is a schematic diagram of the simulation results of the seam network complexity of designs 1, 2, and 3 of the block B of the present invention.
Fig. 23 is a schematic diagram showing the simulation results of the seam network complexity of the designs 4, 5, and 6 of the block B of the present invention.
FIG. 24 is a bar graph of the measured data values after simulation in all the designs of the block B of the present invention.
Detailed Description
The technical solutions of the present invention will be clearly described below with reference to the accompanying drawings, and it is obvious that the described embodiments are not all embodiments of the present invention, and all other embodiments obtained by a person skilled in the art without making any inventive effort are within the scope of protection of the present invention.
As shown in fig. 1, the invention provides a three-dimensional development integral fracturing sequence optimization design method, which comprises the following steps:
s1, constructing a geological engineering integrated model, and continuously drilling and completing a well along with a target area, further encrypting the existing geological engineering integrated model according to real-time logging data to achieve the purpose of high-definition modeling;
as shown in fig. 2, S1 specifically includes the following steps:
s101, collecting original data of a target block, wherein the original data comprise oil reservoir construction information, original seismic data and a seismic reference plane; horizon and fault interpretation schemes and reservoir mechanics inversion results; the block is known as well position coordinates, geological stratification and conventional logging data and oil test and production test data; block known well imaging log data, rock mechanics experimental data, microseism data; existing fracturing construction curves and fracturing design data; the data of the related research results of the research area and the adjacent area are collected as much as possible, and the data are not limited to the above data;
s102, constructing a three-dimensional geological model of the target block by adopting a modeling method combining certainty and randomness, wherein the three-dimensional geological model comprises a construction model, a permeability model, a porosity model and an oil saturation model, and specifically obtaining the model by adopting fault modeling, layered quality control, construction modeling, grid design, lithofacies simulation and attribute simulation.
The model size is about 2800m×1800m×70m, and the grid precision is 20m×20m; the target area has a relatively gentle structure and a structural drop of about 70 m; building four layers of construction models, wherein the development continuity of 2 and 3 layers of oil layers is good; in the three-dimensional geological model, the average porosity is 8.6%, the average permeability is 1.84mD, and the average oil saturation is 67%;
s103, constructing a three-dimensional ground stress and rock mechanics model of the target block on the basis of the geological model; the ground stress model specifically comprises a vertical stress model, a pore pressure model, a horizontal minimum principal stress model and a horizontal maximum principal stress model; the rock mechanics model specifically includes a young's modulus model, a poisson's ratio model, a shear modulus model, and a bulk modulus model.
And constructing a three-dimensional ground stress and rock mechanics model (shown in figure 3) of the target block by using geological engineering integrated software, wherein the grid size of the stress model is the same as that of the three-dimensional geological model, in the three-dimensional ground stress model, the average Young modulus model is 20.05GPa, the average Poisson's ratio model is 0.28, the average pore pressure is 54.44MPa, the average minimum horizontal main stress model is 71.11MPa, the average maximum horizontal main stress model is 82.75MPa, and the pressure of an overburden layer is 78.47MPa.
S2, designing a fracturing sequence scheme, wherein the well layout principle comprises the following points: (1) By combining with the characteristics of the actual bottom layer of the field, adopting the W-shaped stereoscopic staggered well arrangement (shown in figure 4) and adopting the W-shaped stereoscopic staggered well arrangement, the purposes of avoiding vertical adverse interference among horizontal wells, strengthening vertical powerful interference, achieving integral fracturing and enhancing the complexity of a seam net can be achieved; (2) Respectively designing according to the different numbers of longitudinal layers and horizontal wells in the region, wherein the number of layers is designed into multiple layers, and the specific number of layers and the number of horizontal wells in each layer are set according to actual needs; (3) Follow the principle that is favorable for site construction and high-efficiency operation; the fracturing design scheme specifically comprises the following categories:
(1) The method for sequentially fracturing the same layer comprises the following steps:
the first step is to perform fracturing of the first layer of horizontal wells, and single-section well-by-well fracturing or multi-section well-by-well fracturing is performed according to the sequence from left to right (or right to left).
And secondly, carrying out single-stage interlayer fracturing or multi-stage interlayer fracturing from the last fractured horizontal well of the first layer to the first fractured horizontal well of the next adjacent layer.
And thirdly, carrying out fracturing of the next adjacent layer of horizontal well, and carrying out fracturing of the second layer of horizontal well by adopting a single-section well-by-well fracturing or multi-section well-by-well fracturing mode sequentially from the second layer of horizontal well.
And fourthly, circulating the methods of the second step and the third step until the fracturing of the horizontal wells in all layers is completed.
The first layer can be the nearest layer from the ground or the farthest layer from the ground, so that the same-layer sequential fracturing can be performed from top to bottom or from bottom to top.
The following details will be described by taking two layers and three layers as examples:
in the two-layer situation (shown in fig. 5), when the two-layer situation is up to down, from well 1, single-section well-by-well fracturing is firstly adopted for horizontal wells (or described as wells) 1-4 (1- (K-1)), single-section interlayer fracturing is adopted for wells 5-6 (K- (K+1)), single-section well-by-well fracturing is adopted for wells 7-9 ((K+2) -N), and when the two-layer situation is down to up, the two-layer situation is achieved through inversion.
Under the three-layer condition (shown in fig. 6), starting from well 1, carrying out single-section well-by-well fracturing from well 1-2 (1- (K-1)), carrying out single-section interlayer fracturing from well 3-4 (K- (K+1)), carrying out single-section well-by-well fracturing from well 5-6 ((K+2) - (N-1)), carrying out single-section interlayer fracturing from well 7-8 (N- (N+1)), carrying out single-section well-by-well fracturing from well 9-10 ((N+2) -Y), and reversing when the sequence is down.
Under the condition of three layers and multiple sections (shown in fig. 7), starting from a No. 1 well, performing multi-section well-by-well fracturing from a No. 1-2 (1- (K-1)) well, performing multi-section inter-layer fracturing from a No. 3-4 (K- (K+1)) well, performing multi-section well-by-well fracturing from a No. 5-6 ((K+2) - (N-1)) well, performing multi-section inter-layer fracturing from a No. 7-8 (N- (N+1)) well, performing multi-section well-by-well fracturing from a No. 9-10 ((N+2) -Y) well, and reversing when the sequence is down.
Wherein K represents the number of the first layer of horizontal wells, N-K represents the number of the second layer of horizontal wells, and Y-N represents the number of the third layer of horizontal wells; and the multi-stage fracturing operation is controlled to press three stages at most once, and the following steps are the same.
(2) The staggered zipper fracturing is performed by adopting single-section well-by-well alternate fracturing or multi-section well-by-well alternate fracturing among the multi-layer horizontal wells, and the two-layer and three-layer are taken as examples for the detailed explanation below.
In the case of two-layer single-stage (as shown in fig. 8), starting from well 1, and performing single-stage well-by-well alternate fracturing on the upper layer and the lower layer according to the sequence of well 1 to well 6, well 2 and the like (1 ～ 6 ～ 2 ～ 7 ～ 3 ～ 8 ～ 4 ～ 9 ～ 5-K), and reversing the steps when the sequence is from top to bottom.
In the case of three-layer single-stage (as shown in fig. 9), starting from well 1, and performing single-stage well-by-well alternate fracturing on the upper, middle and lower layers according to the sequence of well 1 to well 4, well 8 and the like (1 ～ 4 ～ 8 ～ 5 ～ 2 ～ 6 ～ 9 ～ 10 ～ 7 ～ 3), and reversing the steps when the sequence is from top to bottom.
In the case of three layers and multiple sections (as shown in fig. 10), starting from well 1, and performing multi-section well-by-well alternate fracturing (controlling at most one time for three sections) on the upper, middle and lower layers according to the sequence of well 1 to well 4 to well 8 and the like (1 ～ 4 ～ 8 ～ 5 ～ 2 ～ 6 ～ 9 ～ 10 ～ 7 ～ 3), and reversing the steps when the sequence is from top to bottom.
(3) The edge-center zipper fracturing is carried out by adopting a single-section well-by-well alternative fracturing or multi-section well-by-well alternative fracturing mode from the left side and the right side, and the two-layer and three-layer are taken as examples for the detailed description below.
In the two-layer single-stage case (as shown in fig. 11), in the top-down order, the left side edge starts from well No. 1, while the right side edge starts from well No. 5, fracturing proceeds in the order in the second class, the left side is subjected to 1 ～ 6 ～ 2 ～ 7 ～ 3 sequence, the right side is subjected to 5 ～ 9 ～ 4 ～ 8 ～ 3 sequence, single-section well-by-well alternating fracturing is performed, and the steps are reversed when the sequence is from bottom to top.
In the case of three-layer single-stage (as shown in fig. 12), when the sequence is up to down, the left side edge part starts from the No. 1 well, the right side edge part starts from the No. 3 well, the fracturing is performed according to the sequence in the second class, the left side is in the sequence of 1 ～ 4 ～ 8 ～ 5 ～ 2, the right side is subjected to 3 ～ 7 ～ 10 ～ 9 ～ 6-2, the single-stage well-by-well alternating fracturing is performed according to the sequence, and when the sequence is down to up, the steps are reversed.
In the case of three-layer multi-section (as shown in fig. 13), when the sequence is up to down, the left side edge part starts from the No. 1 well, the right side edge part starts from the No. 3 well, the fracturing is performed according to the sequence in the second class, the left side is in the sequence of 1 ～ 4 ～ 8 ～ 5 ～ 2, the right side is subjected to 3 ～ 7 ～ 10 ～ 9 ～ 6-2, the multi-section well-by-well alternating fracturing is performed (the maximum number of the three sections are pressed for one time) and when the sequence is down to up, the steps are reversed.
Taking the design of two horizontal wells as an example, as shown in fig. 14, when the number of block wells is 3, 6 schemes can be designed: 1 staggered zipper fracturing, 2 edge-center fracturing, 3 same-layer sequential fracturing, 4 multi-section (five-section) staggered zipper fracturing, 5 multi-section (five-section) edge-center fracturing and 6 multi-section (five-section) same-layer sequential fracturing.
As shown in fig. 15, when the number of block wells=4 ports, 10 schemes are designed (fig. 13): 1 staggered zipper fracturing, 2 edge-center fracturing, 3 upper and lower same-layer sequential fracturing, 4 center-edge fracturing, 5 lower and upper same-layer sequential fracturing, 6 (five-segment) staggered zipper fracturing, 7 (five-segment) edge-center fracturing, 8 (five-segment) upper and lower same-layer sequential fracturing, 9 (five-segment) center-edge fracturing and 10 lower and upper same-layer sequential fracturing.
When the number of wells in the block is 5, 6 and 7 …, the design rule is the same as the design rule, and the design of the multi-layer horizontal well is the same as the design rule of two layers.
S3, comprehensively judging the various fracturing schemes designed in the S2, and carrying out scheme optimization by adopting a numerical model method based on the finally formed seam network complexity, the platform reservoir transformation volume SRV, the stress field change condition and the 10-year yield, wherein the judging steps are as follows:
s301, other conditions are the same, only the fracturing sequence is changed, single well crack expansion balance and platform crack expansion balance are calculated, the smaller the balance factor (Z) is, the better the balance factor (Z) is, and a fracturing sequence scheme with uneven 1/3 crack expansion is eliminated;
single well fracture propagation equalization factor Z N (N is the number of stereo development wells), and the block crack expansion balance factor is Z, Z N The method is obtained by the following formula:
in the above, Z N Equalization factors (dimensionless) for single well fracture propagation; n is the number of stereo development wells (ports); n is the number of cracks (bars); x is x n Represents the length (m) of the nth slit;the average length (m) of n cracks is shown.
The block crack propagation equalization factor Z is calculated by the following formula:
Z N and the evaluation result table of Z refers to the following table 1:
table 1 table for evaluating seam net equalization factor
| Z N | Z |
| Not less than 2 (very poor) | Not less than 2 (very poor) |
| 1-2 (Difference) | 1-2 (Difference) |
| 0.5-1 (pass) | 0.5-1 (pass) |
| 0.1-0.5 (good) | 0.1-0.5 (good) |
| < 0.1 (Excellent) | < 0.1 (Excellent) |
S302, calculating a reservoir reconstruction volume SRV of a platform by using a Slab volumetric method for the rest fracturing sequence schemes eliminated in the step S301, wherein the larger the SRV value is, the better the fracturing sequence schemes with the smaller SRV values of 1/2 number are eliminated;
s303, predicting the yield of the rest fracturing sequence schemes subjected to the step S302 for 10 years, calculating the total oil yield of the single well oil yield EUR and the platform oil yield EUR, wherein the total oil yield of the platform oil yield EUR is equal to the total oil yield of all the single wells of the platform, the bigger the total oil yield of the platform is, the better the bigger the total oil yield of the platform is, and the fracturing sequence scheme with the largest total oil yield of the platform is selected as the optimal scheme.
The selection of the best scheme is carried out by the specific number of layers and the number of horizontal wells:
as shown in fig. 16, the method in S2 is adopted to perform W-type double-layer well arrangement, wherein 5 wells are arranged on the upper layer, 3 wells are arranged on the lower layer, the black well is the old well of the area, and the block well is divided into A, B small blocks for operation.
(1) For the a block 4 well:
firstly, designing 10 medium-pressure fracturing schemes for the A block 4 well according to the method in the step S2, wherein the schemes are respectively 1 staggered zipper fracturing, 2 edge-center fracturing, 3 upper and lower same-layer sequential fracturing, 4 center-edge fracturing, 5 lower and upper same-layer sequential fracturing, 6 multi-section staggered zipper fracturing, 7 multi-section edge-center fracturing, 8 multi-section upper and lower same-layer sequential fracturing, 9 multi-section center-edge fracturing and 10 multi-section lower and upper same-layer sequential fracturing.
And secondly, respectively and comprehensively judging 10 design schemes in the first step according to the method in the step S3, wherein the simulation result of the seaming complexity of each scheme is shown in figures 17, 18, 19 and 20.
As shown in fig. 17, the lower upper slotted network of scheme 1 is limited in expansion, the laminated slotted networks of scheme 2 are relatively good in expansion, and the lower well slotted network of scheme 3 is limited in expansion; therefore, scheme 2> scheme 1> scheme 3; as shown in fig. 18, the lower well fracture network of scheme 4 is limited in expansion, and the upper layer fracture of scheme 5 is too long in expansion, so that interference is easy to form; overall schemes 4, 5 do not have good crack propagation; as shown in fig. 19, the middle seam network expansion is limited in scheme 6; scheme 7 each laminated fracture network expands relatively slightly better, but has no single stage fracturing effect well; scheme 8 upper well screen expansion is poor; as shown in fig. 20, the lower well fracture network of scheme 9 is limited in expansion, and the middle well fracture network of scheme 10 is limited in expansion, and the effect of multi-stage fracturing is lower than that of single-stage fracturing as a whole.
Third, through a calculation formula
And calculate Z in each scheme separately in combination with the values in FIG. 21 N And Z, the results are shown in the following table:
table 3 Block A4 well 10 schemes stitch Net evaluation factorsZ N Is the case of (2)
Fourth, SRV was calculated and the results are shown in the following table:
table 4 case table of block a remaining 6 schemes SRV
| Scheme for the production of a semiconductor device | Reservoir retrofit volume SRV (10) 4 m 3 ) |
| Scheme 2 | 490 |
| Scheme 4 | 350 |
| Scheme 7 | 430 |
| Scheme 8 | 440 |
| Scheme 9 | 290 |
| Scheme 10 | 310 |
Fifth, screening scheme 2, scheme 7 and scheme 8 from the above table, and counting the accumulated oil yield EUR of these four schemes, the results are shown in the following table:
TABLE 5EUR Condition Table
| Series of | EUR(t) | Alternatives to |
| Scheme 2 | 110289 | Selecting and using |
| Scheme 7 | 104024 | Obsolete |
| Scheme 8 | 100890 | Obsolete |
The finally determined scheme is as follows: scheme 2.
(2) 3 wells for B block
Step one, designing a fracturing scheme in step 6 according to the method in the step S2, wherein the fracturing scheme is as follows: 1 staggered zipper fracturing, 2 edge-center fracturing, 3 same-layer sequential fracturing, 4 multi-section staggered zipper fracturing, 5 multi-section edge-center fracturing and 6 multi-section same-layer sequential fracturing.
And secondly, respectively carrying out comprehensive judgment on the 6 schemes according to the method in S3, wherein the simulation result of the seaming complexity degree of each scheme is shown in figures 22 and 23.
Third, through a calculation formula
And combining the values in FIG. 24 to calculate Z in each scheme N And Z, the results are shown in the following table:
table 6 block B3 well 6 schemes stitch net evaluation factor Z N Condition table of (2)
Fourth, SRV was calculated and the results are shown in the following table:
table 7 case table of block B remaining 4 schemes SRV
| Scheme for the production of a semiconductor device | Reservoir retrofit volume SRV (10) 4 m 3 ) |
| Scheme 1 | 454 |
| Scheme 2 | 480 |
| Scheme 3 | 440 |
| Scheme 5 | 380 |
Fifthly, screening scheme 1 and scheme 2 from the above table, and counting the accumulated oil yield EUR of the three schemes, wherein the result is shown in the following table:
TABLE 8EUR Condition Table
| Series of | EUR(t) | Alternatives to |
| Scheme 1 | 82563 | Selecting and using |
| Scheme 2 | 82662 | Selecting and using |
The finally determined scheme is as follows: scheme 1 or scheme 2.
And combining the two condition evaluation results, combining the three factors of expanded equalization of a seam net, transformation of volume SRV and ten-year yield prediction, recommending to adopt the fracturing sequence of combining the edge and the center with staggered zippers, so that the production effect can be improved, and the platform working efficiency can be effectively improved.
The above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to examples, it should be understood by those skilled in the art that modifications and equivalents may be made thereto without departing from the scope of the technical solution of the present invention, which is intended to be covered by the claims of the present invention.
Claims (9)
1. The three-dimensional development integral fracturing sequence optimization design method is characterized by comprising the following steps of:
s1, constructing a geological engineering integrated model;
s2, designing a fracturing sequence scheme, wherein the fracturing sequence scheme comprises same-layer sequential fracturing, staggered zipper fracturing and edge-center zipper fracturing;
s3, comprehensively evaluating the designed fracturing sequence scheme, and evaluating the aspects of fracture network simulation, platform reservoir transformation volume prediction and yield prediction by adopting a numerical simulation method.
2. The method for three-dimensional development of integral fracturing sequence optimization design according to claim 1, wherein the step S3 specifically comprises the following steps:
s301, calculating single well crack expansion balance and platform crack expansion balance, and eliminating a fracturing sequence scheme with uneven crack expansion;
s302, calculating the transformation volume of the reservoir of the platform by using a Slab volumetric method for the rest fracturing sequence scheme eliminated in the step S301, and eliminating part of the fracturing sequence scheme according to the transformation volume of the reservoir of the platform;
s303, predicting the yield of the rest fracturing sequence schemes eliminated in the step S302, calculating the accumulated oil yield of the platform, and screening out the fracturing sequence scheme with the largest accumulated oil yield of the platform as an optimal scheme.
3. The three-dimensional development integral fracturing sequence optimization design method according to claim 2, wherein the single-well crack expansion balance is formed by a single-well crack expansion balance factor Z N Determining Z N Calculated by the following formula:
in the above, Z N A single well crack propagation equalization factor; n is the number of three-dimensional development wells; n is the number of cracks; x is x n Represents the length of the nth slit;the average length of n cracks is shown.
4. The method for optimizing and designing a three-dimensional development overall fracturing sequence according to claim 3, wherein the platform crack propagation balance is determined by a platform crack propagation balance factor Z, and Z is calculated by the following formula:
in the above, Z 1 、Z 2 、Z 3 …Z N The single well crack propagation balance factors of the first well, the second well and the third well … Nth well are respectively represented, and N is the number of three-dimensional development wells.
5. The method for three-dimensional development of integral fracturing sequence optimization design according to claim 1, wherein the well arrangement mode in the step S2 adopts 'W-shaped' three-dimensional staggered well arrangement.
6. The method for optimizing the design of the three-dimensional development integral fracturing sequence according to claim 5, wherein in the same-layer sequential fracturing, fracturing is carried out from one end of a region, single-stage well-by-well fracturing or multi-stage well-by-well fracturing is adopted for horizontal wells of the same layer, and single-stage interlayer fracturing or multi-stage interlayer fracturing is adopted for horizontal wells of different layers.
7. The method for optimizing the design of the three-dimensional development integral fracturing sequence according to claim 5, wherein in the staggered zipper fracturing, the fracturing is performed from one end of a region, and a single-section well-by-well alternating fracturing mode or a multi-section well-by-well alternating fracturing mode is adopted.
8. The method for optimizing the design of the three-dimensional development integral fracturing sequence according to claim 5, wherein in the edge-center zipper fracturing, fracturing is performed simultaneously from two ends of a region, and a single-section well-by-well alternating fracturing mode or a multi-section well-by-well alternating fracturing mode is adopted.
9. The method for three-dimensional development of integral fracturing sequence optimization design according to claim 1, wherein the step S1 specifically comprises the following steps:
s101, collecting original data of a target block;
s102, constructing a three-dimensional geological model of a target block by adopting a modeling method combining certainty and randomness, wherein the three-dimensional geological model comprises a construction model, a permeability model, a porosity model and an oil saturation model;
s103, constructing a three-dimensional ground stress and rock mechanics model of the target block on the basis of the geological model; the ground stress model specifically comprises a vertical stress model, a pore pressure model, a horizontal minimum principal stress model and a horizontal maximum principal stress model; the rock mechanics model specifically includes a young's modulus model, a poisson's ratio model, a shear modulus model, and a bulk modulus model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310632853.6A CN116927740B (en) | 2023-05-31 | 2023-05-31 | Three-dimensional development integral fracturing sequence optimization design method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310632853.6A CN116927740B (en) | 2023-05-31 | 2023-05-31 | Three-dimensional development integral fracturing sequence optimization design method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116927740A true CN116927740A (en) | 2023-10-24 |
| CN116927740B CN116927740B (en) | 2024-06-18 |
Family
ID=88374504
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310632853.6A Active CN116927740B (en) | 2023-05-31 | 2023-05-31 | Three-dimensional development integral fracturing sequence optimization design method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116927740B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117951428A (en) * | 2024-03-26 | 2024-04-30 | 中国石油大学(华东) | Geothermal reservoir fracturing parameter optimization method based on injection well fracturing and heating effect integrated evaluation |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104278980A (en) * | 2013-07-09 | 2015-01-14 | 中国石油天然气股份有限公司 | Method for optimizing compact oil horizontal well seam network parameters by adopting improved volume |
| US20180149000A1 (en) * | 2016-11-29 | 2018-05-31 | Conocophillips Company | Engineered stress state with multi-well completions |
| CN114357766A (en) * | 2022-01-04 | 2022-04-15 | 中国石油大学(北京) | Long straight well section well pattern integral volume fracturing optimization design method |
| CN114511139A (en) * | 2022-01-17 | 2022-05-17 | 中国石油大学(北京) | Sequential optimization method and system for multistage fracturing of shale gas horizontal well |
| CN114542038A (en) * | 2020-11-26 | 2022-05-27 | 中国石油天然气股份有限公司 | Fracturing method for natural fracture development section of horizontal well |
| CN114737943A (en) * | 2022-03-08 | 2022-07-12 | 中国石油大学(北京) | Unconventional reservoir three-dimensional well pattern transformation method, device, equipment and medium |
| GB202217776D0 (en) * | 2022-01-20 | 2023-01-11 | Univ Southwest Petroleum | Horizontal wellbore trajectory optimization method for deep shale gas fracturing to maximize srv |
| CN116066048A (en) * | 2021-11-03 | 2023-05-05 | 中国石油天然气股份有限公司 | Method, device, equipment and medium for well group linkage fracturing of stacked oil reservoir horizontal wells |
-
2023
- 2023-05-31 CN CN202310632853.6A patent/CN116927740B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104278980A (en) * | 2013-07-09 | 2015-01-14 | 中国石油天然气股份有限公司 | Method for optimizing compact oil horizontal well seam network parameters by adopting improved volume |
| US20180149000A1 (en) * | 2016-11-29 | 2018-05-31 | Conocophillips Company | Engineered stress state with multi-well completions |
| CN114542038A (en) * | 2020-11-26 | 2022-05-27 | 中国石油天然气股份有限公司 | Fracturing method for natural fracture development section of horizontal well |
| CN116066048A (en) * | 2021-11-03 | 2023-05-05 | 中国石油天然气股份有限公司 | Method, device, equipment and medium for well group linkage fracturing of stacked oil reservoir horizontal wells |
| CN114357766A (en) * | 2022-01-04 | 2022-04-15 | 中国石油大学(北京) | Long straight well section well pattern integral volume fracturing optimization design method |
| CN114511139A (en) * | 2022-01-17 | 2022-05-17 | 中国石油大学(北京) | Sequential optimization method and system for multistage fracturing of shale gas horizontal well |
| GB202217776D0 (en) * | 2022-01-20 | 2023-01-11 | Univ Southwest Petroleum | Horizontal wellbore trajectory optimization method for deep shale gas fracturing to maximize srv |
| CN114737943A (en) * | 2022-03-08 | 2022-07-12 | 中国石油大学(北京) | Unconventional reservoir three-dimensional well pattern transformation method, device, equipment and medium |
Non-Patent Citations (4)
| Title |
|---|
| RIPUDAMAN MANCHANDA 等: "Fracture sequencing in Multi-Well Pads:Impact of Staggering and Lagging Stages in Zipper Fracturing on Well Productivity", 《SPE》, pages 1 - 20 * |
| 冯其红;李东杰;时贤;王森;徐世乾;秦勇;安杰;: "基于扩展有限元的水平井改进拉链式压裂数值模拟", 中国石油大学学报(自然科学版), no. 02 * |
| 刘洪, 廖如刚, 李小斌, 胡昌权, 肖晖, 黄园园, 张潇宇: "页岩气"井工厂"不同压裂模式下裂缝复杂程度研究", 《天然气工业》, pages 70 - 76 * |
| 张进科, 苟利鹏, 吴文瑞, 杨金峰: "水力压裂缝间距及压裂顺序对裂缝扩展影响研究", 《地下空间与工程学报》, pages 603 - 609 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117951428A (en) * | 2024-03-26 | 2024-04-30 | 中国石油大学(华东) | Geothermal reservoir fracturing parameter optimization method based on injection well fracturing and heating effect integrated evaluation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116927740B (en) | 2024-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105201479B (en) | A kind of horizontal well on shale reservoir stratum staged fracturing perforation cluster method for optimally designing parameters | |
| CN111322050B (en) | Shale horizontal well section internal osculating temporary plugging fracturing construction optimization method | |
| CN110469303B (en) | Volume fracturing parameter optimization design method based on four types of transformation volumes | |
| US20160376885A1 (en) | Method and Apparatus for Performance Prediction of Multi-Layered Oil Reservoirs | |
| CN116927740B (en) | Three-dimensional development integral fracturing sequence optimization design method | |
| CN104750896B (en) | A kind of fractured-cavernous carbonate reservoir method for numerical simulation | |
| CN104278980B (en) | Method for optimizing compact oil horizontal well seam network parameters by adopting improved volume | |
| CN105201484A (en) | Vertical well separate layer fracturing interval optimization and construction parameter optimization designing method | |
| CN111456709B (en) | Horizontal well multistage fracturing segmentation clustering method based on logging curve | |
| CN110838175B (en) | Geological model building method for gas injection development oil reservoir | |
| CN114737943B (en) | Unconventional reservoir three-dimensional well pattern reconstruction method, device, equipment and medium | |
| CN106991236A (en) | A kind of refracturing well and stratum selection method based on four-dimensional crustal stress dynamic change | |
| CN105735961B (en) | A kind of hypotonic compact reservoir separate stratum fracfturing series method for optimizing | |
| CN108280275A (en) | A kind of high prediction technique of tight sand hydraulic fracturing seam | |
| RU2672292C1 (en) | Method for oil deposit development with horizontal wells with multi-stage fracing | |
| CN114427432B (en) | Method for determining development potential of residual gas in gas reservoir | |
| CN116127675A (en) | Prediction method for maximum recoverable reserve of shale oil horizontal well volume fracturing | |
| CN114647957B (en) | Acid fracturing effect evaluation method based on acid fracturing 'transformation area' theory | |
| CN114357766A (en) | Long straight well section well pattern integral volume fracturing optimization design method | |
| CN115705454A (en) | Crack propagation simulation fracturing design optimization method based on phase field method | |
| CN113803041B (en) | Volumetric fracturing design method and system | |
| CN108229713B (en) | Optimization design method for multi-layer commingled production scheme of fault block oil reservoir | |
| CN108829994A (en) | Multilayer shale pseudo-three-dimensional fracture modeling stitches high method for solving | |
| CN104535391A (en) | Physical geography data processing method based on layered geography models | |
| CN114444340A (en) | Evaluation method for fracture distribution rule of mudstone reservoir fracturing |
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 |