CN106951660B - Sea facies clastic rock horizontal well reservoir logging interpretation method and device - Google Patents
Sea facies clastic rock horizontal well reservoir logging interpretation method and device Download PDFInfo
- Publication number
- CN106951660B CN106951660B CN201710216848.1A CN201710216848A CN106951660B CN 106951660 B CN106951660 B CN 106951660B CN 201710216848 A CN201710216848 A CN 201710216848A CN 106951660 B CN106951660 B CN 106951660B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- lithologic
- horizontal well
- region
- water
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 191
- 208000035126 Facies Diseases 0.000 title claims abstract description 90
- 238000000034 method Methods 0.000 title claims abstract description 69
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 250
- 238000011160 research Methods 0.000 claims abstract description 83
- 230000035699 permeability Effects 0.000 claims abstract description 54
- 238000012937 correction Methods 0.000 claims description 61
- 230000004044 response Effects 0.000 claims description 44
- 230000015572 biosynthetic process Effects 0.000 claims description 21
- 239000008398 formation water Substances 0.000 claims description 20
- 230000002159 abnormal effect Effects 0.000 claims description 19
- 238000012360 testing method Methods 0.000 claims description 19
- 230000035772 mutation Effects 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 13
- 238000010219 correlation analysis Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 8
- 235000020681 well water Nutrition 0.000 claims description 4
- 239000002349 well water Substances 0.000 claims description 4
- 238000011161 development Methods 0.000 abstract description 10
- 239000007789 gas Substances 0.000 description 17
- 230000007547 defect Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 15
- 238000004519 manufacturing process Methods 0.000 description 10
- 230000018109 developmental process Effects 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 238000005553 drilling Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 239000012530 fluid Substances 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000003860 storage Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 230000003247 decreasing effect Effects 0.000 description 3
- 238000009933 burial Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- -1 and the like Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000033558 biomineral tissue development Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Images
Classifications
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a sea facies clastic rock horizontal well reservoir logging interpretation method and device, and relates to the technical field of clastic rock oil and gas development. The method comprises the following steps: dividing a research area into a plurality of lithologic areas with different lithologic differences, and dividing reservoir beds of the lithologic areas into an undisturbed reservoir bed and a water flooded reservoir bed; establishing a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the water flooded reservoir of each lithologic region; determining the position of a horizontal well track of a reservoir in a lithologic region in a marine facies clastic rock reservoir; correcting a natural gamma curve, a sound wave time difference curve and a stratum deep resistivity curve of a horizontal well in the same reservoir section; determining the water flooding condition of a reservoir section of a horizontal well of a lithologic region; and selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the lithologic region horizontal well according to the water logging condition of the reservoir section of the lithologic region horizontal well, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the lithologic region horizontal well.
Description
Technical Field
The invention relates to the technical field of clastic rock oil and gas development, in particular to a sea-facies clastic rock horizontal well reservoir logging interpretation method and device.
Background
At present, marine clastic rock is an important area for global oil and gas reserves to grow, and can form a huge-scale oil and gas field. The 20 th century is 40 to 50 th, and is the time of large discovery of the marine facies clastic rock oil and gas field in the world, the marine facies clastic rock oil and gas field is mainly concentrated in the middle east Persian bay area and the middle Asia area, the new Tetis structural domain of the middle and new generations is taken as the main part, most of oil and gas reservoirs are mainly giant-oversize structural oil and gas reservoirs, the oil and gas storage scale is huge, and the maximum recoverable storage amount is up to 120 hundred million tons. The oil and gas reserves of the clastic rock of the marine facies in China are mainly distributed in Tarim basins, Ordos basins, Sichuan basins and the like. Compared with the land clastic rock, the marine clastic rock has weak heterogeneity, the well-to-well contrast of sand bodies is strong, the physical property of reservoir layers is good, and a horizontal well is generally adopted for efficient development.
The horizontal well oil recovery is an advanced technology applied to oil and gas development, the horizontal well refers to a well with a well inclination angle reaching or approaching 90 degrees, a well body drills a well with a certain length along the horizontal direction, the oil and gas productivity can be improved by increasing the contact area of a well hole and a stratum and penetrating a zone with good permeability, the drilling trend is controlled according to the oil and gas water distribution so as to reduce the water coning and the gas coning, the oil and gas recovery ratio of a single well can be improved, and the exploration and development efficiency is improved. With the continuous development of the horizontal well drilling technology, the horizontal well oil production technology is widely applied, reservoir interpretation and evaluation play an important role in subsequent work such as a horizontal well completion mode, and the like, so that the horizontal well logging interpretation precision needs to be continuously improved. At present, horizontal well logging interpretation basically follows the idea and method of a straight well, and logging data are used for reservoir interpretation after being simply corrected. However, horizontal well bores extend nearly horizontally in the stratum, and the spatial positions of the horizontal well bores are different from those of the vertical well, so that the classical interpretation model is poor in applicability. Meanwhile, in the prior art, the whole horizontal section is interpreted as the same reservoir, the non-reservoir section is not identified, the position of the well track of the horizontal well in the oil reservoir cannot be determined, and finally, reasonable interpretation parameters and interpretation conclusions cannot be obtained easily, and the water flooding interpretation of the horizontal well is lacked. Therefore, the result of well logging interpretation of the marine clastic rock horizontal well reservoir in the prior art is inaccurate.
Disclosure of Invention
The embodiment of the invention provides a sea-facies clastic rock horizontal well reservoir logging interpretation method and device, and aims to solve the problem that in the prior art, the result of the sea-facies clastic rock horizontal well reservoir logging interpretation is inaccurate.
In order to achieve the purpose, the invention adopts the following technical scheme:
a well logging interpretation method for a marine clastic rock horizontal well reservoir comprises the following steps:
dividing a research area into a plurality of lithological areas with different lithological differences according to a reservoir rock core test result of a marine facies clastic rock vertical well in the research area and a conventional logging curve of the vertical well, and dividing the reservoir in each lithological area into an original reservoir and a water flooded reservoir according to the reservoir water flooding condition of each lithological area;
according to the core data and the logging data of the research area, a core scale logging method is adopted to establish a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the water flooded reservoir of each lithologic area;
determining the position of a horizontal well borehole trajectory of a reservoir of a lithologic region in a marine facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic wave time difference curve and resistivity data of the reservoir of the lithologic region and by combining a vertical well stratum dividing result of the reservoir of the lithologic region;
according to the position of the well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, acquiring a natural gamma curve, a sonic time difference curve and a stratum deep resistivity curve of the vertical well and the horizontal well in the same reservoir section, and correcting the natural gamma curve, the sonic time difference curve and the stratum deep resistivity curve of the horizontal well in the same reservoir section;
determining the water flooding condition of the reservoir section of the horizontal well in the lithologic region according to a preset water flooding qualitative interpretation rule of the horizontal well, a corrected natural gamma curve, a corrected sound wave time difference curve and a corrected stratum deep resistivity curve of the horizontal well;
and selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the lithologic region horizontal well according to the water logging condition of the reservoir section of the lithologic region horizontal well, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the lithologic region horizontal well.
Specifically, according to a reservoir core test result of a marine facies clastic rock vertical well in a research area and a conventional logging curve of the vertical well, the research area is divided into a plurality of lithological areas with different lithological differences, and according to a reservoir flooding condition of each lithological area, the reservoir of each lithological area is divided into an original reservoir and a flooding reservoir, and the method comprises the following steps:
determining the heterogeneity of the marine clastic rock reservoir in the research area according to the reservoir core test result of the marine clastic rock vertical well in the research area and the conventional logging curve of the vertical well;
dividing a research area into a plurality of lithologic areas with different lithologic differences according to the heterogeneity of marine facies clastic rock reservoirs in the research area;
and dividing the reservoir stratum of each lithologic area into an original reservoir stratum and a water flooding reservoir stratum according to the water flooding condition of the reservoir stratum of each lithologic area.
Specifically, according to the core data and the logging data of the research area, a core scale logging method is adopted to establish a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the water flooded reservoir of each lithologic area, and the logging interpretation model comprises the following steps:
dividing a sample layer according to the marine facies clastic rock core well;
establishing a shale content logging interpretation model of an undisturbed reservoir and a water flooded reservoir of each lithologic area; the mud content well logging interpretation model is
Wherein, VshThe shale content of an undisturbed reservoir and a water flooded reservoir in a lithologic region; GCUR is an empirical coefficient; the delta GR is the natural gamma relative value of the undisturbed reservoir and the water flooded reservoir in the lithologic region;
GR is the natural gamma value of an undisturbed reservoir and a water flooded reservoir in a lithologic area; GRminNatural gamma value of pure sandstone; GRmaxNatural gamma value of pure mudstone;
performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the rocky region and the acoustic time difference logging curve data of the undisturbed reservoir and the water flooded reservoir in the rocky region, and establishing an acoustic time difference relative porosity model; the acoustic time difference relative porosity model is phi xAC-y; wherein, phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; AC is the acoustic wave time difference value of the undisturbed reservoir and the water flooded reservoir in the rock region; x and y are model parameters;
performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the lithologic region and the permeability of the undisturbed reservoir and the water flooded reservoir in the lithologic region, and establishing a porosity relative permeability model; the porosity relative permeability model is Perm & ltxe & gtyΦ(ii) a Perm is the permeability of an undisturbed reservoir and a water flooded reservoir in a lithologic region; e is a natural constant; phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; x and y are model parameters;
establishing a water saturation model of an undisturbed reservoir and a water flooded reservoir in a lithologic region according to an Archie formula:
wherein a is a lithology coefficient related to lithology; b is a constant related to lithology; m is a cementation index; n is a saturation index; rw is the resistivity of the formation water of the undisturbed reservoir and the flooded reservoir in the rock region; rt is the stratum resistivity of an undisturbed reservoir and a flooded reservoir in the lithologic region; phi is the porosity of an undisturbed reservoir and a water flooded reservoir in the lithologic region; sw is the water saturation of an undisturbed reservoir and a water flooded reservoir in the lithologic region; the lithology coefficient a related to the lithology, the constant b related to the lithology, the cementation index m, the saturation index n, the formation water resistivity Rw of the undisturbed reservoir and the flooded reservoir of the lithology area are obtained from a preset corresponding relation table; the corresponding relation table comprises lithology coefficients, constants, cementation indexes, saturation indexes and formation water resistivity, which are corresponding to the undisturbed reservoir and the water flooded reservoir of each lithologic area and are related to lithology;
establishing the measurement of the oil saturation of the undisturbed reservoir and the flooded reservoir of the lithologic region according to the water saturation modelThe well interpretation model is as follows: so＝1-Sw(ii) a Wherein S isoThe oil saturation of the reservoir is the original reservoir in the rock region and the oil saturation of the reservoir flooded with water.
Specifically, according to a horizontal well natural gamma curve, a sonic time difference curve and resistivity data of a reservoir stratum of a rock region, determining the position of a horizontal well borehole trajectory of the reservoir stratum of the rock region in a marine facies clastic rock reservoir by combining a vertical well stratum division result of the reservoir stratum of the rock region, and the determining method comprises the following steps of:
identifying the mutation points of the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region according to the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region;
determining an interface of an overlying stratum and a reservoir stratum of the rock region according to mutation points of a horizontal well natural gamma curve and an acoustic wave time difference curve of the reservoir stratum of the rock region;
determining logging response abnormal points inside marine facies clastic rocks of the lithological zone according to the horizontal well resistivity data of the reservoir of the lithological zone;
and determining the position of the horizontal well track of the reservoir of the lithologic region in the marine facies clastic rock reservoir according to the interface of the overlying stratum and the reservoir of the lithologic region and the logging response abnormal point inside the marine facies clastic rock of the lithologic region by combining the vertical well stratum dividing result of the reservoir of the lithologic region.
Specifically, according to the position of a well track of a horizontal well of a reservoir in the lithologic region in the marine facies clastic rock reservoir, acquiring a natural gamma curve, an acoustic time difference curve and a stratum deep resistivity curve of a vertical well and the horizontal well in the same reservoir section, and correcting the natural gamma curve, the acoustic time difference curve and the stratum deep resistivity curve of the horizontal well in the same reservoir section;
according to the formula: ACCorrection ofAnd (3) correcting the acoustic wave time difference curve of a horizontal well of the same reservoir section by α AC + β, wherein AC isCorrection ofThe corrected horizontal well sound wave time difference value, AC is the original sound wave time difference value of the horizontal well, α and β are correction parameters;
according to the formula: rt correction＝αRt+ β ground for horizontal well of same reservoir sectionCorrecting a layer depth resistivity curve; wherein R ist correctionThe horizontal well stratum resistivity value after correction is obtained; rtWhich is the original formation resistivity value of the horizontal well, α and β are the calibration parameters.
Specifically, the preset horizontal well flooding qualitative interpretation rule comprises each flooding condition, a natural gamma value range and a characteristic value, a sound wave time difference value range and a characteristic value, and a resistivity value range and a characteristic value corresponding to each flooding condition;
the method for determining the water logging condition of the reservoir section of the horizontal well in the lithologic region according to the preset water logging qualitative interpretation rule of the horizontal well, the corrected natural gamma curve, the corrected sound wave time difference curve and the corrected stratum deep resistivity curve of the horizontal well comprises the following steps:
and inquiring the corrected water flooding conditions corresponding to the natural gamma curve, the acoustic wave time difference curve and the stratum deep resistivity curve of the horizontal well in the preset water flooding qualitative interpretation rule of the horizontal well, and identifying that the water flooding condition of the reservoir section of the horizontal well in the lithologic region is an oil layer or a reservoir layer of a poor oil layer.
A sea facies clastic rock horizontal well reservoir logging interpretation device comprises:
the device comprises a research area dividing unit, a data processing unit and a data processing unit, wherein the research area dividing unit is used for dividing a research area into a plurality of lithological areas with different lithological differences according to a reservoir rock core test result of a marine facies clastic rock vertical well in the research area and a conventional logging curve of the vertical well, and dividing the reservoir in each lithological area into an original reservoir and a water logging reservoir according to the reservoir water logging condition of each lithological area;
the well logging interpretation model establishing unit is used for establishing a well logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of an undisturbed reservoir and a flooded reservoir of each lithologic area by adopting a core scale well logging method according to core data and well logging data of a research area;
the horizontal well track position determining unit is used for determining the position of the horizontal well track of the reservoir of the rock region in the sea facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic wave time difference curve and resistivity data of the reservoir of the rock region and by combining a vertical well stratum dividing result of the reservoir of the rock region;
the curve data correction unit is used for acquiring natural gamma curves, acoustic time difference curves and stratum deep resistivity curves of a vertical well and a horizontal well in the same reservoir interval according to the position of a well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, and correcting the natural gamma curves, the acoustic time difference curves and the stratum deep resistivity curves of the horizontal well in the same reservoir interval;
the water flooding condition determining unit is used for determining the water flooding condition of the reservoir section of the horizontal well in the lithologic region according to the preset water flooding qualitative interpretation rule of the horizontal well, the corrected natural gamma curve, the corrected sound wave time difference curve and the corrected stratum deep resistivity curve of the horizontal well;
and the horizontal well logging interpretation unit is used for selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the horizontal well of the lithologic region according to the water logging condition of the reservoir section of the horizontal well of the lithologic region, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the horizontal well of the lithologic region.
Furthermore, the study region dividing unit is specifically configured to:
determining the heterogeneity of the marine clastic rock reservoir in the research area according to the reservoir core test result of the marine clastic rock vertical well in the research area and the conventional logging curve of the vertical well;
dividing a research area into a plurality of lithologic areas with different lithologic differences according to the heterogeneity of marine facies clastic rock reservoirs in the research area;
and dividing the reservoir stratum of each lithologic area into an original reservoir stratum and a water flooding reservoir stratum according to the water flooding condition of the reservoir stratum of each lithologic area.
In addition, the well logging interpretation model building unit is specifically configured to:
dividing a sample layer according to the marine facies clastic rock core well;
establishing a shale content logging interpretation model of an undisturbed reservoir and a water flooded reservoir of each lithologic area; the mud content well logging interpretation model is
Wherein, VshThe shale content of an undisturbed reservoir and a water flooded reservoir in a lithologic region; GCUR is an empirical coefficient; the delta GR is the natural gamma relative value of the undisturbed reservoir and the water flooded reservoir in the lithologic region;
GR is the natural gamma value of an undisturbed reservoir and a water flooded reservoir in a lithologic area; GRminNatural gamma value of pure sandstone; GRmaxNatural gamma value of pure mudstone;
performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the rocky region and the acoustic time difference logging curve data of the undisturbed reservoir and the water flooded reservoir in the rocky region, and establishing an acoustic time difference relative porosity model; the acoustic time difference relative porosity model is phi xAC-y; wherein, phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; AC is the acoustic wave time difference value of the undisturbed reservoir and the water flooded reservoir in the rock region; x and y are model parameters;
performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the lithologic region and the permeability of the undisturbed reservoir and the water flooded reservoir in the lithologic region, and establishing a porosity relative permeability model; the porosity relative permeability model is Perm & ltxe & gtyΦ(ii) a Perm is the permeability of an undisturbed reservoir and a water flooded reservoir in a lithologic region; e is a natural constant; phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; x and y are model parameters;
establishing a water saturation model of an undisturbed reservoir and a water flooded reservoir in a lithologic region according to an Archie formula:
wherein a is a lithology coefficient related to lithology; b is a constant related to lithology; m is a cementation index; n is a saturation index; rw is the resistivity of the formation water of the undisturbed reservoir and the flooded reservoir in the rock region; rt is the stratum resistivity of an undisturbed reservoir and a flooded reservoir in the lithologic region; phi is the porosity of an undisturbed reservoir and a water flooded reservoir in the lithologic region; sw is undisturbed reservoir and water in the lithologic regionFlooding the water saturation of the reservoir; the lithology coefficient a related to the lithology, the constant b related to the lithology, the cementation index m, the saturation index n, the formation water resistivity Rw of the undisturbed reservoir and the flooded reservoir of the lithology area are obtained from a preset corresponding relation table; the corresponding relation table comprises lithology coefficients, constants, cementation indexes, saturation indexes and formation water resistivity, which are corresponding to the undisturbed reservoir and the water flooded reservoir of each lithologic area and are related to lithology;
and establishing a logging interpretation model of the oil saturation of the undisturbed reservoir and the flooded reservoir in the lithologic region according to the water saturation model: so＝1-Sw(ii) a Wherein S isoThe oil saturation of the reservoir is the original reservoir in the rock region and the oil saturation of the reservoir flooded with water.
In addition, the horizontal well borehole trajectory position determination unit is specifically configured to:
identifying the mutation points of the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region according to the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region;
determining an interface of an overlying stratum and a reservoir stratum of the rock region according to mutation points of a horizontal well natural gamma curve and an acoustic wave time difference curve of the reservoir stratum of the rock region;
determining logging response abnormal points inside marine facies clastic rocks of the lithological zone according to the horizontal well resistivity data of the reservoir of the lithological zone;
and determining the position of the horizontal well track of the reservoir of the lithologic region in the marine facies clastic rock reservoir according to the interface of the overlying stratum and the reservoir of the lithologic region and the logging response abnormal point inside the marine facies clastic rock of the lithologic region by combining the vertical well stratum dividing result of the reservoir of the lithologic region.
Furthermore, the curve data correction unit is specifically configured to;
according to the formula: ACCorrection ofAnd (3) correcting the acoustic wave time difference curve of a horizontal well of the same reservoir section by α AC + β, wherein AC isCorrection ofThe corrected horizontal well sound wave time difference value, AC is the original sound wave time difference value of the horizontal well, α and β are correction parameters;
according to the formula: rt correction＝αRt+ β the formation deep resistivity curve of the horizontal well in the same reservoir section is corrected, wherein Rt correctionThe horizontal well stratum resistivity value after correction is obtained; rtWhich is the original formation resistivity value of the horizontal well, α and β are the calibration parameters.
In addition, the preset horizontal well flooding qualitative interpretation rule in the flooding condition determination unit comprises natural gamma value ranges and characteristic values, acoustic wave time difference value ranges and characteristic values, and resistivity value ranges and characteristic values corresponding to the flooding conditions and the flooding conditions;
the flooding condition determining unit is specifically configured to:
and inquiring the corrected water flooding conditions corresponding to the natural gamma curve, the acoustic wave time difference curve and the stratum deep resistivity curve of the horizontal well in the preset water flooding qualitative interpretation rule of the horizontal well, and identifying that the water flooding condition of the reservoir section of the horizontal well in the lithologic region is an oil layer or a reservoir layer of a poor oil layer.
According to the sea facies clastic rock horizontal well reservoir logging interpretation method and device provided by the embodiment of the invention, firstly, a research area is divided into a plurality of lithological areas with different lithological differences according to a reservoir core test result of a sea facies clastic rock vertical well in the research area and a conventional logging curve of the vertical well, and the reservoir in each lithological area is divided into an original reservoir and a water logging reservoir according to the reservoir water logging condition of each lithological area; then, according to the core data and the logging data of the research area, a core scale logging method is adopted to establish a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the water flooded reservoir of each lithologic area; determining the position of a horizontal well borehole trajectory of a reservoir of a lithologic region in a marine facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic wave time difference curve and resistivity data of the reservoir of the lithologic region and by combining a vertical well stratum dividing result of the reservoir of the lithologic region; according to the position of the well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, acquiring a natural gamma curve, a sonic time difference curve and a stratum deep resistivity curve of the vertical well and the horizontal well in the same reservoir section, and correcting the natural gamma curve, the sonic time difference curve and the stratum deep resistivity curve of the horizontal well in the same reservoir section; determining the water flooding condition of the reservoir section of the horizontal well in the lithologic region according to a preset water flooding qualitative interpretation rule of the horizontal well, a corrected natural gamma curve, a corrected sound wave time difference curve and a corrected stratum deep resistivity curve of the horizontal well; and selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the lithologic region horizontal well according to the water logging condition of the reservoir section of the lithologic region horizontal well, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the lithologic region horizontal well. Therefore, the method can overcome the defects that the horizontal well reservoir is difficult to interpret and no marine clastic rock geological constraint is generated to correct the horizontal well logging curve in the prior art, and can make up the defect of lack of flooding interpretation in the horizontal well reservoir interpretation, so that the result of the marine clastic rock horizontal well reservoir logging interpretation is more accurate.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
Fig. 1 is a first flowchart of a well logging interpretation method for a marine clastic rock horizontal well reservoir, provided by an embodiment of the present invention;
FIG. 2 is a second flowchart of a well logging interpretation method for a marine clastic rock horizontal well reservoir, provided by an embodiment of the present invention;
FIG. 3 is a graph of the position relationship between the HD11-8H wellbore trajectories and the reservoir in an embodiment of the invention;
FIG. 4 is a graph comparing the natural gamma values of a vertical well reservoir section and a horizontal well reservoir section of a Hardson oil field east river sandstone in a Tarim basin in an embodiment of the invention;
FIG. 5 is a schematic diagram showing a comparison of acoustic time difference curves of a straight well and a horizontal well in northwest region of the Hardson oil field in the Tarim basin in the embodiment of the invention;
FIG. 6 is a schematic diagram showing a comparison of resistivity curves of a straight well and a horizontal well in the northwest region of the Hadison oil field in the Tarim basin in accordance with an embodiment of the present invention;
FIG. 7 is a diagram illustrating the well logging results of horizontal wells HD11-8H reservoir in an embodiment of the invention;
fig. 8 is a schematic structural diagram of a logging interpretation device for a marine clastic rock horizontal well reservoir, provided by an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the process of implementing the embodiment of the present invention, the inventor finds the development and research of the prior art as follows:
the development and research of the prior art mainly relate to the two aspects of correcting the position of a horizontal well borehole path in an oil reservoir and the influence of a horizontal well logging response.
The method comprises the steps of searching a vertical well section corresponding to a horizontal well track in an oil reservoir according to the horizontal well logging curve characteristics, namely identifying a part similar to the vertical well logging curve characteristics on the horizontal well curve, and combining the geological structure characteristics of the oil reservoir to clarify the formation structure form drilled by the horizontal well so as to complete the research on the relation between the horizontal well track and the oil reservoir (ZhouCan, and the like, geophysical progress 2006; Zhao Xiaoqing Jilin university doctor thesis, 2010; Lanwenjian, and the like, drilling and mining process, 2012, red rock, and the like, petroleum and natural gas science, 2012). However, because the well body of the horizontal well is in a horizontal state, the measured curve form is greatly different from that of the vertical well. The curve amplitude changes obviously when a straight well passes through different reservoirs, and the stratum interface is easy to judge, in a horizontal well, only a vertical section and a deflecting section can pass through different strata, the horizontal section is drilled along a certain target reservoir, most of well bores of the horizontal well pass through the same layer, the curve is smooth on the whole, the change is not obvious, and the well logging curve abnormity can occur only at a non-reservoir section, so that the stratum interface is difficult to judge. In the process of drilling a horizontal well, a borehole trajectory is not a straight line, fluctuates up and down, plus lithology and fluid differences, and is not combined with geological research, and the formation interface is determined to have large artificial subjectivity and uncertainty only by using horizontal well logging response.
In addition, in the study of horizontal well logging response influence correction, numerical simulation from a large inclination to the horizontal direction is carried out on a logging instrument, and corresponding logging response is studied. Meanwhile, sound wave and resistivity anisotropy experiments are carried out on the same core samples in the vertical and horizontal directions, then the anisotropy of a reservoir is analyzed, the logging curve difference of the rock in the horizontal and vertical directions is researched, and horizontal well reservoir logging interpretation is carried out on the basis (Luh P C.seg Expanded Abstracts, 1992; Wang super, et al, Tuha oil and gas, 2011; Chen wood silver, et al, overseas logging technology, 2013). In addition, a numerical simulation calculation is carried out on the well logging curve of the horizontal well by adopting a three-dimensional element method, the well logging curve of the horizontal well is corrected by adopting a rapid deconvolution method, the influence of the layer thickness-surrounding rock on the well logging curve under the stratum level condition is researched by adopting a three-dimensional finite element method, and the well logging curve of the horizontal well is corrected based on influence factors (Liyongjie, et al. science and technology and engineering, 2014; Tantang, Wujiyuan, et al. geophysical science and technology, 2012; Liwutiezhu. foreign well logging technology, 2016; Wujiyuan, et al. petroleum geology and engineering, 2016; Chengqing, et al. unconventional oil gas, 2016). However, the response characteristics of horizontal well logging are different from those of a vertical well, and the stratum in the direction perpendicular to the instrument axis is not an isotropic homogeneous body any more but an anisotropic heterogeneous body in most cases. The main factors influencing the logging response of the horizontal well comprise a borehole, well deviation, layer thickness, surrounding rock, drilling fluid invasion and the like, and further comprise instrument eccentricity influence, drilling cutting layer influence and gas containing influence. The geological and engineering factors can influence the logging response of the horizontal well, wherein each influencing factor is very difficult to correct on a straight well, and the horizontal well is more difficult to correct; and the problems of well logging interpretation research of a vertical well reservoir and dislocation of a horizontal well exist, namely the research result of the vertical well cannot be applied to the horizontal well and cannot play a guiding role. In the process of horizontal well reservoir well logging interpretation, the horizontal well interpretation and geological study are neglected and an effective method for sea facies clastic rock horizontal well reservoir interpretation is lacked, wherein the interpretation is generally based on response and correction study of a logging instrument.
Therefore, based on the development and research of the prior art, the invention aims to overcome the defects that in the prior art scheme, horizontal well reservoir logging interpretation is difficult, no marine clastic rock geological constraint exists to correct a horizontal well logging curve, and meanwhile, the defect that the flooding interpretation in the horizontal well reservoir interpretation is lack is overcome, and provides the method for rapidly interpreting the marine clastic rock horizontal well reservoir by using rock cores, conventional logging and horizontal well logging information, so that the accuracy of the marine clastic rock horizontal well reservoir interpretation is improved.
In order to overcome the problems in the prior art and achieve the above object, as shown in fig. 1, an embodiment of the present invention provides a method for explaining well logging in a marine clastic rock horizontal well reservoir, including:
And step 102, according to the core data and the logging data of the research area, adopting a core scale logging method to establish a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the water flooded reservoir of each lithologic area.
103, determining the position of a horizontal well borehole trajectory of the reservoir of the lithologic region in the marine facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic wave time difference curve and resistivity data of the reservoir of the lithologic region and by combining a vertical well stratum dividing result of the reservoir of the lithologic region.
And step 104, acquiring natural gamma curves, acoustic time difference curves and formation deep resistivity curves of a vertical well and a horizontal well in the same reservoir interval according to the position of the well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, and correcting the natural gamma curves, the acoustic time difference curves and the formation deep resistivity curves of the horizontal well in the same reservoir interval.
And 105, determining the water flooding condition of the reservoir section of the horizontal well in the lithologic region according to a preset water flooding qualitative interpretation rule of the horizontal well, a corrected natural gamma curve, a corrected sound wave time difference curve and a corrected stratum deep resistivity curve of the horizontal well.
And 106, selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the lithologic region horizontal well according to the water logging condition of the reservoir section of the lithologic region horizontal well, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the lithologic region horizontal well.
According to the sea facies clastic rock horizontal well reservoir logging interpretation method provided by the embodiment of the invention, firstly, a research area is divided into a plurality of lithological areas with different lithological differences according to a reservoir core test result of a sea facies clastic rock vertical well in the research area and a conventional logging curve of the vertical well, and the reservoir in each lithological area is divided into an original reservoir and a water logging reservoir according to the reservoir water logging condition of each lithological area; then, according to the core data and the logging data of the research area, a core scale logging method is adopted to establish a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the water flooded reservoir of each lithologic area; determining the position of a horizontal well borehole trajectory of a reservoir of a lithologic region in a marine facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic wave time difference curve and resistivity data of the reservoir of the lithologic region and by combining a vertical well stratum dividing result of the reservoir of the lithologic region; according to the position of the well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, acquiring a natural gamma curve, a sonic time difference curve and a stratum deep resistivity curve of the vertical well and the horizontal well in the same reservoir section, and correcting the natural gamma curve, the sonic time difference curve and the stratum deep resistivity curve of the horizontal well in the same reservoir section; determining the water flooding condition of the reservoir section of the horizontal well in the lithologic region according to a preset water flooding qualitative interpretation rule of the horizontal well, a corrected natural gamma curve, a corrected sound wave time difference curve and a corrected stratum deep resistivity curve of the horizontal well; and selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the lithologic region horizontal well according to the water logging condition of the reservoir section of the lithologic region horizontal well, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the lithologic region horizontal well. Therefore, the method can overcome the defects that the horizontal well reservoir is difficult to interpret and no marine clastic rock geological constraint is generated to correct the horizontal well logging curve in the prior art, and can make up the defect of lack of flooding interpretation in the horizontal well reservoir interpretation, so that the result of the marine clastic rock horizontal well reservoir logging interpretation is more accurate.
In order that those skilled in the art will better understand the present invention, a more detailed example is set forth below. It is to be understood that the drawings in the following description are merely exemplary embodiments of the invention, and that the embodiments described are only a few, and not all, embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the scope of the present invention.
As shown in fig. 2, an embodiment of the present invention provides a well logging interpretation method for a marine clastic rock horizontal well reservoir, including:
Here, the closer the core test and conventional logging response of the reservoir are, the weaker the heterogeneity of the reservoir.
Here, the core of the research area and the reservoir with similar conventional logging response (i.e. weak heterogeneity of the reservoir) can be classified into one type, so that the research area can be divided into different lithologic areas.
And 203, dividing the reservoir stratum of each lithologic area into an original reservoir stratum and a water flooding reservoir stratum according to the water flooding condition of the reservoir stratum of each lithologic area.
And 204, dividing a sample layer according to the marine facies clastic rock core well by adopting a core scale logging method according to the core data and the logging data of the research area.
Here, since the longitudinal resolutions of the core data and the logging data are different, in order to reduce errors caused by different resolutions, a method of dividing the readings of the sample layers is required to realize resolution matching between the core data and the logging data. According to the modeling idea of core scale logging, a sample layer can be divided by using a core well, so that the subsequent establishment of a logging interpretation model of an undisturbed reservoir and a water flooded reservoir is completed.
And step 205, establishing shale content logging interpretation models of the undisturbed reservoir and the flooded reservoir of each lithologic area.
Wherein, the mud content well logging interpretation model is as shown in formula (2):
wherein, VshThe unit of the shale content of an undisturbed reservoir and a water flooded reservoir in the rocky region is; GCUR is an empirical coefficient, and in general, an old stratum with larger burial depth can be 2, and a new stratum with shallower burial depth can be 3.7; the delta GR is the natural gamma relative value of the undisturbed reservoir and the water flooded reservoir in the lithologic region; formula (1):
GR is the natural gamma value of the undisturbed reservoir and the flooded reservoir in the lithologic region, and the unit is API; GRminNatural gamma value of pure sandstone; GRmaxIs the natural gamma value of pure mudstone, and has the unit of API.
And step 206, performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the lithologic region and the acoustic time difference logging curve data of the undisturbed reservoir and the water flooded reservoir in the lithologic region, and establishing an acoustic time difference relative porosity model.
The acoustic time difference relative porosity model is phi xAC-y; wherein, phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; AC is the acoustic wave time difference value of the undisturbed reservoir stratum and the water flooded reservoir stratum in the lithologic region, and the unit is mus/ft; x and y are model parameters.
And step 207, performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the lithologic region and the permeability of the undisturbed reservoir and the water flooded reservoir in the lithologic region, and establishing a porosity relative permeability model.
Wherein the porosity relative permeability model is Perm & ltxe & gtyΦ(ii) a Perm is the permeability of an undisturbed reservoir and a flooded reservoir in a lithologic region, and the unit is mD; e is a natural constant; phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; x and y are model parameters.
And step 208, establishing a water saturation model of the undisturbed reservoir and the water flooded reservoir of the lithologic region according to the Archie formula.
Wherein, the water saturation models of the undisturbed reservoir and the flooded reservoir of the lithologic region established by the Archie formula are formula (3):
wherein a is a lithology coefficient related to lithology; b is a constant related to lithology; m is a cementation index; n is a saturation index; rw is the resistivity of the formation water of the undisturbed reservoir and the flooded reservoir in the lithologic region, and the unit is omega.m; rt is the stratum resistivity of an undisturbed reservoir and a water flooded reservoir in a lithologic region, and the unit is omega.m; phi is the porosity of an undisturbed reservoir and a water flooded reservoir in the lithologic region, and the unit is; sw is the water saturation of an undisturbed reservoir and a water flooded reservoir in the lithologic region, and the unit is; the lithology coefficient a related to the lithology, the constant b related to the lithology, the cementation index m, the saturation index n, the formation water resistivity Rw of the undisturbed reservoir and the flooded reservoir of the lithology area are obtained from a preset corresponding relation table; the correspondence table includes the original state of each lithologic regionLithology coefficients related to lithology, constants related to lithology, cementation indexes, saturation indexes and formation water resistivity corresponding to the reservoir and the water flooded reservoir.
And 209, establishing a logging interpretation model of the oil saturation of the undisturbed reservoir and the flooded reservoir of the lithologic region according to the water saturation model.
Wherein, the well logging interpretation model of the oil saturation is shown as a formula (4): so＝1-Sw(ii) a Wherein S isoThe oil saturation of the reservoir is the original reservoir in the rock region and the oil saturation of the reservoir in the water logging region, and the unit is percent.
And step 210, identifying the mutation points of the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir stratum in the rock region according to the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir stratum in the rock region.
And step 211, determining an interface of the overburden and the reservoir of the rock region according to the mutation points of the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir of the rock region.
And 212, determining logging response abnormal points inside the marine facies clastic rocks of the lithological zone according to the horizontal well resistivity data of the reservoir of the lithological zone.
Here, because the difference between the lithology, physical property and oil content of the marine clastic rock of the reservoirs in the overburden and the lithology region is large, the horizontal well logging curve has a mutation point from a high value to a low value at the interface of the overburden and the lithology region. In the marine facies clastic rock, the horizontal well logging curve of the reservoir section is relatively stable, but the logging curve of the non-reservoir section obviously fluctuates, and the position relation between the well track of the horizontal well and the oil reservoir can be determined by searching the mutation point of the stratum interface and the logging response abnormal point in the marine facies clastic rock and combining with the division result of the vertical well stratum.
And 214, acquiring natural gamma curves, acoustic time difference curves and formation deep resistivity curves of the vertical well and the horizontal well in the same reservoir interval according to the position of the well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, and correcting the natural gamma curves, the acoustic time difference curves and the formation deep resistivity curves of the horizontal well in the same reservoir interval.
The sea-phase clastic rock is obvious in wave elutriation, uniform in lithology and physical property, weak in reservoir heterogeneity, and high in correlation between the logging response characteristic values of the vertical well and the logging response characteristic values of the horizontal well in the same reservoir section. Identifying the section with the most homogeneous logging response on the horizontal well section and the vertical well, extracting the average value of the logging response of the section, analyzing the difference of the section and the vertical well in natural gamma, acoustic wave time difference and resistivity curves, establishing a correction template by utilizing the characteristic values of the same reservoir section of the horizontal well and the vertical well, and completing the correction of the logging curve of the horizontal well. The abnormal values of logging response caused by irregular well walls and non-reservoir sections which do not need to be processed need to be removed, so that the reservoir parameters of the horizontal well can be calculated more accurately.
Here, the following formula may be specifically used: ACCorrection ofAnd (3) correcting the acoustic wave time difference curve of a horizontal well of the same reservoir section by α AC + β, wherein AC isCorrection ofThe unit of the corrected horizontal well acoustic wave time difference value is mu s/ft, the unit of the AC is the original acoustic wave time difference value of the horizontal well and is mu s/ft, and α and β are correction parameters.
According to the formula: rt correction＝αRt+ β the formation deep resistivity curve of the horizontal well in the same reservoir section is corrected, wherein Rt correctionThe unit of the horizontal well stratum resistivity value after correction is omega m; rtThe resistivity values of the original stratum of the horizontal well are shown in omega m, and α and β are correction parameters.
Here, the preset horizontal well flooding qualitative interpretation rule includes each flooding condition, and a natural gamma value range and a characteristic value, a sonic time difference value range and a characteristic value, and a resistivity value range and a characteristic value corresponding to each flooding condition.
And step 216, selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the lithologic region horizontal well according to the water logging condition of the reservoir section of the lithologic region horizontal well, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the lithologic region horizontal well.
And comprehensively and qualitatively explaining, quantitatively calculating parameters and a horizontal well track and oil reservoir relation graph, and determining the comprehensive explanation of the horizontal well section of the marine facies clastic rock in the research area by taking actual production data as verification.
In order to make the above steps 201 to 216 clearer, a specific application example of the above steps 201 to 216 is listed below:
reservoir characteristic analysis is carried out on the Hardson oil field east river sandstone in the Tarim basin of China through the steps 201 to 216, high-precision reservoir logging interpretation models of an undisturbed reservoir and a water flooded reservoir in different lithologic areas are established, and the position of the horizontal well HD11-8H borehole trajectory in the reservoir is researched. Meanwhile, taking a lithologic area (such as a northwest area) as an example, horizontal well curve correction is carried out, a horizontal well water logging qualitative interpretation standard is established, and horizontal well reservoir interpretation is carried out on HD 11-8H.
The Hadison east river sandstone in the oilfield is a typical marine facies clastic rock reservoir which is a barrier-free wave-controlled sandy shore sediment, the lithology of the east river sandstone section is gray fine-grain quartz sandstone, the component maturity and the structural maturity are high, the lithology and physical properties are uniform, and the reservoirs drilled by a vertical well and a horizontal well have strong similarity, so that the research result of the well logging explanation of the vertical well reservoir can be fully used in the well logging explanation of the horizontal well.
The Hadison oilfield depositional features are compartmentalized, appearing as differences in lithology. By integrating the rock core and the conventional logging response characteristics, the east river sandstone reservoir area can be divided into 4 lithologic areas, which are respectively: northwest, central, southwest and southeast. Meanwhile, according to whether the reservoir in the lithologic area is flooded with water or not, the reservoir in the lithologic area can be divided into an undisturbed reservoir and a flooded reservoir.
Secondly, by a core scale logging method, comprehensively researching core data and logging data of the area, a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of an undisturbed reservoir and a water flooded reservoir of the four rock areas can be established.
The natural gamma value GR of the target interval of the research area can accurately reflect the change of the shale content, so the shale content is calculated by adopting a natural gamma curve, and the calculation formulas are the formula (1) and the formula (2).
Dividing sample layers by using a core well, testing the porosity according to the acoustic time difference value and the corresponding core, and establishing an acoustic time difference relative porosity calculation model, namely that phi is xAC-y, specifically to an original reservoir and a water flooded reservoir of each rock region, the following steps are performed:
northwest undisturbed reservoir: phi 0.7577AC-40.555 equation (5)
Water flooded reservoir in northwest region: phi 0.3991AC-14.335 formula (6)
Undisturbed reservoir in southeast east region: phi 0.1709AC +4.1951 equation (7)
Water flooded reservoir in southeast region: phi 0.9381AC-49.627 equation (8)
Central zone undisturbed reservoir: phi 0.3233AC-7.8854 formula (9)
Water flooded reservoir in central area: phi 0.1265AC +6.698 equation (10)
Southwest undisturbed reservoir: phi 0.7516AC-37.104 formula (11)
Water flooded reservoir in southwest area: phi 0.8103AC-41.352 equation (12)
According to the single correlation analysis of main influence factors such as core analysis porosity, shale content, carbonate content and the like and permeability, the single correlation between the permeability and the porosity is the best, and a relative permeability calculation model of the reservoir and the porosity is established, namely Perm & ltxe & gtyΦSpecifically, the original reservoir and the water flooded reservoir in each lithologic region are as follows:
northwest undisturbed reservoir: perm 0.034e0.491ΦFormula (13)
Water flooded reservoir in northwest region: perm＝0.025e0.473ΦFormula (14)
Undisturbed reservoir in southeast east region: perm ═ 0.119e0.443ΦFormula (15)
Water flooded reservoir in southeast region: perm ═ 0.041e0.449ΦFormula (16)
Central zone undisturbed reservoir: perm ═ 0.052e0.445ΦFormula (17)
Water flooded reservoir in central area: perm ═ 0.057e0.435ΦFormula (18)
Southwest undisturbed reservoir: perm ═ 0.089e0.43ΦFormula (19)
Water flooded reservoir in southwest area: perm 0.184e0.399ΦFormula (20)
And (3) establishing a water saturation model of the undisturbed reservoir and the water flooded reservoir in the rock region by adopting an Archie formula, wherein the formula is the formula (3) and the formula (4), and parameters a, b, m, n and Rw of each rock region can be shown in the following table 1. Table 1: a corresponding relation table:
3, a set of stable middle mudstone sections (overlying strata) are deposited on the upper section of the east river sandstone in the research area, the characteristics that the natural gamma value is greater than 80API, and the acoustic wave time difference value is greater than 80 mus/ft are shown, the difference with the east river sandstone (the natural gamma GR is 33-38 API, the acoustic wave time difference AC is 60-72 mus/ft, and the resistivity is 1-10 omega.m) is large, and the natural gamma value and the acoustic wave time difference value of a horizontal well logging curve on the interface of the two sections have obvious mutation points from a high value to a low value. In addition, the existence of the non-reservoir section causes the abnormal logging response of the horizontal well section, the non-reservoir section has a large resistivity value (more than 10 omega-m) and obvious curve fluctuation on a horizontal well logging curve, and therefore the position relation between the well track of the horizontal well and the oil reservoir can be determined by searching the mutation point of the stratum interface and the logging response abnormal point inside the marine facies clastic rock.
FIG. 3 is a graph of the position relationship between the HD11-8H wellbore trajectory and the reservoir. The upper part of the Donghe sandstone is a middle mudstone section, the lower part of the Donghe sandstone is a reservoir formation, wherein the HD11-8 is a straight well section of the horizontal well HD11-8H, a very obvious GR abrupt change point from a high value to a low value can be seen at 5200m of the HD11-8H, a vertical line of a well track is made from the abrupt change point, the intersection point of the well track and the well track is an interface of the middle mudstone section and the Donghe sandstone, and the position of the HD11-8H well track in an oil reservoir can be determined by combining the Donghe sandstone stratum division results of the HD11-8 and the HD 112. The Donghe sandstone accounts for 5 small layers, the HD11-8H is wholly positioned in 1 small layer, and meanwhile, the logging curves of the HD11-8H at sections 5284-5349 m, 5382-5450 m and 5461-5492 m are relatively stable and are expressed as the characteristics of a reservoir section; the curve fluctuation is obvious in the sections 5200 to 5284m, 5349 to 5382m and 5450 to 5461m, and the curve fluctuation is characterized by a non-reservoir section.
And 4, respectively counting the characteristic values of natural gamma curves, acoustic wave time difference curves and deep resistivity curves of the same reservoir section of the vertical well and the horizontal well, wherein the extracted characteristic values are the average values of the well logging curves of the reservoir section with the most homogeneous horizontal well and the vertical well, such as the extracted average values of the well logging curves of 5409-5421 m parts of HD11-8H in the northwest region and the extracted average values of the well logging curves of 5103.4-5103.9 m parts of HD11-8 corresponding to the same reservoir section.
Fig. 4 is a graph comparing the natural gamma values of a vertical well reservoir section and a horizontal well reservoir section of the hartson oil field hardson in the talimu basin. Because natural gamma is a radioactive logging method, the influence of the borehole direction on a logging curve is small, and meanwhile, the fact that the value of the natural gamma curve of the horizontal well is basically close to that of the vertical well can be found without correction.
The acoustic wave time difference value of the horizontal well is generally higher than that of the vertical well, and curve correction is needed. Fig. 5 is a schematic diagram comparing acoustic time difference curves of a straight well and a horizontal well in northwest region of east river sandstone of Hadison oil field in Tarim basin, and the correction formula is as follows:
ACcorrection ofα AC + β, here, then:
ACcorrection of0.6639AC +32.077 equation (21)
Horizontal wells generally have higher resistivity values than vertical wells and need to be corrected by curves. FIG. 6 is a schematic diagram showing the comparison of resistivity curves of a straight well and a horizontal well in northwest region of east river sandstone in Hadison oil field in Tarim basin, and the correction formula is Rt correction＝αRt+ β, here, then:
Rt correction＝0.9879Rt+0.2515 formula (22)
The horizontal well logging curve correction is completed through the correction template, non-reservoir sections which do not need to be processed and logging response abnormal values caused by well wall irregularity, such as 5200-5284 m, 5349-5382 m and 5450-5461 m sections of HD11-8H in northwest, including 5276-5284 m sound wave time difference abnormal response sections, need to be removed in the correction process.
And 5, counting main value intervals and characteristic values of curves such as natural gamma, acoustic wave time difference, resistivity and the like of an undisturbed reservoir of the horizontal well and different water flooding levels on the basis of the actual production condition, namely forming the preset qualitative interpretation rule of the water flooding of the horizontal well, and displaying the interpretation rule in a table form, wherein the table 2 shows that whether the horizontal well is flooded or not and the water flooding level can be directly judged by reading the logging response value of the horizontal well.
Table 2 qualitative interpretation criteria for horizontal well flooding:
comparing the oil layer, the poor oil layer and each water flooded reservoir, finding that the natural gamma value of the poor oil layer has larger fluctuation and is mainly influenced by the development condition of the surrounding non-reservoir, and the natural gamma value after water flooding is in a decreasing trend. Generally, the acoustic wave time difference curve of a reservoir after flooding is in an increasing trend, and the resistivity is in a significantly decreasing trend, because the Hadamard area adopts sewage reinjection, the degree of mineralization of injected water is similar to that of formation water, and the resistivity of an oil layer is decreased and gradually approaches to a water layer along with the increase of the flooding degree.
Sixthly, FIG. 7 is a well logging interpretation result diagram of a horizontal well HD11-8H reservoir. On the basis of completing the horizontal well logging curve correction, performing water-logging qualitative analysis on HD11-8H according to the table 2, and finding that 5284-5302 m and 5409-5450 m are oil layer sections, 5302-5349 m, 5382-5409 m and 5461-5492 m are poor oil layer sections, the rest 5200-5284 m, 5349-5382 m and 5450-5461 m are non-reservoir section, wherein the water-logging layer section does not exist, the whole horizontal section reservoir is 174m and accounts for 59.6%, and the high-quality reservoir (namely oil layer) section is only 69m and accounts for 23.6%.
Reservoir parameters are calculated in reservoir sections (oil layers and poor oil layers), the shale content is calculated by using formulas (1) and (2), the porosity is calculated by using a formula (5), the permeability is calculated by using a formula (13), the oil saturation is calculated by using formulas (3) and (4), wherein a, b, m, n and Rw adopt the numerical values of original reservoirs in northwest regions, and thus the well logging interpretation of the horizontal well reservoir is completed.
In the production characteristics, the initial yield of HD11-8H is 42t/d, while the average initial yield of a horizontal well in a research area is higher than 100t/d, which is caused by the low drilling rate of a high-quality reservoir (namely an oil layer) of the well; meanwhile, the well does not contain water in the initial production stage, which means that the well does not drill a water flooded layer or a water layer, and the interpretation conclusion of oil and water is consistent, so that the interpretation method for the sea-facies clastic rock horizontal well reservoir logging provided by the invention is high in precision and can meet the production requirement.
Through the specific examples, the advantages of the invention compared with the prior art can be found:
the idea of horizontal well reservoir interpretation based on marine facies clastic rock geological features and vertical well reservoir research is creatively provided. The method is based on the core and the conventional logging curve, the basic characteristics and classification of the reservoir are researched, and the defect that the reservoir interpretation of the horizontal well is restricted due to lack of geological knowledge in the existing method is overcome; by adopting a core scale logging method, high-precision reservoir logging interpretation models of an undisturbed reservoir and a flooded reservoir in different lithologic regions are established, and the defect that only one interpretation model is used for all reservoirs in the whole research region in the existing method is overcome; the horizontal well geological interface logging response characteristic research is carried out, the response characteristics of a stratum interface and a reservoir section on a horizontal well curve are determined, the passing condition of a well track in a single sand body can be revealed, and the defect that the position relation between a well track of a horizontal well and an oil reservoir is fuzzy in the existing method is overcome. The invention carries out horizontal well reservoir interpretation based on geological constraints and can meet the production requirement.
The method for correcting the logging curve of the horizontal well is creatively provided based on marine facies clastic rock geological constraints and the position relation between the horizontal well and the oil reservoir to extract logging response characteristic values of the same reservoir section of the horizontal well and the vertical well, and establish a chart to finish the correction of the logging curve of the horizontal well. The method has the advantages that the vertical well sections corresponding to the horizontal sections of the same reservoir are picked up, the influence of non-reservoir sections on curves is eliminated, and the defect that the whole horizontal section is interpreted as the same reservoir in the existing method is overcome; identifying a section of the horizontal well section and the vertical well which is corresponding to the most homogeneous logging, extracting the average value of logging response of the section, and comprehensively researching the response difference of the section of the logging response and the average value of logging response of the section of the horizontal well section and the vertical well to lay a foundation for correcting the horizontal well curve of the vertical well; and a correction template is established, the logging curve is corrected, logging response abnormal values caused by irregular well walls are eliminated, and the defect that the horizontal well logging curve is difficult to correct due to large difference with a vertical well environment is overcome. The invention provides a set of horizontal well reservoir interpretation method, which improves the accuracy of horizontal well reservoir interpretation.
The method creatively establishes a marine facies clastic rock horizontal well water logging qualitative interpretation standard, performs horizontal well reservoir interpretation, and completes horizontal well reservoir evaluation and fluid identification. Logging response characteristic values of different flooding levels are extracted, a horizontal well flooding qualitative interpretation standard is established, whether the horizontal well is flooded or not and the flooding level can be directly judged by reading the horizontal well logging response values, and the defect of lack of flooding interpretation in horizontal well reservoir interpretation is overcome; quantitative parameters such as shale content, porosity, permeability and oil saturation can be calculated by using a logging interpretation model of a vertical well, and the defect that reservoir interpretation and fluid identification are difficult to perform due to the fact that a logging interpretation model cannot be established due to the lack of core data of a horizontal well is overcome.
The embodiment of the invention can be combined with geological research, horizontal well reservoir interpretation is carried out based on the geological research and the position relation between the horizontal well and the oil reservoir, the multi-solution property of the horizontal well reservoir logging interpretation is reduced, the horizontal well reservoir evaluation and the fluid identification can be completed by using the corrected horizontal segment logging curve, and the horizontal well reservoir evaluation and the fluid identification can be verified through actual production data of the horizontal well, so that the interpretation precision is high, and the method can be used in production practice.
Corresponding to the method embodiment described in fig. 1 and fig. 2, as shown in fig. 8, an embodiment of the present invention provides a device for explaining well logging in a marine clastic rock horizontal well reservoir, including:
the research area dividing unit 31 is configured to divide the research area into a plurality of lithological areas with different lithological differences according to a reservoir core test result of the marine facies clastic rock vertical well in the research area and a conventional logging curve of the vertical well, and divide the reservoir in each lithological area into an original reservoir and a water flooded reservoir according to a reservoir water flooding condition in each lithological area.
And the logging interpretation model establishing unit 32 is used for establishing a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the flooded reservoir of each lithologic area by adopting a core scale logging method according to the core data and the logging data of the research area.
And the horizontal well track position determining unit 33 is used for determining the position of the horizontal well track of the reservoir of the rock region in the sea facies clastic rock reservoir according to the horizontal well natural gamma curve, the acoustic wave time difference curve and the resistivity data of the reservoir of the rock region and by combining the vertical well stratum dividing result of the reservoir of the rock region.
And the curve data correction unit 34 is configured to acquire a natural gamma curve, a sonic time difference curve and a formation deep resistivity curve of the vertical well and the horizontal well in the same reservoir interval according to the position of the well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, and correct the natural gamma curve, the sonic time difference curve and the formation deep resistivity curve of the horizontal well in the same reservoir interval.
And the flooding condition determining unit 35 is configured to determine a flooding condition of the reservoir section of the horizontal well in the lithologic region according to a preset horizontal well flooding qualitative interpretation rule, the corrected natural gamma curve, the corrected acoustic wave time difference curve, and the corrected stratum depth resistivity curve.
And the horizontal well logging interpretation unit 36 is configured to select, according to the flooding condition of the reservoir section of the horizontal well in the lithologic region, a logging interpretation model corresponding to the flooding condition of the reservoir section of the horizontal well in the lithologic region, and determine the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the horizontal well in the lithologic region.
Furthermore, the study region dividing unit 31 is specifically configured to:
and determining the heterogeneity of the marine clastic rock reservoir in the research area according to the reservoir core test result of the marine clastic rock vertical well in the research area and the conventional logging curve of the vertical well.
According to the heterogeneity of the marine clastic rock reservoir in the research area, the research area is divided into a plurality of lithological areas with different lithological differences.
And dividing the reservoir stratum of each lithologic area into an original reservoir stratum and a water flooding reservoir stratum according to the water flooding condition of the reservoir stratum of each lithologic area.
In addition, the well logging interpretation model building unit 32 is specifically configured to:
and dividing a sample layer according to the marine clastic rock core well.
Establishing a shale content logging interpretation model of an undisturbed reservoir and a water flooded reservoir of each lithologic area; the mud content well logging interpretation model is
Wherein, VshThe shale content of an undisturbed reservoir and a water flooded reservoir in a lithologic region; GCUR is an empirical coefficient; the delta GR is the natural gamma relative value of the undisturbed reservoir and the water flooded reservoir in the lithologic region;
GR is the natural gamma value of an undisturbed reservoir and a water flooded reservoir in a lithologic area; GRminNatural gamma value of pure sandstone; GRmaxIs the natural gamma value of pure mudstone.
Performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the rocky region and the acoustic time difference logging curve data of the undisturbed reservoir and the water flooded reservoir in the rocky region, and establishing an acoustic time difference relative porosity model; the acoustic time difference relative porosity model is phi xAC-y; wherein, phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; AC is the acoustic wave time difference value of the undisturbed reservoir and the water flooded reservoir in the rock region; x and y are model parameters.
Performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the lithologic region and the permeability of the undisturbed reservoir and the water flooded reservoir in the lithologic region, and establishing a porosity relative permeability model; the porosity relative permeability model is Perm & ltxe & gtyΦ(ii) a Perm is the permeability of an undisturbed reservoir and a water flooded reservoir in a lithologic region; e is a natural constant; phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; x and y are model parameters.
Establishing a water saturation model of an undisturbed reservoir and a water flooded reservoir in a lithologic region according to an Archie formula:
wherein a is a lithology coefficient related to lithology; b is a constant related to lithology; m is a cementation index; n is a saturation index; rw is the resistivity of the formation water of the undisturbed reservoir and the flooded reservoir in the rock region; rt is the stratum resistivity of an undisturbed reservoir and a flooded reservoir in the lithologic region; phi is the porosity of an undisturbed reservoir and a water flooded reservoir in the lithologic region; sw is the water saturation of an undisturbed reservoir and a water flooded reservoir in the lithologic region; the lithology coefficient a related to the lithology, the constant b related to the lithology, the cementation index m, the saturation index n, the formation water resistivity Rw of the undisturbed reservoir and the flooded reservoir of the lithology area are obtained from a preset corresponding relation table; the corresponding relation table comprises lithology coefficients, constants, cementation indexes, saturation indexes and formation water resistivity, which are corresponding to the undisturbed reservoir and the water flooded reservoir of each lithologic area and are related to lithology, and the constants, the cementation indexes, the saturation indexes and the formation water resistivity are related to the lithology.
And establishing a logging interpretation model of the oil saturation of the undisturbed reservoir and the flooded reservoir in the lithologic region according to the water saturation model: so＝1-Sw(ii) a Wherein S isoThe oil saturation of the reservoir is the original reservoir in the rock region and the oil saturation of the reservoir flooded with water.
In addition, the horizontal well trajectory position determination unit 33 is specifically configured to:
and identifying the mutation points of the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region according to the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region.
And determining the interface of the overlying stratum and the reservoir stratum of the rock region according to the mutation points of the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir stratum of the rock region.
And determining logging response abnormal points inside the marine facies clastic rocks of the lithological zone according to the resistivity data of the horizontal wells of the reservoir of the lithological zone.
And determining the position of the horizontal well track of the reservoir of the lithologic region in the marine facies clastic rock reservoir according to the interface of the overlying stratum and the reservoir of the lithologic region and the logging response abnormal point inside the marine facies clastic rock of the lithologic region by combining the vertical well stratum dividing result of the reservoir of the lithologic region.
Further, the curve data correcting unit 34 is specifically configured to;
according to the formula: ACCorrection ofAnd (3) correcting the acoustic wave time difference curve of a horizontal well of the same reservoir section by α AC + β, wherein AC isCorrection ofThe corrected horizontal well sound wave time difference value, AC is the original sound wave time difference value of the horizontal well, and α and β are correction parameters.
According to the formula: rt correction＝αRt+ β the formation deep resistivity curve of the horizontal well in the same reservoir section is corrected, wherein Rt correctionThe horizontal well stratum resistivity value after correction is obtained; rtWhich is the original formation resistivity value of the horizontal well, α and β are the calibration parameters.
In addition, the preset horizontal well flooding qualitative interpretation rule in the flooding condition determining unit 35 includes a natural gamma value range and a characteristic value, a sound wave time difference value range and a characteristic value, and a resistivity value range and a characteristic value corresponding to each flooding condition and each flooding condition.
The flooding condition determining unit 35 is specifically configured to: and inquiring the corrected water flooding conditions corresponding to the natural gamma curve, the acoustic wave time difference curve and the stratum deep resistivity curve of the horizontal well in the preset water flooding qualitative interpretation rule of the horizontal well, and identifying that the water flooding condition of the reservoir section of the horizontal well in the lithologic region is an oil layer or a reservoir layer of a poor oil layer.
According to the device for explaining the logging of the reservoir of the marine clastic rock horizontal well, firstly, a research area is divided into a plurality of lithological areas with different lithological differences according to a reservoir core test result of a marine clastic rock vertical well in the research area and a conventional logging curve of the vertical well, and the reservoir in each lithological area is divided into an original reservoir and a water-flooded reservoir according to the reservoir water-flooded condition of each lithological area; then, according to the core data and the logging data of the research area, a core scale logging method is adopted to establish a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the water flooded reservoir of each lithologic area; determining the position of a horizontal well borehole trajectory of a reservoir of a lithologic region in a marine facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic wave time difference curve and resistivity data of the reservoir of the lithologic region and by combining a vertical well stratum dividing result of the reservoir of the lithologic region; according to the position of the well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, acquiring a natural gamma curve, a sonic time difference curve and a stratum deep resistivity curve of the vertical well and the horizontal well in the same reservoir section, and correcting the natural gamma curve, the sonic time difference curve and the stratum deep resistivity curve of the horizontal well in the same reservoir section; determining the water flooding condition of the reservoir section of the horizontal well in the lithologic region according to a preset water flooding qualitative interpretation rule of the horizontal well, a corrected natural gamma curve, a corrected sound wave time difference curve and a corrected stratum deep resistivity curve of the horizontal well; and selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the lithologic region horizontal well according to the water logging condition of the reservoir section of the lithologic region horizontal well, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the lithologic region horizontal well. Therefore, the method can overcome the defects that the horizontal well reservoir is difficult to interpret and no marine clastic rock geological constraint is generated to correct the horizontal well logging curve in the prior art, and can make up the defect of lack of flooding interpretation in the horizontal well reservoir interpretation, so that the result of the marine clastic rock horizontal well reservoir logging interpretation is more accurate.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
Claims (10)
1. A well logging interpretation method for a marine clastic rock horizontal well reservoir is characterized by comprising the following steps:
dividing a research area into a plurality of lithological areas with different lithological differences according to a reservoir rock core test result of a marine facies clastic rock vertical well in the research area and a conventional logging curve of the vertical well, and dividing the reservoir in each lithological area into an original reservoir and a water flooded reservoir according to the reservoir water flooding condition of each lithological area;
according to the core data and the logging data of the research area, a core scale logging method is adopted to establish a logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of the undisturbed reservoir and the water flooded reservoir of each lithologic area;
determining the position of a horizontal well borehole trajectory of a reservoir of a lithologic region in a marine facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic wave time difference curve and resistivity data of the reservoir of the lithologic region and by combining a vertical well stratum dividing result of the reservoir of the lithologic region;
according to the position of the well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, acquiring a natural gamma curve, a sonic time difference curve and a stratum deep resistivity curve of the vertical well and the horizontal well in the same reservoir section, and correcting the natural gamma curve, the sonic time difference curve and the stratum deep resistivity curve of the horizontal well in the same reservoir section;
determining the water flooding condition of the reservoir section of the horizontal well in the lithologic region according to a preset water flooding qualitative interpretation rule of the horizontal well, a corrected natural gamma curve, a corrected sound wave time difference curve and a corrected stratum deep resistivity curve of the horizontal well;
according to the water flooding condition of the reservoir section of the horizontal well of the lithologic region, selecting a logging interpretation model corresponding to the water flooding condition of the reservoir section of the horizontal well of the lithologic region, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the horizontal well of the lithologic region;
determining the position of a horizontal well borehole trajectory of a reservoir in a lithologic region in a marine facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic moveout curve and resistivity data of the reservoir in the lithologic region and by combining a vertical well stratigraphic division result of the reservoir in the lithologic region, wherein the determining comprises the following steps:
identifying the mutation points of the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region according to the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region;
determining an interface of an overlying stratum and a reservoir stratum of the rock region according to mutation points of a horizontal well natural gamma curve and an acoustic wave time difference curve of the reservoir stratum of the rock region;
determining logging response abnormal points inside marine facies clastic rocks of the lithological zone according to the horizontal well resistivity data of the reservoir of the lithological zone;
and determining the position of the horizontal well track of the reservoir of the lithologic region in the marine facies clastic rock reservoir according to the interface of the overlying stratum and the reservoir of the lithologic region and the logging response abnormal point inside the marine facies clastic rock of the lithologic region by combining the vertical well stratum dividing result of the reservoir of the lithologic region.
2. The method for explaining well logging of the sea facies clastic rock horizontal well reservoir according to claim 1, wherein the research area is divided into a plurality of lithological areas with different lithological differences according to the reservoir core test result of the sea facies clastic rock vertical well in the research area and the conventional logging curve of the vertical well, and the reservoir in each lithological area is divided into an original reservoir and a water flooded reservoir according to the reservoir water flooded condition in each lithological area, comprising:
determining the heterogeneity of the marine clastic rock reservoir in the research area according to the reservoir core test result of the marine clastic rock vertical well in the research area and the conventional logging curve of the vertical well;
dividing a research area into a plurality of lithologic areas with different lithologic differences according to the heterogeneity of marine facies clastic rock reservoirs in the research area;
and dividing the reservoir stratum of each lithologic area into an original reservoir stratum and a water flooding reservoir stratum according to the water flooding condition of the reservoir stratum of each lithologic area.
3. The method for explaining logging of the sea facies clastic rock horizontal well reservoir according to claim 1, wherein a logging explanation model of the shale content, the porosity, the permeability and the oil saturation of an undisturbed reservoir and a flooded reservoir of each lithologic region is established by adopting a core scale logging method according to core data and logging data of a research region, and comprises the following steps:
dividing a sample layer according to the marine facies clastic rock core well;
establishing a shale content logging interpretation model of an undisturbed reservoir and a water flooded reservoir of each lithologic area; the mud content well logging interpretation model is
Wherein, VshThe shale content of an undisturbed reservoir and a water flooded reservoir in a lithologic region; GCUR is an empirical coefficient; the delta GR is the natural gamma relative value of the undisturbed reservoir and the water flooded reservoir in the lithologic region;
GR is the natural gamma value of an undisturbed reservoir and a water flooded reservoir in a lithologic area; GRminNatural gamma value of pure sandstone; GRmaxNatural gamma value of pure mudstone;
performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the rocky region and the acoustic time difference logging curve data of the undisturbed reservoir and the water flooded reservoir in the rocky region, and establishing an acoustic time difference relative porosity model; the acoustic time difference relative porosity model is phi xAC-y; wherein, phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; AC is the acoustic wave time difference value of the undisturbed reservoir and the water flooded reservoir in the rock region; x and y are model parameters;
analyzing the porosity of the undisturbed reservoir and the core of the water flooded reservoir of the lithologic area and the original reservoir and the water flooded reservoir of the lithologic areaPerforming single correlation analysis on the permeability, and establishing a porosity relative permeability model; the porosity relative permeability model is Perm & ltxe & gtyΦ(ii) a Perm is the permeability of an undisturbed reservoir and a water flooded reservoir in a lithologic region; e is a natural constant; phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; x and y are model parameters;
establishing a water saturation model of an undisturbed reservoir and a water flooded reservoir in a lithologic region according to an Archie formula:
wherein a is a lithology coefficient related to lithology; b is a constant related to lithology; m is a cementation index; n is a saturation index; rw is the resistivity of the formation water of the undisturbed reservoir and the flooded reservoir in the rock region; rt is the stratum resistivity of an undisturbed reservoir and a flooded reservoir in the lithologic region; phi is the porosity of an undisturbed reservoir and a water flooded reservoir in the lithologic region; sw is the water saturation of an undisturbed reservoir and a water flooded reservoir in the lithologic region; the lithology coefficient a related to the lithology, the constant b related to the lithology, the cementation index m, the saturation index n, the formation water resistivity Rw of the undisturbed reservoir and the flooded reservoir of the lithology area are obtained from a preset corresponding relation table; the corresponding relation table comprises lithology coefficients, constants, cementation indexes, saturation indexes and formation water resistivity, which are corresponding to the undisturbed reservoir and the water flooded reservoir of each lithologic area and are related to lithology;
and establishing a logging interpretation model of the oil saturation of the undisturbed reservoir and the flooded reservoir in the lithologic region according to the water saturation model: so＝1-Sw(ii) a Wherein S isoThe oil saturation of the reservoir is the original reservoir in the rock region and the oil saturation of the reservoir flooded with water.
4. The sea facies clastic rock horizontal well reservoir logging interpretation method is characterized in that according to the position of a horizontal well borehole trajectory of a reservoir of the lithologic region in a sea facies clastic rock reservoir, natural gamma curves, acoustic time difference curves and formation depth resistivity curves of a vertical well and a horizontal well in the same reservoir section are obtained, and the natural gamma curve, the acoustic time difference curve and the formation depth resistivity curve of the horizontal well in the same reservoir section are corrected;
according to the formula: ACCorrection ofAnd (3) correcting the acoustic wave time difference curve of a horizontal well of the same reservoir section by α AC + β, wherein AC isCorrection ofThe corrected horizontal well sound wave time difference value, AC is the original sound wave time difference value of the horizontal well, α and β are correction parameters;
according to the formula: rt correction＝αRt+ β the formation deep resistivity curve of the horizontal well in the same reservoir section is corrected, wherein Rt correctionThe horizontal well stratum resistivity value after correction is obtained; rtWhich is the original formation resistivity value of the horizontal well, α and β are the calibration parameters.
5. The sea facies clastic rock horizontal well reservoir logging interpretation method of claim 4, wherein the preset horizontal well water flooding qualitative interpretation rules include natural gamma value ranges and characteristic values, acoustic wave time difference value ranges and characteristic values, and resistivity value ranges and characteristic values corresponding to each water flooding condition and each water flooding condition;
the method for determining the water logging condition of the reservoir section of the horizontal well in the lithologic region according to the preset water logging qualitative interpretation rule of the horizontal well, the corrected natural gamma curve, the corrected sound wave time difference curve and the corrected stratum deep resistivity curve of the horizontal well comprises the following steps:
and inquiring the corrected water flooding conditions corresponding to the natural gamma curve, the acoustic wave time difference curve and the stratum deep resistivity curve of the horizontal well in the preset water flooding qualitative interpretation rule of the horizontal well, and identifying that the water flooding condition of the reservoir section of the horizontal well in the lithologic region is an oil layer or a reservoir layer of a poor oil layer.
6. The utility model provides a sea facies clastic rock horizontal well reservoir logging interpretation device which characterized in that includes:
the device comprises a research area dividing unit, a data processing unit and a data processing unit, wherein the research area dividing unit is used for dividing a research area into a plurality of lithological areas with different lithological differences according to a reservoir rock core test result of a marine facies clastic rock vertical well in the research area and a conventional logging curve of the vertical well, and dividing the reservoir in each lithological area into an original reservoir and a water logging reservoir according to the reservoir water logging condition of each lithological area;
the well logging interpretation model establishing unit is used for establishing a well logging interpretation model of the shale content, the porosity, the permeability and the oil saturation of an undisturbed reservoir and a flooded reservoir of each lithologic area by adopting a core scale well logging method according to core data and well logging data of a research area;
the horizontal well track position determining unit is used for determining the position of the horizontal well track of the reservoir of the rock region in the sea facies clastic rock reservoir according to a horizontal well natural gamma curve, an acoustic wave time difference curve and resistivity data of the reservoir of the rock region and by combining a vertical well stratum dividing result of the reservoir of the rock region;
the curve data correction unit is used for acquiring natural gamma curves, acoustic time difference curves and stratum deep resistivity curves of a vertical well and a horizontal well in the same reservoir interval according to the position of a well track of the horizontal well of the reservoir in the lithologic region in the marine facies clastic rock reservoir, and correcting the natural gamma curves, the acoustic time difference curves and the stratum deep resistivity curves of the horizontal well in the same reservoir interval;
the water flooding condition determining unit is used for determining the water flooding condition of the reservoir section of the horizontal well in the lithologic region according to the preset water flooding qualitative interpretation rule of the horizontal well, the corrected natural gamma curve, the corrected sound wave time difference curve and the corrected stratum deep resistivity curve of the horizontal well;
the horizontal well logging interpretation unit is used for selecting a logging interpretation model corresponding to the water logging condition of the reservoir section of the horizontal well of the lithologic region according to the water logging condition of the reservoir section of the horizontal well of the lithologic region, and determining the shale content, the porosity, the permeability and the oil saturation of the reservoir section of the horizontal well of the lithologic region;
the horizontal well borehole trajectory position determination unit is specifically configured to:
identifying the mutation points of the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region according to the horizontal well natural gamma curve and the acoustic wave time difference curve of the reservoir in the lithologic region;
determining an interface of an overlying stratum and a reservoir stratum of the rock region according to mutation points of a horizontal well natural gamma curve and an acoustic wave time difference curve of the reservoir stratum of the rock region;
determining logging response abnormal points inside marine facies clastic rocks of the lithological zone according to the horizontal well resistivity data of the reservoir of the lithological zone;
and determining the position of the horizontal well track of the reservoir of the lithologic region in the marine facies clastic rock reservoir according to the interface of the overlying stratum and the reservoir of the lithologic region and the logging response abnormal point inside the marine facies clastic rock of the lithologic region by combining the vertical well stratum dividing result of the reservoir of the lithologic region.
7. The device for explaining well logging of the marine clastic rock horizontal well reservoir according to claim 6, wherein the study region dividing unit is specifically configured to:
determining the heterogeneity of the marine clastic rock reservoir in the research area according to the reservoir core test result of the marine clastic rock vertical well in the research area and the conventional logging curve of the vertical well;
dividing a research area into a plurality of lithologic areas with different lithologic differences according to the heterogeneity of marine facies clastic rock reservoirs in the research area;
and dividing the reservoir stratum of each lithologic area into an original reservoir stratum and a water flooding reservoir stratum according to the water flooding condition of the reservoir stratum of each lithologic area.
8. The device for explaining logging of the marine clastic rock horizontal well reservoir according to claim 6, wherein the logging explanation model establishing unit is specifically configured to:
dividing a sample layer according to the marine facies clastic rock core well;
establishing a shale content logging interpretation model of an undisturbed reservoir and a water flooded reservoir of each lithologic area; the mud content well logging interpretation model is
Wherein, VshThe shale content of an undisturbed reservoir and a water flooded reservoir in a lithologic region; GCUR is experienceA coefficient; the delta GR is the natural gamma relative value of the undisturbed reservoir and the water flooded reservoir in the lithologic region;
GR is the natural gamma value of an undisturbed reservoir and a water flooded reservoir in a lithologic area; GRminNatural gamma value of pure sandstone; GRmaxNatural gamma value of pure mudstone;
performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the rocky region and the acoustic time difference logging curve data of the undisturbed reservoir and the water flooded reservoir in the rocky region, and establishing an acoustic time difference relative porosity model; the acoustic time difference relative porosity model is phi xAC-y; wherein, phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; AC is the acoustic wave time difference value of the undisturbed reservoir and the water flooded reservoir in the rock region; x and y are model parameters;
performing single correlation analysis on the core analysis porosity of the undisturbed reservoir and the water flooded reservoir in the lithologic region and the permeability of the undisturbed reservoir and the water flooded reservoir in the lithologic region, and establishing a porosity relative permeability model; the porosity relative permeability model is Perm & ltxe & gtyΦ(ii) a Perm is the permeability of an undisturbed reservoir and a water flooded reservoir in a lithologic region; e is a natural constant; phi is the porosity of an undisturbed reservoir and a water flooded reservoir of the lithologic region; x and y are model parameters;
establishing a water saturation model of an undisturbed reservoir and a water flooded reservoir in a lithologic region according to an Archie formula:
wherein a is a lithology coefficient related to lithology; b is a constant related to lithology; m is a cementation index; n is a saturation index; rw is the resistivity of the formation water of the undisturbed reservoir and the flooded reservoir in the rock region; rt is the stratum resistivity of an undisturbed reservoir and a flooded reservoir in the lithologic region; phi is the porosity of an undisturbed reservoir and a water flooded reservoir in the lithologic region; sw is the water saturation of an undisturbed reservoir and a water flooded reservoir in the lithologic region; the lithology coefficient a related to lithology, the constant b related to lithology, the cementation index m, the saturation index n, and the rockAcquiring the formation water resistivity Rw of the undisturbed reservoir and the flooded reservoir of the zone from a preset corresponding relation table; the corresponding relation table comprises lithology coefficients, constants, cementation indexes, saturation indexes and formation water resistivity, which are corresponding to the undisturbed reservoir and the water flooded reservoir of each lithologic area and are related to lithology;
and establishing a logging interpretation model of the oil saturation of the undisturbed reservoir and the flooded reservoir in the lithologic region according to the water saturation model: so＝1-Sw(ii) a Wherein S isoThe oil saturation of the reservoir is the original reservoir in the rock region and the oil saturation of the reservoir flooded with water.
9. The device for explaining logging of a marine clastic rock horizontal well reservoir as claimed in claim 6, wherein the curve data correction unit is specifically configured to;
according to the formula: ACCorrection ofAnd (3) correcting the acoustic wave time difference curve of a horizontal well of the same reservoir section by α AC + β, wherein AC isCorrection ofThe corrected horizontal well sound wave time difference value, AC is the original sound wave time difference value of the horizontal well, α and β are correction parameters;
according to the formula: rt correction＝αRt+ β the formation deep resistivity curve of the horizontal well in the same reservoir section is corrected, wherein Rt correctionThe horizontal well stratum resistivity value after correction is obtained; rtWhich is the original formation resistivity value of the horizontal well, α and β are the calibration parameters.
10. The marine facies clastic rock horizontal well reservoir logging interpretation device of claim 9, wherein the preset horizontal well flooding qualitative interpretation rules in the flooding condition determination unit include natural gamma value ranges and eigenvalues, acoustic wave time difference value ranges and eigenvalues, and resistivity value ranges and eigenvalues corresponding to each flooding condition and each flooding condition;
the flooding condition determining unit is specifically configured to:
and inquiring the corrected water flooding conditions corresponding to the natural gamma curve, the acoustic wave time difference curve and the stratum deep resistivity curve of the horizontal well in the preset water flooding qualitative interpretation rule of the horizontal well, and identifying that the water flooding condition of the reservoir section of the horizontal well in the lithologic region is an oil layer or a reservoir layer of a poor oil layer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710216848.1A CN106951660B (en) | 2017-04-05 | 2017-04-05 | Sea facies clastic rock horizontal well reservoir logging interpretation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710216848.1A CN106951660B (en) | 2017-04-05 | 2017-04-05 | Sea facies clastic rock horizontal well reservoir logging interpretation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106951660A CN106951660A (en) | 2017-07-14 |
| CN106951660B true CN106951660B (en) | 2020-08-11 |
Family
ID=59475263
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710216848.1A Active CN106951660B (en) | 2017-04-05 | 2017-04-05 | Sea facies clastic rock horizontal well reservoir logging interpretation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106951660B (en) |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107703561A (en) * | 2017-10-23 | 2018-02-16 | 中石化石油工程技术服务有限公司 | The logging method of organic carbon content is calculated in a kind of shale gas horizontal well |
| CN110306978A (en) * | 2018-03-27 | 2019-10-08 | 中国石油化工股份有限公司 | A kind of heterogeneous reservoir reservoir physical parameter prediction technique |
| CN110320562B (en) * | 2018-03-30 | 2020-11-13 | 中国石油化工股份有限公司 | Method for correcting acoustic time difference in horizontal well of shale gas reservoir |
| CN110412660A (en) * | 2018-04-26 | 2019-11-05 | 中国石油大学(北京) | Reservoir Classification method and apparatus |
| CN109766615B (en) * | 2018-12-29 | 2023-05-30 | 中国石油天然气集团有限公司 | Reservoir water production rate calculation method based on apparent reservoir resistivity curve |
| CN110046414A (en) * | 2019-04-03 | 2019-07-23 | 长江大学 | The Reservoir Parameter Models construction method and equipment of advantage seepage flow section for identification |
| CN110399649B (en) * | 2019-07-03 | 2023-05-30 | 中国石油天然气集团有限公司 | Metamorphic rock quantitative identification method based on diagenetic indicating element |
| CN110656930B (en) * | 2019-10-18 | 2022-11-08 | 中国海洋石油集团有限公司 | Clastic rock stratum pressure calculation method based on lithology-physical property classification |
| CN110688781B (en) * | 2019-10-31 | 2022-10-28 | 西安石油大学 | Well logging interpretation method for low-permeability heterogeneous gas reservoir |
| CN110910267B (en) * | 2019-11-18 | 2023-11-28 | 中国石油天然气股份有限公司 | Method and device for determining sea-phase carbonate rock interlayer |
| CN111079283B (en) * | 2019-12-13 | 2023-06-09 | 四川新网银行股份有限公司 | Method for processing information saturation imbalance data |
| CN111008482B (en) * | 2019-12-18 | 2023-04-07 | 长江大学 | Metamorphic rock reservoir water saturation calculation method and device |
| CN112989551B (en) * | 2019-12-18 | 2022-11-04 | 中国石油天然气股份有限公司 | Metamorphic rock reservoir porosity calculation method and device and computer-readable storage medium |
| CN111025409B (en) * | 2019-12-23 | 2021-03-02 | 中国石油大学(北京) | Flooded layer evaluation method and device and storage medium |
| CN112228037A (en) * | 2020-09-04 | 2021-01-15 | 长江大学 | Porosity explanation method and device based on sound wave propagation theory |
| CN112268923B (en) * | 2020-10-22 | 2022-05-20 | 西安交通大学 | Method for acquiring formation thermal conductivity based on logging curve |
| CN112730808B (en) * | 2020-12-28 | 2023-04-07 | 长江大学 | Method for determining number of continental lake basin sea invasion action periods under non-dense sampling condition |
| CN113107464B (en) * | 2021-05-11 | 2024-05-07 | 中国石油天然气集团有限公司 | Horizontal well stepping type water flooded layer identification logging method |
| CN113356834B (en) * | 2021-07-14 | 2023-04-18 | 中国石油天然气股份有限公司 | Borehole acoustic time difference correction method, correction device and electronic equipment |
| CN116699697B (en) * | 2022-02-25 | 2024-05-07 | 中国石油化工股份有限公司 | Sandstone logging interpretation method and device based on glauconite mineral double-structure model |
| CN116084929B (en) * | 2023-04-10 | 2023-06-16 | 西北大学 | Oil-water interface determining method |
| CN117332301B (en) * | 2023-10-17 | 2024-07-12 | 大庆油田有限责任公司 | Flooding layer interpretation method for reservoir classification evaluation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8939217B2 (en) * | 2010-11-30 | 2015-01-27 | Tempress Technologies, Inc. | Hydraulic pulse valve with improved pulse control |
| CN105317431A (en) * | 2014-07-30 | 2016-02-10 | 中国石油化工股份有限公司 | Method for explaining and evaluating logging parameters of horizontal well |
| CN105631753A (en) * | 2015-12-14 | 2016-06-01 | 中国石油天然气集团公司 | Attitude of stratum-based horizontal well oil reservoir profile modeling method |
| CN105701319A (en) * | 2016-03-08 | 2016-06-22 | 中国石油大学(华东) | Sedimentary microfacies modeling method under horizontal well rule development well net |
-
2017
- 2017-04-05 CN CN201710216848.1A patent/CN106951660B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8939217B2 (en) * | 2010-11-30 | 2015-01-27 | Tempress Technologies, Inc. | Hydraulic pulse valve with improved pulse control |
| CN105317431A (en) * | 2014-07-30 | 2016-02-10 | 中国石油化工股份有限公司 | Method for explaining and evaluating logging parameters of horizontal well |
| CN105631753A (en) * | 2015-12-14 | 2016-06-01 | 中国石油天然气集团公司 | Attitude of stratum-based horizontal well oil reservoir profile modeling method |
| CN105701319A (en) * | 2016-03-08 | 2016-06-22 | 中国石油大学(华东) | Sedimentary microfacies modeling method under horizontal well rule development well net |
Non-Patent Citations (1)
| Title |
|---|
| "哈得逊油田东河砂岩储层水淹变化机理研究";刘强 等;《中国矿业大学学报》;20170404;第47卷(第6期);150-159 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106951660A (en) | 2017-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106951660B (en) | Sea facies clastic rock horizontal well reservoir logging interpretation method and device | |
| CN104047597B (en) | Oily mud shale stratum log standardized method | |
| Sen et al. | Petrophysical heterogeneity of the early Cretaceous Alamein dolomite reservoir from North Razzak oil field, Egypt integrating well logs, core measurements, and machine learning approach | |
| CN104698497B (en) | A kind of fine calibration of thin sand-mud interbed layer position | |
| CN109143399B (en) | Method for identifying carbonate rock sequence interface | |
| CN111856566B (en) | Method and device for predicting thin reservoir in sand body of lake-phase beach dam | |
| CN108629459B (en) | Method and device for detecting hydrocarbon-containing pore of reservoir | |
| Yao et al. | Status and prospects of exploration and exploitation key technologies of the deep petroleum resources in onshore China | |
| CN105842733A (en) | Shale reservoir earthquake identification method | |
| CN107831540A (en) | The direct new method for extracting of reservoir physical parameter | |
| Wimmers et al. | Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir | |
| CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
| CN115857047B (en) | Comprehensive prediction method for earthquake reservoir | |
| CN113376692B (en) | Method and device for optimizing fracturing modification scheme of tight sandstone gas horizontal well | |
| Taheri et al. | Quantifying fracture density within the Asmari reservoir: an integrated analysis of borehole images, cores, and mud loss data to assess fracture-induced effects on oil production in the Southwestern Iranian Region | |
| CN115387785A (en) | Sea-facies carbonate-cuttings limestone reservoir high-permeability strip identification method and device | |
| CN107831537B (en) | Method for describing high sand-to-land ratio nepheloid rock reservoir | |
| Li et al. | Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China | |
| Johansen et al. | Use of Advanced Wireline Logs To Reduce Uncertainties in a Complex Reservoir: A Case Study From the Ivar Aasen Oilfield in the Norwegian Central North Sea | |
| CN113495293B (en) | Reservoir fluid prediction method and device | |
| CN110927819B (en) | Crack development degree characterization method | |
| CN111708100B (en) | Deep thin layer oil and gas reservoir determination method and related device | |
| CN112147698B (en) | Crack development zone identification and feature determination method and system | |
| CN113093275B (en) | Method and device for improving drilling success rate of curved-flow river and curved-flow river delta oilfield | |
| Zhiwei | Research on Constrained Sedimentary Facies Reservoir by Well-seismic Inversion |
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 |