CN106443254B - The electrical parameter test data processing method of porous media containing hydrate - Google Patents
The electrical parameter test data processing method of porous media containing hydrate Download PDFInfo
- Publication number
- CN106443254B CN106443254B CN201610842507.0A CN201610842507A CN106443254B CN 106443254 B CN106443254 B CN 106443254B CN 201610842507 A CN201610842507 A CN 201610842507A CN 106443254 B CN106443254 B CN 106443254B
- Authority
- CN
- China
- Prior art keywords
- hydrate
- porous media
- complex resistivity
- unit
- concentration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
The present invention relates to a kind of electrical parameter test data processing methods of porous media containing hydrate, using the impedance of the porous media containing hydrate as test data, to complex resistivity amplitude, phase angle, four parameters of real and imaginary parts are analyzed and are handled, it is fitted the relationship between each parameter of complex resistivity and hydrate concentration by different modes, can not only the Overall Acquisition electrology characteristic of porous media containing hydrate useful information, and it can realize quantitative assessment to hydrate concentration, and it quickly and efficiently completes from electrical parameter test data to this process of hydrate concentration evaluation result.
Description
Technical field
The invention belongs to oil-gas exploration technical fields, are related to gas hydrates reservoir property Detection Techniques, specifically,
Relate to a kind of electrical parameter test data processing method of porous media containing hydrate.
Background technology
Gas hydrates are a kind of energy resources with great potential, are distributed mainly on land permafrost band and big
In the bottom sediment of shelf margin, have the characteristics that distribution is wide, large storage capacity, energy density are high, cleaning.Gas hydrates store up
Layer is a kind of reservoir of specific type, is needed for the new theory of the Research on Characteristics of gas hydrates reservoir, technology and methods,
More effectively to carry out qualitative recognition to gas hydrates reservoir, to carry out quantitative assessment to hydrate concentration.Natural
In the work of gas hydrate resource investigation, main technological means is carried out by exploitation methods such as geophysics and geochemistry
Exploration, wherein geophysical log is to carry out qualitative and quantitative assessment important means to gas hydrates reservoir, utilizes ground
Ball physical logging technology can detect reservoir physical property under high pressure low temperature environment in situ, and data reliability is higher,
Determining hydrate preservation position, estimation hydrate stock number etc. with irreplaceable role.
Now, the solution of traditional resistivity logging tool and log response is relied primarily on to the evaluation of hydrate concentration
Release result.But, it is contemplated that the particularity and complexity of hydrate reservoir physical property, this electrical parameter of traditional resistivity
The electrology characteristic containing hydrate reservoir cannot be fully portrayed, causes log response information content insufficient, to the solution to log response
The difficulty brought and be difficult to overcome is released, such as uncertainty is high, multi-solution.Complex resistivity includes real part resistivity and imaginary part resistance
Rate, can not only portray to Comprehensive the electrical properties of measured medium (contain hydrate reservoir), and Primary Study also table
The complex resistivity and hydrate concentration of the bright porous media containing hydrate are closely related.It can be seen that complex resistivity is capable of providing
Largely for the electrology characteristic information of analysis and utilization, but there are following urgent problems：(1) how to complex resistivity
Test data carries out effectively analyzing processing, to excavate and realize determining for hydrate reservoir using these electrology characteristic information
Property and quantitative assessment；(2) analysis and the Processing Algorithm of test data how to be realized, to quickly and efficiently complete from electrical parameter
Test data is to this process of hydrate concentration evaluation result.The open report to solve the above problems is not yet found at present.
Invention content
The present invention is directed to the above-mentioned deficiency of the prior art, provides a kind of electrical parameter test data of porous media containing hydrate
Processing method, this method effectively can be analyzed and be handled to complex resistivity measurement data, be portrayed comprehensively to excavate
The useful information of the electrology characteristic of porous media containing hydrate realizes the quantitative assessment to hydrate concentration, and quickly and efficiently
It completes from electrical parameter test data to this process of hydrate concentration evaluation result.
In order to achieve the above object, the present invention provides a kind of electrical parameter test data of porous media containing hydrate processing
Method contains following steps：
The temperature of porous media and the pressure of gas in reaction kettle are obtained during decomposition of hydrate, and according to temperature and pressure
The hydrate concentration in different decomposition stage is calculated by the calculation formula of hydrate concentration for power；
The resistance of the porous media containing hydrate under different test frequency pumping signals effect during acquisition decomposition of hydrate
It is anti-, and sent a telegram in reply accordingly according to the geometric dimension of electric sensor and distance and the impedance computation of the porous media containing hydrate
Resistance rate, and further obtain real part, imaginary part, amplitude and the phase angle of complex resistivity；
The impedance spectrum Nyquist figure that the porous media containing hydrate is drawn according to impedance, by impedance spectrum Nyquist figures obtain with
The basic structure of its corresponding equivalent-circuit model, the influence according to hydrate concentration to equivalent circuit model parameter value
Rule, the relationship being fitted between equivalent circuit model parameter and hydrate concentration, establishes MODELING FOR COMPLEX RESISTIVITY MODEL, the complex resistance
The output parameter of rate model is hydrate concentration, and input parameter is complex resistivity and test frequency；
According to the MODELING FOR COMPLEX RESISTIVITY MODEL, the electric sensor geometric dimension actually measured and distance and test frequency, profit
With solution by iterative method hydrate concentration.
As the preferred design of the above-mentioned processing method of the present invention, the hydrate concentration is carried out according to temperature and pressure
The calculation formula of calculating is expressed as：
In formula, ShFor the hydrate concentration of the porous media containing hydrate；MhFor hydrate molal weight, unit：kg/
mol；ρhFor the density of hydrate, unit：kg/m3；VPFor the volume of porous media hole, unit：m3；VGTo react gas reactor
Volume, unit：m3；T1(i.e. Hydrate in Porous Medium does not have started generation or for porous media when to be free of hydrate
All decompose finish) temperature, unit：K；T2For during decomposition of hydrate when containing hydrate porous media temperature, it is single
Position：K；P1To be free of hydrate in porous media, (i.e. Hydrate in Porous Medium does not have started generation or all decomposes
Finish) when react gas reactor pressure, unit：Pa；P2It is single to react the pressure of gas reactor during decomposition of hydrate
Position：Pa；Z1To be free of hydrate in porous media, (i.e. Hydrate in Porous Medium does not have started generation or all decomposes
Finish) when react gas reactor compressibility factor；Z2To react the compressibility factor of gas reactor during decomposition of hydrate；R is
Mol gas constant, unit：J/(mol·K).
As the preferred design of the above-mentioned processing method of the present invention, the impedance spectrum Nyquist figures are approximately one and coordinate
Axis is in the straight line at 45 degree of angles, and the basic structure of equivalent-circuit model is by with ω RWCWOne resistance R of=1 relationshipWWith one
Capacitance CWIt is composed in series, resistance RWWith capacitance CWRespectively with ω-1/2It is in a linear relationship, it is expressed as：
RW=a ω-1/2+b (2)
CW=c ω-1/2+d (3)
In formula, a, b, c, d are equivalent circuit model parameter；
The expression formula of the equivalent-circuit model is：
In formula, ρ is complex resistivity；S is the area of electric sensor electrode slice, unit：m2；L is electric sensor electrode
The distance between piece, unit：m；ω is angular frequency, unit：The π of rad/s, ω=2 f, f are test frequency, unit：Hz；
It is further represented as：
According to hydrate concentration to the affecting laws of equivalent circuit model parameter value, it is fitted equivalent circuit model parameter
A, the relationship between b, c, d and hydrate concentration, equivalent circuit model parameter a, b, c, d are respectively between hydrate concentration
Relationship be expressed as：
B=40.401Sh+30.732 R2=0.9882 (7)
D=0.00002 (9)
In formula, R2For degree of fitting, i.e. matched curve and for the degree of correlation between the data of fitting；
And then the relationship between the complex resistivity and hydrate concentration of the porous media containing hydrate is obtained, establish complex resistance
Rate model, the MODELING FOR COMPLEX RESISTIVITY MODEL are expressed as：
In formula, ShFor the hydrate concentration of the porous media containing hydrate.
In order to achieve the above object, at the present invention also provides a kind of electrical parameter test data of porous media containing hydrate
Reason method, contains following steps：
The temperature of porous media and the pressure of gas in reaction kettle are obtained during decomposition of hydrate, and according to temperature and pressure
The hydrate concentration in different decomposition stage is calculated by the calculation formula of hydrate concentration for power；
The resistance of the porous media containing hydrate under different test frequency pumping signals effect during acquisition decomposition of hydrate
It is anti-, and sent a telegram in reply accordingly according to the geometric dimension of electric sensor and distance and the impedance computation of the porous media containing hydrate
Resistance rate, and further obtain real part, imaginary part, amplitude and the phase angle of complex resistivity；
Respectively under linear coordinate system and non-linear scale system, using test frequency as reference axis horizontal axis, respectively with complex resistance
Rate real part, imaginary part, amplitude and phase angle are the reference axis longitudinal axis, draw the Dispersion of complex resistivity real part, imaginary part, amplitude and phase angle
Curve acts on lower complex resistivity real part, imaginary part, amplitude and phase according to Dispersion curve acquisition difference test frequency pumping signal
The characteristic parameter of the Dispersion curve of four, angle parameter；
Under linear coordinate system, using hydrate concentration as reference axis horizontal axis, with complex resistivity real part, imaginary part, amplitude and
The characteristic parameter of phase angle Dispersion curve is the reference axis longitudinal axis, draws the frequency dispersion of complex resistivity real part, imaginary part, amplitude and phase angle
Relational graph between characteristic characteristic parameter and hydrate concentration, according to the characteristic parameter of Dispersion curve and hydration
Relationship between object saturation degree establishes hydrate concentration computation model, for determining hydrate concentration.
In order to achieve the above object, at invention further provides a kind of electrical parameter test data of porous media containing hydrate
Reason method, contains following steps：
The temperature of porous media and the pressure of gas in reaction kettle are obtained during decomposition of hydrate, and according to temperature and pressure
The hydrate concentration in different decomposition stage is calculated by the calculation formula of hydrate concentration for power；
The resistance of the porous media containing hydrate under different test frequency pumping signals effect during acquisition decomposition of hydrate
It is anti-, and sent a telegram in reply accordingly according to the geometric dimension of electric sensor and distance and the impedance computation of the porous media containing hydrate
Resistance rate, and further obtain real part, imaginary part, amplitude and the phase angle of complex resistivity；
Under linear coordinate system, using hydrate concentration as reference axis horizontal axis, with complex resistivity real part, imaginary part, amplitude and
The frequency dispersion degree of phase angle is the reference axis longitudinal axis, and the frequency dispersion degree and hydrate for drawing complex resistivity real part, imaginary part, amplitude and phase angle are saturated
Relational graph between degree, according to the pass between the frequency dispersion degree and hydrate concentration of complex resistivity real part, imaginary part, amplitude and phase angle
System, establishes hydrate concentration computation model, for determining hydrate concentration.
The present invention also provides a kind of software realizations of the electrical parameter of porous media containing hydrate test data processing method
Method contains following steps：
Using Matlab as development platform, exploitation can carry out MODELING FOR COMPLEX RESISTIVITY MODEL the code of numerical solution, utilize
MODELING FOR COMPLEX RESISTIVITY MODEL solution code is converted into the C++ codes of standard by the coder tools of Matlab；
Using Visual Studio as development platform, by creating engineering, writing source file, setting engineering attribute, compiling four
A step is handled above-mentioned through the transformed C++ codes for carrying out numerical solution to MODELING FOR COMPLEX RESISTIVITY MODEL of Matlab,
And then generate dynamic link library (DLL, Dynamic Link Library) file；
Using LabVIEW as development platform, the subprogram (i.e. sub- VI) for realizing following functions is developed：The reality of impedance data
When obtain calculated with pretreatment, complex resistivity, real-time display and the preservation of impedance data and complex resistivity, wherein pretreatment master
To include filtering and being averaged, it includes real part, imaginary part, amplitude and the phase angle for calculating complex resistivity, impedance number that complex resistivity, which calculates,
When according to being preserved in real time with complex resistivity data, hard disk is stored in text or binary format；Utilize the library function in LabVIEW
Node directly invokes the above-mentioned dll file for carrying out numerical solution to MODELING FOR COMPLEX RESISTIVITY MODEL, and MODELING FOR COMPLEX RESISTIVITY MODEL is exported
Hydrate concentration data carry out real-time display and preservation；
Using the program issuing function of LabVIEW, executable file is first generated, then installation journey is generated by executable file
Sequence, and then install and be applied to computer.
Compared with prior art, the beneficial effects of the present invention are：
(1) present invention is analyzed and is handled by the complex resistivity measurement data to the porous media containing hydrate, can
The electrology characteristic of Overall Acquisition porous media containing hydrate obtains the relationship between complex resistivity and hydrate concentration, establishes
MODELING FOR COMPLEX RESISTIVITY MODEL realizes the quantitative assessment to hydrate concentration, and can quickly and efficiently complete to test from electrical parameter
Data are to this process of hydrate concentration evaluation result.
(2) processing method through the invention is capable of the electrology characteristic of Overall Acquisition porous media containing hydrate, obtains electricity
Learn the quantitative relationship between characteristic and hydrate concentration, the i.e. feature and saturation of the Dispersion curve of each parameter of complex resistivity
It is full to establish corresponding hydrate for each parameter frequency dispersion degree of relationship, complex resistivity between degree and the relationship between hydrate concentration
With degree computation model, the quantitative assessment to hydrate concentration is realized, and can quickly and efficiently complete to test from electrical parameter
Data are to this process of hydrate concentration evaluation result.
(3) MODELING FOR COMPLEX RESISTIVITY MODEL that the present invention is established takes full advantage of the electrical properties of the porous media containing hydrate, tool
It is that resistivity and dielectric constant are utilized simultaneously for body, existing resistance characteristic has capacitance characteristic again, to hydrate concentration
It is calculated so that definitely, the accuracy higher of model, the scope of application is wider for model physical significance.
(4) the present invention also provides a kind of software realizations of the electrical parameter of porous media containing hydrate test data processing method
Method, using the method for cross-platform hybrid programming so that model solution is more quickly and accurate.
Description of the drawings
Fig. 1 is that electrical parameter of the embodiment of the present invention tests system schematic.
Fig. 2 is the flow of the electrical parameter test data processing method of porous media containing hydrate in the embodiment of the present invention one
Figure.
Fig. 3 is the electrical parameter test data processing method software realization of porous media containing hydrate in the embodiment of the present invention one
Method flow diagram.
Fig. 4 is the impedance spectrum Nyquist figures of the porous media containing hydrate in the embodiment of the present invention one.
Fig. 5 is the equivalent-circuit model figure of the porous media containing hydrate in the embodiment of the present invention one.
Fig. 6 is resistance R in the equivalent-circuit model of the porous media containing hydrate in the embodiment of the present invention oneWWith capacitance CWPoint
Not and ω-1/2Between relational graph.
Fig. 7 is the equivalent circuit model parameter a and hydrate of the porous media containing hydrate in the specific embodiment of the invention one
Relational graph between saturation degree.
Fig. 8 is the equivalent circuit model parameter b and hydrate of the porous media containing hydrate in the specific embodiment of the invention one
Relational graph between saturation degree.
Fig. 9 is the equivalent circuit model parameter c and hydrate of the porous media containing hydrate in the specific embodiment of the invention one
Relational graph between saturation degree.
Figure 10 is the equivalent circuit model parameter d and hydrate of the porous media containing hydrate in the specific embodiment of the invention one
Relational graph between saturation degree.
Figure 11 is the Dispersion of each parameter of complex resistivity of the porous media containing hydrate in the specific embodiment of the invention two
Figure.
Figure 12 is to be included in 0.1Hz-10Hz test frequency ranges under log-log coordinate system in the specific embodiment of the invention two
The Dispersion figure of each parameter of complex resistivity of hydrate porous media.
Figure 13 is the Dispersion of each parameter of complex resistivity of the porous media containing hydrate in the specific embodiment of the invention two
Relational graph between the absolute value and hydrate concentration of slope of curve k.
The frequency dispersion degree of complex resistivity amplitude when Figure 14 is methane hydrate saturation degree difference in the specific embodiment of the invention three
Variation diagram.
The frequency dispersion degree of complex resistivity real part when Figure 15 is methane hydrate saturation degree difference in the specific embodiment of the invention three
Variation diagram.
In figure, 1, reaction kettle, 2, Low Constant Temperature Water Baths, 3, gas bomb, 4, booster pump, 5, electric sensor outlet hole, 6,
Impedance analyzer, 7, temperature sensor outlet hole, 8, pressure sensor outlet hole, 9, data collector, 10, computer.
Specific implementation mode
In the following, the present invention is specifically described by illustrative embodiment.It should be appreciated, however, that not into one
In the case of step narration, element, structure and features in an embodiment can also be advantageously incorporated into other embodiment
In.
Below by being illustrated for containing the processing of methane hydrate porous media electrical parameter test data.Pass through electricity
Parameter test system tests pressure, temperature and the impedance in the reaction kettle in methane hydrate decomposable process, obtains
Corresponding pressure, temperature and testing impedance data.See that Fig. 1, the structure composition of the electrical parameter test system include：Reaction
Kettle, for generating and methane decomposition hydrate；The Low Constant Temperature Water Baths being connect with reaction kettle provide constant low temperature for reaction kettle
Environment；Gas bomb provides required gas to generate methane hydrate；Booster pump is connected between gas bomb and reaction kettle, right
Gas in gas bomb carries out pressurized treatment, is then charged in reaction kettle, can either provide the gas generated needed for methane hydrate
Body is capable of providing hyperbaric environment again；Electric sensor in reaction kettle, for measuring the porous media containing methane hydrate
Impedance；The impedance analyzer being connect with electric sensor keeps electric sensor complete for providing pumping signal for electric sensor
At the function of measuring the impedance of porous media containing methane hydrate, meanwhile, the impedance of measurement is sent to computer by impedance analyzer；
Temperature sensor in reaction kettle, for measuring the temperature in decomposition of hydrate or generating process；In reaction kettle
Pressure sensor, for measuring the pressure in decomposition of hydrate or generating process；Data collector, for acquiring the temperature measured
And pressure, and the temperature and pressure of acquisition is sent to computer.
Referring to Fig. 1, the reaction kettle is cylinder, and lower space holds the aqueous sodium chloride that mass fraction is 3.5%
The natural sea sand of liquid saturation, upper space are filled with high pressure gas, and the upper surface of natural sea sand is contacted with high pressure gas and upper surface
It is horizontal.
Electric sensor described above includes two metal electrode films for being relatively parallel placement, and two sheet metals are vertically extra large
Sand surface is fully inserted into sea sand, and the distance of two electrode slices is 7cm, and electrode slice is square, length of side 5cm.
The low temperature environment refers to temperature and is not higher than 5 DEG C.
Embodiment one：Referring to Fig. 2, a kind of electrical parameter test data processing method of porous media containing hydrate is described to contain
Hydrate porous media is contained by what the sodium-chloride water solution that methane hydrate, natural sea sand and mass fraction are 3.5% formed
Methane hydrate porous media, it is assumed that the chemical formula of methane hydrate is CH4·nH2O, n are hydrate index, are to n values
5.89 carrying out correlation computations.For ease of description, porous media containing hydrate, the processing method is claimed to contain following step below
Suddenly：
Step 1：Obtain during decomposition of hydrate the temperature of porous media and the pressure of gas in reaction kettle, and according to
The hydrate concentration in different decomposition stage is calculated by the calculation formula of hydrate concentration for temperature and pressure；The water
The calculation formula for closing object saturation degree is expressed as：
In formula, ShFor the hydrate concentration of the porous media containing hydrate；MhFor hydrate molal weight, unit：kg/
mol；ρhFor the density of hydrate, unit：kg/m3；VPFor the volume of porous media hole, unit：m3；VGTo react gas reactor
Volume, unit：m3；T1(i.e. Hydrate in Porous Medium does not have started generation or for porous media when to be free of hydrate
All decompose finish) temperature, unit：K；T2For during decomposition of hydrate when containing hydrate porous media temperature, it is single
Position：K；P1To be free of hydrate in porous media, (i.e. Hydrate in Porous Medium does not have started generation or all decomposes
Finish) when react gas reactor pressure, unit：Pa；P2It is single to react the pressure of gas reactor during decomposition of hydrate
Position：Pa；Z1To be free of hydrate in porous media, (i.e. Hydrate in Porous Medium does not have started generation or all decomposes
Finish) when react gas reactor compressibility factor；Z2To react the compressibility factor of gas reactor during decomposition of hydrate；R is
Mol gas constant, unit：J/(mol·K).
Step 2：Obtain lower porous Jie containing hydrate of different test frequency pumping signals effect during decomposition of hydrate
The impedance of matter, and obtained accordingly according to the impedance computation of the geometric dimension of electric sensor and distance and the porous media containing hydrate
Complex resistivity, and further obtain complex resistivity real part, imaginary part, amplitude and phase angle；Wherein test frequency is in 0.1Hz-
Value within the scope of 1MHz.
The amplitude of complex resistivity is calculated according to formula (2), the expression formula of formula (2) is：
In formula, | ρ | it is complex resistivity amplitude, unit：Ω·m；| Z | it is impedance magnitude, unit：Ω；S is the face of electrode slice
Product, unit：m2；L is the distance between electrode slice, unit：m.
In the present embodiment, the distance between electrode slice is 7cm, and the area of electrode slice is 25cm2。
The real and imaginary parts of complex resistivity are calculated according to formula (3), the expression formula of formula (3) is：
In formula, ρ ' is the real part of complex resistivity, unit：Ω·m；ρ " is the imaginary part of complex resistivity, unit：Ω·m；θ is
Phase angle.
Step 3：Since impedance real part Z ' and imaginary impedance Z " the impedance spectrum Nyquist figures constituted are common impedances
It is basic can to obtain corresponding equivalent-circuit model according to the curve in impedance spectrum Nyquist figures for modal data representation
Structure.The impedance spectrum Nyquist figures that the porous media containing hydrate is drawn according to the impedance that step 2 measures, by impedance spectrum
Nyquist figures obtain the basic structure of equivalent-circuit model corresponding thereto, according to hydrate concentration to equivalent circuit mould
The affecting laws of shape parameter value, the relationship being fitted between equivalent circuit model parameter and hydrate concentration, establish complex resistance
The output parameter of rate model, the MODELING FOR COMPLEX RESISTIVITY MODEL is hydrate concentration, and input parameter is complex resistivity and test frequency.
Referring to Fig. 4, the impedance spectrum Nyquist figures are a straight line with reference axis in 45 degree of angles, are methane hydrates
Impedance spectrum Nyquist figures when saturation degree is 30.51%, referring to Fig. 5, the basic structure of equivalent-circuit model is by with ω
RWCWOne resistance R of=1 relationshipWWith a capacitance CWIt is composed in series, and resistance RWWith capacitance CWWith ω-1/2It is linear to close
System, resistance R is obtained by the real and imaginary parts of complex resistivityWWith capacitance CW.Referring to Fig. 6, resistance RWWith capacitance CWWith ω-1/2Line
Sexual intercourse is expressed as：
RW=a ω-1/2+b (4)
CW=c ω-1/2+d (5)
In formula, a, b, c, d are equivalent circuit model parameter；
The expression formula of the equivalent-circuit model is：
In formula, ρ is complex resistivity；ω is angular frequency, unit：The π of rad/s, ω=2 f, f are test frequency, unit：Hz；
It is further represented as：
According to hydrate concentration to the affecting laws of equivalent circuit model parameter value, it is fitted equivalent circuit model parameter
A, the relationship between b, c, d and hydrate concentration, referring to Fig. 7-10, wherein Fig. 7 is between hydrate concentration and parameter a
Relational graph, relational graphs of the Fig. 8 between hydrate concentration and parameter b, Fig. 9 is between hydrate concentration and parameter c
Relational graph can be obtained by Figure 10, and d approximations are taken as 0.00002.Equivalent circuit model parameter a, b, c, d respectively with hydrate concentration
Between relationship be expressed as：
B=40.401Sh+30.732R2=0.9882 (9)
D=0.00002 (11)
In formula, R2For degree of fitting, i.e. matched curve and for the degree of correlation between the data of fitting, value is closer to 1
Indicate that fitting degree is higher, fitting effect is better；
And then the relationship between the complex resistivity and hydrate concentration of the porous media containing hydrate is obtained, establish complex resistance
Rate model, the MODELING FOR COMPLEX RESISTIVITY MODEL are expressed as：
In formula, ShFor the hydrate concentration of the porous media containing hydrate.
Step 4：According to the MODELING FOR COMPLEX RESISTIVITY MODEL, the electric sensor geometric dimension actually measured and distance and test
Frequency utilizes solution by iterative method hydrate concentration.
Referring to Fig. 3, the step of above-mentioned data processing method software implementation method of the invention, is：
Step 1：It is development platform with Matlab 2012, exploitation can carry out numerical solution to above-mentioned MODELING FOR COMPLEX RESISTIVITY MODEL
Code, model solution code is converted into the C++ codes of standard using the coder tools of Matlab.In view of formula (12)
The complexity of shown MODELING FOR COMPLEX RESISTIVITY MODEL, i.e., non-linear complex number equation, in the present embodiment using Fixed point method or it is general repeatedly
It for method or Newton iteration method, and works out Matlab codes and carries out algorithm realization, code saves as .m files, by required parameter
Type inputs, such as parameter type double indicates that double precision, selection need to generate the standard code of C Plus Plus form, to
It completes Matlab code conversions to be transplantable standard C++ codes.
Step 2：It is development platform with Visual Studio 2010, by creating engineering, writing source file, setting work
Four steps such as journey attribute, compiling are by above-mentioned through the transformed C+ for carrying out numerical solution to MODELING FOR COMPLEX RESISTIVITY MODEL of Matlab
+ code is handled, and then generates dynamic link library (DLL, Dynamic LinkLibrary) file, and file name suffix is
.dll。
Step 3：It is development platform with LabVIEW 2012, develops the subprogram (i.e. sub- VI) for realizing following functions：
The real-time acquisition of impedance data and pretreatment, complex resistivity calculating, real-time display and the preservation of impedance data and complex resistivity,
In, pretreatment includes mainly filtering and is averaged, complex resistivity calculate include the calculating real part of complex resistivity, imaginary part, amplitude and
Phase angle is stored in hard disk when impedance data preserves in real time with complex resistivity data with text or binary format；
It is directly invoked using the library function node in LabVIEW above-mentioned for carrying out numerical solution to MODELING FOR COMPLEX RESISTIVITY MODEL
Dll file, the Selection and call library function node in LabVIEW, and configure required defeated in hydrate concentration calculating process
Enter parameter, i.e., the data type of electric sensor geometric dimension and distance, test frequency, complex resistivity and each parameter is (such as
Double, i.e. double precision), and the hydrate concentration data of MODELING FOR COMPLEX RESISTIVITY MODEL output are subjected to real-time display and preservation.
Step 4：Using the program issuing function of LabVIEW, executable file is first generated, then give birth to by executable file
At installation procedure, installs and be applied to computer, without installing LabVIEW and Matlab.
The operating procedure for generating executable file is as follows：
An engineering project is created first in LabVIEW, addition in the project generates the material needed for executable file,
These materials include mainly in the present embodiment：The LabVIEW programs of data acquisition and procession software are calculated in hydrate concentration
When need the dll file that calls and the library file msvcr120.dll of 2010 compilers of Visual Studio acquiescence；Its
In, the LabVIEW programs of the data acquisition and procession software be described in above-mentioned steps three subprogram (i.e. sub- VI,
It is exactly .vi files)
The click right in Program Generating specification creates an application program .exe, inputs and make by oneself in first item information
The .vi files of data acquisition and procession software are placed in source file and start inside vi, by other by the application name of justice
Dll file be placed on the position for including always；
Linguistic property when selecting the icon of application program, running clicks generation and can be obtained required .exe files later,
That is executable file.
The operating procedure for generating installation procedure is as follows：
Continue to create in above-mentioned engineering first, right button option program generates specification, the source document in installation procedure attribute
Above-mentioned generated .exe files are placed into the target view on the right by part part；By LabVIEW in additional setup procedures
Runtime engine cover in installation procedure, Setup Type selection all；It clicks to generate later and can be obtained required installation journey
Sequence file includes mainly setup.exe files, setup.ini files, bin file folder, license files in file
With supportfiles files.
Embodiment two：A kind of electrical parameter test data processing method of porous media containing hydrate is described more containing hydrate
Hole medium is hydrated by what the sodium-chloride water solution that methane hydrate, natural sea sand and mass fraction are 3.5% formed containing methane
Object porous media, it is assumed that the chemical formula of methane hydrate is CH4·nH2O, n are hydrate index, are 5.89 progress phases to n values
It closes and calculates.For ease of description, porous media containing hydrate, the processing method is claimed to contain following steps below：
Step 1：Obtain during decomposition of hydrate the temperature of porous media and the pressure of gas in reaction kettle, and according to
The hydrate concentration in different decomposition stage is calculated by the calculation formula of hydrate concentration for temperature and pressure；The water
The calculation formula for closing object saturation degree is expressed as：
In formula, ShFor the hydrate concentration of the porous media containing hydrate；MhFor hydrate molal weight, unit：kg/
mol；ρhFor the density of hydrate, unit：kg/m3；VPFor the volume of porous media hole, unit：m3；VGTo react gas reactor
Volume, unit：m3；T1(i.e. Hydrate in Porous Medium does not have started generation or for porous media when to be free of hydrate
All decompose finish) temperature, unit：K；T2For during decomposition of hydrate when containing hydrate porous media temperature, it is single
Position：K；P1To be free of hydrate in porous media, (i.e. Hydrate in Porous Medium does not have started generation or all decomposes
Finish) when react gas reactor pressure, unit：Pa；P2It is single to react the pressure of gas reactor during decomposition of hydrate
Position：Pa；Z1To be free of hydrate in porous media, (i.e. Hydrate in Porous Medium does not have started generation or all decomposes
Finish) when react gas reactor compressibility factor；Z2To react the compressibility factor of gas reactor during decomposition of hydrate；R is
Mol gas constant, unit：J/(mol·K).
Step 2：Obtain lower porous Jie containing hydrate of different test frequency pumping signals effect during decomposition of hydrate
The impedance of matter, and obtained accordingly according to the impedance computation of the geometric dimension of electric sensor and distance and the porous media containing hydrate
Complex resistivity, and further obtain complex resistivity real part, imaginary part, amplitude and phase angle；Wherein test frequency is in 0.1Hz-
Value within the scope of 1MHz.
The amplitude of complex resistivity is calculated according to formula (2), the expression formula of formula (2) is：
In formula, | ρ | it is complex resistivity amplitude, unit：Ω·m；| Z | it is impedance magnitude, unit：Ω；S is the face of electrode slice
Product, unit：m2；L is the distance between electrode slice, unit：m.
In the present embodiment, the distance between electrode slice is 7cm, and the area of electrode slice is 25cm2。
The real and imaginary parts of complex resistivity are calculated according to formula (3), the expression formula of formula (3) is：
In formula, ρ ' is the real part of complex resistivity, unit：Ω·m；ρ " is the imaginary part of complex resistivity, unit：Ω·m；θ is
Phase angle.
Step 3：Respectively under linear coordinate system and non-linear scale system, using test frequency as reference axis horizontal axis, respectively with
Complex resistivity real part, imaginary part, four parameters of amplitude and phase angle be the reference axis longitudinal axis, draw complex resistivity real part, imaginary part absolute value,
The Dispersion curve of amplitude and phase angle absolute value, under being acted on according to Dispersion curve acquisition difference test frequency pumping signal
Complex resistivity real part, four imaginary part, amplitude and phase angle parameters Dispersion curve characteristic parameter；
Referring to Figure 11, Figure 11 give methane hydrate saturation degree be 30.51% when, the amplitude of complex resistivity, phase angle are exhausted
To the Dispersion curve of value, real part and imaginary part absolute value within the scope of test frequency 0.1Hz-1MHz.Song in analysis chart 11
Line can obtain：With the increase of test frequency, the amplitude of complex resistivity is decreased to 1.69 Ω m, real part by 29.04 Ω m
It is decreased to 1.68 Ω m by 20.90 Ω m, imaginary part absolute value is decreased to increase to 0.11 again after 0 Ω m by 20.16 Ω m
Ω m, influence of the test frequency to complex resistivity amplitude, real part and imaginary part are gradually reduced.In test frequency 0.1Hz-100Hz
The amplitude of complex resistivity, real part and imaginary part dispersion phenomenon become apparent in range.With the increase of test frequency, complex resistivity
Phase angle absolute value first reduces to be increased afterwards, and 4.15 ° are increased to again after being decreased to 0.01 ° by 43.96 °.Test frequency 10Hz~
100Hz range internal phase angle absolute values reduce speed；After the process slowly reduced by one section, when test frequency is close
When 200kHz, the phase angle absolute value of complex resistivity has the tendency that significantly increasing；The imaginary part absolute value and phase angle of complex resistivity are absolute
The extreme point for being worth Dispersion curve is 200kHz, it is believed that extreme point is the boundary for swashing electrical effect and electromagnetic induction effect
Point, (i.e. frequency is less than extreme value dot frequency 200kHz) effect of induced polarization is major influence factors before extreme point, extreme point it
(i.e. frequency is higher than extreme value dot frequency 200kHz) electromagnetic induction effect is gradually reinforced afterwards.
Referring to Figure 12, Figure 12 give methane hydrate saturation degree be 30.51% when, 0.1Hz~10Hz test frequency models
Enclose complex resistivity real part under interior log-log coordinate system, four imaginary part, amplitude and phase angle parameters Dispersion curve.Complex resistivity
Amplitude, real and imaginary parts absolute value and test frequency at linear approximate relationship, be expressed as：
Lgy=klgx+b (13)
In formula, y is amplitude, real part or the imaginary part absolute value of complex resistivity；X is test frequency；K is Dispersion curve
Slope indicates that the significance degree of each parameter frequency dispersion, the absolute value of the more notable then k of frequency dispersion are bigger；B is when test frequency x is 1Hz
When, the amplitude of complex resistivity, the logarithm of real part or imaginary part absolute value.
Linear fit is carried out based on the principle of least square, obtains amplitude, real part, imaginary part absolute value and the frequency of complex resistivity
Between the slope of fitting a straight line be respectively k1=-0.4606 (R2=0.9997), k2=-0.4276 (R2=0.9983), k3=-
0.5123(R2=0.9992).
Step 4：Under linear coordinate system, using hydrate concentration as reference axis horizontal axis, with complex resistivity real part, imaginary part
Absolute value, amplitude Dispersion slope of a curve be the reference axis longitudinal axis, draw complex resistivity real part, imaginary part absolute value, amplitude
Dispersion slope of a curve absolute value and hydrate concentration between relational graph, it is exhausted according to Dispersion slope of a curve
Relationship between value and hydrate concentration establishes hydrate concentration computation model, for determining hydrate concentration.
The amplitude of complex resistivity, the frequency dispersion of real part and imaginary part absolute value are sought under the conditions of different methane hydrate saturation degrees
Characteristic curve k, amplitude, real part and the imaginary part that hydrate concentration and complex resistivity are given referring to Figure 13, Figure 13 are absolute
It is worth the relationship of the slope k absolute value under Dispersion curve log-log coordinate system.As seen from Figure 13, full with methane hydrate
Increase with degree, the absolute value reduction of the amplitude, real part and imaginary part absolute value Dispersion slope of a curve k of complex resistivity, telegram in reply
The frequency dispersion degree of the amplitude of resistance rate, real part and imaginary part absolute value gradually weakens.Under same saturation degree, complex resistivity imaginary part is absolute
The frequency dispersion degree of value is most strong, and complex resistivity amplitude is taken second place, and complex resistivity real part is most weak.It is bent to hydrate concentration and Dispersion
Line slope k absolute values are fitted, and obtain the relationship between hydrate concentration and slope k absolute value, establish hydrate saturation
Computation model is spent, for determining hydrate concentration.
Embodiment three：A kind of electrical parameter test data processing method of porous media containing hydrate is described more containing hydrate
Hole medium is hydrated by what the sodium-chloride water solution that methane hydrate, natural sea sand and mass fraction are 3.5% formed containing methane
Object porous media, it is assumed that the chemical formula of methane hydrate is CH4·nH2O, n are hydrate index, are 5.89 progress phases to n values
It closes and calculates.For ease of description, porous media containing hydrate, the processing method is claimed to contain following steps below：
Step 1：Obtain during decomposition of hydrate the temperature of porous media and the pressure of gas in reaction kettle, and according to
The hydrate concentration in different decomposition stage is calculated by the calculation formula of hydrate concentration for temperature and pressure；The water
The calculation formula for closing object saturation degree is expressed as：
In formula, ShFor the hydrate concentration of the porous media containing hydrate；MhFor hydrate molal weight, unit：kg/
mol；ρhFor the density of hydrate, unit：kg/m3；VPFor the volume of porous media hole, unit：m3；VGTo react gas reactor
Volume, unit：m3；T1(i.e. Hydrate in Porous Medium does not have started generation or for porous media when to be free of hydrate
All decompose finish) temperature, unit：K；T2For during decomposition of hydrate when containing hydrate porous media temperature, it is single
Position：K；P1To be free of hydrate in porous media, (i.e. Hydrate in Porous Medium does not have started generation or all decomposes
Finish) when react gas reactor pressure, unit：Pa；P2It is single to react the pressure of gas reactor during decomposition of hydrate
Position：Pa；Z1To be free of hydrate in porous media, (i.e. Hydrate in Porous Medium does not have started generation or all decomposes
Finish) when react gas reactor compressibility factor；Z2To react the compressibility factor of gas reactor during decomposition of hydrate；R is
Mol gas constant, unit：J/(mol·K).
Step 2：Obtain lower porous Jie containing hydrate of different test frequency pumping signals effect during decomposition of hydrate
The impedance of matter, and obtained accordingly according to the impedance computation of the geometric dimension of electric sensor and distance and the porous media containing hydrate
Complex resistivity, and further obtain complex resistivity real part, imaginary part, amplitude and phase angle；Wherein test frequency is in 0.1Hz-
Value within the scope of 1MHz.
The amplitude of complex resistivity is calculated according to formula (2), the expression formula of formula (2) is：
In formula, | ρ | it is complex resistivity amplitude, unit：Ω·m；| Z | it is impedance magnitude, unit：Ω；S is the face of electrode slice
Product, unit：m2；L is the distance between electrode slice, unit：m.
In the present embodiment, the distance between electrode slice is 7cm, and the area of electrode slice is 25cm2。
The real and imaginary parts of complex resistivity are calculated according to formula (3), the expression formula of formula (3) is：
In formula, ρ ' is the real part of complex resistivity, unit：Ω·m；ρ " is the imaginary part of complex resistivity, unit：Ω·m；θ is
Phase angle.
Step 3：Under linear coordinate system, using hydrate concentration as reference axis horizontal axis, with complex resistivity real part, imaginary part,
The frequency dispersion degree of amplitude and phase angle is the reference axis longitudinal axis, draws frequency dispersion degree and the hydration of complex resistivity real part, imaginary part, amplitude and phase angle
Relational graph between object saturation degree, according to the frequency dispersion degree and hydrate concentration of complex resistivity real part, imaginary part, amplitude and phase angle it
Between relationship, hydrate concentration computation model is established, for determining hydrate concentration.
Frequency dispersion degree is bigger to indicate degree of dependence of the parameters such as complex resistivity real part, imaginary part, amplitude and phase angle to test frequency
It is stronger, i.e., it is more sensitive to the variation of test frequency.
The frequency dispersion degree of the complex resistivity amplitude is expressed as：
In formula, ZPFEFor the frequency dispersion degree of complex resistivity amplitude；ZLFor the amplitude of low frequency complex resistivity；ZHFor high frequency complex resistivity
Amplitude.
The frequency dispersion degree of the complex resistivity real part is expressed as：
In formula, RPFEFor the frequency dispersion degree of complex resistivity real part；RLFor low frequency complex resistivity real part；RHIt is real for high frequency complex resistivity
Portion.
The frequency dispersion degree of the complex resistivity imaginary part is expressed as：
In formula, XPFEFor the frequency dispersion degree of complex resistivity imaginary part；XLFor low frequency complex resistivity imaginary part；XHIt is empty for high frequency complex resistivity
Portion.
The frequency dispersion degree of the complex resistivity phase angle is expressed as：
In formula, θPFEFor the frequency dispersion degree of complex resistivity phase angle；θLFor low frequency complex resistivity phase angle；θHFor high frequency complex resistivity phase
Angle.
In the present embodiment, when calculating frequency dispersion is spent, selection 10Hz is low frequency point, 1MHz is high frequency points, in actually calculating
The high and low visual concrete condition of frequency point is suitably chosen.Complex resistance is calculated separately according to calculating formula (14), (15), (16), (17)
The frequency of the frequency dispersion degree of rate amplitude, the frequency dispersion degree of complex resistivity real part, the frequency dispersion degree of complex resistivity imaginary part and complex resistivity phase angle
Divergence.
When giving methane hydrate saturation degree difference referring to Figure 14, Figure 15, Figure 14, the telegram in reply of the porous media containing hydrate
The frequency dispersion degree situation of change of resistance rate amplitude；When Figure 15 gives methane hydrate saturation degree difference, the porous media containing hydrate
The frequency dispersion degree situation of change of complex resistivity real part.By Figure 14, Figure 15 it is found that within the scope of 10Hz~1MHz test frequencies, with
The frequency dispersion degree of the increase of hydrate concentration, complex resistivity amplitude and complex resistivity real part is gradually reduced.It is linear fitting it is found that
The frequency dispersion degree of complex resistivity amplitude, the frequency dispersion degree of complex resistivity real part and hydrate concentration are substantially in a linear relationship.
Embodiment provided above only with illustrating the present invention for convenience, and it is not intended to limit the protection scope of the present invention,
Technical solution scope of the present invention, person of ordinary skill in the field make various simple deformations and modification, should all include
In the above claim.
Claims (1)
1. a kind of electrical parameter test data processing method of porous media containing hydrate, which is characterized in that contain following steps：
Step 1: the temperature of porous media and the pressure of gas in reaction kettle during obtaining decomposition of hydrate, and according to temperature
The hydrate concentration in different decomposition stage is calculated by the calculation formula of hydrate concentration with pressure, specifically：
The calculation formula calculated the hydrate concentration according to temperature and pressure is expressed as：
In formula, ShFor the hydrate concentration of the porous media containing hydrate；MhFor hydrate molal weight, unit：kg/mol；ρh
For the density of hydrate, unit：kg/m3；VPFor the volume of porous media hole, unit：m3；VGTo react the body of gas reactor
Product, unit：m3；T1The temperature of porous media, unit when for without hydrate：K；T2To contain hydration during decomposition of hydrate
The temperature of porous media, unit when object：K；P1To be free of the pressure for reacting gas reactor when hydrate, unit in porous media：
Pa；P2To react the pressure of gas reactor, unit during decomposition of hydrate：Pa；Z1It is anti-when to be free of hydrate in porous media
Answer the compressibility factor of gas reactor；Z2To react the compressibility factor of gas reactor during decomposition of hydrate；R is that moles of gas is normal
Number, unit：J/(mol·K)；
Step 2: the porous media containing hydrate under different test frequency pumping signals effect during obtaining decomposition of hydrate
Impedance, and obtained accordingly again according to the impedance computation of the geometric dimension of electric sensor and distance and the porous media containing hydrate
Resistivity, and further obtain real part, imaginary part, amplitude and the phase angle of complex resistivity；
Step 3: drawing the impedance spectrum Nyquist figures of the porous media containing hydrate according to impedance, schemed by impedance spectrum Nyquist
To the basic structure of equivalent-circuit model corresponding thereto, according to hydrate concentration to equivalent circuit model parameter value
Affecting laws, the relationship being fitted between equivalent circuit model parameter and hydrate concentration, establish MODELING FOR COMPLEX RESISTIVITY MODEL, described multiple
The output parameter of resistivity models be hydrate concentration, input parameter be complex resistivity and test frequency, specifically：
The impedance spectrum Nyquist figures are a straight line with reference axis in 45 degree of angles, the basic structure of equivalent-circuit model
By with ω RWCWOne resistance R of=1 relationshipWWith a capacitance CWIt is composed in series, resistance RWWith capacitance CWRespectively with ω-1/2It is in
Linear relationship is expressed as：
RW=a ω-1/2+b (2)
CW=c ω-1/2+d (3)
In formula, a, b, c, d are equivalent circuit model parameter；
The expression formula of the equivalent-circuit model is：
In formula, ρ is complex resistivity；S is the area of electric sensor electrode slice, unit：m2；L is between electric sensor electrode slice
Distance, unit：m；ω is angular frequency, unit：The π of rad/s, ω=2 f, f are test frequency, unit：Hz；
It is further represented as：
According to hydrate concentration to the affecting laws of equivalent circuit model parameter value, fitting equivalent circuit model parameter a, b,
C, the relationship between d and hydrate concentration, the pass between hydrate concentration respectively equivalent circuit model parameter a, b, c, d
System is expressed as：
A=782.78Sh 2-177.1Sh+406.23 R2=0.9583 (6)
B=40.401Sh+30.732 R2=0.9882 (7)
C=0.0023Sh 0.0344 R2=0.9207 (8)
D=0.00002 (9)
In formula, R2For degree of fitting, i.e. matched curve and for the degree of correlation between the data of fitting；
And then the relationship between the complex resistivity and hydrate concentration of the porous media containing hydrate is obtained, establish complex resistivity mould
Type, the MODELING FOR COMPLEX RESISTIVITY MODEL are expressed as：
In formula, ShFor the hydrate concentration of the porous media containing hydrate；
Step 4: according to the MODELING FOR COMPLEX RESISTIVITY MODEL, the electric sensor geometric dimension actually measured and distance and test frequency,
Utilize solution by iterative method hydrate concentration.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610842507.0A CN106443254B (en) | 2016-09-22 | 2016-09-22 | The electrical parameter test data processing method of porous media containing hydrate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610842507.0A CN106443254B (en) | 2016-09-22 | 2016-09-22 | The electrical parameter test data processing method of porous media containing hydrate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106443254A CN106443254A (en) | 2017-02-22 |
| CN106443254B true CN106443254B (en) | 2018-08-28 |
Family
ID=58166747
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610842507.0A Active CN106443254B (en) | 2016-09-22 | 2016-09-22 | The electrical parameter test data processing method of porous media containing hydrate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106443254B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115015086B (en) * | 2022-07-26 | 2024-01-26 | 中国石油大学(华东) | Hydrate stratum permeability on-site in-situ quantitative evaluation method based on complex conductivity |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105334546A (en) * | 2015-09-23 | 2016-02-17 | 中国石油大学(华东) | Simulated experiment testing method of gas hydrate in porous medium |
| CN105547964A (en) * | 2016-01-15 | 2016-05-04 | 中国石油大学(华东) | One-way flowing determination method for air-water relative permeability under state of different hydrate saturations |
| CN105675441A (en) * | 2016-01-15 | 2016-06-15 | 中国石油大学(华东) | Gas-water relative permeability measurement method under radial flow condition at different hydrate saturations |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| BR112015004403A2 (en) * | 2012-11-16 | 2017-07-04 | Chevron Usa Inc | methods and systems for determining minimum porosity for clathrate in sediment |
-
2016
- 2016-09-22 CN CN201610842507.0A patent/CN106443254B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105334546A (en) * | 2015-09-23 | 2016-02-17 | 中国石油大学(华东) | Simulated experiment testing method of gas hydrate in porous medium |
| CN105547964A (en) * | 2016-01-15 | 2016-05-04 | 中国石油大学(华东) | One-way flowing determination method for air-water relative permeability under state of different hydrate saturations |
| CN105675441A (en) * | 2016-01-15 | 2016-06-15 | 中国石油大学(华东) | Gas-water relative permeability measurement method under radial flow condition at different hydrate saturations |
Non-Patent Citations (1)
| Title |
|---|
| 多孔介质中水合物的热物理参数测量;刁少波等;《岩矿测试》;20080630;第27卷(第3期);第165-168页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106443254A (en) | 2017-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | A finite-element time-domain forward solver for electromagnetic methods with complex-shaped loop sources | |
| Streich | 3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: Direct solution and optimization for high accuracy | |
| Vinje et al. | Traveltime and amplitude estimation using wavefront construction | |
| Key et al. | Adaptive finite-element modeling using unstructured grids: The 2D magnetotelluric example | |
| Cross et al. | Construction and application of a stratigraphic inverse model | |
| US10788598B2 (en) | Method of predicting parameters of a geological formation | |
| Kretz et al. | An integrated reservoir characterization study matching production data and 4D seismic | |
| CN113887078B (en) | Deep-ultra-deep sea-facies carbonate reservoir occurrence depth lower limit prediction method | |
| CN106443254B (en) | The electrical parameter test data processing method of porous media containing hydrate | |
| Keller et al. | SHEMAT-Suite: An open-source code for simulating flow, heat and species transport in porous media | |
| CN111897012A (en) | Shale gas single-well geological comprehensive evaluation method | |
| CN109577969B (en) | Method for calculating pore pressure of carbonate rock stratum based on rock compression coefficient | |
| GUO et al. | Forward-modeling and feasibility study of non-repeating acquired time-lapse seismic exploration | |
| CN109582923A (en) | A kind of method and processing terminal based on two-parameter calculated hydration object saturation degree | |
| Moczo et al. | Numerical wave propagation simulation | |
| CN114428286A (en) | Gas saturation prediction method based on pre-stack seismic data | |
| Zhao et al. | Reliability experiment on the calculation of natural gas hydrate saturation using acoustic wave data | |
| Li et al. | Discrimination of effective source rocks and evaluation of the hydrocarbon resource potential in Marsel, Kazakhstan | |
| Riva et al. | Exact sensitivity matrix and influence of the number of pilot points in the geostatistical inversion of moment equations of groundwater flow | |
| Jiang et al. | Determination of vertical hydraulic conductivity of aquitards in a multilayered leaky system using water-level signals in adjacent aquifers | |
| Wilkinson et al. | Uncertainty in regional estimates of capacity for carbon capture and storage | |
| Gasperikova et al. | Kimberlina 1.2 ccus geophysical models and synthetic data sets | |
| Qi et al. | Analytical solution for stress distribution around arbitrary stopes using evolutionary complex variable methods | |
| Liu et al. | Performance investigations of auxiliary‐space Maxwell solver preconditioned iterative algorithm for controlled‐source electromagnetic induction problems with electrical anisotropy | |
| US12019038B2 (en) | Evaluation of source rock samples from subterranean reservoirs |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |