CN107831537B - Method for describing high sand-to-land ratio nepheloid rock reservoir - Google Patents
Method for describing high sand-to-land ratio nepheloid rock reservoir Download PDFInfo
- Publication number
- CN107831537B CN107831537B CN201710006504.8A CN201710006504A CN107831537B CN 107831537 B CN107831537 B CN 107831537B CN 201710006504 A CN201710006504 A CN 201710006504A CN 107831537 B CN107831537 B CN 107831537B
- Authority
- CN
- China
- Prior art keywords
- inversion
- well
- sand
- reservoir
- participating
- 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 24
- 238000000034 method Methods 0.000 title claims abstract description 20
- 239000004576 sand Substances 0.000 claims abstract description 45
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 4
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 4
- 230000000704 physical effect Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 230000008021 deposition Effects 0.000 abstract description 4
- 238000011161 development Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 2
- CDBYLPFSWZWCQE-UHFFFAOYSA-L Sodium Carbonate Chemical compound [Na+].[Na+].[O-]C([O-])=O CDBYLPFSWZWCQE-UHFFFAOYSA-L 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6163—Electromagnetic
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6167—Nuclear
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a high sand ratio nepheloid sedimentary rock reservoir description method, which comprises the following steps: step 1, fine synthesis record calibration is carried out, and the corresponding relation between lithology and seismic reflection characteristics is determined; step 2, analyzing a multi-well intersection map, and preferably selecting a well logging curve capable of effectively distinguishing sand mudstones; step 3, performing reservoir inversion under waveform indication by using the effective logging curve; step 4, starting from the wells, performing mutual comparison and careful analysis on the wells participating in inversion, and determining an inversion chromatographic block; and 5, comparing and verifying the real drilling sand bodies participating in the inversion well and the non-participating in the inversion well with the well-side inversion channel, and accurately tracking and describing the sand bodies by using inversion and drilling data. The method for describing the high sand ratio nepheloid sedimentary rock reservoir can accurately and quickly describe the sand body form matched with the actual deposition characteristics, and can accurately reflect the development and distribution characteristics of nepheloid sedimentary rock.
Description
Technical Field
The invention relates to the technical field of oil exploration, in particular to a method for describing a high sand-to-land ratio nephelometric rock reservoir.
Background
The method is formed for accurately and finely describing the reservoir.
In recent years, the exploration of the nepheloid accumulated rock of the oil field is greatly successful, and the nepheloid accumulated rock reservoir becomes an important storage increasing place in the near term. However, the nephelometric sedimentary rock sedimentary body developed by the continental lake basin is complex in development, the lithofacies changes rapidly, and the phenomenon of mutual superposition exists. Along with the deep exploration, the turbid integrated bodies which are large in scale and easy to identify in earthquake are basically identified; some overlapped sands have thin mud-rock interlayer between adjacent sands, and the high sand turbidity deposition rocks are not easy to be identified on the seismic section, so the distribution form and the deposition characteristics of the sands are difficult to be cleaned directly from the explanation of seismic data and the section of the conventional inversion result. The existing geophysical technical method (such as conventional seismic description, attribute extraction and the like) which is used for searching large-scale and is easy to identify in earthquake is not suitable for increasingly complicated geological trap types of the nephritic rocks, and the overall exploration deployment and reserve upgrading evaluation of the nephritic rocks are influenced. Therefore, a description method for the turbidimetric rock reservoir in the high sand area is invented, and the problems are solved.
Disclosure of Invention
The invention aims to provide a high sand-land ratio nepheloid rock reservoir description method for accurately implementing sand body boundaries and clearing sand body distribution under the condition of high sand-land ratio.
The object of the invention can be achieved by the following technical measures: a high sand ratio nepheloid rock reservoir description method comprises the following steps: step 1, fine synthesis record calibration is carried out, and the corresponding relation between lithology and seismic reflection characteristics is determined; step 2, analyzing a multi-well intersection map, and preferably selecting a well logging curve capable of effectively distinguishing sand mudstones; step 3, performing reservoir inversion under waveform indication by using the effective logging curve; step 4, starting from the wells, performing mutual comparison and careful analysis on the wells participating in inversion, and determining an inversion chromatographic block; and 5, comparing and verifying the real drilling sand bodies participating in the inversion well and the non-participating in the inversion well with the well-side inversion channel, and accurately tracking and describing the sand bodies by using inversion and drilling data.
The object of the invention can also be achieved by the following technical measures:
in step 1, calibrating through fine synthetic records in a seismic workstation, determining the corresponding relation between lithology and seismic reflection characteristics, and analyzing the describable probability of the high sand area turbidimetric rock reservoir on a conventional seismic section.
In step 2, the well logging curves are normalized, and well logging curves which can reflect stratum lithologic changes and are sensitive are found out through multi-well intersection map analysis, wherein the well logging curves comprise natural potential and natural gamma curves.
In step 3, reservoir inversion under waveform indication is carried out by using the optimized logging curve and by means of Schmei inversion software, the principle is a process of optimizing reflection coefficient combination under the guidance of seismic waveform characteristics, a waveform phase control inversion idea is adopted, and the seismic waveform characteristics are used for replacing a variation function to analyze reservoir space structure change, so that the fine reservoir prediction capability is improved.
In step 4, starting from the wells, performing mutual comparison and careful analysis on the wells participating in inversion, and preliminarily determining a chromatographic block.
In step 5, comparing and verifying the real drilling sand bodies participating in the inversion well and the non-participating in the inversion well with the well-side inversion channel, and extrapolating the sand body boundary by utilizing waveform similarity and inversion color change rules.
The method for describing the high sand-to-land ratio nepheloid sedimentary rock reservoir can accurately describe the nepheloid sedimentary rock reservoir with smaller mudstone interval between the vertically adjacent sand bodies, can accurately and quickly describe the sand body form matched with the actual deposition characteristic, and can accurately reflect the development and distribution characteristics of nepheloid sedimentary rock.
Drawings
FIG. 1 is a flow chart depicting one embodiment of a method of the present invention for characterizing a high sand turbidimetric reservoir;
FIG. 2 is a schematic illustration of a well lithology and seismic reflection signature match relationship in an embodiment of the present invention;
FIG. 3 is a schematic representation of a preferred lithology sensitive log in an embodiment of the present invention;
FIG. 4 is a schematic representation of another preferred litho-sensitive log in an embodiment of the present invention;
FIG. 5 is a schematic diagram of the inversion principle of waveform indication according to an embodiment of the present invention;
FIG. 6 is a waveform indicating inversion flow in an embodiment of the present invention;
FIG. 7 is a schematic diagram illustrating the selection principle of the inversion color patch block according to an embodiment of the present invention;
FIG. 8 is a schematic diagram of inversion results and well matching in an embodiment of the invention.
Detailed Description
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
As shown in fig. 1, fig. 1 is a flow chart depicting a method of the present invention for high sand turbidimetric reservoirs.
Step 101, fine synthesis record calibration is carried out, and the corresponding relation between lithology and seismic reflection characteristics is determined. And calibrating through fine synthetic records in a seismic workstation, determining the corresponding relation between lithology and seismic reflection characteristics, and analyzing the describable probability of the high sand area turbidimetric rock reservoir on a conventional seismic section. In one embodiment, the frequency band of the seismic data of the three-dimensional earthquake in the target interval is narrow, the main frequency is low, and the signal-to-noise ratio is low. Meanwhile, the sand-ground ratio of the turbid accumulated rock of the Dongying 4 sand groups in the region is nearly 50%, the sand body cannot be traced on a seismic section through fine synthetic record calibration, the strong axis and the weak axis can represent the reflection characteristics of the turbid accumulated sand body, and the boundary of the sand body cannot be accurately traced, as shown in figure 2.
And 102, analyzing a multi-well intersection map, and preferably selecting a well logging curve capable of effectively distinguishing sand mudstones. And carrying out normalization processing on the logging curve, and finding out the logging curve which can reflect the stratum lithology change and is sensitive, such as a natural potential and a natural gamma curve, through analysis of a multi-well intersection map. In one embodiment, the logging curve is first subjected to unified correction and processing, including mudstone baseline correction and normalization processing; through analysis of a multi-well intersection diagram, a natural gamma curve is found to be sensitive to lithology changes, and sand shale can be well distinguished, as shown in fig. 3 and 4.
And 103, performing reservoir inversion under the waveform indication by using the effective logging curve. And performing reservoir inversion under the waveform indication by using the optimized logging curve and by means of Schmei inversion software. The principle is a process of optimizing the combination of reflection coefficients under the guidance of seismic waveform characteristics, as shown in fig. 5. The method breaks through the limitation of the traditional wave impedance inversion based on a spatial domain interpolation algorithm, is not influenced by well point distribution, adopts the wave phase control inversion idea, utilizes seismic waveform characteristics to replace a variation function to analyze reservoir space structure change, and improves the fine prediction capability of a reservoir. In one embodiment, the reservoir inversion under the waveform indication is performed by using the natural gamma curve sensitive to lithology, which is preferably selected in step 102, through a Simie inversion software, and the core is a process of optimizing the reflection coefficient combination under the guidance of seismic waveform characteristics. The specific flow is shown in fig. 6.
And step 104, starting from the wells, performing mutual comparison and careful analysis on the wells participating in inversion, and determining an inversion chromatographic block. Effective sandstone with good physical properties is reflected in black and white. And finally, extrapolating the sand body boundary by utilizing the waveform similarity and the inversion color change rule. In one embodiment, after the preliminary generation of the inverted data volume in step 103, sand description rules are first formulated. And starting from the wells, performing mutual comparison and careful analysis on the wells participating in inversion, and preliminarily determining a chromatographic block. Statistics shows that the shale content is in positive correlation with the GR value, and sandstone with the shale content of more than 40% has poor physical property and is a dry layer; in order to identify effective sand bodies, invalid sandstones (light ash and deep ash) are determined by setting the argillaceous content to be more than 40% and the GR to be more than 100 API; the effective sandstone having good physical properties is reflected in black and white, as shown in fig. 7.
And 105, comparing and verifying the real drilling sand bodies participating in the inversion well and the non-participating in the inversion well with the well-side inversion channel, and accurately tracking and describing the sand bodies by fully utilizing data such as inversion, drilling and the like. In one embodiment, after step 104 is completed, the inverted profile of the through-well is cut and the through-well solid borehole sand is compared to the off-well inversion trace. By contrast, the drilling sand body is well matched with the inversion section, the main oil-containing sand body is high in coincidence rate, the oil-water relation is not contradictory, and the total coincidence degree is more than 80%, as shown in fig. 8. The method has the advantages that the longitudinal and transverse resolution of the inversion section is obviously higher than that of the conventional seismic section, the sand body sharp point vanishing point is clear, the problem that the seismic section cannot be subjected to sand body description is solved, and the inversion body is utilized to accurately track and describe the sand body in combination with well drilling and other data.
Claims (1)
1. The method for describing the high sand ratio nepheloid rock reservoir is characterized by comprising the following steps of:
step 101, fine synthesis record calibration is carried out, and the corresponding relation between lithology and seismic reflection characteristics is determined; calibrating through fine synthetic records in a seismic workstation, determining the corresponding relation between lithology and seismic reflection characteristics, and analyzing the describable probability of the high sand area nephelometric rock reservoir on a conventional seismic section;
102, analyzing a multi-well intersection map, and preferably selecting a well logging curve capable of effectively distinguishing sand mudstones; carrying out normalization processing on the logging curve, and finding out a natural gamma curve which can reflect that the lithological change of the stratum is sensitive through analysis of a multi-well intersection map;
103, performing reservoir inversion under waveform indication by using the effective logging curve; reservoir inversion under the waveform indication is carried out by using the optimized logging curve and by means of Schmei inversion software, and for the process of optimizing the reflection coefficient combination under the guidance of seismic waveform characteristics, the waveform phase control inversion idea is adopted, and the seismic waveform characteristics are used for replacing variation functions to analyze reservoir space structure changes, so that the fine prediction capability of the reservoir is improved;
104, after the inversion data volume is preliminarily generated in the step 103, firstly making a sand body description principle; starting from a well, performing mutual comparison and analysis on the wells participating in inversion, and preliminarily determining a chromatographic block, wherein the sandstone with the shale content of more than 40% has poor physical property and is a dry layer; the shale content exceeds 40 percent, and the GR is more than 100API, so that the sandstone is invalid; effective sandstone with good physical properties is reflected by black and white, and finally, the sand body boundary is extrapolated by utilizing waveform similarity and an inversion color change rule;
and 105, cutting the inversion section of the well after the completion of the step 104, comparing and verifying the real drilling sand bodies participating in the inversion well and not participating in the inversion well with the well-side inversion channel, and accurately tracking and describing the sand bodies by using inversion and drilling data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710006504.8A CN107831537B (en) | 2017-01-05 | 2017-01-05 | Method for describing high sand-to-land ratio nepheloid rock reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710006504.8A CN107831537B (en) | 2017-01-05 | 2017-01-05 | Method for describing high sand-to-land ratio nepheloid rock reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107831537A CN107831537A (en) | 2018-03-23 |
| CN107831537B true CN107831537B (en) | 2020-07-14 |
Family
ID=61643059
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710006504.8A Active CN107831537B (en) | 2017-01-05 | 2017-01-05 | Method for describing high sand-to-land ratio nepheloid rock reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107831537B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112882092B (en) * | 2021-01-11 | 2022-06-21 | 中国海洋大学 | Well-seismic joint inversion method, system, medium and application based on multi-attribute CNN |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102841378A (en) * | 2011-06-23 | 2012-12-26 | 中国石油天然气股份有限公司 | Method for predicting reservoir stratum by seismic inversion data |
| CN104698496A (en) * | 2013-12-05 | 2015-06-10 | 中国石油化工股份有限公司 | Small sand body boundary identification method and small sand body space quantitative description method |
| CN104898161A (en) * | 2014-03-05 | 2015-09-09 | 中国石油集团东方地球物理勘探有限责任公司 | Effective sandstone predicting method based on logging response simulator |
| CN105652323A (en) * | 2015-08-11 | 2016-06-08 | 中国石油化工股份有限公司 | Reservoir stratum prediction method |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2859383B1 (en) * | 2012-06-11 | 2018-10-17 | Saudi Arabian Oil Company | System for displaying of petrophysical property height profile for horizontal or slant wellbores |
| CN102759748B (en) * | 2012-07-04 | 2016-05-25 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | The cross analysis Seismic Reservoir Prediction method that based target is analyzed |
| CN103176211B (en) * | 2013-03-07 | 2016-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | Based on gas-bearing reservoir prediction method and the device of many sensibility elasticities parameter |
| CN103643949B (en) * | 2013-12-20 | 2016-06-01 | 中国石油天然气集团公司 | A kind of reservoir contains quantitative forecast method and the device of oil gas |
| CN104502966B (en) * | 2014-12-23 | 2017-06-13 | 中国石油天然气集团公司 | A kind of Forecasting Methodology and system of thin reservoir |
| WO2016110298A1 (en) * | 2015-01-06 | 2016-07-14 | Entreprise Nationale De Geophysique | Method for estimating the volume of clay (v clay) on the basis of the formation density, of the porosity and of the acoustic velocities in clay-sandstone reservoirs |
| CN105807320A (en) * | 2016-03-17 | 2016-07-27 | 成都创源油气技术开发有限公司 | Shale gas reservoir seismic prediction method |
-
2017
- 2017-01-05 CN CN201710006504.8A patent/CN107831537B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102841378A (en) * | 2011-06-23 | 2012-12-26 | 中国石油天然气股份有限公司 | Method for predicting reservoir stratum by seismic inversion data |
| CN104698496A (en) * | 2013-12-05 | 2015-06-10 | 中国石油化工股份有限公司 | Small sand body boundary identification method and small sand body space quantitative description method |
| CN104898161A (en) * | 2014-03-05 | 2015-09-09 | 中国石油集团东方地球物理勘探有限责任公司 | Effective sandstone predicting method based on logging response simulator |
| CN105652323A (en) * | 2015-08-11 | 2016-06-08 | 中国石油化工股份有限公司 | Reservoir stratum prediction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107831537A (en) | 2018-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106951660B (en) | Sea facies clastic rock horizontal well reservoir logging interpretation method and device | |
| CN109061765B (en) | Trap evaluation method for heterogeneous thin sandstone interbed reservoir | |
| CN105759312A (en) | Reservoir characteristic method well-seismic calibration method | |
| CN106842307A (en) | It is a kind of based on just drilling the subtle hydrocarbon reservoir method that the lower waveform separation of constraint is retrieved again | |
| CN109541685B (en) | River channel sand body identification method | |
| CN105093306A (en) | Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration | |
| CN109143399B (en) | Method for identifying carbonate rock sequence interface | |
| CN105527653A (en) | Virtual well construction method based on geological information | |
| Euzen et al. | Well log cluster analysis: an innovative tool for unconventional exploration | |
| CN105240006A (en) | Oil and water layer recognition method suitable for volcanic reservoir | |
| CN108508489B (en) | Seismic inversion method based on waveform micro-variation matching | |
| CN112711068B (en) | Method and device for predicting effective reservoir of oil gas in sandstone | |
| CN117607995A (en) | Quasi-impedance seismic lithology inversion method utilizing fusion of natural gamma and wave impedance | |
| CN107831537B (en) | Method for describing high sand-to-land ratio nepheloid rock reservoir | |
| CN115857047B (en) | Comprehensive prediction method for earthquake reservoir | |
| CN111257933A (en) | Novel method for predicting oil and gas reservoir based on low-frequency shadow phenomenon | |
| CN107329168B (en) | Sedimentary body lithology identification method and system | |
| CN113514884A (en) | Compact sandstone reservoir prediction method | |
| CN110389381B (en) | Sand reservoir prediction method and device based on seismic attributes | |
| CN110297264B (en) | Low-permeability gas reservoir thin reservoir dessert earthquake prediction method | |
| CN113109875A (en) | Inversion method of carbonate rock reservoir under full waveform velocity field constraint | |
| Zhao et al. | An automatical infill shot method for uniform imaging of target layer | |
| CN111239823A (en) | Method for identifying distribution of invaded rocks | |
| CN113495293B (en) | Reservoir fluid prediction method and device | |
| CN115932967B (en) | Crack identification method based on ant body algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |