阅读说明：本技术 一种沉积微相识别方法 (A kind of sedimentary facies identification method ) 是由 徐立恒 梁宇 孙辉 黄勇 朱权 张秀丽 何秋丽 杨会东 于 2018-05-24 设计创作，主要内容包括：本发明公开了一种沉积微相识别方法,属于油田储层预测领域。该方法包括以下步骤：对目标工区进行地震采集,获取叠前地震数据,并在目标工区内钻井和测井,获取钻井数据和测井数据。利用沉积微相数据、叠前地震数据,以及沉积微相的变差函数曲线,在微相-纵波阻抗关系图版和微相-横波阻抗关系图版的约束下进行叠前地质统计学反演,建立沉积微相三维模型。根据沉积微相三维模型提取目标层位的微相平面图和微相剖面图,完成对目标工区沉积微相的识别。本发明提供的沉积微相识别方法降低了不同的微相在单一波阻抗上存在多解性的可能,提高了储层沉积微相识别精度,对沉积微相结构的识别误差较小。(The invention discloses a kind of sedimentary facies identification methods, belong to oilfield reservoir prediction field.Method includes the following steps: carrying out earthquake-capturing to target work area, earthquake data before superposition, and drilling well and well logging in target work area are obtained, well data and log data are obtained.Using the variogram curve of sedimentary micro data, earthquake data before superposition and sedimentary micro, prestack geostatistical inversion is carried out under microfacies-p-wave impedance relationship plate and the constraint of microfacies-S-wave impedance relationship plate, establishes sedimentary micro threedimensional model.The microfacies plan view and microfacies sectional view of destination layer position are extracted according to sedimentary micro threedimensional model, complete the identification to target work area sedimentary micro.Sedimentary facies identification method provided by the invention reduces different microfacies, and there are the possibility of multi-solution in single wave impedance, improve Sedimentary Micro Facies accuracy of identification, smaller to the identification error of sedimentary micro structure.)

1. a kind of sedimentary facies identification method, which is characterized in that the described method includes: carrying out earthquake-capturing to target work area, obtain Earthquake data before superposition, and drilling well and well logging in the target work area are taken, well data and log data are obtained；

Using the well data and the log data, sedimentary micro data are obtained；

P-wave impedance data and S-wave impedance data are obtained using the log data, and count the sedimentary micro data, institute P-wave impedance data and the S-wave impedance data are stated, microfacies-p-wave impedance relationship plate and microfacies-S-wave impedance relationship are obtained Plate；

Using the variogram curve of the sedimentary micro data, the earthquake data before superposition and the sedimentary micro, in institute State progress prestack geostatistics under microfacies-p-wave impedance relationship plate and the constraint of the microfacies-S-wave impedance relationship plate Sedimentary micro threedimensional model is established in inverting；

The microfacies plan view and microfacies sectional view of destination layer position are extracted according to the sedimentary micro threedimensional model, are completed to the mesh Mark the identification of work area sedimentary micro.

2. the method according to claim 1, wherein the well data includes: core data, the well logging number According to including: gamma, microelectrode and microinverse；

It is described to utilize the well data and the log data, obtain sedimentary micro data, comprising:

Using the gamma, the microelectrode and the microinverse, gamma curve, Micro-electrode Log Curves and microinverse curve are obtained；

According to the gamma curve, the Micro-electrode Log Curves and the microinverse curve, with reference to the core data, explanation is obtained Sedimentary micro data.

3. the method according to claim 1, wherein the log data further include: velocity of longitudinal wave data, shear wave Speed data harmony wave density data；

It is described to obtain p-wave impedance data and S-wave impedance data using the log data, comprising: to utilize the longitudinal wave sound wave Speed data and the sound wave density data obtain p-wave impedance data；

Using the shear wave velocity data and the sound wave density data, S-wave impedance data are obtained.

4. according to the method described in claim 2, it is characterized in that, described utilize the sedimentary micro data, the prestack The variogram curve for shaking data and the sedimentary micro, in the microfacies-p-wave impedance relationship plate and the microfacies- Prestack geostatistical inversion is carried out under the constraint of S-wave impedance relationship plate, establishes sedimentary micro threedimensional model, comprising:

The deposition for going out the target work area according to the gamma curve, the Micro-electrode Log Curves and the microinverse curve interpretation is micro- Phase, the type of the sedimentary micro include: channel sand, interchannel sand, interchannel mud；

Count the variogram curve of the channel sand, the interchannel sand, the interchannel mud；

Utilize the sedimentary micro data, the earthquake data before superposition and the channel sand, the interchannel sand, the interchannel The variogram curve of mud, in the microfacies-p-wave impedance relationship plate and the constraint of the microfacies-S-wave impedance relationship plate Lower progress prestack geostatistical inversion, establishes sedimentary micro threedimensional model.

Technical field

The present invention relates to oilfield reservoirs to predict field, in particular to a kind of sedimentary facies identification method.

Background technique

In In Oil Field Exploration And Development, most of oil field is in high water cut stage, the remaining oil dispersion in oil field, well spacing density Greatly, exploitation oil gas is difficult.And identify the distribution characteristics of oilfield reservoir sedimentary micro, it easily facilitates and oil gas is excavated, improve The rate of oil and gas recovery.Therefore it provides a kind of sedimentary facies identification method is highly important.

Pertinent literature " research of STF method for predicting reservoir and application based on sedimentary micro " and " fusion earthquake and well logging letter The sedimentary microfacies of delta stochastic modeling of breath is studied " one kind is disclosed based on post-stack seismic data and p-wave impedance and microfacies relationship Sedimentary facies identification method under plate control.

Inventor find the relevant technologies the prior art has at least the following problems:

Summary of the invention

The embodiment of the invention provides a kind of sedimentary facies identification methods, can solve above-mentioned technical problem.Particular technique side Case is as follows:

The embodiment of the invention provides a kind of sedimentary facies identification methods, which comprises carries out ground to target work area Shake acquisition obtains earthquake data before superposition, and drilling well and well logging in the target work area, obtains well data and log data；

Using the well data and the log data, sedimentary micro data are obtained；

P-wave impedance data and S-wave impedance data are obtained using the log data, and count the sedimentary micro number According to, p-wave impedance data and the S-wave impedance data, microfacies-p-wave impedance relationship plate and microfacies-shear wave resistance are obtained Anti- relationship plate；

Using the variogram curve of the sedimentary micro data, the earthquake data before superposition and the sedimentary micro, Prestack geology system is carried out under the microfacies-p-wave impedance relationship plate and the constraint of the microfacies-S-wave impedance relationship plate Meter learns inverting, establishes sedimentary micro threedimensional model；

The microfacies plan view and microfacies sectional view of destination layer position are extracted according to the sedimentary micro threedimensional model, are completed to institute State the identification of target work area sedimentary micro.

In a kind of possible design, the well data includes: core data, and the log data includes: gamma, micro- Electrode and microinverse；

It is described to utilize the well data and the log data, obtain sedimentary micro data, comprising:

Using the gamma, the microelectrode and the microinverse, it is bent to obtain gamma curve, Micro-electrode Log Curves and microinverse Line；

It is explained according to the gamma curve, the Micro-electrode Log Curves and the microinverse curve with reference to the core data Obtain sedimentary micro data.

In a kind of possible design, the log data further include: velocity of longitudinal wave data, shear wave velocity data and sound wave Density data；

It is described to obtain p-wave impedance data and S-wave impedance data using the log data, comprising: to utilize the longitudinal wave Acoustic speed data and the sound wave density data obtain p-wave impedance data；

Using the shear wave velocity data and the sound wave density data, S-wave impedance data are obtained.

It is described to utilize the sedimentary micro data, the earthquake data before superposition and described in a kind of possible design The variogram curve of sedimentary micro, in the microfacies-p-wave impedance relationship plate and the microfacies-S-wave impedance relationship plate Constraint under carry out prestack geostatistical inversion, establish sedimentary micro threedimensional model, comprising:

Go out the heavy of the target work area according to the gamma curve, the Micro-electrode Log Curves and the microinverse curve interpretation Product microfacies, the type of the sedimentary micro includes: channel sand, interchannel sand, interchannel mud；

Count the variogram curve of the channel sand, the interchannel sand, the interchannel mud；

Utilize the sedimentary micro data, the earthquake data before superposition and the channel sand, the interchannel sand, described The variogram curve of interchannel mud, in the microfacies-p-wave impedance relationship plate and the microfacies-S-wave impedance relationship plate Constraint is lower to carry out prestack geostatistical inversion, establishes sedimentary micro threedimensional model.

Technical solution provided in an embodiment of the present invention has the benefit that

Sedimentary facies identification method provided in an embodiment of the present invention, by obtaining deposition using well data and log data Microfacies data, and p-wave impedance data and S-wave impedance data are obtained using log data, statistics well point sedimentary micro data are indulged Acoustic Impedance Data and S-wave impedance data obtain microfacies-impedance relationship plate and microfacies-S-wave impedance relationship plate, recycle The variogram curve of sedimentary micro data, earthquake data before superposition and sedimentary micro, in microfacies-p-wave impedance relationship plate With progress prestack geostatistical inversion under the constraint of microfacies-S-wave impedance relationship plate, sedimentary micro three-dimensional mould can establish Type.As it can be seen that sedimentary facies identification method provided in an embodiment of the present invention is carried out based on earthquake data before superposition, there is shear wave resistance Anti- information can sufficiently react Sedimentary Micro Facies information.Also, the sedimentary facies identification method can obtain microfacies-longitudinal wave Impedance relationship plate and microfacies-S-wave impedance relationship plate, reducing different microfacies, there are multi-solutions in single wave impedance Possibility, improve Sedimentary Micro Facies accuracy of identification, it is smaller to the identification error of sedimentary micro structure.

Detailed description of the invention

To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment Attached drawing is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other Attached drawing.

Fig. 1 is three kinds of sedimentary micro form schematic diagrames in 3rd area grand celebration Yuan Bei provided in an embodiment of the present invention；

Fig. 2 is the microfacies-p-wave impedance relationship plate and microfacies-shear wave in 3rd area grand celebration Yuan Bei provided in an embodiment of the present invention Impedance relationship plate；

Fig. 3 is the variogram schematic diagram of the channel sand in 3rd area grand celebration Yuan Bei provided in an embodiment of the present invention；

Fig. 4 is the variogram schematic diagram of the interchannel sand in 3rd area grand celebration Yuan Bei provided in an embodiment of the present invention；

Fig. 5 is the variogram schematic diagram of the interchannel mud in 3rd area grand celebration Yuan Bei provided in an embodiment of the present invention；

Fig. 6 is the sedimentary micro threedimensional model schematic diagram in 3rd area grand celebration Yuan Bei provided in an embodiment of the present invention；

Fig. 7 is the sedimentary micro sectional view in 3rd area grand celebration Yuan Bei provided in an embodiment of the present invention；

Fig. 8 is the sedimentary micro plan view in 3rd area grand celebration Yuan Bei provided in an embodiment of the present invention.

Specific embodiment

To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached drawing to embodiment party of the present invention Formula is described in further detail.

The embodiment of the invention provides a kind of sedimentary facies identification methods, method includes the following steps:

Step 1 carries out earthquake-capturing to target work area, obtains earthquake data before superposition, and drilling well and survey in target work area Well obtains well data and log data.

Step 2 utilizes well data and log data, acquisition sedimentary micro data.

Step 3 obtains p-wave impedance data and S-wave impedance data using log data, and counts well point sedimentary micro number According to, p-wave impedance data and S-wave impedance data, microfacies-p-wave impedance relationship plate and microfacies-S-wave impedance relational graph are obtained Version.

Step 4, using the variogram of sedimentary micro data, earthquake data before superposition and sedimentary micro, it is vertical in microfacies- Prestack geostatistical inversion is carried out under wave impedance relationship plate and the constraint of microfacies-S-wave impedance relationship plate, establishes deposition Microfacies threedimensional model.

Step 5, the microfacies plan view and microfacies sectional view that destination layer position is extracted according to sedimentary micro threedimensional model, completion pair The identification of target work area sedimentary micro.

Sedimentary facies identification method provided in an embodiment of the present invention, by obtaining deposition using well data and log data Microfacies data, and p-wave impedance data and S-wave impedance data are obtained using log data, statistics well point sedimentary micro data are indulged Acoustic Impedance Data and S-wave impedance data obtain microfacies-impedance relationship plate and microfacies-S-wave impedance relationship plate, recycle The variogram curve of sedimentary micro data, earthquake data before superposition and sedimentary micro, in microfacies-p-wave impedance relationship plate With progress prestack geostatistical inversion under the constraint of microfacies-S-wave impedance relationship plate, sedimentary micro three-dimensional mould can establish Type.As it can be seen that sedimentary facies identification method provided in an embodiment of the present invention is carried out based on earthquake data before superposition, there is shear wave resistance Anti- information can sufficiently react Sedimentary Micro Facies information.Also, the sedimentary facies identification method can obtain microfacies-longitudinal wave Impedance relationship plate and microfacies-S-wave impedance relationship plate, reducing different microfacies, there are multi-solutions in single wave impedance Possibility, improve Sedimentary Micro Facies accuracy of identification, it is smaller to the identification error of sedimentary micro structure.

Wherein, " microfacies-p-wave impedance relationship plate, microfacies-S-wave impedance relationship plate " is it is to be understood that p-wave impedance Relationship plate, the relationship plate of S-wave impedance and sedimentary micro with sedimentary micro, it can be seen that sedimentary micro from the plate Distribution situation.

Just above-mentioned each step is illustrated respectively below:

In embodiments of the present invention, well data includes: core data, and log data includes: gamma, microelectrode and micro- ladder Degree.

Using well data and log data, sedimentary micro data are obtained, comprising:

Step 201 utilizes gamma, microelectrode and microinverse, acquisition gamma curve, Micro-electrode Log Curves and microinverse curve.

Step 202 is explained and is deposited with reference to core data according to gamma curve, Micro-electrode Log Curves and microinverse curve Microfacies data.

Through the above steps, sedimentary micro data can be accurately obtained, establish sedimentary micro threedimensional model to be subsequent It lays the foundation.

" gamma curve, Micro-electrode Log Curves and microinverse curve " in above-mentioned steps 201 well logging common for this field Curve, can by logging operation obtain gamma, microelectrode and microinverse data, obtained after sorting-out in statistics the gamma curve, Micro-electrode Log Curves and microinverse curve.

" core data " in above-mentioned steps 202 it is to be understood that underground core lithology.

In order to accurately obtain p-wave impedance data, microfacies-p-wave impedance relationship plate and microfacies-are established so as to subsequent S-wave impedance relationship plate, log data further include: velocity of longitudinal wave data, shear wave velocity data harmony wave density data；

P-wave impedance data and S-wave impedance data are obtained using log data, comprising: utilize longitudinal wave acoustic speed data Harmony wave density data obtains p-wave impedance data；Using shear wave acoustic speed data harmony wave density data, shear wave resistance is obtained Anti- data.

Wherein, " p-wave impedance data " it also will be understood that are as follows: p-wave impedance data are obtained vertical by using logging operation The multiplication of wave velocity data harmony wave density data can be obtained." S-wave impedance data " it also will be understood that are as follows: S-wave impedance data, The shear wave velocity data harmony wave density data multiplication obtained by using logging operation can be obtained.

It is more accurate to the identification of target work area sedimentary micro in order to guarantee, utilize sedimentary micro data, pre-stack seismic number According to and sedimentary micro variogram curve, in microfacies-p-wave impedance relationship plate and microfacies-S-wave impedance relationship plate Constraint under carry out prestack geostatistical inversion, establish sedimentary micro threedimensional model, comprising:

Step 401, the sedimentary micro for going out target work area according to gamma curve, Micro-electrode Log Curves and microinverse curve interpretation, The type of sedimentary micro includes: channel sand, interchannel sand, interchannel mud；

Step 402, the variogram curve for counting channel sand, interchannel sand, interchannel mud；

Step 403 utilizes the variation of sedimentary micro data, earthquake data before superposition and channel sand, interchannel sand, interchannel mud Function curve carries out prestack geology system under microfacies-p-wave impedance relationship plate and the constraint of microfacies-S-wave impedance relationship plate Meter learns inverting, establishes sedimentary micro threedimensional model.

Wherein " pre-stack seismic " it also will be understood that are as follows: the seismic data before being superimposed through common reflection point, by target work area Shake acquisition obtains, and compared with conventional poststack seismic data, reservoir information is more abundant.

The type of " sedimentary micro " in above-mentioned steps 401 is by gamma curve, Micro-electrode Log Curves and microinverse curve Morphological feature explain.

The variogram curve of channel sand, interchannel sand, interchannel mud in above-mentioned steps 402, is explained based on well point The variation value of different distance h point is calculated using lower column variation value formula for sedimentary micro, then using it is multiple point variation be worth into Row fitting obtains variogram curve, to reflect different sedimentary micro spatial distribution characteristics.

Wherein, h is the distance between well and well, and unit is rice；(x_i,y_i)、(x_i+h,y_i) it is well point coordinate；Z(x_i,y_i)、 Z(x_i+h,y_i) it is well point microfacies；N (h) is to meet the well point number that distance is h；γ (h) is variation value.

The specific implementation process of the model of sedimentary micro is established in above-mentioned steps 403 are as follows: firstly, being based on well point sedimentary micro Data obtain embryo deposit microfacies model using prestack geostatistical inversion algorithm by variogram curve.Then, it ties Microfacies-p-wave impedance relationship plate and microfacies-S-wave impedance relationship plate are closed, is longitudinal wave resistance by embryo deposit microfacies model conversation Anti- model and S-wave impedance model, by p-wave impedance model, S-wave impedance model compared with seismic data cube carries out correlation, phase Closing property, which is met the requirements, is set to final result output for embryo deposit microfacies model, otherwise recalculates embryo deposit microfacies mould Type, until meeting correlation requirement.

Hereinafter the present invention will be further described through by specific embodiment.