CN111722277B - Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration - Google Patents
Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration Download PDFInfo
- Publication number
- CN111722277B CN111722277B CN201910222303.0A CN201910222303A CN111722277B CN 111722277 B CN111722277 B CN 111722277B CN 201910222303 A CN201910222303 A CN 201910222303A CN 111722277 B CN111722277 B CN 111722277B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- solution
- type
- threshold value
- trap
- 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 39
- 238000011161 development Methods 0.000 claims abstract description 19
- 230000018109 developmental process Effects 0.000 claims abstract description 19
- 238000011160 research Methods 0.000 claims abstract description 16
- 230000011218 segmentation Effects 0.000 claims abstract description 11
- 238000012795 verification Methods 0.000 claims abstract description 11
- 230000002452 interceptive effect Effects 0.000 claims abstract description 10
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 238000005553 drilling Methods 0.000 claims description 18
- 230000004044 response Effects 0.000 claims description 14
- 239000011324 bead Substances 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000007476 Maximum Likelihood Methods 0.000 claims description 7
- 230000001427 coherent effect Effects 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000005211 surface analysis Methods 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 abstract description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 239000011435 rock Substances 0.000 description 8
- 239000002904 solvent Substances 0.000 description 7
- 239000012530 fluid Substances 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 4
- 230000000903 blocking effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 229930195733 hydrocarbon Natural products 0.000 description 4
- 150000002430 hydrocarbons Chemical class 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001125 extrusion Methods 0.000 description 2
- 239000013505 freshwater Substances 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 241001481833 Coryphaena hippurus Species 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a solution-breaking trap description method for solution-breaking reservoir exploration, which comprises the following steps: developing forward simulation research aiming at the structural style of the solution breaking body, establishing an earthquake recognition mode of the solution breaking body, and determining the reservoir type in the current solution breaking body; analyzing by combining a logging curve and the emptying loss characteristic of a drilled well in an area or an adjacent area where the current dissolved aggregate is located, and determining a reservoir porosity threshold value and a reservoir threshold value of the current dissolved aggregate; determining the boundary of the current solution-breaking trap after interactive verification in a three-dimensional space; and obtaining the reservoir volume, the reservoir effective volume and the effective reservoir thickness of the current solution. The invention takes the separability of the fracture zone and the classification of the fracture-hole-seam aggregate reservoir into consideration, and enables the description of the fracture solution trapping to be refined. By means of fracture zone segmentation description, classified reservoir prediction and volume estimation and combined with three-dimensional visual image display, the transverse segmentation, longitudinal development depth and internal structure difference of the solution-breaking trap can be described in space.
Description
Technical Field
The invention relates to the technical field of exploration, in particular to a solution breaking trap description method and device for solution breaking reservoir exploration.
Background
The fractured-dissolved body is a carbonate rock fracture-cave type reservoir body which develops along a fractured zone and is formed by the action of various factors such as the fractured zone, atmospheric fresh water, deep heat flow and the like as the name suggests.
In recent years, a series of works are carried out by scholars aiming at the aspects of carbonate reservoir development forming mechanism, reservoir body classification, reservoir body prediction, fracture comprehensive research and the like, however, two major factors of a fracture zone and a fracture cave which are important components of a fractured-solute trap are fractured and researched for a long time, the research work of how to integrally describe the fractured-solute as a trap is relatively little, great influence is brought to subsequent resource amount estimation and reserve calculation, and the resource comprehensive evaluation, exploration, development and deployment in the fractured-solute development area are restricted.
Accordingly, the present invention provides a method and apparatus for solution trap description for solution reservoir exploration.
Disclosure of Invention
To solve the above problems, the present invention provides a method for describing a dissolved trap for use in the exploration of a dissolved hydrocarbon reservoir, the method comprising the steps of:
developing forward simulation research aiming at the structural style of the solution breaking body, establishing an earthquake recognition mode of the solution breaking body, and determining the reservoir type in the current solution breaking body;
based on the reservoir type, analyzing by using the emptying loss characteristics of the drilled well in the area where the current solution breaking body is located or the adjacent area and combining a logging curve, and determining a reservoir porosity threshold value and a reservoir threshold value of the current solution breaking body;
determining the boundary of the current cut solution trap after interactive verification on a three-dimensional space based on the reservoir threshold value;
and obtaining the reservoir volume, the reservoir effective volume and the effective reservoir thickness of the current solution breaking according to the determined reservoir type, the reservoir porosity threshold value, the reservoir threshold value and the boundary of the current solution breaking trap.
According to one embodiment of the invention, the reservoir type in the current solution is determined by:
performing seismic response characteristic analysis on the current solution breaking body through the solution breaking body seismic identification mode, and judging the type of the reservoir of the current solution breaking body to be a first type of reservoir when bead reflection occurs; when the disordered weak reflection occurs, judging that the reservoir type of the current solution breaking reservoir is a second reservoir type; and when the linear weak reflection occurs, judging the reservoir type of the current solution breaking reservoir to be a third reservoir type.
According to an embodiment of the present invention, the reservoir thresholds include a first type of reservoir threshold, a second type of reservoir threshold, and a third type of reservoir threshold, wherein:
setting a preset multiple of the energy threshold value as the interface energy of the area or the adjacent area as the first type of reservoir threshold value;
setting the automatic fracture extraction value to be larger than a first preset value as the second type reservoir threshold value;
and setting the automatic fracture extraction value larger than the second preset value and smaller than a third preset value as the third type reservoir threshold value.
According to an embodiment of the present invention, the step of determining the boundary of the current disconnected fluid trap after interactive verification on the three-dimensional space based on the reservoir threshold value further includes the following steps:
determining the type of the fractured zone according to the characteristics of the section and the plane pattern difference segmentation of the fractured zone so as to segment the fractured zone, wherein the segmented fractured zones are respectively used as the boundaries of upward traps of the current fractured solution traps in the same fractured zone;
dividing boundaries of reservoir extension development at two sides of a fracture zone by using the amplitude change rate and the reservoir threshold value of the drilling and logging statistics;
and interactively verifying the broken solution contour determined by the tensor body, the broken fracture zone range detected by the maximum likelihood body and the reservoir plane boundary described by the amplitude change rate on a three-dimensional space, and determining the boundary of the current broken solution trap.
According to an embodiment of the present invention, the classification of the fracture zone is determined according to the trend surface analysis and the coherent attribute in combination with the fracture response characteristics of the profile, wherein the classification of the fracture zone comprises: compression, translation, and tension.
According to one embodiment of the invention, the reservoir volume of the current dilacerate is obtained by:
and carrying out three-dimensional carving on the reservoir in the current solution according to the total energy volume and the enhanced coherent body, and calculating the volume of the reservoir in the current solution.
According to one embodiment of the invention, the reservoir effective volume is obtained by:
and carrying out volume correction by utilizing the porosity counted by the well-drilled well logging curve to obtain the effective volume of the reservoir of the current solution-breaking trap.
According to one embodiment of the invention, the effective reservoir thickness is obtained by:
and obtaining the trap geological reserves by using the oil or gas saturation of the current block and referring to the volume coefficient of the adjacent reserve submitting area, and obtaining the effective reservoir thickness according to the reservoir volume and the effective reservoir volume.
According to one embodiment of the invention, the effective reservoir thickness is calculated by:
and dividing the reservoir volume and the reservoir effective volume to obtain the effective reservoir thickness.
According to another aspect of the present invention there is also provided a solution trap description apparatus for solution reservoir exploration, the apparatus comprising:
the determining module is used for developing forward simulation research aiming at the structural style of the solution breaking body, establishing an earthquake recognition mode of the solution breaking body and determining the reservoir type in the current solution breaking body;
the parameter module is used for analyzing by combining a logging curve by utilizing the emptying loss characteristic of a drilled well in an area where the current solution breaking body is located or an adjacent area based on the reservoir type, and determining a reservoir porosity threshold value and a reservoir threshold value of the current solution breaking body;
the boundary module is used for determining the boundary of the current solution breaking trap after interactive verification in a three-dimensional space based on the reservoir threshold value;
and the calculation module is used for obtaining the reservoir volume, the reservoir effective volume and the effective reservoir thickness of the current solution breaking body according to the determined reservoir type, the reservoir porosity threshold value, the reservoir threshold value and the boundary of the current solution breaking body trap.
The method and the device for describing the broken solution trapping for exploring the broken solution hydrocarbon reservoir, provided by the invention, take the separability of a broken zone and the classification of a broken-cave-seam aggregate reservoir into consideration, so that the description of the broken solution trapping is more refined. According to the invention, through fracture zone segmentation description, classified reservoir prediction and volume estimation, and combined with three-dimensional visual image display, the characteristics of the horizontal segmentation, the longitudinal development depth, the internal structure difference and the like of the fractured solvent trap can be spatially described. The method can describe the solution trap in the actual exploration process, provides reference for drilling the solution trap with subsequent optimal reserve abundance meeting economic benefits, and provides reliable reference basis for improving the drilling success rate.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by the practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
FIG. 1 shows a flow diagram of an interrupted solution trap description method for exploration of an interrupted solution reservoir according to one embodiment of the present invention;
FIG. 2 shows a flow diagram for determining a boundary of a current broken fluid trap in a broken fluid trap description method for use in broken fluid reservoir exploration, according to an embodiment of the present invention;
FIG. 3 shows a flow diagram of an interrupted solution trap description method for exploration of an interrupted solution reservoir according to another embodiment of the invention;
FIG. 4 illustrates an exemplary geological model and forward modeling results according to one embodiment of the present invention;
FIG. 5 is a schematic diagram illustrating a first type of reservoir threshold determination according to an embodiment of the present invention;
FIG. 6 illustrates a second type and third type reservoir threshold determination according to one embodiment of the present invention;
FIG. 7 shows a schematic view of a fracture zone segment according to one embodiment of the present invention;
FIG. 8 illustrates a schematic diagram of the cross-sectional solvent trap boundary mutually verified by planar attributes and maximum likelihood, according to an embodiment of the present invention;
FIG. 9 shows a schematic of the estimation of effective reservoir thickness and different types of reservoir volumes of the solution breaker according to one embodiment of the invention; and
FIG. 10 shows a block diagram of an apparatus for solution trap description for solution reservoir exploration, according to an embodiment of the present invention.
Detailed Description
To make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention are described in detail below with reference to the accompanying drawings.
The fractured sol is a carbonate rock fracture-cave type reservoir body which develops along a fractured zone and is formed by the action of various factors such as the fractured zone, atmospheric fresh water, deep heat flow and the like as the name suggests. After the carbonate rock is extruded by a plurality of stages of structures, a broken zone with a certain scale develops along the deep broken zone, and after a plurality of stages of karst water seeps down along the broken zone or local hot liquid upwelles, the broken zone and cracks are corroded and reformed to form columnar corrosion-dissolving holes and hole reservoir bodies, and under the covering layer plugging of covering marl rock, mudstone and the like and the shielding of lateral compact marl rock, a trap type formed by irregular fracture control karst crack hole bodies is formed, and the trap type is called as a broken solution trap for short. The trap forms a special kind of dissolved hydrocarbon reservoir after the oil gas moves vertically along the deep fracture and fills into the reservoir in the later period.
The carbonate trap can be subdivided into an epibiotic fracture-cave trap, an inner curtain fracture-cave trap, a fracture-control fracture-cave trap, a composite trap, a submarine trap and the like according to a reservoir type, a master control factor and a cover condition, and the karst fracture-cave trap is one of karst fracture-cave traps and belongs to a fracture-control fracture-cave trap.
The characteristics and the cause of the fracture solution determine that a fracture system and a fracture-cave system formed by corrosion are the core of the trap description, so compared with the zonal control carbonate karst fracture-cave trap, the fracture does not appear only as a blocking condition or a dredging condition, a fracture zone is also a storage space, and the fracture-cave and the fracture need to be respectively described as two reservoir types.
In recent years, a series of works are carried out by scholars aiming at the aspects of carbonate rock reservoir development formation mechanism, reservoir body classification, reservoir body prediction, fracture comprehensive research and the like, and a large amount of literature investigation shows that the research progresses and results related to the fracture solution trapping mainly focus on the aspects of fracture solution reservoir characteristics, fracture solution prediction, fracture solution contour detection, fracture solution internal structure representation, fracture solution carving method research and the like.
In general, the description of the soluble trap in the prior art has the following problems:
(1) The identification description of the fracture zone and the prediction of the fracture-cavity reservoir stratum are carried out separately, the linear characteristics of the fracture zone and the prediction of the fractures in the adjacent area are emphasized, the boundary of the fracture-cavity is marked by taking a layer system and an abnormal boundary as limits, and the characteristics of a 'body' on a fracture-cavity space structure and the complexity of an internal structure are reflected.
(2) In the early carbonate rock trapping research, attention is paid to the fact where a fracture extends to, a sliding fracture zone has a dolphin effect and a ribbon effect, structural patterns on sections and planes of different parts of the same fracture zone are large in difference, stress backgrounds are different, reservoir development strength is also different, a cut-off karst trapping boundary is not applicable to the cut-off karst trapping boundary by a layer-control karst fracture-cave type trapping boundary cutting method, except that the cut-off karst trapping boundary is controlled by reservoir development and is also influenced by the fracture zone in a segmented mode.
(3) The method for detecting the sliding fractured zone (coherence, ant body, discontinuous property and the like) and predicting the fracture-cavity reservoir (amplitude change rate, wave impedance inversion and the like) is more, and a trap description method and an idea of a system are not formed.
Ignoring the description of the compartmentalization of the fractured zone and the difficulty in characterizing the "fractured solution" spatial volume for conventional carbonate trap description, fig. 1 shows a flow chart of a fractured solution trap description method for fractured solution reservoir exploration according to one embodiment of the invention.
As shown in fig. 1, in step S101, forward simulation research is performed on the structure pattern of the solution fracture, an solution fracture seismic recognition mode is established, and the reservoir type in the current solution fracture is determined.
Preferably, the reservoir type in the current solution is determined by:
performing seismic response characteristic analysis on the current solution breaking body through a solution breaking body seismic identification mode, and judging the type of the reservoir of the current solution breaking body to be a first type of reservoir when bead reflection occurs; when the mixed and disorderly weak reflection occurs, judging that the type of the reservoir with the current solution is a second type of reservoir; and when the linear weak reflection occurs, judging the reservoir type of the current solution breaking reservoir to be a third reservoir type.
In step S102, based on the reservoir type, the reservoir porosity threshold value and the reservoir threshold value of the current solution breaking are determined by analyzing the emptying loss characteristics of the drilled well in the area or the adjacent area where the current solution breaking is located in combination with the well logging curve.
Preferably, the reservoir thresholds include a first type reservoir threshold, a second type reservoir threshold, and a third type reservoir threshold, wherein:
setting a preset multiple of the energy threshold value as the interface energy of the area or the adjacent area as a first type reservoir threshold value;
setting the automatic fracture extraction value larger than the first preset value as a second type reservoir threshold value;
and setting the automatic fracture extraction value larger than the second preset value and smaller than a third preset value as a third reservoir threshold value.
In step S103, based on the reservoir threshold, after interactive verification in the three-dimensional space, the boundary of the current solution trap is determined.
Preferably, the boundary of the current disconnected solvent trap is determined by a method as shown in fig. 2:
first, in step S201, according to the characteristics of the section and the plane pattern difference segmentation of the fractured zone, the category of the fractured zone is determined to segment the fractured zone, and the segmented fractured zones are respectively used as the boundary where the current fractured solution is trapped in the same fractured zone and is trapped upward.
Preferably, the category of the fracture zone is determined according to the trend surface analysis and the coherent attribute in combination with the fracture response characteristics of the profile, wherein the category of the fracture zone includes: compression, translation, and tension.
Then, in step S202, the boundary of reservoir extension development on both sides of the fracture zone is divided by using the amplitude change rate in combination with the reservoir threshold value of the drilling and logging statistics.
Finally, in step S203, the fractured-solvent profile determined by the tensor body, the fractured-zone range detected by the maximum likelihood body, and the reservoir plane boundary delineated by the amplitude change rate are interactively verified in a three-dimensional space, and the boundary of the current fractured-solvent trap is determined.
In step S104, a reservoir volume, a reservoir effective volume, and an effective reservoir thickness of the current solution breaker are obtained according to the determined reservoir type, the reservoir porosity threshold, the reservoir threshold, and the boundary of the current solution breaker trap.
Preferably, the reservoir volume of the current dilysed solution is obtained by: and carrying out three-dimensional carving on the reservoir in the current solution according to the total energy volume and the enhanced coherent body, and calculating the volume of the reservoir in the current solution.
Preferably, the reservoir effective volume is obtained by: and (5) carrying out volume correction by using the porosity counted by the well-drilled logging curve to obtain the effective volume of the reservoir with the current solution trap.
Preferably, the effective reservoir thickness is obtained by: and obtaining the trap geological reserve by utilizing the oil or gas saturation of the current block and referring to the volume coefficient of the adjacent reserve submitting area, and obtaining the effective reservoir thickness according to the reservoir volume and the effective reservoir volume.
Preferably, the effective reservoir thickness is calculated by: and dividing the volume of the reservoir and the effective volume of the reservoir to obtain the effective reservoir thickness.
The method shown in fig. 1 is characterized in that firstly, on the basis of geological pattern research, a solution fracture seismic identification mode is established through a forward modeling analysis technology, the trap boundary is determined by carrying out prediction on different types of reservoirs through post-stack seismic data, the three-dimensional carving volume and the effective reservoir thickness of the different types of reservoirs are obtained by determining different types of reservoirs and porosity threshold values through well drilling data and well logging data, the trap resource quantity is estimated, and solution fracture trap description is completed.
FIG. 3 shows a flow diagram of a method for solution trap description for solution reservoir exploration, according to another embodiment of the present invention. As shown in fig. 3, the geological research results are analyzed for three-dimensional seismic data, well drilling data, and well logging data in a predetermined area.
In step S301, forward simulation analysis is performed to define a solution seismic recognition mode, and different reservoir recognition prediction methods are established. In particular, forward simulation research is carried out on the fractured solutions with different structural styles, and the seismic response characteristics of two most important factors, namely 'sliding fracture zone' and 'fracture hole' in the fractured solutions are determined.
In step S302, the logging data and the loss information are used to determine reservoir and other key thresholds. Specifically, a reservoir threshold value and a reservoir porosity threshold value are determined, and the reservoir porosity threshold value are determined by combining the drilling emptying loss characteristic of the area where the current broken solution is located or the drilled emptying loss characteristic of the adjacent area with the well logging curve analysis.
In step S303, on the basis of the attribute analysis, the solution trap boundary is drawn by checking the plane and the three-dimensional space against each other.
Preferably, the fracture zone is segmented firstly, the sliding fracture zone has the characteristic of section and plane pattern difference segmentation, different segments have different breaking degrees and different reservoir development degrees due to different stress backgrounds, the segmentation is carried out by combining technologies such as trend surface and coherent body with the seismic section fracture pattern, and the different segments are used as the boundary of the fracture solution trap on the same fracture zone and the upward trap.
And then predicting a reservoir development area, wherein a fracture solution contains fracture-cavity reservoirs with different erosion degrees, the fracture solution presents response characteristics such as disorder, abnormal strong reflection, blank reflection and the like in earthquake, and the boundaries of the extension development of the reservoirs at two sides of the fracture zone are divided by utilizing attributes such as amplitude change rate and the like in combination with reservoir threshold values of drilling and logging statistics.
And finally, carrying out spatial verification on the trap boundary, wherein the solution fracture is strong in longitudinal and transverse heterogeneity, the range of a fracture zone is detected by utilizing a tensor body profile and a maximum likelihood, and the reservoir development boundary which is described with the amplitude change rate attribute is interactively verified in a three-dimensional space, so that the solution fracture trap boundary is comprehensively determined.
In step S304, the different reservoir space types are three-dimensionally sculpted, and the volume and effective thickness are estimated in combination with the key threshold parameters. Specifically, the solution-breaking trapping can be divided into three types of fracture, cave and hole, parameters such as porosity and the like are different, the volume of the reservoir is calculated through classified volume carving, effective reservoir thicknesses of different types of reservoirs are obtained through threshold value parameter correction, and trapping resource quantity estimation is completed through other parameters such as porosity and the like.
The Tahe oil field is a typical karst fracture-cave type oil and gas reservoir, a trapping target is distributed along a paleo-water system and fractures, in the process that the oil field is expanded to the periphery exploration, the exploration target continuously extends to the deep part, the surface karst effect is weakened, a sliding fracture zone develops, a broken zone with multi-stage activities forms a fracture-cave type reservoir through fluid corrosion transformation, and a series of sliding fracture-cave type oil and gas reservoirs distributed along a deep large sliding fracture zone are formed after oil and gas filling. The solution-breaking trap has the characteristics of large longitudinal buried depth, strong transverse heterogeneity, no structural control, no uniform oil-water interface, mudstone top seal of an overlying area and side seal of compact limestone.
The existing carbonate rock trap description is mainly based on reservoir development control factor classification, layer system division and blocking to carry out trap boundary drawing and thickness description, and structural blocking, fracture blocking, sensitive attribute boundary determination, energy body description reservoir thickness and the like are considered more. These series-controlled, structure-controlled traps are difficult to describe and characterize in detail such as a cut-off solution.
In one embodiment, the broken soluble body trap in the northward region is described by the broken soluble body trap description method provided by the present invention, taking a certain three-dimensional work area (hereinafter referred to as W three-dimensional) in the northward region as an example.
Firstly, by analyzing the seismic response characteristics of the reservoir through well drilling emptying, leakage sections and well logging interpretation in the northward region, the string bead reflection, weak reflection and chaotic reflection can be the response of a good reservoir, meanwhile, statistical classification of drilling oil and gas results also proves that the success rate of 'clear fracture + abnormal reflection' is high, and the weak reflection in a fracture zone can also drill a large-scale reservoir meeting high yield.
FIG. 4 is a typical geologic model of an amputated solution and forward modeling results. Referring to fig. 4, the left half of the translational section, the left half of the extrusion section and the left half of the tension section are forward modeling models established according to typical geological models and seismic profile characteristics, and the right half of the left half are simulation results obtained by performing forward modeling on the forward modeling models.
An interrupted solution earthquake recognition mode is established by combining drilling geological response characteristic analysis and forward modeling results: fracture zone development scale seam holes (first type reservoirs) are expressed as 'bead string' reflection, and fracture zone development more uniform seam holes (second type reservoirs) are expressed as disorder weak reflection; when the fracture zone is mainly of translational nature, the fracture zone is narrow and has uniform physical properties (reservoir of the third category), it shows a "linear" weak reflection.
From three key factors of the broken solution, namely fracture, cave and hole, the step which is critical to the identification and description of the broken solution trap is also the prediction of the three reservoirs of different types, so that the threshold value of the reservoir is more critical. The large-scale caverns (reservoirs of the first type) are mainly characterized as 'beads' in earthquake, W three-dimensional drilled wells are few, different three-dimensional seismic data are different in energy, a uniform threshold value cannot be used, the statistics on lost wells with beads drilled in adjacent regions is considered, the energy of the beads is proportional to the stable T74 interface energy of Ordovician carbonate, the ratio of the bead threshold value to the T74 interface energy is 0.41 time, and the energy threshold value of the beads is 1298 (as shown in figure 5) when the ratio is used in W three-dimensional.
For fracture zones and hole reservoirs related to fractures, AFE (full-name automatic fracture extraction, which is a fracture and fracture prediction means in the field of petroleum and geophysical exploration) is a better prediction means, main fracture zones (second type of reservoirs) are used as good reservoir spaces, emptying loss occurs in the drilling process of a drilling well, the AFE values of the emptying loss section of the drilling well are counted, the average value is 201.2, and the AFE value larger than 201.2 can be used as a type of fracture (main fracture) reservoir threshold value (as shown in fig. 6).
For the third type of reservoir (fine fractures), counting the values of the AFE section fracture response positions in the W2 well side drilling process, wherein the small fracture emptying loss characteristics are not obvious, the average value of the AFE attribute is 28, and the AFE value which is larger than 28 and smaller than 201 can be used as the predicted threshold value of the W three-dimensional second type of fracture (fine fractures) reservoir (as shown in FIG. 6).
When a broken solvent trap boundary is determined, as the sliding fracture zone has subsection difference along the trend, the trap depiction relates to boundary division on the fracture zone trend, the fracture zone is divided into three sections according to extrusion, translation and tension by adopting trend surface analysis and coherent attribute combined with section fracture response characteristics, and different section conversion parts can be used as boundary basis of the trap along the fracture trend (figure 7);
when determining the boundary range of a reservoir in which the fractured fluid is eroded along the fractured zone, utilizing the abnormal amplitude change rate attribute in a certain time window to perform carving; because the broken solution is irregular in the longitudinal direction, the tensor body and the maximum likelihood body reflect the internal structure information of the broken solution with two different scales, the broken solution contour determined by the tensor body, the broken zone range detected by the maximum likelihood body and the reservoir plane boundary described by the amplitude change rate attribute are subjected to interactive verification on a three-dimensional space (figure 8), and finally the boundary of the broken solution trap is determined.
When the trap volume and the effective thickness of the solution breaking body are estimated, the solution breaking body comprises three different reservoirs (reservoirs) including beads (first type reservoirs), a main fracture zone (first type reservoirs) and fine cracks (first type reservoirs), and the total energy body is obtained
In the formula, x i Representing the amplitude value of the ith sample point in a time window, N is the number of sample points in the time window), and AFE body (further enhancement processing is carried out on the basis of coherence, enhanced coherence for short)And (3) performing three-dimensional carving on the reservoir threshold value, calculating the volumes of three types of reservoirs, and performing volume correction by using the porosity counted by the well-drilled well logging curve to finally obtain the effective volume of the solution-breaking trap reservoir (figure 9). And obtaining the trap geological reserve by using the oil or gas saturation of the block where the W three-dimension is located and referring to the volume coefficient of the adjacent reserve submitting area, and dividing the carved reservoir volume and the effective reservoir volume to obtain the effective reservoir thickness.
FIG. 10 shows a block diagram of an apparatus for solution trap description for solution reservoir exploration, according to an embodiment of the present invention. As shown in fig. 10, the description apparatus 1000 includes: a determination module 1001, a parameter module 1002, a boundary module 1003, and a calculation module 1004.
The determining module 1001 is configured to perform forward modeling research on a structural style of the solution breaker, establish an earthquake recognition mode of the solution breaker, and determine a reservoir type in the current solution breaker.
The parameter module 1002 is configured to determine a reservoir porosity threshold value and a reservoir threshold value of a current solution breaker by analyzing, based on a reservoir type, a well-drilled evacuation loss characteristic of an area where the current solution breaker is located or an adjacent area in combination with a well logging curve.
The boundary module 1003 is configured to determine a boundary of the current solution trap after interactive verification in a three-dimensional space based on the reservoir threshold.
The calculating module 1004 is configured to obtain a reservoir volume, a reservoir effective volume, and an effective reservoir thickness of the current solution breaker according to the determined reservoir type, the reservoir porosity threshold, the reservoir threshold, and the boundary of the current solution breaker trap.
In conclusion, the solution fracture trapping description method and device for solution fracture hydrocarbon reservoir exploration, provided by the invention, take fracture zone separability and fracture-cave-seam aggregate reservoir classification into consideration, so that solution fracture trapping description is more refined. According to the invention, through fracture zone segmentation description, classified reservoir prediction and volume estimation, and combined with three-dimensional visual image display, the characteristics of the horizontal segmentation, the longitudinal development depth, the internal structure difference and the like of the fractured solvent trap can be spatially described. The method can describe the solution trap in the actual exploration process, provides reference for drilling the solution trap with subsequent optimal reserve abundance meeting economic benefits, and provides reliable reference basis for improving the drilling success rate.
It is to be understood that the disclosed embodiments of the invention are not limited to the particular structures, process steps, or materials disclosed herein but are extended to equivalents thereof as would be understood by those ordinarily skilled in the relevant arts. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting.
Reference in the specification to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. Thus, the appearances of the phrase "one embodiment" or "an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment.
Although the embodiments of the present invention have been described above, the above description is only for the convenience of understanding the present invention, and is not intended to limit the present invention. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (8)
1. A method of solution trap description for solution reservoir exploration, the method comprising the steps of:
carrying out forward simulation research aiming at the structural style of the solution breaking body, establishing a solution breaking body seismic recognition mode, and determining the reservoir type in the current solution breaking body, wherein the seismic response characteristic analysis is carried out on the current solution breaking body through the solution breaking body seismic recognition mode, and the reservoir type of the current solution breaking body is judged to be a first-type reservoir, a second-type reservoir or a third-type reservoir according to the reflection condition;
based on the reservoir type, analyzing by utilizing the emptying loss characteristics of the drilled region where the current solution breaking body is located or the drilled adjacent region and combining a logging curve, determining a reservoir porosity threshold value and a reservoir threshold value of the current solution breaking body, wherein the reservoir threshold value comprises a first type of reservoir threshold value, a second type of reservoir threshold value and a third type of reservoir threshold value, and setting the preset multiple of the energy threshold value as the interface energy of the region where the current solution breaking body is located or the adjacent region as the first type of reservoir threshold value; setting the automatic fracture extraction value larger than a first preset value as the second type reservoir threshold value; setting the automatic fracture extraction value to be larger than a second preset value and smaller than a third preset value as the third type reservoir threshold value;
determining the boundary of the current broken solution trap after interactive verification in a three-dimensional space based on the reservoir threshold value;
obtaining the reservoir volume, the reservoir effective volume and the effective reservoir thickness of the current solution breaking according to the determined reservoir type, the reservoir porosity threshold value, the reservoir threshold value and the boundary of the current solution breaking trap;
based on the reservoir threshold value, the step of determining the boundary of the current solution trap after interactive verification in the three-dimensional space further comprises the following steps:
according to the characteristics of the section and the plane pattern difference segmentation of the fracture zone, determining the category of the fracture zone to segment the fracture zone, wherein the segmented fracture zone is respectively used as the boundary of the current broken solution trap in the same fracture zone and trap upwards;
dividing boundaries of reservoir extension development at two sides of a fracture zone by using the amplitude change rate and the reservoir threshold value of the drilling and logging statistics;
and interactively verifying the broken solution contour determined by the tensor body, the broken fracture zone range detected by the maximum likelihood body and the reservoir plane boundary described by the amplitude change rate on a three-dimensional space, and determining the boundary of the current broken solution trap.
2. The method of claim 1, wherein the type of reservoir in the current solution is determined by:
performing seismic response characteristic analysis on the current solution breaking body through the solution breaking body seismic identification mode, and judging the type of the reservoir of the current solution breaking body to be a first type of reservoir when bead reflection occurs; when the mixed and disorderly weak reflection occurs, judging that the type of the reservoir with the current solution is a second type of reservoir; and when the linear weak reflection occurs, judging that the reservoir type of the current solution breaking reservoir is a third reservoir type.
3. The method of claim 1, wherein determining a fracture zone classification based on trend surface analysis and coherence properties in combination with profile fracture response characteristics, wherein fracture zone classification comprises: compression, translation, and tension.
4. The method of any one of claims 1-3, wherein the reservoir volume for current dilission is obtained by:
and carrying out three-dimensional carving on the reservoir in the current solution according to the total energy volume and the enhanced coherent body, and calculating the volume of the reservoir in the current solution.
5. The method of claim 4, wherein the reservoir effective volume is obtained by:
and carrying out volume correction by utilizing the porosity counted by the well-drilled well logging curve to obtain the effective volume of the reservoir of the current solution-breaking trap.
6. The method of claim 5, wherein the effective reservoir thickness is obtained by:
and obtaining the trap geological reserves by using the oil or gas saturation of the current block and referring to the volume coefficient of the adjacent reserve submitting area, and obtaining the effective reservoir thickness according to the reservoir volume and the effective reservoir volume.
7. The method of claim 6, wherein the effective reservoir thickness is calculated by:
and dividing the reservoir volume and the reservoir effective volume to obtain the effective reservoir thickness.
8. An interrupted solution trap description device for exploration of an interrupted solution reservoir, characterized in that the method according to any of claims 1-7 is performed, the device comprising:
the determining module is used for developing forward simulation research aiming at the structural style of the solution breaking body, establishing an earthquake recognition mode of the solution breaking body and determining the reservoir type in the current solution breaking body;
the parameter module is used for analyzing by combining a logging curve by utilizing the emptying loss characteristic of a drilled well in an area where the current solution breaking body is located or an adjacent area based on the reservoir type, and determining a reservoir porosity threshold value and a reservoir threshold value of the current solution breaking body;
the boundary module is used for determining the boundary of the current solution trapping after interactive verification on a three-dimensional space based on the reservoir threshold value;
and the calculation module is used for obtaining the reservoir volume, the reservoir effective volume and the effective reservoir thickness of the current solution breaking body according to the determined reservoir type, the reservoir porosity threshold value, the reservoir threshold value and the boundary of the current solution breaking body trap.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910222303.0A CN111722277B (en) | 2019-03-22 | 2019-03-22 | Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910222303.0A CN111722277B (en) | 2019-03-22 | 2019-03-22 | Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111722277A CN111722277A (en) | 2020-09-29 |
| CN111722277B true CN111722277B (en) | 2023-03-14 |
Family
ID=72562806
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910222303.0A Active CN111722277B (en) | 2019-03-22 | 2019-03-22 | Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111722277B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114594529B (en) * | 2020-12-07 | 2023-12-01 | 中国石油化工股份有限公司 | Method for predicting reservoir of buried hill based on azimuth seismic attribute change rate |
| CN113139302B (en) * | 2021-05-20 | 2023-03-10 | 电子科技大学 | Area growth-based solution breaking identification method |
| CN114427458A (en) * | 2021-12-17 | 2022-05-03 | 中国石油化工股份有限公司 | Sea-facies carbonatite fracture control fracture-cave type oil and gas reservoir exploration method |
| CN114487023B (en) * | 2021-12-24 | 2024-03-12 | 常州大学 | EIT experimental device and method for monitoring production profile dynamics of multi-branch broken solution oil-gas reservoir combined production reservoir |
| CN115822563B (en) * | 2022-12-30 | 2023-06-09 | 中国石油化工股份有限公司 | Ultra-deep fracture-control hole type oil and gas reservoir high-yield well design method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105319585A (en) * | 2015-10-14 | 2016-02-10 | 中国石油化工股份有限公司 | Method for utilizing thin-layer interference amplitude recovery to identify oil and gas reservoir |
| CN107272065A (en) * | 2017-08-04 | 2017-10-20 | 中国石油化工股份有限公司 | The disconnected solution profile testing method of carbonate rock |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6370491B1 (en) * | 2000-04-04 | 2002-04-09 | Conoco, Inc. | Method of modeling of faulting and fracturing in the earth |
| CN105626058B (en) * | 2015-12-30 | 2018-11-16 | 中国石油天然气股份有限公司 | Method and device for determining development degree of reservoir karst |
| CN105572739B (en) * | 2016-02-19 | 2019-03-12 | 中国石油集团川庆钻探工程有限公司 | Method for judging carbonatite hole crack development characteristics |
-
2019
- 2019-03-22 CN CN201910222303.0A patent/CN111722277B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105319585A (en) * | 2015-10-14 | 2016-02-10 | 中国石油化工股份有限公司 | Method for utilizing thin-layer interference amplitude recovery to identify oil and gas reservoir |
| CN107272065A (en) * | 2017-08-04 | 2017-10-20 | 中国石油化工股份有限公司 | The disconnected solution profile testing method of carbonate rock |
Non-Patent Citations (2)
| Title |
|---|
| New insights into the carbonate karstic fault system and reservoir formation in the Southern Tahe area of the Tarim Basin;Xinbian Lu et al.;《Marine and Petroleum Geology》;20170617;第587-605页 * |
| 塔河油田碳酸盐岩断溶体刻画技术研究与应用;王震等;《石油物探》;20190131;第58卷(第1期);第149-154页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111722277A (en) | 2020-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111722277B (en) | Method and device for describing broken solution trap for broken solution oil and gas reservoir exploration | |
| US11156743B2 (en) | Lateral placement and completion design for improved well performance of unconventional reservoirs | |
| CN110308487B (en) | Quantitative characterization method for solution-fractured oil reservoir | |
| US9081918B2 (en) | Methods and systems regarding models of underground formations | |
| CN112698392A (en) | Multilayer system reservoir three-dimensional drawing and oil-gas spatial distribution and constant volume method | |
| CN110632652B (en) | Method and device for determining oil-gas migration mode of fault block oil-gas reservoir | |
| CN110308497B (en) | Method for fusing different types of reservoir body models of fracture-cavity oil reservoir | |
| US20170096880A1 (en) | Reservoir permeability upscaling | |
| CN112489208A (en) | Crack sheet extraction method based on ant algorithm and three-dimensional geological model construction method | |
| CN111830558B (en) | Fracture zone engraving method | |
| CN116449423A (en) | Method, device, equipment and storage medium for analyzing single sand body configuration of braided river | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| CN106948812B (en) | Method and device for determining permeability lower limit value of hypertonic zone | |
| CN111680386B (en) | Method and device for analyzing sealing performance of brittle cover layer of broken block trap | |
| CN114114452B (en) | Characterization method and system for broken solution type carbonate reservoir | |
| CN112394416B (en) | Heterogeneous fracture control reservoir prediction method and device | |
| CN112505761B (en) | Reservoir gas content detection method and device | |
| CN116559938A (en) | Method for establishing fracture model of oil-gas reservoir and electronic equipment | |
| CN112649862B (en) | Broken solution identification method and device based on stratum structure information separation | |
| CN113534248A (en) | Method, device and system for quantitatively analyzing closure of trap fault | |
| CN112285774A (en) | Method for rapidly evaluating walk-slip fault sealing capability based on multiple geological parameters | |
| CN114114411A (en) | Quantitative determination method and device for three-dimensional fault conductibility | |
| US10460051B2 (en) | Computationally-efficient modeling of viscous fingering effect for enhanced oil recovery (EOR) agent injected at multiple injection concentrations | |
| CN112965108B (en) | Method and system for determining vertical closure of trap cap layer | |
| CN115932967B (en) | Crack identification method based on ant body algorithm |
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 |