CN103163553B - Earthquake hydrocarbon detection method and device based on multiple pore medium model - Google Patents

Abstract

The invention relates to the technical field of geophysical technology, in particular to a method and a device for detecting seismic hydrocarbons based on a multi-pore medium model. The method comprises the following steps: obtaining a data set relating to subsurface reservoir rock and deriving reservoir information therefrom, the reservoir information comprising: reservoir environment and physical properties, petrophysical properties, and micro-pore structure information; obtaining model parameters required by a multi-pore medium model according to reservoir information, and establishing a rock physical template of the fluid saturated rock; carrying out amplitude preservation processing on the seismic data, extracting an angle gather, and carrying out synchronous inversion on the angle gather before stacking to obtain elastic parameters before stacking; projecting the prestack elastic parameters onto the rock physical template to form a plurality of projection data points, and calculating the porosity and/or saturation of the measured medium; the hydrocarbon detection is carried out by judging the oil-gas distribution of the underground reservoir according to the porosity and/or the saturation of the detected medium, the oil-gas distribution of an unconventional reservoir can be directly predicted, and the multi-solution property of the hydrocarbon detection is effectively reduced.

Description

Based on earthquake detecting method of hydrocarbon and the device of multiple pore medium model

Technical field

The present invention relates to technical field of physical geography, particularly relate to the earthquake detecting method of hydrocarbon based on multiple pore medium model and device.

Background technology

In the past few decades, seismic exploration technique there occurs earth-shaking change, seismic interpretation also experienced by multiple developing stage such as structure elucidation, formation lithology explanation, production seismics explanation, Rock physical analysis simultaneously, and explains future development towards EARTHQUAKE QUANTIFICATION.But the target of seismic prospecting is all oil gas all the time, the pursuit of this constancy continues to promote the progress of hydrocarbon indication theory and technology.

Hydrocarbon indication technology refers to the method utilizing seismic reflection or refracting feature look for oil and gas to hide.In earthquake reflected wave exploration, the seismic signal that underground returns not only comprises the tectonic information of reflection subsurface interface, also containing the amplitude information reflecting formation lithology and fluid.In seismic interpretation in the past, people are confined to structure elucidation, and along with modern real amplitude acquisition and processing is greatly strengthened, seismic amplitude has become the Main Basis identifying potential hydrocarbon reservoir.

Early 1980s, first propose to utilize Amplitudeversusangle identification " bright spot " type gas reservoir, his working mark uses the appearance of the change of amplitude offset distance (Amplitude Versus Offset, AVO) technology really.AVO technology makes seismic amplitude explain and turns to prestack gradually from poststack, directly can carry out hydrocarbon indication and predicting reservoir hydrocarbon occurrence.Subsequently, many scholars propose different approximate expressions to Zuo Pulizi (Zoeppritz) equation that reflection interfacial energy is distributed.Can extract corresponding AVO attribute by these formula, and the combination of these attributes derives a series of detecting factor relevant with fluid, as P*G, LMR parameter, K-μ fluid factor etc.In addition, the prestack elastic parameter inversion technology being representative with elastic impedance inverting and prestack Simultaneous Retrieving improves the ability of hydrocarbon indication further.

Another essence that seismic event is propagated in rock is exactly the attenuation by absorption of energy, causes waveform to change.In general, hydrocarbon-bearing pool has lower quality factor, and seismic reflection energy attenuation may be obvious, and therefore people develop some seismic attenuation technology to detect the oil-gas possibility of reservoir.

In the last few years, the fast development of earthquake rock physics technology facilitated hydrocarbon indication and was explained, particularly with the EARTHQUAKE QUANTIFICATION forecasting techniques that rock physics template is representative by former qualitative description quantification of progressively marching toward.First first this concept of rock physics template proposed in 2003 by Ф Degaar and Avseth, develops various forms of quantitative interpretation template subsequently, comprise AI-V _p/ V _stemplate, PGT template etc.Geology and earthquake are closely connected by rock physics mould plate technique, are a kind of very important instruments, greatly reduce the risk in seismic prospecting and near-gar tesistance.

Along with socio-economic development is to the sustainable growth of hydrocarbon resources demand, conventional gas and oil resource is day by day exhausted, people have invested aboundresources more sight and the huge unconventional petroleum resources of Exploration Potential, as compact sandstone gas, shale gas, fine and close oil, shale wet goods.Ooze owing to having low hole, the strong feature such as nonuniformity, complex hydrocarbon water relation, the seismic response that these unconventionaloil pool are hidden is more weak, and disguise is very strong, and current hydrocarbon indication technology also exists deficiency to a certain extent.

AVO technology, prestack elastic parameter inversion, seismic attenuation etc. with seismic properties be the hydrocarbon indication technology that instructs can only qualitative recognition oily abnormal, because between unconventionaloil pool Tibetan with country rock, elastic parameter difference is less, these technology are difficult to effectively differentiate to there is obvious multi-solution.For rock physics mould plate technique, its core is petrophysical model, but the unconventional reservoir oozed for low hole at present does not form comparatively perfect rock physics theoretical model, if high hole sand shale model is carried out EARTHQUAKE QUANTIFICATION prediction and can be there is larger risk in rashly using.

Summary of the invention

The embodiment of the present invention provides earthquake detecting method of hydrocarbon based on multiple pore medium model and device, for solving the exploitation and research of hiding unconventionaloil pool in prior art.

A kind of earthquake detecting method of hydrocarbon based on multiple pore medium model in the embodiment of the present invention, wherein, described method comprises:

Obtain the data set relevant to subsurface reservoir rock, and obtain reservoir information from described data centralization, described reservoir information comprises: reservoir environment transitivity, petrophysical property and micropore structure information;

Obtain the model parameter needed for multiple pore medium model according to described reservoir information, set up the rock physics template of fluid saturated rock based on described multiple pore medium model;

Extract angle gathers after relative amplitude preserved processing is carried out to geological data, by described angle gathers by prestack Simultaneous Retrieving, obtain prestack elastic parameter;

Described prestack elastic parameter is projected in described rock physics template and form multiple data for projection point, and calculate factor of porosity and/or the saturation degree of measured medium;

Judge that subsurface reservoir hydrocarbon occurrence is to carry out hydrocarbon indication by the factor of porosity of described measured medium and/or saturation degree.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, described data set comprises: geologic report, landwaste record, core data, logging trace.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, described rock physics template is the combination of any two described physical quantitys in elastic parameter; Described elastic parameter at least comprises following physical quantity: velocity of longitudinal wave V _p, shear wave velocity V _s, p-wave impedance Z _p, S-wave impedance Z _s, elastic impedance EI, transformed wave elastic impedance PSEI, P-S wave velocity ratio V _p/ V _s, Poisson ratio v, bulk modulus K, modulus of shearing μ, Young modulus E, compressional wave modulus P and Lame's constant λ.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, factor of porosity and/or the saturation degree of described calculating measured medium comprise: search the template net lattice point that each described data for projection point is nearest, and then obtain factor of porosity and/or the saturation degree of corresponding described measured medium.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, described reservoir environment transitivity comprises: temperature, load pressure, pore pressure, layer position, lithology, mineralogical composition, mineral content, porosity ranges, saturation degree scope, pore fluid type.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, described petrophysical property comprises: the velocity of longitudinal wave V of ultrasonic measurement _p, shear wave velocity V _sand the corresponding relation between rock porosity, rock saturation degree.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, the described rock physics template setting up fluid saturated rock based on multiple pore medium model comprises: according to the scale value of described porosity ranges and saturation degree range set factor of porosity and saturation degree, described multiple pore medium model is utilized to calculate any two described physical quantitys in the described elastic parameter of each group factor of porosity and fluid saturated rock corresponding to the scale value of saturation degree, draw the X plot of any two described physical quantitys in described elastic parameter and indicate the scale value of described factor of porosity and saturation degree, form described rock physics template.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, the described elastic parameter of the fluid saturated rock that the described scale value utilizing described multiple pore medium model to calculate each group factor of porosity and saturation degree is corresponding specifically comprises:

Obtain three parameter of pore structure and three saturation parameters; Three described parameter of pore structure are respectively: aspect ratio α, scale factor x and percent continuity ξ; Three described saturation parameters are respectively: total water saturation S _w, intercommunicating pore water saturation S _cwand isolated hole water saturation S _iw; Pass between three described saturation parameters is: S _w=S _cwξ+S _iw(1-ξ);

According to formula: φ _iso=φ (1-ξ), calculates the factor of porosity φ of solid matrix _iso, wherein, φ is the total porosity of rock; Described solid matrix comprises: isolated hole and mineral grain;

Again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{\max }=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}})K_{\min }+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_{\mathrm{hc}}{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_w{P}_i^{\′\′}}{(1-{\mathrm{\φ}}_{\mathrm{iso}})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′\′}}\\ {\widehat{\mathrm{\μ}}}_{\mathrm{mat}}=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}}){\mathrm{\μ}}_{\min }}{(1-{\mathrm{\φ}}_{\mathrm{iso}})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′\′}}\end{array}\right.$

Calculate the elastic modulus of solid matrix, wherein: K _minand μ _minbe respectively the average external volume modulus of described mineral grain and average modulus of shearing, K _hcand K _wbe respectively the bulk modulus of hydro carbons and water, P with Q is the geometric factor relevant with pore morpholohy;

According to formula: φ _con=φ ξ, calculates the factor of porosity φ of dry skeleton _con; Described dry skeleton comprises: isolated hole, intercommunicating pore and mineral grain;

Again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{\mathrm{dry}}=\frac{(1-\mathrm{\φ})K_{\min }+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_{\mathrm{hc}}{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_w{P}_i^{\′\′}}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′\′}+{\mathrm{\φ}}_{\mathrm{con}}\underset{k=1}{\overset{n}{\mathrm{\Σ}}}{x}_k{\tilde{P}}_k}\\ {\widehat{\mathrm{\μ}}}_{\mathrm{dry}}=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}}){\mathrm{\μ}}_{\min }}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′\′}+{\mathrm{\φ}}_{\mathrm{con}}\underset{k=1}{\overset{n}{\mathrm{\Σ}}}{x}_k{\tilde{Q}}_k}\end{array}\right.$

Calculate the elastic modulus of described dry skeleton; Wherein,

According to the elastic modulus of described solid matrix and the elastic modulus of dry skeleton, calculate the overall elasticity modulus of described fluid saturated rock:

$\left\{\begin{array}{c}{K}_{\mathrm{sat}}={\hat{K}}_{\mathrm{dry}}+\frac{{(1-{\hat{K}}_{\mathrm{dry}}/{\hat{K}}_{\mathrm{mat}})}^2}{{\mathrm{\φ}}_{\mathrm{con}}/K_{\mathrm{fl}}^{*}+(1-{\mathrm{\φ}}_{\mathrm{com}})/{\hat{K}}_{\mathrm{mat}}-{\hat{K}}_{\mathrm{dry}}/{\hat{K}}_{\mathrm{mat}}^2}\\ {\mathrm{\μ}}_{\mathrm{sat}}={\widehat{\mathrm{\μ}}}_{\mathrm{dry}}\end{array}\right.$

Wherein,

$\frac{1}{K_{\mathrm{fl}}^{*}}=\frac{1-S_{\mathrm{cw}}}{K_{\mathrm{hc}}}+\frac{S_{\mathrm{cw}}}{K_w}$

Again according to formula: ρ _sat=(1-φ) ρ _min+ φ (1-S _w) ρ _hc+ φ S _wρ _wcalculate the density p of described fluid saturated rock _sat, wherein: ρ _minfor the average density of described mineral grain, ρ _hcand ρ _wbe respectively the density of hydro carbons and water.

Physical quantity described in other of described fluid saturated rock is by the elastic modulus K of described solid matrix _sat, dry skeleton elastic modulus μ _satwith the density p of described fluid saturated rock _satthese three gain of parameters.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, described model parameter comprises: mineral grain parameter, fluid parameter, factor of porosity parameter, parameter of pore structure, saturation parameters; Wherein, mineral grain parameter comprises: average external volume modulus K _min, average shear modulus μ _min, average density ρ _min; Fluid parameter comprises: hydro carbons bulk modulus K _hc, hydro carbons density p _hc, water volume modulus K _w, water-mass density ρ _w.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, described mineral grain parameter calculates according to Hill model.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, described fluid parameter obtains according to Flag program computation.

The above-mentioned earthquake detecting method of hydrocarbon based on multiple pore medium model, wherein, the described aspect ratio α in described parameter of pore structure and described scale factor x obtains according to pore texture characterization technique; Described percent continuity ξ is according to velocity of longitudinal wave V _pand the relation between saturation degree is estimated to obtain.

The embodiment of the present invention also provides a kind of earthquake hydrocarbon indication device based on multiple pore medium model, and wherein, described device comprises:

Data capture unit, for obtaining the data set relevant to subsurface reservoir rock, and obtain reservoir information from described data centralization, described reservoir information comprises: reservoir environment transitivity, petrophysical property and micropore structure information;

Template sets up unit, for obtaining the model parameter needed for multiple pore medium model according to described reservoir information, sets up the rock physics template of fluid saturated rock based on described multiple pore medium model;

Data processing unit, extracts angle gathers after carrying out relative amplitude preserved processing to geological data, by described angle gathers by prestack Simultaneous Retrieving, obtains prestack elastic parameter;

Projecting cell, forms multiple data for projection point for described prestack elastic parameter being projected in described rock physics template, and calculates factor of porosity and/or the saturation degree of measured medium;

By the factor of porosity of described measured medium and/or saturation degree, judging unit, for judging that subsurface reservoir hydrocarbon occurrence is to carry out hydrocarbon indication.

The above-mentioned earthquake hydrocarbon indication device based on multiple pore medium model, wherein, described projecting cell specifically for searching the nearest template net lattice point of each described data for projection point, and then obtains factor of porosity and/or the saturation degree of corresponding described measured medium.

The above-mentioned earthquake hydrocarbon indication device based on multiple pore medium model, wherein, described template sets up unit specifically for the scale value according to porosity ranges and saturation degree range set factor of porosity and saturation degree, described multiple pore medium model is utilized to calculate any two described physical quantitys in the elastic parameter of each group factor of porosity and fluid saturated rock corresponding to the scale value of saturation degree, draw the X plot of any two described physical quantitys in described elastic parameter and indicate the scale value of described factor of porosity and saturation degree, forming described rock physics template.

The above-mentioned earthquake hydrocarbon indication device based on multiple pore medium model, wherein, described template is set up unit and is also comprised: computing unit, for obtaining three parameter of pore structure and three saturation parameters; Three described parameter of pore structure are respectively: aspect ratio α, scale factor x and percent continuity ξ; Three described saturation parameters are respectively: total water saturation S _w, intercommunicating pore water saturation S _cwand isolated hole water saturation S _iw; Pass between three described saturation parameters is: S _w=S _cwξ+S _iw(1-ξ);

Described computing unit is according to formula: φ _iso=φ (1-ξ), calculates the factor of porosity φ of solid matrix _iso, wherein, φ is the total porosity of rock; Described solid matrix comprises: isolated hole and mineral grain;

Described computing unit is again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{\max }=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}})K_{\min }+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_{\mathrm{hc}}{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_w{P}_i^{\′\′}}{(1-{\mathrm{\φ}}_{\mathrm{iso}})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′\′}}\\ {\widehat{\mathrm{\μ}}}_{\mathrm{mat}}=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}}){\mathrm{\μ}}_{\min }}{(1-{\mathrm{\φ}}_{\mathrm{iso}})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′\′}}\end{array}\right.$

Calculate the elastic modulus of solid matrix, wherein: K _minand μ _minbe respectively the average external volume modulus of described mineral grain and average modulus of shearing, K _hcand K _wbe respectively the bulk modulus of hydro carbons and water, P with Q is the geometric factor relevant with pore morpholohy;

Described computing unit is also according to formula: φ _con=φ ξ, calculates the factor of porosity φ of dry skeleton _con; Described dry skeleton comprises: isolated hole, intercommunicating pore and mineral grain;

Described computing unit is according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{\mathrm{dry}}=\frac{(1-\mathrm{\φ})K_{\min }+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_{\mathrm{hc}}{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_w{P}_i^{\′\′}}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′\′}+{\mathrm{\φ}}_{\mathrm{con}}\underset{k=1}{\overset{n}{\mathrm{\Σ}}}{x}_k{\tilde{P}}_k}\\ {\widehat{\mathrm{\μ}}}_{\mathrm{dry}}=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}}){\mathrm{\μ}}_{\min }}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′\′}+{\mathrm{\φ}}_{\mathrm{con}}\underset{k=1}{\overset{n}{\mathrm{\Σ}}}{x}_k{\tilde{Q}}_k}\end{array}\right.$

Calculate the elastic modulus of described dry skeleton;

Described computing unit also according to the elastic modulus of described solid matrix and the elastic modulus of dry skeleton, calculates the overall elasticity modulus of described fluid saturated rock:

$\left\{\begin{array}{c}{K}_{\mathrm{sat}}={\hat{K}}_{\mathrm{dry}}+\frac{{(1-{\hat{K}}_{\mathrm{dry}}/{\hat{K}}_{\mathrm{mat}})}^2}{{\mathrm{\φ}}_{\mathrm{con}}/K_{\mathrm{fl}}^{*}+(1-{\mathrm{\φ}}_{\mathrm{com}})/{\hat{K}}_{\mathrm{mat}}-{\hat{K}}_{\mathrm{dry}}/{\hat{K}}_{\mathrm{mat}}^2}\\ {\mathrm{\μ}}_{\mathrm{sat}}={\widehat{\mathrm{\μ}}}_{\mathrm{dry}}\end{array}\right.$

Wherein,

$\frac{1}{K_{\mathrm{fl}}^{*}}=\frac{1-S_{\mathrm{cw}}}{K_{\mathrm{hc}}}+\frac{S_{\mathrm{cw}}}{K_w}$

Described computing unit is again according to formula: ρ _sat=(1-φ) ρ _min+ φ (1-S _w) ρ _hc+ φ S _wρ _wcalculate the density p of described fluid saturated rock _sat, wherein: ρ _minfor the average density of described mineral grain, ρ _hcand ρ _wbe respectively the density of hydro carbons and water.

Physical quantity described in other of described fluid saturated rock is by the elastic modulus K of described solid matrix _sat, dry skeleton elastic modulus μ _satwith the density p of described fluid saturated rock _satthese three gain of parameters.

The present invention contrasts prior art and has following beneficial effect: the multiple pore medium model of proposition embodies Complicated Pore Structures and microfluid uneven distribution to the impact of rock elasticity feature, be more suitable for describing unconventional reservoir particularly compact reservoir, the inversion result therefore based on this model is more reliable.Multiple pore medium model is combined with rock physics mould plate technique, and adopt Template Map method to be converted into by prestack seismic attributes more can the factor of porosity of inverting reservoir oil-gas possibility and saturation parameters, the hydrocarbon occurrence of direct prediction unconventional reservoir, effectively reduces the multi-solution of hydrocarbon indication.Template Map method counting yield is very high, is more suitable for industrial applications.

Accompanying drawing explanation

Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:

The process flow diagram of the earthquake detecting method of hydrocarbon based on multiple pore medium model that Fig. 1 provides for the embodiment of the present invention;

The structural representation of the earthquake hydrocarbon indication device based on multiple pore medium model that Fig. 2 provides for the embodiment of the present invention;

A kind of tight sand micropore structure nanometer CT scan schematic diagram that Fig. 3 provides for the embodiment of the present invention;

The parameter of pore structure schematic diagram of the experimentally data analysis microcosmic that Fig. 4 provides for the embodiment of the present invention;

The Z set up based on multiple pore medium model that Fig. 5 provides for the embodiment of the present invention _p-V _p/ V _stemplate;

The p-wave impedance Z that angle gathers is obtained by prestack Simultaneous Retrieving that Fig. 6 provides for the embodiment of the present invention _pwith P-S wave velocity ratio V _p/ V _sdiagrammatic cross-section;

The prestack elastic parameter that Fig. 7 provides for the embodiment of the present invention is at Z _p-V _p/ V _sprojection in template;

The earthquake hydrocarbon indication diagrammatic cross-section that Fig. 8 provides for the embodiment of the present invention.

Embodiment

For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with accompanying drawing, the embodiment of the present invention is described in further detail.At this, schematic description and description of the present invention is for explaining the present invention, but not as a limitation of the invention.

Fig. 1 gives the process flow diagram of the earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides, and described method specifically comprises:

Step 101, obtain the data set relevant to subsurface reservoir rock, and obtain reservoir information from described data centralization, described reservoir information comprises: reservoir environment transitivity, petrophysical property and micropore structure information; Preferably, described data set comprises: geologic report, landwaste record, core data, logging trace etc.Be can obtain the available data relevant from existing subsurface reservoir rock at this data set, the present invention is the associated reservoir information extracting data centralization, and for setting up the rock physics template of fluid saturated rock.

Step 102, obtains the model parameter needed for multiple pore medium model according to described reservoir information, sets up the rock physics template of fluid saturated rock based on described multiple pore medium model; Concrete, this step sets up the rock physics template of rock, and the model parameter needed for multiple pore medium model wherein can calculate or estimate to obtain from reservoir information.

Step 103, extracts angle gathers after carrying out relative amplitude preserved processing to geological data, by described angle gathers by prestack Simultaneous Retrieving, obtains prestack elastic parameter; Concrete, the known geological data obtained is carried out the method that relative amplitude preserved processing can adopt prior art, obtains angle gathers with this, and adopted by angle gathers existing prestack Simultaneous Retrieving technology to obtain prestack elastic parameter.

Step 104, projects to described prestack elastic parameter in described rock physics template and forms multiple data for projection point, and calculate factor of porosity and/or the saturation degree of measured medium;

By the factor of porosity of described measured medium and/or saturation degree, step 105, judges that subsurface reservoir hydrocarbon occurrence is to carry out hydrocarbon indication.Concrete, when judging subsurface reservoir hydrocarbon occurrence, factor of porosity or the saturation degree of measured medium can be observed, obtain subsurface reservoir hydrocarbon occurrence thus; Preferably, in order to reduce the simple multi-solution utilizing the factor of porosity of measured medium or saturation degree to judge, factor of porosity can also be utilized to be multiplied by saturation degree and to obtain air content to judge, make the reservoir hydrocarbons distribution results that obtains more accurate.

At this, the multiple pore medium model that the embodiment of the present invention proposes embodies Complicated Pore Structures and microfluid uneven distribution to the impact of rock elasticity feature, be more suitable for describing unconventional reservoir particularly compact reservoir, the inversion result therefore based on this model is more reliable.Multiple pore medium model is combined with rock physics mould plate technique, and the factor of porosity and saturation parameters that adopt Template Map method to be converted into by prestack seismic attributes more to reflect reservoir oil-gas possibility, the hydrocarbon occurrence of direct prediction unconventional reservoir, effectively reduces the multi-solution of hydrocarbon indication.Template Map method counting yield is very high, is more suitable for industrial applications.

The earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides, preferably, the model parameter needed for multiple pore medium model comprises: mineral grain parameter, fluid parameter, factor of porosity parameter, parameter of pore structure, saturation parameters; Wherein, mineral grain parameter comprises: average external volume modulus K _min, average shear modulus μ _min, average density ρ _min; Fluid parameter comprises: hydro carbons bulk modulus K _hc, hydro carbons density p _hc, water volume modulus K _w, water-mass density ρ _w.Above-mentioned parameter can by calculating or estimating to obtain; Concrete, according to the difference of media type, adopt computing method to the parameter conveniently calculated, adopt method of estimation to obtain to calculating the parameter being difficult to obtain, concrete calculating or method of estimation can adopt the state of the art.

Preferably, described mineral grain parameter calculates according to Hill model; In another preferred embodiment, described fluid parameter obtains according to Flag program computation.Concrete, Hill model and Flag program (UH & CSM, Fluid/DHI association) they are the prior aries that this area is commonly used, and do not repeat them here.

The earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides preferably, is the combination of any two described physical quantitys in elastic parameter according to described rock physics template; Described elastic parameter at least comprises following physical quantity: velocity of longitudinal wave V _p, shear wave velocity V _s, p-wave impedance Z _p, S-wave impedance Z _s, elastic impedance EI, transformed wave elastic impedance PSEI, P-S wave velocity ratio V _p/ V _s, Poisson ratio v, bulk modulus K, modulus of shearing μ, Young modulus E, compressional wave modulus P and Lame's constant λ.Concrete, it is generally acknowledged that fluid saturated rock is isotropic, only have any two physical quantitys to be independently in elastic parameter, such as velocity of longitudinal wave V _pwith shear wave velocity V _s, or bulk modulus K and modulus of shearing μ, and other physical quantitys in elastic parameter just can be transformed by other two physical quantitys arbitrary in this elastic parameter.

The earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides, preferably, the described rock physics template setting up fluid saturated rock based on multiple pore medium model comprises: according to the scale value of described porosity ranges and saturation degree range set factor of porosity and saturation degree, described multiple pore medium model is utilized to calculate any two described physical quantitys in the described elastic parameter of each group factor of porosity and fluid saturated rock corresponding to the scale value of saturation degree, draw the X plot of any two described physical quantitys in described elastic parameter and indicate the scale value of described factor of porosity and saturation degree, form described rock physics template.

The earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides, preferably, factor of porosity and/or the saturation degree of described calculating measured medium comprise: search the template net lattice point that each described data for projection point is nearest, and then obtain factor of porosity and/or the saturation degree of corresponding described measured medium.Concrete, the value of several prestack elastic parameters obtained to be projected contraposition one by one according to the content of the transverse and longitudinal coordinate in rock physics template and scale, form multiple data for projection point with this, the immediate factor of porosity of template net lattice point corresponding to data for projection point and the scale value of saturation degree are each factor of porosity of measured medium and the value of saturation degree.

The earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides, preferably, described reservoir environment transitivity comprises: temperature, load pressure, pore pressure, layer position, lithology, mineralogical composition, mineral content, porosity ranges, saturation degree scope, pore fluid type.

The earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides, preferably, described petrophysical property comprises: the velocity of longitudinal wave V of ultrasonic measurement _p, shear wave velocity V _sand the corresponding relation between rock porosity, rock saturation degree.I.e. physical quantity velocity of longitudinal wave V in elastic parameter _p, shear wave velocity V _sand the corresponding relation between rock porosity, rock saturation degree.

The earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides, preferably, the described elastic parameter of the fluid saturated rock that the described scale value utilizing described multiple pore medium model to calculate each group factor of porosity and saturation degree is corresponding specifically comprises:

Obtain three parameter of pore structure and three saturation parameters; Preferably, first suppose: 1, blowhole is by isolating hole and intercommunicating pore two parts form; 2, the connective good hole of the preferential filling of the hydro carbons of postchronous migration.Obtain three described parameter of pore structure according to above-mentioned hypothesis to be respectively: aspect ratio α, scale factor x and percent continuity ξ; Three described saturation parameters are respectively: total water saturation S _w, intercommunicating pore water saturation S _cwand isolated hole water saturation S _iw

Pass so between three described saturation parameters is: S _w=S _cwξ+S _iw(1-ξ);

According to formula: φ _iso=φ (1-ξ), calculates the factor of porosity φ of solid matrix _iso, wherein, φ is the total porosity of rock; Described solid matrix comprises: isolated hole and mineral grain;

Again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{\max }=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}})K_{\min }+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_{\mathrm{hc}}{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_w{P}_i^{\′\′}}{(1-{\mathrm{\φ}}_{\mathrm{iso}})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′\′}}\\ {\widehat{\mathrm{\μ}}}_{\mathrm{mat}}=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}}){\mathrm{\μ}}_{\min }}{(1-{\mathrm{\φ}}_{\mathrm{iso}})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′\′}}\end{array}\right.$

Calculate the elastic modulus of solid matrix, wherein: K _minand μ _minbe respectively the average external volume modulus of described mineral grain and average modulus of shearing, K _hcand K _wbe respectively the bulk modulus of hydro carbons and water, P with Q is the geometric factor relevant with pore morpholohy; At this, the value of P with Q is relevant with pore morpholohy, and it can be obtained by prior art, no longer calculates at this.

According to formula: φ _con=φ ξ, calculates the factor of porosity φ of dry skeleton _con; Described dry skeleton comprises: isolated hole, intercommunicating pore and mineral grain;

Again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{\mathrm{dry}}=\frac{(1-\mathrm{\φ})K_{\min }+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_{\mathrm{hc}}{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_w{P}_i^{\′\′}}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′\′}+{\mathrm{\φ}}_{\mathrm{con}}\underset{k=1}{\overset{n}{\mathrm{\Σ}}}{x}_k{\tilde{P}}_k}\\ {\widehat{\mathrm{\μ}}}_{\mathrm{dry}}=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}}){\mathrm{\μ}}_{\min }}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′\′}+{\mathrm{\φ}}_{\mathrm{con}}\underset{k=1}{\overset{n}{\mathrm{\Σ}}}{x}_k{\tilde{Q}}_k}\end{array}\right.$

Calculate the elastic modulus of described dry skeleton; Wherein,

According to the elastic modulus of described solid matrix and the elastic modulus of dry skeleton, Gassmann equation is utilized to calculate the overall elasticity modulus of described fluid saturated rock:

$\left\{\begin{array}{c}{K}_{\mathrm{sat}}={\hat{K}}_{\mathrm{dry}}+\frac{{(1-{\hat{K}}_{\mathrm{dry}}/{\hat{K}}_{\mathrm{mat}})}^2}{{\mathrm{\φ}}_{\mathrm{con}}/K_{\mathrm{fl}}^{*}+(1-{\mathrm{\φ}}_{\mathrm{com}})/{\hat{K}}_{\mathrm{mat}}-{\hat{K}}_{\mathrm{dry}}/{\hat{K}}_{\mathrm{mat}}^2}\\ {\mathrm{\μ}}_{\mathrm{sat}}={\widehat{\mathrm{\μ}}}_{\mathrm{dry}}\end{array}\right.$

Wherein,

$\frac{1}{K_{\mathrm{fl}}^{*}}=\frac{1-S_{\mathrm{cw}}}{K_{\mathrm{hc}}}+\frac{S_{\mathrm{cw}}}{K_w}$

Again according to formula: ρ _sat=(1-φ) ρ _min+ φ (1-S _w) ρ _hc+ φ S _wρ _wcalculate the density of described fluid saturated rock, wherein: ρ _minfor the average density of described mineral grain, ρ _hcand ρ _wbe respectively the density of hydro carbons and water.

Physical quantity described in other of described fluid saturated rock is by the elastic modulus K of described solid matrix _sat, dry skeleton elastic modulus μ _satwith the density p of described fluid saturated rock _satthese three gain of parameters, such as $V_P/V_s=\sqrt{K_{\mathrm{sat}}/{\mathrm{\μ}}_{\mathrm{sat}}+4/3}.$

The earthquake detecting method of hydrocarbon based on multiple pore medium model that the embodiment of the present invention provides, preferably, described aspect ratio α in described parameter of pore structure and described scale factor x obtains according to pore texture characterization technique, namely observes being accurate to the trickle phenomenon of Nano grade; Percent continuity ξ is according to velocity of longitudinal wave V _pand the relation between saturation degree is estimated to obtain, and concrete method of estimation can adopt method of the prior art.

The rock physics template adopting the multiple pore medium model that provides of the above embodiment of the present invention to set up fluid saturated rock embodies Complicated Pore Structures and microfluid uneven distribution to the impact of rock elasticity feature, be more suitable for describing unconventional reservoir particularly compact reservoir, the inversion result therefore based on this model is more reliable.Multiple pore medium model is combined with rock physics mould plate technique, and adopt Template Map method to be converted into by prestack seismic attributes more can the factor of porosity of inverting reservoir oil-gas possibility and saturation parameters, the hydrocarbon occurrence of direct prediction unconventional reservoir, effectively reduces the multi-solution of hydrocarbon indication.Template Map method counting yield is very high, is more suitable for industrial applications.

The embodiment of the present invention also provides a kind of earthquake hydrocarbon indication device based on multiple pore medium model, and wherein as shown in Figure 2, described device comprises:

Data capture unit 201, for obtaining the data set relevant to subsurface reservoir rock, and obtain reservoir information from described data centralization, described reservoir information comprises: reservoir environment transitivity, petrophysical property and micropore structure information; Be can obtain the available data relevant from existing subsurface reservoir rock at this data set, the present invention is the associated reservoir information extracting data centralization, and for setting up the rock physics template of fluid saturated rock.

Template sets up unit 202, for obtaining the model parameter needed for multiple pore medium model according to described reservoir information, sets up the rock physics template of fluid saturated rock based on described multiple pore medium model; Concrete, this step sets up the rock physics template of rock, and the model parameter needed for multiple pore medium model wherein can calculate or estimate to obtain from reservoir information.

Data processing unit 203, extracts angle gathers after carrying out relative amplitude preserved processing to geological data, by described angle gathers by prestack Simultaneous Retrieving, obtains prestack elastic parameter; Concrete, the known geological data obtained is carried out the method that relative amplitude preserved processing can adopt prior art, obtains angle gathers with this, and adopted by angle gathers existing prestack Simultaneous Retrieving technology to obtain prestack elastic parameter.

Projecting cell 204, forms multiple data for projection point for described prestack elastic parameter being projected in described rock physics template, and calculates factor of porosity and/or the saturation degree of measured medium;

By the factor of porosity of described measured medium and/or saturation degree, judging unit 205, for judging that subsurface reservoir hydrocarbon occurrence is to carry out hydrocarbon indication.Concrete, when judging subsurface reservoir hydrocarbon occurrence, factor of porosity or the saturation degree of measured medium can be observed, obtain subsurface reservoir hydrocarbon occurrence thus; Preferably, in order to reduce the simple multi-solution utilizing the factor of porosity of measured medium or saturation degree to judge, factor of porosity can also be utilized to be multiplied by saturation degree and to obtain air content to judge, make the reservoir hydrocarbons distribution results that obtains more accurate.

Earthquake hydrocarbon indication device based on multiple pore medium model provided by the present invention, in one preferred embodiment, the model parameter needed for multiple pore medium model comprises: mineral grain parameter, fluid parameter, factor of porosity parameter, parameter of pore structure, saturation parameters; Wherein, mineral grain parameter comprises: average external volume modulus K _min, average shear modulus μ _min, average density ρ _min; Fluid parameter comprises: hydro carbons bulk modulus K _hc, hydro carbons density p _hc, water volume modulus K _w, water-mass density ρ _w.Above-mentioned parameter can by calculating or estimating to obtain; Concrete, according to the difference of media type, adopt computing method to the parameter conveniently calculated, adopt method of estimation to obtain to calculating the parameter being difficult to obtain, concrete calculating or method of estimation can adopt the state of the art.

In one preferred embodiment, described mineral grain parameter calculates according to Hill model; In another preferred embodiment, described fluid parameter obtains according to Flag program computation.Concrete, Hill model and Flag program are the prior aries that this area is commonly used, and do not repeat them here.

Earthquake hydrocarbon indication device based on multiple pore medium model of the present invention, in one preferred embodiment, is the combination of any two described physical quantitys in elastic parameter according to described rock physics template; Described elastic parameter at least comprises following physical quantity: velocity of longitudinal wave V _p, shear wave velocity V _s, p-wave impedance Z _p, S-wave impedance Z _s, elastic impedance EI, transformed wave elastic impedance PSEI, P-S wave velocity ratio V _p/ V _s, Poisson ratio v, bulk modulus K, modulus of shearing μ, Young modulus E, compressional wave modulus P and Lame's constant λ.Concrete, it is generally acknowledged that fluid saturated rock is isotropic, only have any two physical quantitys to be independently in elastic parameter, such as velocity of longitudinal wave V _pwith shear wave velocity V _s, or bulk modulus K and modulus of shearing μ, and other physical quantitys in elastic parameter just can be transformed by other two physical quantitys arbitrary in this elastic parameter.

A kind of earthquake hydrocarbon indication device based on multiple pore medium model that the embodiment of the present invention provides, preferably, described projecting cell specifically for searching the nearest template net lattice point of each described data for projection point, and then obtains factor of porosity and/or the saturation degree of corresponding described measured medium.Concrete, the value of several prestack elastic parameters obtained to be projected contraposition one by one according to the content of the transverse and longitudinal coordinate in rock physics template and scale, form multiple data for projection point with this, the immediate factor of porosity of template net lattice point corresponding to data for projection point and the scale value of saturation degree are each factor of porosity of measured medium and the value of saturation degree.

A kind of earthquake hydrocarbon indication device based on multiple pore medium model that the embodiment of the present invention provides, preferably, described template sets up unit specifically for the scale value according to porosity ranges and saturation degree range set factor of porosity and saturation degree, described multiple pore medium model is utilized to calculate any two described physical quantitys in the elastic parameter of each group factor of porosity and fluid saturated rock corresponding to the scale value of saturation degree, draw the X plot of any two described physical quantitys in described elastic parameter and indicate the scale value of described factor of porosity and saturation degree, form described rock physics template.

Earthquake hydrocarbon indication device based on multiple pore medium model of the present invention, in one preferred embodiment, described reservoir environment transitivity comprises: temperature, load pressure, pore pressure, layer position, lithology, mineralogical composition, mineral content, porosity ranges, saturation degree scope, pore fluid type.

Earthquake hydrocarbon indication device based on multiple pore medium model of the present invention, in a preferred embodiment, described petrophysical property comprises: the velocity of longitudinal wave V of ultrasonic measurement _p, shear wave velocity V _sand the corresponding relation between rock porosity, rock saturation degree.I.e. physical quantity velocity of longitudinal wave V in elastic parameter _p, shear wave velocity V _sand the corresponding relation between rock porosity, rock saturation degree.

A kind of earthquake hydrocarbon indication device based on multiple pore medium model that the embodiment of the present invention provides, preferably, described template is set up unit and is also comprised: computing unit, for obtaining three parameter of pore structure and three saturation parameters; Three described parameter of pore structure are respectively: aspect ratio α, scale factor x and percent continuity ξ; Three described saturation parameters are respectively: total water saturation S _w, intercommunicating pore water saturation S _cwand isolated hole water saturation S _iw; Pass between three described saturation parameters is: S _w=S _cwξ+S _iw(1-ξ);

Preferably, described computing unit is according to formula: φ _iso=φ (1-ξ), calculates the factor of porosity φ of solid matrix _iso, wherein, φ is the total porosity of rock; Described solid matrix comprises: isolated hole and mineral grain;

Described computing unit is again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{\max }=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}})K_{\min }+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_{\mathrm{hc}}{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_w{P}_i^{\′\′}}{(1-{\mathrm{\φ}}_{\mathrm{iso}})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′\′}}\\ {\widehat{\mathrm{\μ}}}_{\mathrm{mat}}=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}}){\mathrm{\μ}}_{\min }}{(1-{\mathrm{\φ}}_{\mathrm{iso}})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′\′}}\end{array}\right.$

Calculate the elastic modulus of solid matrix, wherein: K _minand μ _minbe respectively the average external volume modulus of described mineral grain and average modulus of shearing, K _hcand K _wbe respectively the bulk modulus of hydro carbons and water, P with Q is the geometric factor relevant with pore morpholohy; At this, the value of P with Q is relevant with pore morpholohy, and it can be obtained by prior art, no longer calculates at this.

Described computing unit is also according to formula: φ _con=φ ξ, calculates the factor of porosity φ of dry skeleton _con; Described dry skeleton comprises: isolated hole, intercommunicating pore and mineral grain;

Described computing unit is according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{\mathrm{dry}}=\frac{(1-\mathrm{\φ})K_{\min }+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_{\mathrm{hc}}{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{K}_w{P}_i^{\′\′}}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{P}_i^{\′\′}+{\mathrm{\φ}}_{\mathrm{con}}\underset{k=1}{\overset{n}{\mathrm{\Σ}}}{x}_k{\tilde{P}}_k}\\ {\widehat{\mathrm{\μ}}}_{\mathrm{dry}}=\frac{(1-{\mathrm{\φ}}_{\mathrm{iso}}){\mathrm{\μ}}_{\min }}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{\mathrm{iso}}(1-S_{\mathrm{iw}})\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{\mathrm{iso}}{S}_{\mathrm{iw}}\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}_i{Q}_i^{\′\′}+{\mathrm{\φ}}_{\mathrm{con}}\underset{k=1}{\overset{n}{\mathrm{\Σ}}}{x}_k{\tilde{Q}}_k}\end{array}\right.$

Calculate the elastic modulus of described dry skeleton;

Described computing unit, also according to the elastic modulus of described solid matrix and the elastic modulus of dry skeleton, utilizes Gassmann equation to calculate the overall elasticity modulus of described fluid saturated rock:

$\left\{\begin{array}{c}{K}_{\mathrm{sat}}={\hat{K}}_{\mathrm{dry}}+\frac{{(1-{\hat{K}}_{\mathrm{dry}}/{\hat{K}}_{\mathrm{mat}})}^2}{{\mathrm{\φ}}_{\mathrm{con}}/K_{\mathrm{fl}}^{*}+(1-{\mathrm{\φ}}_{\mathrm{com}})/{\hat{K}}_{\mathrm{mat}}-{\hat{K}}_{\mathrm{dry}}/{\hat{K}}_{\mathrm{mat}}^2}\\ {\mathrm{\μ}}_{\mathrm{sat}}={\widehat{\mathrm{\μ}}}_{\mathrm{dry}}\end{array}\right.$

Wherein,

$\frac{1}{K_{\mathrm{fl}}^{*}}=\frac{1-S_{\mathrm{cw}}}{K_{\mathrm{hc}}}+\frac{S_{\mathrm{cw}}}{K_w}$

Described computing unit is again according to formula: ρ _sat=(1-φ) ρ _min+ φ (1-S _w) ρ _hc+ φ S _wρ _wcalculate the density of described fluid saturated rock, wherein: ρ _minfor the average density of described mineral grain, ρ _hcand ρ _wbe respectively the density of hydro carbons and water.

Physical quantity described in other of described fluid saturated rock is by the elastic modulus K of described solid matrix _sat, dry skeleton elastic modulus μ _satwith the density p of described fluid saturated rock _satthese three gain of parameters, such as $V_P/V_s=\sqrt{K_{\mathrm{sat}}/{\mathrm{\μ}}_{\mathrm{sat}}+4/3}.$

Earthquake hydrocarbon indication device based on multiple pore medium model of the present invention, in one preferred embodiment, described aspect ratio α in described parameter of pore structure and described scale factor x obtains according to pore texture characterization technique, namely observes being accurate to the trickle phenomenon of Nano grade; Described percent continuity ξ is according to velocity of longitudinal wave V _pand the relation between saturation degree is estimated to obtain, and concrete method of estimation can adopt method of the prior art.

Below with certain basin DAMAGE OF TIGHT SAND GAS RESERVOIRS for example, the implementation procedure of the earthquake detecting method of hydrocarbon based on multiple pore medium model is described.

As shown in Figure 3, the micropore structure nanometer CT scan schematic diagram of certain the tight sand sample provided for the embodiment of the present invention.The skeleton that black part is divided into mineral grain to form, light-colored part represents microscopic void.In this embodiment, the micropore structure of tight sand is very complicated, the hole containing different scale, form and connectedness.From micropore structure image, the morphological feature of each hole can be analyzed, and then obtain corresponding parameter of pore structure, as aspect ratio α, scale factor x and percent continuity ξ.

As shown in Figure 4, be the embodiment of the present invention experimentally data analysis parameter of pore structure.Horizontal ordinate is water saturation, and unit is %; Compressional wave coordinate is velocity of longitudinal wave, and unit is km/s.Diamond spot is sample II-01 (factor of porosity φ=12.1%, permeability κ=0.054mD) measured data of experiment, dotted line is that Patchy saturated model predicts the outcome, the predict of multiple pore medium model when other four lines represent that connective coefficient ξ is 0,0.4,0.8,1 respectively.The permeability of sample II-01 is very low, and the non-constant of its pore connectivity is described, its measured data of experiment also fully confirms this point, and being located substantially on percent continuity is near the curve of 0.

As shown in Figure 5, be p-wave impedance Z that the embodiment of the present invention is set up based on multiple pore medium model _pwith P-S wave velocity ratio V _p/ V _sbetween Z _p-V _p/ V _stemplate.Horizontal ordinate is p-wave impedance, and unit is g/cm ³* km/s; Ordinate is P-S wave velocity ratio.The porosity ranges that template provides is 2% ~ 14%, and saturation degree scope is the corresponding certain rock porosity of circle points in 0% ~ 100%, figure and rock saturation degree.

As shown in Figure 6, be the prestack elastic parameter that angle gathers is obtained by prestack Simultaneous Retrieving that the embodiment of the present invention provides: actual p-wave impedance and actual P-S wave velocity ratio V _p1/ V _s1diagrammatic cross-section.Horizontal ordinate is earthquake Taoist monastic name, and ordinate is whilst on tour, unit millimeter.Utilize log data to set up initial model, then according to AVO inversion theory, prestack Simultaneous Retrieving is carried out to prestack angle gathers, obtain the actual p-wave impedance in prestack elastic parameter and actual P-S wave velocity ratio V _p1/ V _s1, wherein as seen from Figure 6, two wells can be comprised, Y1 well and Y2 well.

As shown in Figure 7, be the embodiment of the present invention, the actual p-wave impedance of prestack elastic parameter and actual P-S wave velocity ratio V _p1/ V _s1at Z _p-V _p/ V _sprojection in template.Horizontal ordinate is p-wave impedance, and unit is g/cm ³* km/s; Ordinate is P-S wave velocity ratio, and the grizzled particle in figure is data for projection point.P-S wave velocity ratio part based on sandstone, substantially cover by template scale mark, clearly can be reflected factor of porosity and the saturation degree of reservoir rock by the scale in rock physics template.On this basis, utilize Template Map method can further inverted porosity and saturation degree.

As shown in Figure 8, be the embodiment of the present invention, according to the earthquake hydrocarbon indication diagrammatic cross-section that Template Map method calculates.Upper figure is factor of porosity section, and middle figure is saturation degree section, and figure below is factor of porosity × saturation degree section.Horizontal ordinate is earthquake Taoist monastic name, and ordinate is whilst on tour, unit millimeter.The confined pressure gas saturation of inverting is all 0 substantially, can not disturb reservoir gas-bearing property analysis, so the gas bearing anomaly in saturation degree section just represents gas reservoir, most high saturation can reach 90%.The section of the product of factor of porosity and saturation degree further can reflect the gas-bearing property of reservoir.Two gassiness Favorable Zones can be identified from these hydrocarbon indication sections, be respectively I district and II district, near Y1 well, remarkable gas bearing anomaly detected in 1040 to 1060ms scope.Y1 well is prolific well, and output is 11.43 ten thousand sides/sky, and this confirms the correctness of hydrocarbon indication result.

Above-described embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only the specific embodiment of the present invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (16)

1. based on an earthquake detecting method of hydrocarbon for multiple pore medium model, it is characterized in that, described method comprises:

Obtain the data set relevant to subsurface reservoir rock, and obtain reservoir information from described data centralization, described reservoir information comprises: reservoir environment transitivity, petrophysical property and micropore structure information;

Obtain the model parameter needed for multiple pore medium model according to described reservoir information, set up the rock physics template of fluid saturated rock based on described multiple pore medium model;

Extract angle gathers after relative amplitude preserved processing is carried out to geological data, by described angle gathers by prestack Simultaneous Retrieving, obtain prestack elastic parameter;

Described prestack elastic parameter is projected in described rock physics template and form multiple data for projection point, and calculate factor of porosity and/or the saturation degree of measured medium;

Judge that subsurface reservoir hydrocarbon occurrence is to carry out hydrocarbon indication by the factor of porosity of described measured medium and/or saturation degree.

2. according to claim 1 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, described data set comprises: geologic report, landwaste record, core data, logging trace.

3. according to claim 1 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, described rock physics template is the combination of any two physical quantitys in elastic parameter; Described elastic parameter at least comprises following physical quantity: velocity of longitudinal wave V _p, shear wave velocity V _s, p-wave impedance Z _p, S-wave impedance Z _s, elastic impedance EI, transformed wave elastic impedance PSEI, P-S wave velocity ratio V _p/ V _s, Poisson ratio ν, bulk modulus K, modulus of shearing μ, Young modulus E, compressional wave modulus P and Lame's constant λ.

4. according to claim 1 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, factor of porosity and/or the saturation degree of described calculating measured medium comprise: search the template net lattice point that each described data for projection point is nearest, and then obtain factor of porosity and/or the saturation degree of corresponding described measured medium.

5. according to claim 3 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, described reservoir environment transitivity comprises: temperature, load pressure, pore pressure, layer position, lithology, mineralogical composition, mineral content, porosity ranges, saturation degree scope, pore fluid type.

6. according to claim 1 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, described petrophysical property comprises: the velocity of longitudinal wave V of ultrasonic measurement _p, shear wave velocity V _sand the corresponding relation between rock porosity, rock saturation degree.

7. according to claim 5 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, the described rock physics template setting up fluid saturated rock based on multiple pore medium model comprises: according to the scale value of described porosity ranges and saturation degree range set factor of porosity and saturation degree, described multiple pore medium model is utilized to calculate any two described physical quantitys in the described elastic parameter of each group factor of porosity and fluid saturated rock corresponding to the scale value of saturation degree, draw the X plot of any two described physical quantitys in described elastic parameter and indicate the scale value of described factor of porosity and saturation degree, form described rock physics template.

8. according to claim 7 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, the described elastic parameter of the fluid saturated rock that the described scale value utilizing described multiple pore medium model to calculate each group factor of porosity and saturation degree is corresponding specifically comprises:

Obtain three parameter of pore structure and three saturation parameters; Three described parameter of pore structure are respectively: aspect ratio α, scale factor x and percent continuity ξ; Three described saturation parameters are respectively: total water saturation S _w, intercommunicating pore water saturation S _cwand isolated hole water saturation S _iw; Pass between three described saturation parameters is: S _w=S _cwξ+S _iw(1-ξ);

According to formula: φ _iso=φ (1-ξ), calculates the factor of porosity φ of solid matrix _iso, wherein, φ is the total porosity of rock; Described solid matrix comprises: isolated hole and mineral grain;

Again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{mat}=\frac{(1-{\mathrm{\φ}}_{iso})K_{\min }+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{K}_{hc}{P}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{K}_w{P}_i^{\′\′}}{(1-{\mathrm{\φ}}_{iso})+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{P}_i^{\′\′}}\\ {\widehat{\mathrm{\μ}}}_{mat}=\frac{(1-{\mathrm{\φ}}_{iso}){\mathrm{\μ}}_{\min }}{(1-{\mathrm{\φ}}_{iso})+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{Q}_i^{\′\′}}\end{array}\right.$

Calculate the elastic modulus of solid matrix, wherein: K _minand μ _minbe respectively the average external volume modulus of described mineral grain and average modulus of shearing, K _hcand K _wbe respectively the bulk modulus of hydro carbons and water, P with Q is the geometric factor relevant with pore morpholohy;

According to formula: φ _con=φ ξ, calculates the factor of porosity φ of dry skeleton _con; Described dry skeleton comprises: isolated hole, intercommunicating pore and mineral grain;

Again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{dry}=\frac{(1-\mathrm{\φ})K_{\min }+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{K}_{hc}{P}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{K}_w{P}_i^{\′\′}}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{P}_i^{\′\′}+{\mathrm{\φ}}_{con}\underset{k=1}{\overset{n}{\Σ}}{x}_k{\tilde{P}}_k}\\ {\widehat{\mathrm{\μ}}}_{dry}=\frac{(1-{\mathrm{\φ}}_{iso}){\mathrm{\μ}}_{\min }}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{Q}_i^{\′\′}+{\mathrm{\φ}}_{con}\underset{k=1}{\overset{n}{\Σ}}{x}_k{\tilde{Q}}_k}\end{array}\right.$

Calculate the elastic modulus of described dry skeleton; Wherein,

According to the elastic modulus of described solid matrix and the elastic modulus of dry skeleton, calculate the overall elasticity modulus of described fluid saturated rock:

$\left\{\begin{array}{c}{K}_{sat}={\hat{K}}_{dry}+\frac{{(1-{\hat{K}}_{dry}/{\hat{K}}_{mat})}^2}{{\mathrm{\φ}}_{con}/K_{fl}^{*}+(1-{\mathrm{\φ}}_{con})/{\hat{K}}_{mat}-{\hat{K}}_{dry}/{\hat{K}}_{mat}^2}\\ {\mathrm{\μ}}_{sat}={\widehat{\mathrm{\μ}}}_{dry}\end{array}\right.$

Wherein,

$\frac{1}{K_{fl}^{*}}=\frac{1-S_{cw}}{K_{hc}}+\frac{S_{cw}}{K_w}$

Again according to formula: ρ _sat=(1-φ) ρ _min+ φ (1-S _w) ρ _hc+ φ S _wρ _wcalculate the density p of described fluid saturated rock _sat, wherein: ρ _minfor the average density of described mineral grain, ρ _hcand ρ _wbe respectively the density of hydro carbons and water.

Physical quantity described in other of described fluid saturated rock is by the elastic modulus K of described solid matrix _sat, dry skeleton elastic modulus μ _satwith the density p of described fluid saturated rock _satthese three gain of parameters.

9. according to claim 8 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, described model parameter comprises: mineral grain parameter, fluid parameter, factor of porosity parameter, parameter of pore structure, saturation parameters; Wherein, mineral grain parameter comprises: average external volume modulus K _min, average shear modulus μ _min, average density ρ _min; Fluid parameter comprises: hydro carbons bulk modulus K _hc, hydro carbons density p _hc, water volume modulus K _w, water-mass density ρ _w.

10. according to claim 9 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, it is characterized in that, described mineral grain parameter calculates according to Hill model.

11. according to claim 9 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, and it is characterized in that, described fluid parameter obtains according to Flag program computation.

12., according to claim 9 based on the earthquake detecting method of hydrocarbon of multiple pore medium model, is characterized in that, the described aspect ratio α in described parameter of pore structure and described scale factor x obtains according to pore texture characterization technique; Described percent continuity ξ is according to velocity of longitudinal wave V _pand the relation between saturation degree is estimated to obtain.

13. 1 kinds of earthquake hydrocarbon indication devices based on multiple pore medium model, it is characterized in that, described device comprises:

Data capture unit, for obtaining the data set relevant to subsurface reservoir rock, and obtain reservoir information from described data centralization, described reservoir information comprises: reservoir environment transitivity, petrophysical property and micropore structure information;

Template sets up unit, for obtaining the model parameter needed for multiple pore medium model according to described reservoir information, sets up the rock physics template of fluid saturated rock based on described multiple pore medium model;

Data processing unit, extracts angle gathers after carrying out relative amplitude preserved processing to geological data, by described angle gathers by prestack Simultaneous Retrieving, obtains prestack elastic parameter;

Projecting cell, forms multiple data for projection point for described prestack elastic parameter being projected in described rock physics template, and calculates factor of porosity and/or the saturation degree of measured medium;

By the factor of porosity of described measured medium and/or saturation degree, judging unit, for judging that subsurface reservoir hydrocarbon occurrence is to carry out hydrocarbon indication.

The 14. earthquake hydrocarbon indication devices based on multiple pore medium model according to claim 13, it is characterized in that, described projecting cell specifically for searching the nearest template net lattice point of each described data for projection point, and then obtains factor of porosity and/or the saturation degree of corresponding described measured medium.

The 15. earthquake hydrocarbon indication devices based on multiple pore medium model according to claim 13, it is characterized in that, described template sets up unit specifically for the scale value according to porosity ranges and saturation degree range set factor of porosity and saturation degree, described multiple pore medium model is utilized to calculate any two physical quantitys in the elastic parameter of each group factor of porosity and fluid saturated rock corresponding to the scale value of saturation degree, draw the X plot of any two described physical quantitys in described elastic parameter and indicate the scale value of described factor of porosity and saturation degree, form described rock physics template.

The 16. earthquake hydrocarbon indication devices based on multiple pore medium model according to claim 15, it is characterized in that, described template is set up unit and is also comprised: computing unit, for obtaining three parameter of pore structure and three saturation parameters; Three described parameter of pore structure are respectively: aspect ratio α, scale factor x and percent continuity ξ; Three described saturation parameters are respectively: total water saturation S _w, intercommunicating pore water saturation S _cwand isolated hole water saturation S _iw; Pass between three described saturation parameters is: S _w=S _cwξ+S _iw(1-ξ);

Described computing unit is according to formula: φ _iso=φ (1-ξ), calculates the factor of porosity φ of solid matrix _iso, wherein, φ is the total porosity of rock; Described solid matrix comprises: isolated hole and mineral grain;

Described computing unit is again according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{mat}=\frac{(1-{\mathrm{\φ}}_{iso})K_{\min }+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{K}_{hc}{P}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{K}_w{P}_i^{\′\′}}{(1-{\mathrm{\φ}}_{iso})+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{P}_i^{\′\′}}\\ {\widehat{\mathrm{\μ}}}_{mat}=\frac{(1-{\mathrm{\φ}}_{iso}){\mathrm{\μ}}_{\min }}{(1-{\mathrm{\φ}}_{iso})+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{Q}_i^{\′\′}}\end{array}\right.$

Calculate the elastic modulus of solid matrix, wherein: K _minand μ _minbe respectively the average external volume modulus of described mineral grain and average modulus of shearing, K _hcand K _wbe respectively the bulk modulus of hydro carbons and water, P with Q is the geometric factor relevant with pore morpholohy;

Described computing unit is also according to formula: φ _con=φ ξ, calculates the factor of porosity φ of dry skeleton _con; Described dry skeleton comprises: isolated hole, intercommunicating pore and mineral grain;

Described computing unit is according to formula:

$\left\{\begin{array}{c}{\hat{K}}_{dry}=\frac{(1-\mathrm{\φ})K_{\min }+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{K}_{hc}{P}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{K}_w{P}_i^{\′\′}}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{P}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{P}_i^{\′\′}+{\mathrm{\φ}}_{con}\underset{k=1}{\overset{n}{\Σ}}{x}_k{\tilde{P}}_k}\\ {\widehat{\mathrm{\μ}}}_{dry}=\frac{(1-{\mathrm{\φ}}_{iso}){\mathrm{\μ}}_{\min }}{(1-\mathrm{\φ})+{\mathrm{\φ}}_{iso}(1-S_{iw})\underset{i=1}{\overset{n}{\Σ}}{x}_i{Q}_i^{\′}+{\mathrm{\φ}}_{iso}{S}_{iw}\underset{i=1}{\overset{n}{\Σ}}{x}_i{Q}_i^{\′\′}+{\mathrm{\φ}}_{con}\underset{k=1}{\overset{n}{\Σ}}{x}_k{\tilde{Q}}_k}\end{array}\right.$

Calculate the elastic modulus of described dry skeleton;

Described computing unit also according to the elastic modulus of described solid matrix and the elastic modulus of dry skeleton, calculates the overall elasticity modulus of described fluid saturated rock:

$\left\{\begin{array}{c}{K}_{sat}={\hat{K}}_{dry}+\frac{{(1-{\hat{K}}_{dry}/{\hat{K}}_{mat})}^2}{{\mathrm{\φ}}_{con}/K_{fl}^{*}+(1-{\mathrm{\φ}}_{con})/{\hat{K}}_{mat}-{\hat{K}}_{dry}/{\hat{K}}_{mat}^2}\\ {\mathrm{\μ}}_{sat}={\widehat{\mathrm{\μ}}}_{dry}\end{array}\right.$

Wherein,

$\frac{1}{K_{fl}^{*}}=\frac{1-S_{cw}}{K_{hc}}+\frac{S_{cw}}{K_w}$

Described computing unit is again according to formula: ρ _sat=(1-φ) ρ _min+ φ (1-S _w) ρ _hc+ φ S _wρ _wcalculate the density p of described fluid saturated rock _sat, wherein: ρ _minfor the average density of described mineral grain, ρ _hcand ρ _wbe respectively the density of hydro carbons and water.

Physical quantity described in other of described fluid saturated rock is by the elastic modulus K of described solid matrix _sat, dry skeleton elastic modulus μ _satwith the density p of described fluid saturated rock _satthese three gain of parameters.