CN105589097B - The phased dedusting method of clastic reservoir rock seismic description - Google Patents

The phased dedusting method of clastic reservoir rock seismic description Download PDF

Info

Publication number
CN105589097B
CN105589097B CN201410571913.9A CN201410571913A CN105589097B CN 105589097 B CN105589097 B CN 105589097B CN 201410571913 A CN201410571913 A CN 201410571913A CN 105589097 B CN105589097 B CN 105589097B
Authority
CN
China
Prior art keywords
grey matter
seismic
petrofacies
reservoir
rock
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
Application number
CN201410571913.9A
Other languages
Chinese (zh)
Other versions
CN105589097A (en
Inventor
肖淑明
王建伟
刘志萍
王建华
蔡宏兴
赵海燕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Shengli Oilfield Co Xianhe Oil Production Plant
Original Assignee
China Petroleum and Chemical Corp
Sinopec Shengli Oilfield Co Xianhe Oil Production Plant
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Shengli Oilfield Co Xianhe Oil Production Plant filed Critical China Petroleum and Chemical Corp
Priority to CN201410571913.9A priority Critical patent/CN105589097B/en
Publication of CN105589097A publication Critical patent/CN105589097A/en
Application granted granted Critical
Publication of CN105589097B publication Critical patent/CN105589097B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The present invention provides a kind of phased dedusting method of clastic reservoir rock seismic description, and the phased dedusting method of the clastic reservoir rock seismic description includes:Step 1, grey matter lithofacies distribution statistics is completed using work area drilling data;Step 2, the seismic response forward modeling analysis of work area individual well target zone grey matter rock petrofacies is carried out, obtains grey matter petrofacies and seismic response relational expression;Step 3, asking for grey matter petrofacies influences the seismic properties correction body of reservoir prediction;Step 4, reservoir Application of Logging-constrained Inversion and relative earthquake attributes extraction are carried out;And step 5, data is bored in fact with reference to drilling well, and fine description is carried out using the Application of Logging-constrained Inversion data volume after correction or attribute volume.The phased dedusting method of the clastic reservoir rock seismic description substantially solves the grey matter influence that reservoir prediction is carried out using 3-D seismics, realizes fine description and the evaluation of reservoir, positive effect has been seen in exploration and development application.

Description

The phased dedusting method of clastic reservoir rock seismic description
Technical field
The present invention relates to oil field development technical field, especially relates to a kind of the phased of clastic reservoir rock seismic description and goes Grey method.
Background technology
Realized using three dimensional seismic data highly important in reservoir prediction and description always oil-gas exploration and development field One evaluation content.Seismic reflection character influence of the carbonate rock on clastic reservoir rock is very big, and removing grey matter influences to realize reservoir Accurate Prediction and description are always the field important research content.At present, the technology mainly utilizes individual well dedusting experiment understanding ash Confrontation seismic reservoir response relation, then explain that personnel intuitively distinguish the seismic response features of reservoir, this method is main between well It is inadequate by explanation personnel's subjective judgement, fine degree.In addition, mainly removing grey matter using geophysical method in recent years influences Method for predicting reservoir also achieve greater advance, including the methods of wavelet decomposition, more attributes joint and prestack inversion, but to Drilling data depends on seism processing using limited, is related to that highly professional, method is complicated is difficult to popularization and application. Compared with high matured exploration area, drilling well density is higher, and the cementing situation of carbonate rock is totally clear and definite, and numerous geologists need a kind of fast Prompt method effective, economical and practical, technical costs is low, to serve daily production work.For this, we have invented a kind of new Clastic reservoir rock seismic description phased dedusting method, solve above technical problem.
The content of the invention
Combining closely for drilling well and seismic data is strengthened it is an object of the invention to provide a kind of, meets seismology and statistics General principle is learned, can be widely applied to similar to the clastic reservoir rock seismic description in the oil-gas exploration and development work under the conditions of geology Phased dedusting method.
The purpose of the present invention can be achieved by the following technical measures:The phased dedusting side of clastic reservoir rock seismic description Method, the phased dedusting method of the clastic reservoir rock seismic description include:Step 1, grey matter petrofacies are completed using work area drilling data Distribution statisticses;Step 2, carry out work area individual well target zone grey matter rock petrofacies seismic response forward modeling analysis, obtain grey matter petrofacies with Seismic response relational expression;Step 3, asking for grey matter petrofacies influences the seismic properties correction body of reservoir prediction;Step 4, reservoir is carried out Application of Logging-constrained Inversion and relative earthquake attributes extraction;And step 5, data is bored in fact with reference to drilling well, utilizes the survey after correction Well constrained inversion data volume or attribute volume carry out fine description.
The purpose of the present invention can be also achieved by the following technical measures:
In step 1, target zone grey matter rock is carried out with reference to interval transit time, natural potential, a variety of well-log informations of natural gamma The identification of phase, statistical lithofacies thickness, ultimately form grey matter petrofacies flat distribution map.
In step 2, the seismic response forward modeling analysis of work area individual well target zone grey matter rock petrofacies, the displacement of grey matter petrofacies are completed Seismic properties variable y and grey matter petrofacies thickness x is closely related after normal country rock, more well fitting attribute variables and grey matter petrofacies thickness Variable draws relational expression:Y=289.73e0.0881x, so as to clear and definite grey matter petrofacies and the response relation formula of Seismic reflection character.
In step 3, according to grey matter petrofacies and seismic response relational expression y=289.73e0.0881x, controlled in grey matter petrofacies Grey matter petrofacies are asked in lower computing influences the seismic properties correction body of reservoir prediction, and association attributes school is carried out to three dimensional seismic data Just.
In step 4, reservoir Application of Logging-constrained Inversion and relative earthquake are carried out using the 3-d seismic data set after correction Attributes extraction, 3-d seismic data set can more Accurate Prediction reservoir after correction.
In clastic rock oil basin, carbonate rock has a significant impact to the seismic description of reservoir and prediction, so as to cause The accurate sexual deviation of layer description, probing is caused to be defeated.Present invention aim to address the dedusting problem of such Earthquakes data, So as to improve the earthquake prediction of reservoir and description accuracy.Former technology is confined to individual well dedusting experiment and processing, inter-wells predicting Still sentence knowledge by explanation personnel's subjectivity, quantitative assessment can not be carried out;And simple geophysical method requires that professional technical is strong, bore Well and seismic data combination degree are inadequate, it is impossible to meet the evaluation requirement in higher probing area.The present invention is utilizing more mouthfuls of individual wells just Drill and grey matter content and seismic properties response relation formula are established in experimental basis;NULL area interval of interest grey matter content, it is determined that Grey matter petrofacies planar distribution feature;Using grey matter petrofacies as background control condition, according to grey matter petrofacies and seismic properties response relation Formula, whole removing influences except the grey matter of research area's 3-D seismics association attributes, and then completes layer description and prediction work.In example In, the present invention is obvious to solve the grey matter influence that reservoir prediction is carried out using 3-D seismics, realizes the fine of reservoir and retouches State and evaluate, positive effect has been seen in exploration and development application.
Brief description of the drawings
Fig. 1 phased dedusting methods of clastic reservoir rock seismic description in the embodiment of Shengli Oil Field work area for the present invention Flow chart;
Fig. 2 is in the work area embodiment central gray depth of stratum and seismic amplitude response relation figure using the present invention;
Fig. 3 is in the work area embodiment central gray petrofacies thickness distribution figure using the present invention;
Fig. 4 is to influence seismic reservoir response corrections body figure in the work area embodiment central gray petrofacies using the present invention;
Fig. 5 be using the present invention in the work area embodiment dedusting of inverting section influence before and after comparison diagram;
Fig. 6 be using the present invention in the work area embodiment amplitude attribute dedusting influence before and after comparison diagram.
Embodiment
For enable the present invention above and other objects, features and advantages become apparent, it is cited below particularly go out preferable implementation Example, and coordinate institute's accompanying drawings, it is described in detail below.
Fig. 1 is the flow chart of the present invention phased dedusting method of clastic reservoir rock seismic description in the work area embodiment. In step 101, grey matter lithofacies distribution statistics is completed using all drilling datas in work area.With reference to interval transit time, natural potential, nature A variety of well-log informations such as gamma carry out the identification of target zone grey matter petrofacies, statistical lithofacies thickness, ultimately form grey matter petrofacies and put down Mapping.Flow enters step 102.
Fig. 2 is in the work area embodiment central gray petrofacies thickness distribution figure using the present invention.Wherein, the ash that well logging recognition goes out Matter petrofacies thickness parameter is most sensitive to seismic response.
In step 102, the seismic response forward modeling analysis of work area individual well target zone grey matter rock petrofacies, the displacement of grey matter petrofacies are completed Seismic properties variable (y) and grey matter petrofacies thickness (x) are closely related after normal country rock, more well fitting attribute variables and grey matter petrofacies Gauge variation draws relational expression:Y=289.73e0.0881x, so as to the response relation of clear and definite grey matter petrofacies and Seismic reflection character Formula.Flow enters step 103.
Fig. 3 is in the work area embodiment central gray depth of stratum and seismic amplitude response relation figure using the present invention.Wherein, Amplitude attribute influence relation is fairly obvious after individual well removes grey matter petrofacies, defers to the relational expression fitted.
In step 103, asking for grey matter petrofacies influences the seismic properties correction body of reservoir prediction.According to grey matter petrofacies and earthquake Response relation formula y=289.73e0.0881x, grey matter petrofacies are asked in computing under the control of grey matter petrofacies influences the earthquake of reservoir prediction Attribute corrects body, and association attributes correction is carried out to three dimensional seismic data.Flow enters step 104.
Fig. 4 is to influence seismic reservoir response corrections body figure in the work area embodiment central gray petrofacies using the present invention.Wherein, Bottom figure layer is the seismic response features under the influence of no grey matter petrofacies, and top figure layer is the actual seismic under the influence of grey matter petrofacies Response characteristic, this invention are exactly phased lower removal difference therebetween.
In step 104, reservoir Application of Logging-constrained Inversion and relative earthquake are carried out using the 3-d seismic data set after correction Attributes extraction.Reservoir Application of Logging-constrained Inversion and seismic attribute abstraction presentation notable difference, school are carried out compared to using initial data body Data volume can more Accurate Prediction reservoir after just.Flow enters step 105.
Fig. 5 be using the present invention in the work area embodiment dedusting of inverting section influence before and after comparison diagram.Wherein, upper figure is The Application of Logging-constrained Inversion section of grey matter petrofacies influence is not removed, and the Application of Logging-constrained Inversion that figure below influences to remove grey matter petrofacies cuts open Face, it is coincide with drilling well, positive effect.
Fig. 6 be using the present invention in the work area embodiment amplitude attribute dedusting influence before and after comparison diagram.Wherein, upper figure is The amplitude attribute plan of grey matter petrofacies influence is not removed, figure below is to remove the amplitude attribute plan that grey matter petrofacies influence, with Drilling well coincide, positive effect.
In step 105, data is bored in fact with reference to drilling well, is entered using the Application of Logging-constrained Inversion data volume after correction or attribute volume Row fine description.So far, influence of the grey matter petrofacies to layer description and prediction in three dimensional seismic data is eliminated, reservoir prediction Accuracy is obviously improved, and improves exploration and development benefit.Flow terminates.

Claims (1)

1. the phased dedusting method of clastic reservoir rock seismic description, it is characterised in that the clastic reservoir rock seismic description it is phased Dedusting method includes:
Step 1, grey matter rock lithofacies distribution statistics is completed using work area drilling data;
Step 2, the seismic response forward modeling analysis of work area individual well target zone grey matter rock petrofacies is carried out, obtains grey matter rock petrofacies and earthquake Response relation formula;
Step 3, asking for grey matter rock petrofacies influences the seismic properties correction body of reservoir prediction;
Step 4, reservoir Application of Logging-constrained Inversion and relative earthquake attributes extraction are carried out;And
Step 5, data is bored in fact with reference to drilling well, essence is carried out using the Application of Logging-constrained Inversion data volume after correction or earthquake attribute volume Thin description;
In step 1, target zone grey matter rock petrofacies are carried out with reference to interval transit time, natural potential, a variety of well-log informations of natural gamma Identification, statistical lithofacies thickness, ultimately form grey matter rock petrofacies flat distribution map;
In step 2, the seismic response forward modeling analysis of work area individual well target zone grey matter rock petrofacies is completed, grey matter rock petrofacies are replaced just Seismic properties variable y and grey matter rock petrofacies thickness x is closely related after normal country rock, more well fitting seismic properties variables and grey matter rock rock Phase thickness variable draws relational expression:y=289.73e0.0881x, so as to clear and definite grey matter rock petrofacies and seismic response relational expression;
In step 3, according to grey matter rock petrofacies and seismic response relational expression y=289.73e0.0881x, under the control of grey matter rock petrofacies Grey matter rock petrofacies are asked in computing influences the seismic properties correction body of reservoir prediction, and association attributes school is carried out to three dimensional seismic data Just;
In step 4, reservoir Application of Logging-constrained Inversion and relative earthquake attribute are carried out using the 3-d seismic data set after correction Extraction, 3-d seismic data set can more Accurate Prediction reservoir after correction.
CN201410571913.9A 2014-10-23 2014-10-23 The phased dedusting method of clastic reservoir rock seismic description Active CN105589097B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410571913.9A CN105589097B (en) 2014-10-23 2014-10-23 The phased dedusting method of clastic reservoir rock seismic description

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410571913.9A CN105589097B (en) 2014-10-23 2014-10-23 The phased dedusting method of clastic reservoir rock seismic description

Publications (2)

Publication Number Publication Date
CN105589097A CN105589097A (en) 2016-05-18
CN105589097B true CN105589097B (en) 2018-02-13

Family

ID=55928821

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410571913.9A Active CN105589097B (en) 2014-10-23 2014-10-23 The phased dedusting method of clastic reservoir rock seismic description

Country Status (1)

Country Link
CN (1) CN105589097B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107132574B (en) * 2017-04-11 2020-01-07 广州海洋地质调查局 Prediction method for lithologic combination of marine stratum clastic rocks
CN107103165B (en) * 2017-06-02 2020-12-01 西南石油大学 Method for calculating density of steep dip shear crack line of oil field reservoir
CN109061752B (en) * 2018-06-26 2020-01-17 西南石油大学 Resistivity curve correction method for ash-containing texture layer

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5661697A (en) * 1995-12-18 1997-08-26 Atlantic Richfield Company Method and apparatus for detection of sand formations in amplitude-versus-offset seismic surveys
CN101980053B (en) * 2010-10-18 2013-12-18 中国石油化工股份有限公司 Complicated reef flat reservoir predicting method
CN103513278B (en) * 2012-06-19 2016-02-10 中国石油化工股份有限公司 Seismic wave groups thickness is utilized to carry out the method for reservoir prediction

Also Published As

Publication number Publication date
CN105589097A (en) 2016-05-18

Similar Documents

Publication Publication Date Title
CN102636812B (en) Method for obtaining storage space volume of carbonate reservoir
CN102466815B (en) Triassic clastic rock petroleum reservoir identification method
CN104898161B (en) Effective sandstone predicting method based on logging response simulator
US10386516B2 (en) Time-lapse 4D scattering for imaging hydraulically induced fractures
Refunjol et al. Inversion and attribute-assisted hydraulically induced microseismic fracture characterization in the North Texas Barnett Shale
CN109164487A (en) A kind of method based on model foundation average velocity field and finely at the method for figure
CN104698492A (en) Method for calculating abnormal stratum pressure
Gupta et al. Rock typing in Wolfcamp formation
Price et al. Influence of stratigraphy on barriers to fracture growth and completion optimization in the Meramec Stack Play, Anadarko Basin, Oklahoma
CN105589097B (en) The phased dedusting method of clastic reservoir rock seismic description
Schneider et al. Interpretation of fractured zones using seismic attributes—Case study from Teapot Dome, Wyoming, USA
Willis Upscaling anisotropic geomechanical properties using Backus averaging and petrophysical clusters in the Vaca Muerta Formation
CN106033127A (en) Ground stress orientation seismic prediction method based on shear wave speed variation rate
Marfurt Seismic attributes and the road ahead
Oyeyemi et al. Seismic attributes analysis for reservoir characterization; offshore Niger Delta
Marroqun et al. Seismic attribute-based characterization of coalbed methane reservoirs: An example from the Fruitland Formation, San Juan basin, New Mexico
CN106597547A (en) Method for accurately describing earthquake in thin reservoir
Abe et al. Reservoir characterization and formation evaluation of some part of Niger Delta using 3D seismic and well log data
Adeoti et al. Structural Style and Reservoir Distribution in Deep-Water Niger Delta: A Case Study of" Nanny Field"
Lyu et al. Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China
Fox et al. Induced seismicity study in the Kiskatinaw seismic monitoring and mitigation area, British Columbia
CN104698493A (en) Method for calculating abnormal stratum pressure
CN105093330B (en) Method for searching side-well fracture-cavity reservoir body through multi-seismic-trace accumulated amplitude difference spectrum solution
CN107783192B (en) Carbonate rock ancient stream channel charges property earthquake prediction method and system
CN114063154A (en) Slip fault displacement amount calculation method, slip fault displacement amount calculation device, electronic apparatus, and medium

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant