Advantageous region division method based on shale gas exploration prediction
Technical Field
The invention belongs to the technical field of exploration, and particularly relates to a beneficial zoning method based on shale gas exploration prediction.
Background
Shale gas is used as an unconventional natural gas which is rich in organic matters and is stored in the shale rich in organic matters, and development of the shale gas has very important significance for global energy structure transformation.
In oil exploration, accurately evaluating and dividing the favorable region has important significance for the development of oil and gas fields. For shale gas, the division of the beneficial zone is more important.
Conventional classification utilizes reservoir characteristic analysis, construction and fault evaluation, resource quantity estimation, resource quality evaluation, cost benefit analysis and the like, but the evaluation on shale gas and dense gas is not specific and is not closely related to productivity.
Disclosure of Invention
In order to solve the problems, the invention provides a beneficial zone dividing method based on shale gas exploration prediction, which evaluates shale gas and productivity in close relation to finish accurate division of beneficial zones and accurately guides production of shale gas in later period.
In order to achieve the aim, the technical scheme adopted by the invention is that the beneficial region dividing method based on shale gas exploration prediction comprises the following steps:
s100, acquiring three-dimensional seismic time domain gather data and well data of a work area and ground elevation data;
S201, performing time migration processing on the three-dimensional seismic time domain gather data to obtain time migration data of a time domain;
S202, converting three-dimensional seismic time domain gather data into angle gather data;
s203, obtaining the correlation of the density, the speed or the wave impedance related to the earthquake in the well data, and determining the optimal correlation parameter;
s300, obtaining an elastic parameter body through inversion by utilizing the horizon data of the target layer and combining with the optimal correlation parameters;
S401, interpreting the plane distribution of faults according to the phase characteristics of the horizon data of the target layer to obtain a fault plane distribution diagram, grading to obtain faults of all levels, converting the horizon data of the target layer into depth domain data by the product of an elastic parameter body and a speed body by utilizing the speed body in the elastic parameter body, and adding the depth domain data with ground elevation data to obtain a buried depth result;
S402, obtaining a class I continuous reservoir thickness result by adopting a sampling point statistical method;
S403, obtaining other parameters of the I-type continuous reservoir by using a fitting mode, wherein the other parameters comprise T0C, total air content, porosity and water saturation;
S404, calculating Young modulus and Poisson 'S ratio by using an elastic parameter body, and then obtaining a brittleness index by taking an average value of the Young modulus and Poisson' S ratio;
s500, dividing the reservoir into favorable areas with different levels by using the burial depth result, the class I continuous reservoir thickness result, the T0C, the total air content, the porosity and the water saturation brittleness index result value and the set reservoir division evaluation standard value.
Further, the horizon data of the target layer is picked up according to the wave group characteristics, the wave group relation and the reflection characteristics on the time offset data, and the horizon data of the target layer is obtained.
Further, obtaining the correlation of the density, the velocity or the wave impedance and the longitudinal and transverse wave velocity ratio related to the earthquake in the well data comprises the steps of selecting the well data of the I-type continuous reservoir section, carrying out histogram analysis on the well data, the density, the velocity or the wave impedance and the longitudinal and transverse wave velocity ratio, and determining which parameter or operation result of the parameter of the earthquake can separate the I-type continuous reservoir from the non-reservoir through the histogram analysis to be used as the optimal correlation parameter.
Further, by utilizing the horizon data of the target layer and combining with the optimal correlation parameters, obtaining an elastic parameter body through inversion, and determining that the elastic parameter body directly related to the well and the earthquake comprises density, longitudinal wave speed and transverse wave speed, wherein the longitudinal wave speed and transverse wave speed ratio is the operation result of the longitudinal wave speed or the transverse wave speed.
Further, inverting to obtain a speed body in the elastic parameter body, obtaining a target layer, adding the obtained horizon of the depth domain with the ground elevation to obtain the plane distribution of the interpretation fault of the burial depth result according to the phase characteristics of the target layer of the offset data, and obtaining a fault plane distribution map.
Further, when the sample point statistics method is adopted to obtain the class I continuous reservoir thickness result, the class I continuous reservoir thickness result is obtained by extracting the thickness of the sample point.
Further, other parameters of the type I continuous reservoir including T0C, total gas content, porosity and water saturation are obtained by using a fitting mode, wherein fitting is performed according to the uphole T0C, total gas content, porosity and water saturation and the uphole density or speed, which correlation degree of the uphole T0C, total gas content, porosity and water saturation and the speed or density is the highest is found, a corresponding relation is found, and the T0C, total gas content, porosity and water saturation are calculated by using the relation.
Further, the evaluation criteria of the favorable region is that the burial depth is less than 500 meters, the thickness of the class I continuous reservoir is more than or equal to 10 meters, the TOC is more than 3 percent, the total air content is more than 3 percent, the porosity is more than 4 percent, the water saturation is less than or equal to 40 percent, and the brittleness index is more than 55 percent;
The evaluation criteria of the second kind of favorable area are that the burial depth is less than 500 meters, the thickness of the class I continuous reservoir is less than 10 meters, TOC is more than 2 percent and less than 3 percent, the total air content is more than 2 percent and less than 3 percent, the porosity is more than 3 percent and less than 4 percent, the water saturation is more than 40 percent, and the brittleness index is more than 45 percent and less than 55 percent.
The beneficial effect of adopting this technical scheme is:
The method utilizes the thickness, the water saturation, the TOC, the total air content, the porosity and the brittleness index of the type I continuous reservoir predicted by the earthquake to divide the favorable region and select related parameters, the type I continuous reservoir thickness and the related parameters are closely related to the production energy, the application effect is very good, the comprehensive evaluation is carried out by adopting the high-quality shale predicted by the earthquake in the past, and the evaluation is carried out without adopting the type I continuous reservoir thickness and the related parameters. In addition, the shale gas does not adopt earthquake to predict the water saturation, and the invention increases the water saturation parameter which has a control factor in the shale gas reservoir. In general, the method considers the thickness of the I-type continuous reservoir and related parameters, is more reasonable and reliable, and additionally considers the water saturation parameters to participate in evaluation, so that the method has more practicability and reliability.
Drawings
FIG. 1 is a flow diagram of an advantageous compartmentalization method based on shale gas exploration prediction in accordance with the present invention;
FIG. 2 is a histogram of longitudinal wave velocity, transverse-to-transverse wave velocity ratio and density in an embodiment of the invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings, in order to make the objects, technical solutions and advantages of the present invention more apparent.
In this embodiment, referring to fig. 1, the invention provides a beneficial region division method based on shale gas exploration prediction, which includes the steps of:
s100, acquiring three-dimensional seismic time domain gather data and well data of a work area and ground elevation data;
S201, performing time migration processing on the three-dimensional seismic time domain gather data to obtain time migration data of a time domain;
Specifically, the horizon data of the target layer is picked up on the time offset data according to the wave group characteristics, the wave group relation and the reflection characteristics, and the horizon data of the target layer is obtained.
S202, converting the three-dimensional seismic time domain gather data into angle gather data.
S203, obtaining the correlation of the density, the speed or the wave impedance related to the earthquake in the well data, and determining the optimal correlation parameter;
Specifically, the method for acquiring the correlation of the density, the velocity or the wave impedance and the longitudinal and transverse wave velocity ratio related to the earthquake in the well data comprises the steps of selecting the well data of a type I continuous reservoir section, carrying out histogram analysis on the well data, the density, the velocity or the wave impedance and the longitudinal and transverse wave velocity ratio, and determining which parameter or operation result of the parameter of the well and the earthquake can separate the type I continuous reservoir from a non-reservoir through the histogram analysis to be used as the parameter for determining the optimal correlation.
S300, obtaining an elastic parameter body through inversion by utilizing the horizon data of the target layer and combining with the optimal correlation parameters;
And (3) obtaining an elastic parameter body through inversion by utilizing the horizon data of the target layer and combining with the optimal correlation parameter, and determining that the elastic parameter body directly related to the well and the earthquake comprises density, longitudinal wave speed and transverse wave speed, wherein the longitudinal wave speed and transverse wave speed ratio is the longitudinal wave speed or transverse wave speed operation result.
As shown in fig. 2, it is obvious from the graph that the density can better distinguish the type I continuous reservoir from the non-reservoir section, and the longitudinal wave velocity and the transverse wave velocity are more overlapped than the type I continuous reservoir and the non-reservoir section, and cannot be effectively separated.
S401, interpreting the plane distribution of faults according to the phase characteristics of the horizon data of the target layer to obtain a fault plane distribution diagram, grading to obtain faults of all levels, converting the horizon data of the target layer into depth domain data by using a velocity body in the elastic parameter body through the product of the elastic parameter body and the velocity body, and adding the depth domain data with ground elevation data to obtain a buried depth result.
Specifically, inverting to obtain a velocity body in the elastic parameter body, obtaining a target layer, adding the obtained horizon of the depth domain and the ground elevation to obtain the plane distribution of the interpretation fault of the burial depth result according to the phase characteristics of the target layer of the offset data, and obtaining a fault plane distribution map.
S402, obtaining a class I continuous reservoir thickness result by adopting a sampling point statistical method;
specifically, when the sample point statistics method is adopted to obtain the layer thickness result of the class I continuous reservoir, the class I continuous reservoir thickness result is obtained by extracting the thickness of the sample point and the longitudinal wave velocity body.
S403, obtaining other parameters of the I-type continuous reservoir by using a fitting mode, wherein the other parameters comprise T0C, total air content, porosity and water saturation;
Specifically, other parameters of the I-type continuous reservoir are obtained by using a fitting mode, including T0C, total air content, porosity and water saturation, wherein fitting is performed according to the uphole T0C, total air content, porosity and water saturation and the uphole density or speed, which correlation degree of the uphole T0C, total air content, porosity and water saturation and the speed or density is the highest is found, a corresponding relational expression is found, and the T0C, total air content, porosity and water saturation are calculated by using the relational expression.
S404, calculating Young modulus and Poisson 'S ratio by using an elastic parameter body, and then obtaining a brittleness index by taking an average value of the Young modulus and Poisson' S ratio;
specific:
Young's modulus:
Poisson ratio:
s500, dividing the reservoir into favorable areas with different levels by using the burial depth result, the class I continuous reservoir thickness result, the T0C, the total air content, the porosity and the water saturation brittleness index result value and the set reservoir division evaluation standard value.
Specific:
The evaluation standard of the favorable region is that the burial depth is less than 500 meters, the thickness of the I-type continuous reservoir is more than or equal to 10 meters, the TOC is more than 3 percent, the total air content is more than 3 percent, the porosity is more than 4 percent, the water saturation is less than or equal to 40 percent, and the brittleness index is more than 55 percent;
The evaluation criteria of the second kind of favorable area are that the burial depth is less than 500 meters, the thickness of the class I continuous reservoir is less than 10 meters, TOC is more than 2 percent and less than 3 percent, the total air content is more than 2 percent and less than 3 percent, the porosity is more than 3 percent and less than 4 percent, the water saturation is more than 40 percent, and the brittleness index is more than 45 percent and less than 55 percent.
The foregoing has shown and described the basic principles and main features of the present invention and the advantages of the present invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, and that the above embodiments and descriptions are merely illustrative of the principles of the present invention, and various changes and modifications may be made without departing from the spirit and scope of the invention, which is defined in the appended claims. The scope of the invention is defined by the appended claims and equivalents thereof.