CN107045154A - A kind of method and apparatus of identification attitude of stratum in horizontal well environment - Google Patents

A kind of method and apparatus of identification attitude of stratum in horizontal well environment Download PDF

Info

Publication number
CN107045154A
CN107045154A CN201710069496.1A CN201710069496A CN107045154A CN 107045154 A CN107045154 A CN 107045154A CN 201710069496 A CN201710069496 A CN 201710069496A CN 107045154 A CN107045154 A CN 107045154A
Authority
CN
China
Prior art keywords
stratum
value
horizontal well
well
orientation
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.)
Pending
Application number
CN201710069496.1A
Other languages
Chinese (zh)
Inventor
张伟
梁卫
吴意明
杨玉卿
徐大年
张中庆
罗东红
王瑛
张国华
张宏伟
刘保银
朱佳音
商超
李建
万钧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Oilfield Services Ltd
China National Offshore Oil Corp CNOOC
Original Assignee
China Oilfield Services Ltd
China National Offshore Oil Corp CNOOC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Oilfield Services Ltd, China National Offshore Oil Corp CNOOC filed Critical China Oilfield Services Ltd
Priority to CN201710069496.1A priority Critical patent/CN107045154A/en
Publication of CN107045154A publication Critical patent/CN107045154A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Abstract

The invention discloses a kind of method and apparatus of the identification attitude of stratum in horizontal well environment, this method includes:Stratigraphic dip and formation azimuth angle are extracted according to orientation Natural Gamma-ray Logging Curves;Horizontal well well logging data are pre-processed using stratigraphic dip and formation azimuth angle, curve relevant comparative is carried out to pretreated log data and activity layering obtains bed boundary, stratum initial model is built;The forward response of three-dimensional horizontal well stratum initial model is calculated using finite element method, multi-parameter iterative inversion is carried out with least square method;Well track position of stratum graph of a relation is drawn according to inversion result, stratigraphic sequence relation etc. and inversion result is exported.Pass through the solution of the present invention, accurate formation parameter and anisotropy information, measurement point are obtained to layer frontier distance information, perfect Horizontal Well Log Interpretation data currently, is evaluated, saturation computation and reserves calculate and provide accurate reliable parameter for later stage well log interpretation.

Description

A kind of method and apparatus of identification attitude of stratum in horizontal well environment
Technical field
The present invention relates to a kind of side of the identification attitude of stratum in oil field development technical field, more particularly to horizontal well environment Method and device.
Background technology
As oil-gas exploration and development difficulty is continuously increased, logging while drilling technology research obtains foot with connector for logging while drilling research and development Enough attention, and obtained quick development.With the extensive use of the complicated technology well such as high angle hole, horizontal well, to complexity The research in construction area requires that the relation of tomography and layer interdigit is finer, generally requires to break away from traditional 2D seismic interpretation sides Method.In order to solve the problem of oil development exploration, the new method and new technology how everybody paid close attention to change into the life of reality Force of labor, the practical demand in reality that solves turns into industry focus of attention.International three big oil well logging service companies release one after another With position sensing capability with orientation electromagnetic wave logging instrument is bored, Schlumberger is proposed PeriScope side within 2005 Position resistivity measuring instrument, Baker Hughes Inc in 2006 is proposed with brill orientation electromagnetic wave logger APR, Halliburton in 2007 Company is proposed azimuthal resistivity measuring instrument ADR.Traditional is coaxial with receiving coil with electromagnetic wave logging instrument transmitting coil is bored, The average value that obtained formation signal is formation information is measured, does not possess azimuth characteristic, it is impossible to accurately obtain formation resistivity The information such as anisotropy.Orientation is very different with electromagnetic wave logging instrument is bored with traditional instrument, orientation electromagnetic wave measurement instrument Device is mixed using axially inclined or transverse coil, can preferably provide formation azimuth information, indicates the anisotropy on stratum, And recognize stratigraphic boundary.
Attitude of stratum identification uses horizontal well well logging data parameters and extracted and 3-d inversion, is obtained by 3-d inversion Attitude of stratum parameter must be described.Compared with straight well, horizontal well and high angle hole borehole environment are generally more complicated, and log is removed Outside by mud, well, intrusion effect, while also influenceed by factors such as country rock, layer frontier distance, inclination angle, formation anisotropies, The asymmetry on stratum causes traditional one-dimensional, two-dimension method to be no longer applicable simultaneously.The formation parameter of electric logging response is three-dimensional anti- Drill and belong to nonlinear problem, the multi-solution of the complexity of ground environment, the one-sidedness of single well logging information and inverting adds survey Well inverting, well log interpretation and the difficulty evaluated.Joint inversion provides effective way to solve this problem.Joint inversion refers to Geological model parametric inversion is carried out using two or more log data of different physical mechanisms.Joint inversion is in itself It is by increasing the effective information of particular probe target (increase constraint), to reach the mesh for more accurately reflecting geologic objective body 's.
Homologous physical property is that three-dimensional joint inversion increases the primary condition of the source effective information.Same well difference well logging system Resistivity logging while drilling well logging and the cable resistance rate of row are logged well, and have identical target properties measurement item for identical undisturbed formation Mesh, makes it possible its three-dimensional joint inversion.With bore electric logging and nuclear magnetic resonance log due to logger between the two detection feature and Each the detection feature between difference log has differences, therefore the more single inverting of three-dimensional joint inversion has what is quite enriched Well logging information.
With three-dimensional joint inversion is bored based on damped least square method, electric logging while drilling and Orientation density are logged well, Azimuth gamma well-logging response, well logging moment stratigraphic model parameter organic unity get up to carry out inverting.But, traditional resistor rate Joint inversion carries out one-dimensional and (or) two dimensional inversion primarily directed to rotational symmetry stratum, does not both apply to asymmetricly Rotating fields Horizontal well and high angle hole, the attitude of stratum information such as anisotropy, stratigraphic dip and layer frontier distance can not be provided.Cause This, it would be desirable to be able to the stratum based on three dimensions stratigraphic model of the information such as stratigraphic dip, layer frontier distance and anisotropy is provided Occurrence recognition methods.
The content of the invention
In order to solve the above problems, the present invention proposes the method and dress of the identification attitude of stratum in a kind of horizontal well environment Put, result in accurate formation parameter and anisotropy information, measurement point to layer frontier distance information, improve current level Borehole logging tool interpretation data, is evaluated, saturation computation and reserves calculate and provide accurate reliable parameter for later stage well log interpretation.
In order to achieve the above object, the method that the present invention proposes the identification attitude of stratum in a kind of horizontal well environment, should Method includes:
The orientation Natural Gamma-ray Logging Curves or Orientation density log obtained according to horizontal well well logging extract ground Inclination layer and formation azimuth angle;
Obtain multiple sources away from the horizontal well well logging data under multiple frequencies;Stratigraphic dip and stratum using extraction Azimuth is pre-processed to horizontal well well logging data, the log data standardized;To the log data of standardization Curve relevant comparative is carried out, and stratum AUTOMATIC ZONING is carried out using activity method, bed boundary is obtained;Wherein, horizontal well is surveyed with brill Well data include instrumental azimuth resistivity log response value, density and porosity curve;
The ground stratum boundary divided according to horizontal well well logging data, offset well log data, the stratigraphic dip, AUTOMATIC ZONINGization Face builds stratum initial model, and stratum initial model is initialized and parameter setting;
Using the three-dimensional horizontal well stratum initial model of finite element method calculating with orientation instrument forward response is bored, contrast is real Data and the forward response of stratum initial model are surveyed, object function is set up, multi-parameter iterative inversion is carried out with least square method The model knots modification for solving object function optimal solution obtains inversion result;
According to the distance of inversion result, stratigraphic sequence relation, pendage situation and measurement point to layer border, well is drawn Track-position of stratum graph of a relation, and coloured according to every layer of true formation resistivity inversion result, display attitude of stratum letter directly perceived Breath, and export inversion result.
Alternatively, the orientation Natural Gamma-ray Logging Curves or Orientation density log obtained according to horizontal well well logging Extract stratigraphic dip and formation azimuth angle refers to:According to orientation Natural Gamma-ray Logging Curves or the elevation of Orientation density log Difference and non-linear fitting method extract the stratigraphic dip and formation azimuth angle;
It is described according to the orientation Natural Gamma-ray Logging Curves or the depth displacement and Nonlinear Quasi of Orientation density log Conjunction method, which extracts the stratigraphic dip and formation azimuth angle, to be included:
Every orientation natural gamma in a plurality of orientation Natural Gamma-ray Logging Curves measured orientation gamma-ray tool device Log carries out relevant comparative's analysis, obtains the depth displacement on correspondence stratum, obtains the measurement point on pendage bed boundary;And The expanded view that the borehole wall intersects with pendage bed boundary is represented by the first expression formula;First expression formula includes:Y=A sin(ωx-β)+y0=y0+A cosβsin(ωx)+(-A sinβ)cos(ωx);Wherein, y represents the borehole wall and pendage layer The functional value of interface intersection;A representative function value changes amplitudes;y0Representative function value y average;β representative function values y is met just The initial phase of string function;ω representative function values y meets the cycle of SIN function;
Undetermined coefficient in first expression formula is solved by nonlinear fitting, the first matrix side is obtained using least square method Journey;First matrix equation includes:Wherein, a0=y0, a1=A cos β, a2=-A sin β,a0、a1、a2Represent the system of polynomials of fitting Number;Represent independent variable value during fitting function;
First matrix equation is solved, sinusoids parameter expression formula and inclination angle tendency expression formula is obtained;Sine curve is joined Counting expression formula is:It is inclined to expression formula in inclination angle:Wherein, Dip represents ground Inclination layer, Dir represents tendency, DelRepresent investigation depth, yminRepresent the minimum point on one week that the borehole wall intersects with dipping bed.
Alternatively, pretreatment bag is carried out to horizontal well well logging data using the stratigraphic dip and formation azimuth angle of extraction Include:
Utilize the azimuthal resistivity of stratigraphic dip and formation azimuth angle to the azimuth information in horizontal well well logging data Log data carries out directional correction and the disposal of gentle filter, and carries out borehole environment correction to the data after smothing filtering.
Alternatively, using formation azimuth angle to the azimuthal resistivity log data progress side in horizontal well well logging data Bit correction includes:
The orientation curve in horizontal well well logging is fitted, and is expressed with first orientation curvilinear function;Wherein, first orientation Curvilinear function includes:Y (x)=A sin (x+ φ)+B;X is formation azimuth angle;Y (x) is the function of first orientation curvilinear function Value;B is the y of first orientation curvilinear function average value;A is amplitude when first orientation curve meets SIN function;φ is first Beginning phase;
Multiple characteristic values on first orientation curvilinear function are chosen, and multiple characteristic values are substituted into first orientation curve letter Number, asks for A, B and φ;
A, B and φ for obtaining are substituted into first orientation curvilinear function, the first equation with known coefficient is obtained;
Using the corresponding functional value of the first equation as azimuthal resistivity log data directional correction value.
Alternatively, multiple characteristic values include:(0, y (0)), (1, y (1)), (2, y (2)), (3, y (3));A, B and φ value Respectively:
φ=a tan (y (0)-B)/(y (2)-B);
A=(y (3)-y (1))/(2sin (φ)).
Alternatively, the method for smothing filtering includes:Least square moving average method and weighted moving average.
Alternatively, the data after smothing filtering are carried out with borehole environment correction includes:
According to log data apparent resistivity and borehole environment influence factor, inquire about in the borehole correction charts storehouse pre-established Correction coefficient, according to correction coefficient to log data apparent resistivity carry out borehole environment correction;
Wherein, pre-establishing the borehole correction charts storehouse includes:
Fixed borehole environment influence factor, changes stratum background resistivity, and borehole correction is calculated using forward modeling simulated program Plate;Borehole environment influence factor includes borehole size and drilling fluid resistivity;
The borehole correction charts in the case of different borehole size and drilling fluid resistivity are calculated, and set up borehole correction charts Storehouse.
Alternatively, the log data progress curve relevant comparative of standardization is included:To the horizontal well after smothing filtering with The log for boring different investigation depths in log data is normalized, and the data after normalization are entered under fixed window length Row correlation is contrasted;
Wherein, the relative coefficient in correlation contrast is calculated according to following second relational expressions:
In formula, i represents sampling point position, and n represents that window is long,Contrast bar plot mean in window length is represented,It is right to represent Than bar plot mean, xiRepresent correlation curve ith sample point value, yiRepresent by correlation curve ith sample point curve value, ZiRepresent coefficient correlation;
Carrying out stratum AUTOMATIC ZONING using activity method includes:
The curve activity of calculated level well well logging data;
Whether judgment curves activity is more than default activity threshold value, if curve activity is more than default activity threshold value, The position that will be greater than the extreme point of activity threshold value is defined as layering interfaces.
Alternatively, the log of different investigation depths in the horizontal well well logging data after smothing filtering is returned One change includes:Extreme value normalization, porosity normalization and density normalization;
Wherein, resistivity sampling extreme value method for normalizing is included:
Neutron porosity sampling linear method for normalizing is included:
Neutron density sampling linear method for normalizing is included:
Wherein, xijRepresent the ith sample point value of j-th strip curve, xmin jRepresent the minimum value of j-th strip curve, xmax jTable Show the maximum of j-th strip curve, YijRepresent i-th of value of j-th strip curve after normalization;TNPLjRepresent porosity curve j-th Value, and replace original value with the value after normalization;ALCDjRepresent that neutron density is worth for j-th, and original is replaced with the value after normalizing Initial value.
Alternatively, curve activity is defined using the 3rd following relational expressions:
Wherein, E (d) represents d activity function value, and x (i) represents log measured value,Represent that log exists Average value in interval [d-n, d+n];D represents sampling point position;N represents that sample window is long;I is represented in a sample window length Several sampled points.
Alternatively, using the three-dimensional horizontal well stratum initial model of finite element method calculating with brill orientation instrument forward response Including:
Maxwell's Maxwell equations under given boundary condition are solved, Maxwell is equations turned for wave equation:Wherein, E represents electric-field intensity;μ represents amplitude;ω represents angular frequency;ε represents dielectric Constant;
Wave equation is attributed to field energy functional with reference to boundary condition:
Wherein,Table Show object function;μ0Represent initial amplitude;V represents domain;ω represents angular frequency;ε represents dielectric constant;J represents electric current Density;
Using finite elements subdivision field domain, and corresponding Interpolation-Radix-Function is chosen, spatial spreading is carried out to energy functional, obtained To following discretization functional form:
Wherein, AeRepresentative function It is discretely formed coefficient matrix;BeRepresent equation right-hand vector discrete matrix;CeRepresentative function value discrete matrix;EeRepresent independent variable square Battle array;M represents independent variable number;The transposition of T representing matrixs;
The parameter needed for acquisition inverting is solved to discretization functional.
Alternatively, the model that the multi-parameter iterative inversion solution object function optimal solution is carried out with least square method changes Variable, which obtains inversion result, to be included:
The residual error of actual measured value and the analogue value is solved by least square method:
m≥n;Wherein, r:Rn→RmIt is x nonlinear function;R tables Show independent variable vector;M represents unknown quantity number;N represents forward response curve number;
The Jacobian matrix of the object function is solved by gradient method, Jacobi's system of linear equations is formed, institute is solved State the direction of steepest descent that system of linear equations calculates the object function;Wherein, calculate described according to following 4th relational expression The gradient of object function:J (x) is r (x) Jacobian matrix;
Using golden section approach, determine to change step-length on the direction of steepest descent;
Direction is changed according to the model of determination and changes step-length and determines model knots modification, and mould is changed according to model knots modification Type, completes an iteration inverting;
According to default stopping criterion for iteration, recursive call iterative inversion, until stopping criterion for iteration is met, output knot Really.
Alternatively, according to the distance of inversion result, stratigraphic sequence relation, pendage situation and measurement point to layer border, Drawing well track-position of stratum graph of a relation includes:
By AUTOMATIC ZONING to stratum division result, definitely order relation layer by layer;
The measurement point obtained by inverting is to layer frontier distance, stratigraphic dip inversion result, inverting formation true dip angle, layer Thick, well track passes through bed boundary position and well track and bed boundary relative position relation;
According to the vertical continuity of bed boundary and stratum true dip angle situation of change, the ground outside instrument investigative range is estimated Bed boundary is moved towards and thickness change, and non-inverting subterranean formation zone outside instrument investigative range is continued according to existing inversion result, To ensure the continuity and integrality of attitude of stratum figure;
According to formation resistivity inversion result and Anisotropic inversion result, resistivity value is converted into chromatic value and is filled in In corresponding stratum, and change the change for qualitatively judging resistivity value by color light and shade;
The attitude of stratum two dimensional image ultimately formed shows that in the horizontal Different Strata has different resistivity value, correspondence Different colours, are changed by change in resistance, bed boundary change in location and thickness in color change marker;In the vertical, Fill Color change represents the change occurred in layer with horizontal displacement changes in resistance rate in layer, is used to refer to subterranean formation zone continuous The quality of property.
In order to achieve the above object, the invention also provides a kind of device of the identification attitude of stratum in horizontal well environment, The device includes:Extraction module, processing module, modeling module, computing module and display module;
Extraction module, orientation Natural Gamma-ray Logging Curves or Orientation density for being obtained according to horizontal well well logging are surveyed Well curve extracts stratigraphic dip and formation azimuth angle;
Processing module, for obtain multiple sources away from the horizontal well well logging data under multiple frequencies;Utilize extraction Stratigraphic dip and the formation azimuth angle is pre-processed to the horizontal well well logging data, the well logging number standardized According to;Log data to the standardization carries out curve relevant comparative, and carries out stratum AUTOMATIC ZONING using activity method, obtains ground Bed boundary;Wherein, the horizontal well well logging data are write music including instrumental azimuth resistivity log response value, density and hole Line;
Modeling module, for according to the horizontal well well logging data, offset well log data, the stratigraphic dip, from The bed boundary of dynamic stratification point builds stratum initial model, and the stratum initial model is initialized and parameter is set Put;
Computing module, for using finite element method calculate three-dimensional horizontal well stratum initial model with boring orientation instrument just Response, contrast measured data and the forward response of stratum initial model are drilled, object function is set up, is carried out with least square method many The model knots modification that parameter iteration inverting solves the object function optimal solution obtains inversion result;
Display module, for being arrived according to the inversion result, stratigraphic sequence relation, the pendage situation and measurement point The distance on layer border, draws well track-position of stratum graph of a relation, and carry out according to every layer of true formation resistivity inversion result Color, display attitude of stratum information directly perceived, and export inversion result.
The embodiment of the present invention includes:The orientation Natural Gamma-ray Logging Curves or orientation obtained according to horizontal well well logging are close Spend log and extract stratigraphic dip and formation azimuth angle;Obtain multiple sources away from the horizontal well well logging number under multiple frequencies According to;Horizontal well well logging data are pre-processed using the stratigraphic dip and formation azimuth angle of extraction, standardized Log data;Log data to standardization carries out curve relevant comparative, and carries out stratum AUTOMATIC ZONING using activity method, obtains Bed boundary;Wherein, horizontal well well logging data are write music including instrumental azimuth resistivity log response value, density and hole Line;The bed boundary structure divided according to horizontal well well logging data, offset well log data, the stratigraphic dip, AUTOMATIC ZONINGization Stratum initial model is built, and stratum initial model is initialized and parameter setting;Three-dimensional water is calculated using finite element method Horizontal well stratum initial model with boring orientation instrument forward response, the forward response of contrast measured data and stratum initial model, Object function is set up, the model knots modification that multi-parameter iterative inversion solution object function optimal solution is carried out with least square method is obtained To inversion result;According to the distance of inversion result, stratigraphic sequence relation, pendage situation and measurement point to layer border, draw Well track-position of stratum graph of a relation, and coloured according to every layer of true formation resistivity inversion result, display stratum directly perceived production Shape information, and export inversion result.By the solution of the present invention, obtain accurate formation parameter and anisotropy information, survey Amount point arrives layer frontier distance information, perfect Horizontal Well Log Interpretation data currently, be later stage well log interpretation evaluate, saturation degree Calculate and reserves calculate and provide accurate reliable parameter.
Compared with prior art and product, scheme of the embodiment of the present invention has the following advantages:
The embodiment of the present invention is using attitude of stratum, this method in the three-dimensional horizontal well stratum of orientation resistivity logging while drilling well logging recognition Stratigraphic dip, measurement point are resulted in layer frontier distance, mud using azimuthal resistivity log data inverting attitude of stratum parameter Depth of invasion, invaded zone resistivity, formation resistivity, formation anisotropy, well track-position of stratum relation etc. are starched, it is perfect Existing logging evaluation data, can obtain comprehensive formation information.In later stage interpretation and evaluation, by every layer of formation resistivity it is each to The opposite sex, layer frontier distance inversion result can accurately calculate reservoir oil saturation, mobile oil saturation, Reservoir Thickness etc., be Evaluating reservoir provides reliable parameter.
Brief description of the drawings
The accompanying drawing in the embodiment of the present invention is illustrated below, the accompanying drawing in embodiment is for entering one to the present invention Step understands, is used to explain the present invention together with specification, does not constitute limiting the scope of the invention.
Fig. 1 recognizes the method stream of attitude of stratum to be utilized in the horizontal well environment of the embodiment of the present invention with brill 3 D resistivity Cheng Tu;
Fig. 2 is the horizontal well stratum modeling schematic diagram of the embodiment of the present invention;
Fig. 3 is for orientation under the horizontal well of the embodiment of the present invention with brill apparatus measures with stratigraphic dip change schematic diagram;
Fig. 4 for apparent resistivity under the horizontal well environment of the embodiment of the present invention phase difference resistivity and Amplitude Ratio differential resistivity Value is with layer frontier distance change schematic diagram;
Fig. 5 is apparent resistivity in anisotropic formation in the horizontal well of the embodiment of the present invention with change of pitch angle schematic diagram;
Fig. 6 draws well rail for the Inversion for the parameters of formation result to certain section of measured level well data of the embodiment of the present invention Mark-position of stratum graph of a relation;
Fig. 7 recognizes the device group of attitude of stratum to be utilized in the horizontal well environment of the embodiment of the present invention with brill 3 D resistivity Into block diagram.
Embodiment
For the ease of the understanding of those skilled in the art, the invention will be further described below in conjunction with the accompanying drawings, not It can be used for limiting the scope of the invention.
Scheme of the embodiment of the present invention is entered to solve the joint inversion of traditional resistor rate primarily directed to rotational symmetry stratum The one-dimensional and (or) two dimensional inversion of row, does not both apply to the horizontal well and high angle hole of asymmetricly Rotating fields, can not provide it is each to The problem of attitude of stratum information such as the opposite sex, stratigraphic dip and layer frontier distance, it is proposed that in a kind of new horizontal well environment Recognize the scheme of attitude of stratum.This solution provides can obtain the initial based on three dimensions horizontal well stratum of formation parameter information Model with bore azimuthal resistivity Log Strata occurrence recognition methods, and the program utilize horizontal well in azimuthal resistivity log well Data combine apparent resistivity curve, orientation gamma curve and neutron density porosity to anisotropic fine indicative function Curve, realizes Inversion for the parameters of formation and attitude of stratum identification in three dimensions horizontal well ground environment.
In order to achieve the above object, the method that the present invention proposes the identification attitude of stratum in a kind of horizontal well environment, such as Shown in Fig. 1, this method includes S101-S105:
S101, the orientation Natural Gamma-ray Logging Curves or Orientation density log that are obtained according to horizontal well well logging are carried Take stratigraphic dip and formation azimuth angle.
Alternatively, the orientation Natural Gamma-ray Logging Curves or Orientation density log obtained according to horizontal well well logging Extract stratigraphic dip and formation azimuth angle refers to:According to orientation Natural Gamma-ray Logging Curves or the elevation of Orientation density log Difference and non-linear fitting method extract the stratigraphic dip and formation azimuth angle.
It is described according to the orientation Natural Gamma-ray Logging Curves or the depth displacement and Nonlinear Quasi of Orientation density log Conjunction method, which extracts the stratigraphic dip and formation azimuth angle, includes step S201-S203:
Every orientation nature in S201, a plurality of orientation Natural Gamma-ray Logging Curves measured orientation gamma-ray tool device GR logging curve carries out relevant comparative's analysis, obtains the depth displacement on correspondence stratum, obtains the measurement on pendage bed boundary Point;And the expanded view that the borehole wall intersects with pendage bed boundary is represented by the first expression formula.
Alternatively, first expression formula includes:
Y=A sin (ω x- β)+y0=y0+A cosβsin(ωx)+(-A sinβ)cos(ωx);
Wherein, y represents the functional value that the borehole wall intersects with pendage bed boundary;A representative function value changes amplitudes;y0Represent Functional value y average;β representative function values y meets the initial phase of SIN function;ω representative function values y meets sinusoidal letter Several cycles.
S202, by nonlinear fitting solve the first expression formula in undetermined coefficient, utilize least square method obtain first Matrix equation.
Alternatively, the first matrix equation includes:Wherein, a0= y0, a1=A cos β, a2=-A sin β,a0、a1、a2Represent to intend The multinomial coefficient of conjunction;Represent independent variable value during fitting function.
S203, to the first matrix equation solve, obtain sinusoids parameter expression formula and inclination angle tendency expression formula;It is sinusoidal bent Line parameter expression is:It is inclined to expression formula in inclination angle:Wherein, Dip tables Show stratigraphic dip, Dir represents to be inclined to (or formation azimuth angle), DelRepresent investigation depth, yminRepresent what the borehole wall intersected with dipping bed Minimum point on one week, as shown in Figure 3.
S102, obtain multiple sources away from the horizontal well well logging data under multiple frequencies;Utilize the stratigraphic dip of extraction Horizontal well well logging data are pre-processed with formation azimuth angle, the log data standardized;Survey to standardization Well data carry out curve relevant comparative, and carry out stratum AUTOMATIC ZONING using activity method, obtain bed boundary;Wherein, horizontal well Well logging data include instrumental azimuth resistivity log response value, density and porosity curve.
Alternatively, pretreatment bag is carried out to horizontal well well logging data using the stratigraphic dip and formation azimuth angle of extraction Include:
Using stratigraphic dip and formation azimuth angle to the orientation resistance in the azimuth information in horizontal well well logging data Rate well logging, azimuthal resistivity log data carry out directional correction and the disposal of gentle filter, and the data after smothing filtering are carried out Borehole environment is corrected.
Alternatively, using formation azimuth angle to the azimuthal resistivity well logging in horizontal well well logging data, orientation resistance Rate log data, which carries out directional correction, includes S301-S304:
Orientation curve in S301, fitting horizontal well well logging, and expressed with first orientation curvilinear function;Wherein, One orientation curvilinear function includes:Y (x)=A sin (x+ φ)+B;X is formation azimuth angle;Y (x) is first orientation curvilinear function Functional value;B is the y of first orientation curvilinear function average value;A is amplitude when first orientation curve meets SIN function;φ For initial phase.
S302, the multiple characteristic values chosen on first orientation curvilinear function, and multiple characteristic values are substituted into first orientation song Line function, asks for A, B and φ;
Alternatively, multiple characteristic values include:(0, y (0)), (1, y (1)), (2, y (2)), (3, y (3));A, B and φ value Respectively:
φ=a tan (y (0)-B)/(y (2)-B);
A=(y (3)-y (1))/(2sin (φ)).
S303, by A, B and φ for obtaining substitute into first orientation curvilinear function, obtain with known coefficient the first equation Formula.
S304, using the corresponding functional value of the first equation as azimuthal resistivity log data directional correction value.
In embodiments of the present invention, to orientation resistivity logging, azimuthal resistivity log data carry out directional correction with Afterwards, before relevant comparative and AUTOMATIC ZONING is carried out, the data of measurement need to further be pre-processed, further pre-processed Content is filtered including line smoothing, removes invalid data and missing data is carried out to supplement to greatest extent in zone of reasonableness.
Alternatively, the method for smothing filtering includes:Least square moving average method and weighted moving average.
In embodiments of the present invention, above two smooth filtering method contains a variety of concrete implementation types again, specifically Including:
Linear function is smooth:
Quadratic function is smooth:
Bell shaped function exponential smoothing:
Hamming function exponential smoothing:
Wherein, TiIth sample value is represented, m represents filter window,Represent ith sample point value after filtering.
Alternatively, the data after smothing filtering are carried out with borehole environment correction includes:
According to log data apparent resistivity and borehole environment influence factor, inquire about in the borehole correction charts storehouse pre-established Correction coefficient, according to correction coefficient to log data apparent resistivity carry out borehole environment correction;
Wherein, pre-establishing the borehole correction charts storehouse includes:
Fixed borehole environment influence factor, changes stratum background resistivity, and borehole correction is calculated using forward modeling simulated program Plate;Borehole environment influence factor includes borehole size and drilling fluid resistivity;
The borehole correction charts in the case of different borehole size and drilling fluid resistivity are calculated, and set up borehole correction charts Storehouse.
Alternatively, the log data progress curve relevant comparative of standardization is included:To the horizontal well after smothing filtering with The log for boring different investigation depths in log data is normalized, and the data after normalization are entered under fixed window length Row correlation is contrasted.
In embodiments of the present invention, to smothing filtering and invalid data processing after horizontal well well logging data need into Row normalized, to carry out correlation contrast under fixed window length to the data after normalization.
Alternatively, the log of different investigation depths in the horizontal well well logging data after smothing filtering is returned One change includes:Extreme value normalization, porosity normalization and density normalization;
Wherein, resistivity sampling extreme value method for normalizing is included:
Neutron porosity sampling linear method for normalizing is included:
Neutron density sampling linear method for normalizing is included:
Wherein, xijRepresent the ith sample point value of j-th strip curve, xmin jRepresent the minimum value of j-th strip curve, xmax jTable Show the maximum of j-th strip curve, YijRepresent i-th of value of j-th strip curve after normalization;TNPLjRepresent porosity curve j-th Value, and replace original value with the value after normalization;ALCDjRepresent that neutron density is worth for j-th, and original is replaced with the value after normalizing Initial value.
In embodiments of the present invention, the log of different investigation depths is returned in horizontal well well logging data After one change processing, correlation contrast is further carried out under fixed window length to the data after normalization.
Alternatively, the relative coefficient in correlation contrast is calculated according to following second relational expressions:
Wherein, i represents sampling point position, and n represents that window is long,Contrast bar plot mean in window length is represented,It is right to represent Than bar plot mean, xiRepresent correlation curve ith sample point value, yiRepresent by correlation curve ith sample point curve value, ZiRepresent coefficient correlation.
In embodiments of the present invention, the log of different investigation depths is returned in horizontal well well logging data After one change processing, in addition it is also necessary to further carry out stratum AUTOMATIC ZONING using activity method.
Alternatively, carrying out stratum AUTOMATIC ZONING using activity method includes S401-S402:
S401, calculated level well well logging data curve activity.
Alternatively, curve activity is defined using the 3rd following relational expressions:
Wherein, E (d) represents d activity function value, and x (i) represents log measured value,Represent that log exists Average value in interval [d-n, d+n];D represents sampling point position;N represents that sample window is long;I is represented in a sample window length Several sampled points.
Whether S402, judgment curves activity are more than default activity threshold value, if curve activity is more than default activity threshold Value, then the position that will be greater than the extreme point of activity threshold value is defined as layering interfaces.
In embodiments of the present invention, the step S402 layering interfaces divided are needed to be screened and screened, reject pseudobed, Obtain true bed boundary.Specifically, the present embodiment combination horizontal well bed boundary vertical depth, stratigraphic dip information, are calculated adjacent Stratum boundary meets personally vertical range, and effective interface judgement is carried out by the size of distance, two adjacent interfaces apart from it is too small when regard wherein one It is individual to be considered as pseudobed for inefficient frontier, according to activity value size, interface is screened and screened, reject pseudobed, obtain true layer Interface.
S103, according to horizontal well well logging data, offset well log data, the stratigraphic dip, AUTOMATIC ZONINGization point Bed boundary builds stratum initial model, and stratum initial model is initialized and parameter setting.
In embodiments of the present invention, with reference to AUTOMATIC ZONING produce sequence relation (or AUTOMATIC ZONING result), stratigraphic dip, Interface, apparent resistivity, well track vertical depth-tiltedly depth-horizontal displacement change, set up preliminary 3-dimensional multi-layered geological model, such as Fig. 2 institutes Show.According to apparent resistivity separation characteristic, the parameter such as initialization stratum background resistivity, slurry compounding, anisotropy;With according to well Bore parameters in the Given informations such as footpath, mud, initialization stratum initial model;According to stratigraphic dip, bed boundary, vertical depth initialization Well track and bed boundary location parameter in the initial model of stratum.
S104, calculated using finite element method three-dimensional horizontal well stratum initial model with boring orientation instrument forward response, Measured data and the forward response of stratum initial model are contrasted, object function is set up, carrying out multi-parameter with least square method changes The model knots modification for solving object function optimal solution for inverting obtains inversion result.
Alternatively, using the three-dimensional horizontal well stratum initial model of finite element method calculating with brill orientation instrument forward response Including S501-S504:
Maxwell's Maxwell equations under S501, the given boundary condition of solution, Maxwell is equations turned for fluctuation side Journey.
In embodiments of the present invention, finite element method calculate the stratum initial model response essence set up be attributed to solution to It is under the conditions of deckle circle the problem of Maxwell equations, Maxwell is equations turned for wave equation.
Alternatively, the wave equation can be:Wherein, E represents electric-field intensity; μ represents amplitude;ω represents angular frequency;ε represents dielectric constant.
S502, with reference to boundary condition wave equation is attributed to field energy functional.
Alternatively, the field energy functional can be:
Wherein,Represent object function;μ0Represent initial amplitude;V represents domain;ω represents angular frequency;ε is represented Dielectric constant;J represents current density.
S503, using finite elements subdivision field domain, and choose corresponding Interpolation-Radix-Function, energy functional is carried out space from Dissipate, obtain following discretization functional form:
Wherein, AeRepresentative function It is discretely formed coefficient matrix;BeRepresent equation right-hand vector discrete matrix;CeRepresentative function value discrete matrix;EeRepresent independent variable square Battle array;M represents independent variable number;The transposition of T representing matrixs.
S504, to discretization functional solve obtain inverting needed for parameter.
In embodiments of the present invention, three-dimensional horizontal well stratum introductory die is calculated using finite element method by above-mentioned steps Type with boring after orientation instrument forward response, it is necessary to contrast the forward response of measured data and stratum initial model, set up mesh Scalar functions, the model knots modification that multi-parameter iterative inversion solution object function optimal solution is carried out with least square method obtains inverting As a result.
Alternatively, the model that the multi-parameter iterative inversion solution object function optimal solution is carried out with least square method changes Variable, which obtains inversion result, includes S601-S605:
S601, the residual error by least square method solution actual measured value and the analogue value:
m≥n;Wherein, r:Rn→RmIt is x nonlinear function, Non-linear least square problem sees the special case of unconstrained minimization as;R represents independent variable vector;M represents unknown quantity Number;N represents forward response curve number.
S602, the Jacobian matrix for solving by gradient method the object function, form Jacobi's system of linear equations, ask Solve the direction of steepest descent that the system of linear equations calculates the object function.If J (x) is r (x) Jacobian matrix:
Based on above-mentioned Jacobian matrix, the gradient of the object function is calculated according to following 4th relational expression:
In embodiments of the present invention, the 4th relational expression is made to be equal to 0, you can to solve object function direction of steepest descent.
S603, using golden section approach, on the direction of steepest descent determine change step-length;
S604, change according to the model of determination direction and change step-length and determine model knots modification, and according to model knots modification Change model, complete an iteration inverting;
S605, according to default stopping criterion for iteration, recursive call iterative inversion is defeated until meet stopping criterion for iteration Go out result, so as to obtain structure parameters of soil layer.
S105, according to inversion result, stratigraphic sequence relation, pendage situation and measurement point to the distance on layer border, paint Well track processed-position of stratum graph of a relation, and coloured according to every layer of true formation resistivity inversion result, intuitively show stratum Occurrence information, and export inversion result.
Alternatively, according to the distance of inversion result, stratigraphic sequence relation, pendage situation and measurement point to layer border, Drawing well track-position of stratum graph of a relation includes S701-S705:
S701, by AUTOMATIC ZONING to stratum division result, definitely order relation layer by layer;
S702, the measurement point obtained by inverting turn to layer frontier distance (as shown in Figure 4), stratigraphic dip inversion result Stratum true dip angle, thickness, well track are changed through bed boundary position and well track and bed boundary relative position relation;
S703, the vertical continuity according to bed boundary and stratum true dip angle situation of change, estimate instrument investigative range it Outer bed boundary trend and thickness change, are carried out to non-inverting subterranean formation zone outside instrument investigative range according to existing inversion result Continue, to ensure the continuity and integrality of attitude of stratum figure;
S704, according to formation resistivity inversion result and Anisotropic inversion result, resistivity value is converted into chromatic value It is filled in corresponding stratum, and changes the change for qualitatively judging resistivity value by color light and shade;
S705, the attitude of stratum two dimensional image ultimately formed show that in the horizontal Different Strata has different resistivity value, Correspondence different colours, are changed by change in resistance, bed boundary change in location and thickness in color change marker;In longitudinal direction On, Fill Color change represents the change occurred in layer with horizontal displacement changes in resistance rate in layer, is used to refer to subterranean formation zone Successional quality.
The present embodiment is according to inversion result, to measurement point correspondence position by log interval to inversion result resampling, The square waveform result to Layer by layer inversion carries out non-linear interpolation simultaneously, obtains smooth result curve, is easy to draw map.
Theoretical foundation based on such scheme, is described below several specific embodiments of scheme of the embodiment of the present invention.
Embodiment one
Utilized in the horizontal well environment that the present embodiment is provided with the method for boring 3 D resistivity identification attitude of stratum, including step Rapid S11-S16:
S11, orientation natural gamma and azimuthal resistivity log response value obtained according to horizontal wellbore logging, using depth displacement and Curve matching extracts logger-stratum relative dip angle (i.e. stratigraphic dip) and formation azimuth angle.
S12, the angle of inclination and formation azimuth angle obtained according to step S11, to orientation resistivity logging and orientation resistance Rate log data carries out directional correction.
S13, to carry out directional correction after standardized data, ask for curve activity by relevant comparative methods, use activity Method carries out bed boundary division, and gamma ray curve, azimuthal resistivity log, orientation resistance are considered when bed boundary is divided Rate curve, neutron density curve and porosity curve.
The embodiment of the present invention carries out stratum AUTOMATIC ZONING using relevant comparative and activity method, comprises the following steps S21-S23:
The smooth filtering of S21, log data and invalid data are rejected;
S22, data normalization processing;
S23, relevant comparative and activity method AUTOMATIC ZONING.
The present embodiment need to be pre-processed before relevant comparative and AUTOMATIC ZONING is carried out to the data of measurement, pretreatment Content is filtered including line smoothing, removes invalid data and missing data is carried out to supplement to greatest extent in zone of reasonableness.Often Smooth filtering method has least square moving average method, weighted moving average etc..Wherein above two method is included again A variety of concrete implementation types:
Linear function is smooth:
Quadratic function is smooth:
Bell shaped function exponential smoothing:
Hamming function exponential smoothing:
Wherein, TiIth sample value is represented, m represents filter window,Represent ith sample point value after filtering.
Data after filtering and invalid data processing are normalized the present embodiment, and method for normalizing is:
Extreme value method for normalizing:
Porosity method for normalizing:
Density method for normalizing:
Wherein, xijRepresent the ith sample point value of j-th strip curve, xmin jRepresent the minimum value of j-th strip curve, xmax jTable Show the maximum of j-th strip curve, YijRepresent i-th of value of j-th strip curve after normalization;TNPLjRepresent porosity curve j-th Value, and replace original value with the value after normalization;ALCDjRepresent that neutron density is worth for j-th, and original is replaced with the value after normalizing Initial value.
The embodiment of the present invention carries out relevant comparative and activity method AUTOMATIC ZONING, including step S31-S33 to data:
S31, calculated level well well logging data curve activity.
Alternatively, curve activity is defined using the 3rd following relational expressions:
Wherein, E (d) represents d activity function value, and x (i) represents log measured value,Represent that log exists Average value in interval [d-n, d+n];D represents sampling point position;N represents that sample window is long;I is represented in a sample window length Several sampled points.
S32, with according to the curve activity tried to achieve in S31, setting activity threshold value, layering interfaces being determined according to activity threshold value.
S33, the bed boundary to step S32 divisions are screened and screened, and are rejected pseudobed, are obtained true bed boundary.
The present embodiment combination horizontal well bed boundary vertical depth, stratigraphic dip information, calculate neighboring layer interfaces see it is vertical away from From, effective interface judgement is carried out by the size of distance, two adjacent interfaces apart from it is too small when regard one of them and regarded as inefficient frontier For pseudobed, according to activity value size, interface is screened and screened, reject pseudobed, obtain true bed boundary.
S14, the sequence relation produced with reference to AUTOMATIC ZONING, inclination angle, interface, apparent resistivity, well track vertical depth-tiltedly deep- Horizontal displacement changes, and sets up preliminary three-dimensional horizontal well multilayer geological model.
In embodiments of the present invention, the layering interfaces that abovementioned steps are divided are needed to be screened and screened, reject pseudobed, Obtain true bed boundary.Specifically, the present embodiment combination horizontal well bed boundary vertical depth, stratigraphic dip information, are calculated adjacent Stratum boundary meets personally vertical range, and effective interface judgement is carried out by the size of distance, two adjacent interfaces apart from it is too small when regard wherein one It is individual to be considered as pseudobed for inefficient frontier, according to activity value size, interface is screened and screened, reject pseudobed, obtain true layer Interface.
S15, using finite element method calculate three-dimensional formation initial model forward response, pass through forward modeling response and survey Data comparison, sets up object function, and the model knots modification of object function optimal solution is solved with least square method, is iterated anti- Drill calculating.
The present embodiment calculates three-dimensional horizontal well stratum initial model forward response using finite element method;Utilize least square Method is iterated inverting.
Wherein, calculate three-dimensional formation initial model forward response using finite element method and include S41-S44:
Maxwell's Maxwell equations under S41, the given boundary condition of solution, Maxwell is equations turned for fluctuation side Journey.
In embodiments of the present invention, finite element method calculate the stratum initial model response essence set up be attributed to solution to It is under the conditions of deckle circle the problem of Maxwell equations, Maxwell is equations turned for wave equation.
Alternatively, the wave equation can be:Wherein, E represents electric-field intensity; μ represents amplitude;ω represents angular frequency;ε represents dielectric constant.
S42, with reference to boundary condition wave equation is attributed to field energy functional.
Alternatively, the field energy functional can be:
Wherein,Represent object function;μ0Represent initial amplitude;V represents domain;ω represents angular frequency;ε is represented Dielectric constant;J represents current density.
S43, using finite elements subdivision field domain, and choose corresponding Interpolation-Radix-Function, energy functional is carried out space from Dissipate, obtain following discretization functional form:
Wherein, AeRepresentative function It is discretely formed coefficient matrix;BeRepresent equation right-hand vector discrete matrix;CeRepresentative function value discrete matrix;EeRepresent independent variable square Battle array;M represents independent variable number;The transposition of T representing matrixs.
S44, to discretization functional solve obtain inverting needed for parameter.
In embodiments of the present invention, three-dimensional horizontal well stratum introductory die is calculated using finite element method by above-mentioned steps Type with boring after orientation instrument forward response, it is necessary to contrast the forward response of measured data and stratum initial model, set up mesh Scalar functions, the model knots modification that multi-parameter iterative inversion solution object function optimal solution is carried out with least square method obtains inverting As a result.
Alternatively, the model that the multi-parameter iterative inversion solution object function optimal solution is carried out with least square method changes Variable, which obtains inversion result, includes S51-S55:
S51, the residual error by least square method solution actual measured value and the analogue value:
m≥n;Wherein, r:Rn→RmIt is x nonlinear function, Non-linear least square problem sees the special case of unconstrained minimization as;R represents independent variable vector;M represents unknown quantity Number;N represents forward response curve number.
S52, the Jacobian matrix for solving by gradient method the object function, form Jacobi's system of linear equations, ask Solve the direction of steepest descent that the system of linear equations calculates the object function.If J (x) is r (x) Jacobian matrix:
Based on above-mentioned Jacobian matrix, the gradient of the object function is calculated according to following 4th relational expression:
In embodiments of the present invention, the 4th relational expression is made to be equal to 0, you can to solve object function direction of steepest descent.
S53, using golden section approach, on the direction of steepest descent determine change step-length;
S54, change according to the model of determination direction and change step-length and determine model knots modification, and changed according to model knots modification Varying model, completes an iteration inverting;
S55, according to default stopping criterion for iteration, recursive call iterative inversion, until meeting stopping criterion for iteration, is exported As a result, so as to obtain structure parameters of soil layer.
S16, the distance according to inversion result, stratigraphic sequence relation, pendage situation and measurement point to layer border, draw Well track-position of stratum graph of a relation, and coloured according to every layer of true formation resistivity inversion result, display stratum directly perceived production Shape information, and export inversion result.
Alternatively, according to the distance of inversion result, stratigraphic sequence relation, pendage situation and measurement point to layer border, Drawing well track-position of stratum graph of a relation includes S61-S65:
S61, by AUTOMATIC ZONING to stratum division result, definitely order relation layer by layer;
S62, the measurement point obtained by inverting to layer frontier distance, stratigraphic dip inversion result, inverting formation true dip angle, Thickness, well track pass through bed boundary position and well track and bed boundary relative position relation;
S63, the vertical continuity according to bed boundary and stratum true dip angle situation of change, are estimated outside instrument investigative range Bed boundary trend and thickness change, to non-inverting subterranean formation zone outside instrument investigative range according to inversion result prolong It is continuous, to ensure the continuity and integrality of attitude of stratum figure;
S64, according to formation resistivity inversion result and Anisotropic inversion result, resistivity value is converted into chromatic value and filled out Fill in corresponding stratum, and change the change for qualitatively judging resistivity value by color light and shade;
S65, the attitude of stratum two dimensional image ultimately formed show that in the horizontal Different Strata has different resistivity value, Correspondence different colours, are changed by change in resistance, bed boundary change in location and thickness in color change marker;In longitudinal direction On, Fill Color change represents the change occurred in layer with horizontal displacement changes in resistance rate in layer, is used to refer to subterranean formation zone Successional quality.
The present embodiment is according to inversion result, to measurement point correspondence position by log interval to inversion result resampling, The square waveform result to Layer by layer inversion carries out non-linear interpolation simultaneously, obtains smooth result curve, is easy to draw map.
Embodiment two
The present embodiment is provided can utilize azimuthal resistivity well logging under the bad hole environment such as horizontal well and high angle hole Carry out formation anisotropy and know method for distinguishing;S111-S113 is realized as follows:
S111, using log depth displacement and non-linear fitting method, extract stratigraphic dip and formation azimuth angle.
Alternatively, using log depth displacement and non-linear fitting method, stratigraphic dip and formation azimuth angle bag are extracted Include S211-S213:
S211, the four orientation natural curves measured using ABG loggers (azimuth gamma well-logging instrument) (orientation nature GR logging curve) stratigraphic dip extraction is carried out, relevant comparative's analysis is carried out to four orientation natural curves, respective layer is obtained Six points on depth displacement, as pendage aspect.Because the expanded view that the borehole wall intersects with dipping bed is shown as on image Monocyclic SIN function, meets the first expression formula:
Y=A sin (ω x- β)+y0=y0+A cosβsin(ωx)+(-A sinβ)cos(ωx);
Wherein, y represents the functional value that the borehole wall intersects with pendage bed boundary;A representative function value changes amplitudes;y0Represent Functional value y average;β representative function values y meets the initial phase of SIN function;ω representative function values y meets sinusoidal letter Several cycles.
S212, by nonlinear fitting solve the first expression formula in undetermined coefficient, utilize least square method obtain first Matrix equation.
Alternatively, the first matrix equation includes:Wherein, a0= y0, a1=A cos β, a2=-A sin β,a0、a1、a2Represent to intend The multinomial coefficient of conjunction;Represent independent variable value during fitting function.
S213, to the first matrix equation solve, obtain sinusoids parameter expression formula and inclination angle tendency expression formula;It is sinusoidal bent Line parameter expression is:It is inclined to expression formula in inclination angle:Wherein, Dip tables Show stratigraphic dip, Dir represents to be inclined to (or formation azimuth angle), DelRepresent investigation depth, yminRepresent what the borehole wall intersected with dipping bed Minimum point on one week.
Alternatively, using log depth displacement and non-linear fitting method, stratigraphic dip and formation azimuth angle are extracted also Including S311-S313:
S311, three using other horizontal wellbore logging apparatus measures or four band azimuth information curves (orientation nature gal Horse, azimuthal resistivity well logging, azimuthal resistivity, Orientation density, orientation porosity etc.) stratigraphic dip extraction is carried out, it is bent to orientation Line carries out relevant comparative's analysis, can obtain the depth displacement of respective layer, as three in pendage aspect and more than three Point.Because the expanded view that the borehole wall intersects with dipping bed shows as monocyclic SIN function on image, equation is met:
Y=A sin (ω x- β)+y0Or y=A cos (ω x- β)-y0
S312, with analytic method or non-linear fitting method, solve the undetermined coefficient in equation in S311, utilize a most young waiter in a wineshop or an inn Multiplication obtains matrix equation:
Wherein, a0=y0, a1=A cos β, a2=-A sin β, a0、a1、a2Represent the multinomial coefficient of fitting;Represent that independent variable takes during fitting function Value.
S313, the matrix equation solved in S312, can obtain sinusoids parameter expression formula:
Or
Expression formula is inclined to inclination angle:
Or
Wherein, Dip represents stratigraphic dip, and Dir represents to be inclined to (or formation azimuth angle), DelRepresent investigation depth, yminRepresent Minimum point on the borehole wall intersect with dipping bed one week.
S112, the stratigraphic dip using step S111 acquisitions and formation azimuth angle, to orientation resistivity logging and orientation electricity Resistance rate log data carries out directional correction, and azimuth information meets sinusoidal (or cosine) rule, and the cycle is 2, can set curvilinear equation For:
Y=A sin (x+ φ)+B or y=A cos (x+ φ) 6B;
It is known in x=0,π,Value, by analytic solutions ask A, B andAnalytic solutions are as follows:
Or
Try to achieve after φ,A, B value, take function y=A sin (x+ φ)+B or y=A cos (x + φ)+B in x=0,π,The value at place is worth as after directional correction.
S113, carry out stratum AUTOMATIC ZONING using relevant comparative and activity method, including the smooth filtering of log data and invalid Data are rejected, and are filtered using the curve-fitting method more than 3 points, or include linear interpolation, non-linear using interpolation method Interpolation, area interpolation;Data normalization processing is using index method for normalizing, linear normalization method;Relevant comparative and activity Layering or inflection method layering.
S114, thin layer number is determined, rejects the pseudobed that AUTOMATIC ZONING is produced.
S115, the stratigraphic dip extracted according to step S111, the individual-layer data that step S113 is obtained, with reference to apparent resistivity structure Build stratum initial model.
S115, the azimuthal resistivity well logging forward response for calculating the initial stratum initial model built, with measured curve pair Than carrying out multi-parameter iterative inversion with least square method;The stratum initial model response set up is calculated using finite element method.
Orientation is carried out with the numerical simulation for boring electromagnetic wave instrument using finite element method, have studied the azimuthal resistivity of instrument The relation of log response same level well formation anisotropy, result of study shows that azimuthal resistivity well logging has very well to anisotropy Indicative function, in anisotropic formation, Amplitude Ratio resistivity is less than phase difference resistivity, and increases with anisotropy coefficient, The two difference increase;With change of pitch angle, orientation connector for logging while drilling log response separation, inclination angle is bigger, and curve separation is bigger. Now, comprehensive is effective using formation resistivity and azimuthal resistivity well logging progress formation anisotropy identification.Therefore, this hair The bright azimuthal resistivity well logging for proposing to be obtained using apparatus measures and resistivity curve are carried out formation anisotropy and know method for distinguishing, And invasion profile figure is provided, it is respectively provided with significance in being evaluated in geosteering and after boring.The present invention is not increasing existing well logging Truly believe on the extraction of complete paired parallel horizontal wells formation anisotropy on the premise of technical costs, more effective reaction stratum Breath, perfect existing level borehole logging tool evaluates data;The present invention is used under the horizontal well ground environment that 3-d inversion is obtained simultaneously Anisotropically layer resistivity and invasion profile, more press close to true stratum condition, obtain comprehensive formation parameter information, carry out more It is that authentic and valid formation evaluation and reserves are calculated.
Embodiment three
As illustrated, Fig. 1 produces to be utilized in the horizontal well environment of the embodiment of the present invention with brill 3 D resistivity identification stratum The flow chart of the method specific embodiment of shape;Fig. 2 is TVD tables in the horizontal well stratum modeling schematic diagram of the embodiment of the present invention, figure Show a layer boundary position, RM represents mud resistivity, and RI represents depth of invasion, and RXO represents invaded zone resistivity, and RS represents country rock electricity Resistance rate, Rh-Rv represents purpose horizontal resistivity-vertical resistivity, and DH represents borehole diameter;Fig. 3 is the water of the embodiment of the present invention Orientation resistivity logging while drilling logger direction signal changes with instrument-bed boundary relative inclination under horizontal well, and DIP represents instrument axis Direction and bed boundary normal direction corner dimension;Fig. 4 is the present invention with the phase difference resistivity and amplitude ratio resistance for boring instrument Rate with measurement point to layer frontier distance change, after wherein ARC48PC represents under 2MHz frequencies 48in sources away from borehole environment correction Phase difference apparent resistivity, ARC32PC represents that 32in sources are away from phase difference apparent resistivity after borehole environment correction under 2MHz frequencies, ARC16PC represents that 16in sources are away from phase difference apparent resistivity after borehole environment correction under 2MHz frequencies, and ARC48AC represents 2MHz frequencies Lower 48in sources correct amplitude than apparent resistivity away from borehole environment, and ARC32AC represents that 32in sources are away from borehole environment under 2MHz frequencies Amplitude is corrected than apparent resistivity, ARC16AC represents that apparent resistance is compared in 16in sources away from borehole environment correction amplitude under 2MHz frequencies Rate;Fig. 5 is apparent resistivity in the anisotropic formation of the embodiment of the present invention with change of pitch angle;Fig. 6 is right in the embodiment of the present invention The formation anisotropy inversion result figure of certain horizontal well reality well-log information, wherein UP is reservoir coboundary, and DOWN is that reservoir is following Boundary, TVD is well vertical depth, and BGRC is average natural gamma, and BCGT is 0 ° of orientation gamma, and BCGB is 180 ° of orientation gammas, ZRH48PT be 48in sources away from lower 0 ° of azimuthal resistivity, ZRH48PT is 48in sources away from lower 180 ° of azimuthal resistivities, other resistivity Class of a curve is same, and ALCDLC is neutron density curve, and TNPL is neutron porosity.
The embodiment of the present invention calculate when, fully take into account well in horizontal well, high angle hole three dimensions complex environment, Mud, inclination angle, intrusion, anisotropy, measurement point background resistivity etc. into layer frontier distance, thickness and each layer are multifactor Influence, and combine that other log datas carry out stratum modelings and boundary division, inclination angle are extracted, fully taken into account in calculating process Each influence factor coupling, closer in practical logging ambient conditions, obtained result more presses close to actual value.Solve simultaneously Traditional inversion method of having determined can not carry out the inverting of asymmetric stratum and large dip angle strata;Solving conventional method can not be definitely Layer orientation, also can not accurately calculate measurement point to layer frontier distance, so that the problem of attitude of stratum is described can not be carried out;The opposing party Face, the identification information of formation anisotropy is supplied by 3 D resistivity inverting, horizontal resistivity and vertical resistivity is obtained, in reality There is very important application value in the production application of border.
Above example is only one embodiment of the present invention, and it describes more specific and in detail, but can not therefore and It is interpreted as the limitation to the scope of the claims of the present invention.Its concrete structure and size can be adjusted correspondingly according to actual needs.Should When, it is noted that for the person of ordinary skill of the art, without departing from the inventive concept of the premise, can also make Several modifications and improvements, these belong to protection scope of the present invention.
In order to achieve the above object, the invention also provides a kind of device 1 of the identification attitude of stratum in horizontal well environment, As shown in fig. 7, the device includes:Extraction module 11, processing module 12, modeling module 13, computing module 14 and display module 15;
Extraction module 11, for the orientation Natural Gamma-ray Logging Curves or Orientation density obtained according to horizontal well well logging Log extracts stratigraphic dip and formation azimuth angle;
Processing module 12, for obtain multiple sources away from the horizontal well well logging data under multiple frequencies;Utilize extraction The stratigraphic dip and formation azimuth angle the horizontal well well logging data are pre-processed, the well logging standardized Data;Log data to the standardization carries out curve relevant comparative, and carries out stratum AUTOMATIC ZONING using activity method, obtains Bed boundary;Wherein, the horizontal well well logging data include instrumental azimuth resistivity log response value, density and porosity Curve;
Modeling module 13, for according to the horizontal well well logging data, offset well log data, the stratigraphic dip, The bed boundary of AUTOMATIC ZONINGization point builds stratum initial model, and the stratum initial model is initialized and parameter is set Put;
Computing module 14, for the three-dimensional horizontal well stratum initial model of utilization finite element method calculating with brill orientation instrument Forward response, contrast measured data and the forward response of stratum initial model, set up object function, are carried out with least square method The model knots modification that multi-parameter iterative inversion solves the object function optimal solution obtains inversion result;
Display module 15, for according to the inversion result, stratigraphic sequence relation, the pendage situation and measurement point To the distance on layer border, well track-position of stratum graph of a relation is drawn, and carry out according to every layer of true formation resistivity inversion result Coloring, display attitude of stratum information directly perceived, and export inversion result.
Compared with prior art, the present invention includes:It is bent according to the orientation gamma ray log that horizontal well well logging is obtained Line or Orientation density log extract stratigraphic dip and formation azimuth angle;Obtain multiple sources away from the horizontal well under multiple frequencies Well logging data;Horizontal well well logging data are pre-processed using the stratigraphic dip and formation azimuth angle of extraction, obtained To the log data of standardization;Log data to standardization carries out curve relevant comparative, and carries out stratum certainly using activity method Dynamic layering, obtains bed boundary;Wherein, horizontal well well logging data include instrumental azimuth resistivity log response value, density And porosity curve;Divided according to horizontal well well logging data, offset well log data, the stratigraphic dip, AUTOMATIC ZONINGization Bed boundary builds stratum initial model, and stratum initial model is initialized and parameter setting;Utilize finite element method Calculate three-dimensional horizontal well stratum initial model with boring orientation instrument forward response, contrast measured data and stratum initial model Forward response, sets up object function, and the mould that multi-parameter iterative inversion solves object function optimal solution is carried out with least square method Type knots modification obtains inversion result;According to inversion result, stratigraphic sequence relation, pendage situation and measurement point to layer border Distance, draws well track-position of stratum graph of a relation, and is coloured according to every layer of true formation resistivity inversion result, intuitively Attitude of stratum information is shown, and exports inversion result.By the solution of the present invention, accurate formation parameter is obtained and each to different Property information, measurement point to layer frontier distance information, perfect Horizontal Well Log Interpretation data currently, is that later stage well log interpretation is commented Valency, saturation computation and reserves, which are calculated, provides accurate reliable parameter.
The azimuth gamma well-logging curve that scheme of the embodiment of the present invention is obtained according to well logging first extracts formation azimuth angle And stratigraphic dip;Then multiple sources in obtaining horizontal well of being logged well from resistivity logging while drilling are bent away from the resistivity logging under, multi-frequency Line and direction signal log, pass through data prediction and AUTOMATIC ZONING, the Resistivity log standardized, stratum Layering interfaces data;Result, strata division interface data, resistivity and direction signal are extracted finally according to dip azimuth angle Measured value and offset well log data, build the horizontal well stratigraphic model under three dimensions, and according to log response and model parameter Actual physical meaning, appropriate constraint is carried out to model parameter span;By the three-dimensional Model of Horizontal Well of Inversion Calculation, obtain Obtain attitude of stratum relevant parameter, including bed boundary position, stratigraphic dip, formation resistivity, anisotropy coefficient, measurement point to layer Frontier distance;According to inversion result and stratigraphic sequence relation, pendage situation, measurement point to layer frontier distance, well is drawn Track-position of stratum graph of a relation, and coloured according to every layer of true formation resistivity inversion result, display attitude of stratum letter directly perceived Breath.This method using existing well logging data extract real formation electrical parameter under horizontal well environment (formation resistivity, Resistivity anisotropy information) and sectional parameter (stratigraphic dip, bed boundary, measurement point to layer frontier distance), do not increasing ground In the case of shaking the well logging projects such as well logging, sound wave, still attitude of stratum information can be obtained so that well-log information is able to abundant profit With.Meanwhile, this method results in accurate formation parameter and anisotropy information, measurement point to layer frontier distance information, complete It has been apt to existing Horizontal Well Log Interpretation data, has been evaluated for later stage well log interpretation, saturation computation and reserves calculate and provide accurate Reliable parameter.
It should be noted that embodiment described above be for only for ease of it will be understood by those skilled in the art that, and It is not used in and limits the scope of the invention, on the premise of the inventive concept of the present invention is not departed from, those skilled in the art couple Any obvious replacement and improvement that the present invention is made etc. is within protection scope of the present invention.

Claims (14)

1. the method for the identification attitude of stratum in a kind of horizontal well environment, it is characterised in that methods described includes:
The orientation Natural Gamma-ray Logging Curves or Orientation density log obtained according to horizontal well well logging extract stratum and inclined Angle and formation azimuth angle;
Obtain multiple sources away from the horizontal well well logging data under multiple frequencies;Utilize the stratigraphic dip of extraction and stratum Azimuth is pre-processed to the horizontal well well logging data, the log data standardized;To the standardization Log data carries out curve relevant comparative, and carries out stratum AUTOMATIC ZONING using activity method, obtains bed boundary;Wherein, it is described Horizontal well well logging data include instrumental azimuth resistivity log response value, density and porosity curve;
The ground stratum boundary divided according to the horizontal well well logging data, offset well log data, the stratigraphic dip, AUTOMATIC ZONINGization Face builds stratum initial model, and the stratum initial model is initialized and parameter setting;
Using the three-dimensional horizontal well stratum initial model of finite element method calculating with orientation instrument forward response is bored, number is surveyed in contrast According to the forward response with stratum initial model, object function is set up, multi-parameter iterative inversion solution is carried out with least square method The model knots modification of the object function optimal solution obtains inversion result;
According to the distance of the inversion result, stratigraphic sequence relation, the pendage situation and measurement point to layer border, draw Well track-position of stratum graph of a relation, and coloured according to every layer of true formation resistivity inversion result, display stratum directly perceived production Shape information, and export inversion result.
2. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 1, it is characterised in that described according to water The orientation Natural Gamma-ray Logging Curves or Orientation density log that horizontal well well logging is obtained extract stratigraphic dip and stratum side Parallactic angle refers to:According to depth displacement and the nonlinear fitting side of the orientation Natural Gamma-ray Logging Curves or Orientation density log Method extracts the stratigraphic dip and formation azimuth angle;
Depth displacement and the nonlinear fitting side according to the orientation Natural Gamma-ray Logging Curves or Orientation density log Method, which extracts the stratigraphic dip and formation azimuth angle, to be included:
Every orientation gamma ray log in a plurality of orientation Natural Gamma-ray Logging Curves measured orientation gamma-ray tool device Curve carries out relevant comparative's analysis, obtains the depth displacement on correspondence stratum, obtains the measurement point on pendage bed boundary;And pass through First expression formula represents expanded view that the borehole wall intersects with the pendage bed boundary;First expression formula includes:Y=A sin(ωx-β)+y0=y0+A cosβsin(ωx)+(-A sinβ)cos(ωx);Wherein, y represents the borehole wall and pendage layer The functional value of interface intersection;A representative function value changes amplitudes;y0Representative function value y average;β representative function values y is met just The initial phase of string function;ω representative function values y meets the cycle of SIN function;
Undetermined coefficient in first expression formula is solved by nonlinear fitting, the first matrix side is obtained using least square method Journey;First matrix equation includes:Wherein, a0=y0, a1=A Cos β, a2=-A sin β,a0、a1、a2Represent the multinomial of fitting Formula coefficient;Represent independent variable value during fitting function;
First matrix equation is solved, sinusoids parameter expression formula and inclination angle tendency expression formula is obtained;It is described sinusoidal bent Line parameter expression is:It is inclined to expression formula in the inclination angle:Wherein, Dip represents the stratigraphic dip, and Dir represents tendency, DelRepresent investigation depth, yminRepresent one week that the borehole wall intersects with dipping bed On minimum point.
3. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 1, it is characterised in that described utilize carries The stratigraphic dip and formation azimuth angle taken carries out pretreatment to the horizontal well well logging data to be included:
It is electric to the orientation of the azimuth information in the horizontal well well logging data using the stratigraphic dip and formation azimuth angle Resistance rate log data carries out directional correction and the disposal of gentle filter, and carries out borehole environment school to the data after the smothing filtering Just.
4. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 3, it is characterised in that the utilization institute State formation azimuth angle includes to the azimuthal resistivity log data progress directional correction in the horizontal well well logging data:
The orientation curve in horizontal well well logging is fitted, and is expressed with first orientation curvilinear function;Wherein, the first orientation Curvilinear function includes:Y (x)=A sin (x+ φ)+B;X is the formation azimuth angle;Y (x) is the first orientation curvilinear function Functional value;B is the y of first orientation curvilinear function average value;A is shaking when the first orientation curve meets SIN function Width;φ is initial phase;
Multiple characteristic values on the first orientation curvilinear function are chosen, and the multiple characteristic value is substituted into the first orientation Curvilinear function, asks for described A, B and φ;
Described A, B and φ for obtaining are substituted into the first orientation curvilinear function, the first equation with known coefficient is obtained;
Using the corresponding functional value of first equation as the azimuthal resistivity log data directional correction value.
5. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 4, it is characterised in that the multiple spy Value indicative includes:(0, y (0)), (1, y (1)), (2, y (2)), (3, y (3));Described A, B and φ value is respectively:
φ=atan (y (0)-B)/(y (2)-B);
A=(y (3)-y (1))/(2*sin (φ)).
6. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 3, it is characterised in that the smooth filter The method of ripple includes:Least square moving average method and weighted moving average.
7. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 3, it is characterised in that described to described Data after smothing filtering, which carry out borehole environment correction, to be included:
According to log data apparent resistivity and borehole environment influence factor, the school inquired about in the borehole correction charts storehouse pre-established Positive coefficient, borehole environment correction is carried out according to the correction coefficient to the log data apparent resistivity;
Wherein, pre-establishing the borehole correction charts storehouse includes:
The fixed borehole environment influence factor, changes stratum background resistivity, and borehole correction is calculated using forward modeling simulated program Plate;The borehole environment influence factor includes borehole size and drilling fluid resistivity;
The borehole correction charts in the case of different borehole size and drilling fluid resistivity are calculated, and set up borehole correction charts storehouse.
8. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 1, it is characterised in that described to described The log data of standardization, which carries out curve relevant comparative, to be included:To in the horizontal well well logging data after smothing filtering not Log with investigation depth is normalized, and correlation contrast is carried out under fixed window length to the data after normalization;
Wherein, the relative coefficient in the correlation contrast is calculated according to following second relational expressions:
2
In formula, i represents sampling point position, and n represents that window is long,Contrast bar plot mean in window length is represented,Represent by contrast bar Plot mean, xiRepresent correlation curve ith sample point value, yiRepresent by correlation curve ith sample point curve value, ZiTable Show coefficient correlation;
It is described to be included using activity method progress stratum AUTOMATIC ZONING:
Calculate the curve activity of the horizontal well well logging data;
Judge whether the curve activity is more than default activity threshold value, if the curve activity is more than the default activity Threshold value, the then position that will be greater than the extreme point of the activity threshold value is defined as layering interfaces.
9. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 8, it is characterised in that described to smooth In the filtered horizontal well well logging data log of different investigation depths be normalized including:Extreme value normalizing Change, porosity normalization and density are normalized;
Wherein, resistivity sampling extreme value method for normalizing is included:
Neutron porosity sampling linear method for normalizing is included:
Neutron density sampling linear method for normalizing is included:
Wherein, xijRepresent the ith sample point value of j-th strip curve, xminjRepresent the minimum value of j-th strip curve, xmaxjRepresent jth The maximum of bar curve, YijRepresent i-th of value of j-th strip curve after normalization;TNPLjRepresent that porosity curve is worth for j-th, and Original value is replaced with the value after normalization;ALCDjRepresent that neutron density is worth for j-th, and original value is replaced with the value after normalizing.
10. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 8, it is characterised in that the curve Activity is defined using the 3rd following relational expressions:
Wherein, E (d) represents d activity function value, and x (i) represents log measured value,Represent log in interval Average value in [d-n, d+n];D represents sampling point position;N represents that sample window is long;I represents which is individual in a sample window length Sampled point.
11. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 8, it is characterised in that the utilization Finite element method calculates including with brill orientation instrument forward response for three-dimensional horizontal well stratum initial model:
Maxwell's Maxwell equations under given boundary condition are solved, the Maxwell is equations turned for wave equation:Wherein, E represents electric-field intensity;μ represents amplitude;ω represents angular frequency;ε represents dielectric Constant;
The wave equation is attributed to field energy functional with reference to the boundary condition:
Wherein,Represent Object function;μ0Represent initial amplitude;V represents domain;ω represents angular frequency;ε represents dielectric constant;J represents that electric current is close Degree;
Using finite elements subdivision field domain, and corresponding Interpolation-Radix-Function is chosen, spatial spreading is carried out to the energy functional, obtained To following discretization functional form:
Wherein, AeRepresentative function is discrete The efficiency of formation matrix;BeRepresent equation right-hand vector discrete matrix;CeRepresentative function value discrete matrix;EeRepresent independent variable matrix;M Represent independent variable number;The transposition of T representing matrixs;
The parameter needed for acquisition inverting is solved to the discretization functional.
12. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 8, it is characterised in that described with most The model knots modification that small least square method carries out the multi-parameter iterative inversion solution object function optimal solution obtains inversion result bag Include:
The residual error of actual measured value and the analogue value is solved by least square method:
Wherein, r:Rn→RmIt is x nonlinear function;R is represented Independent variable vector;M represents unknown quantity number;N represents forward response curve number;
The Jacobian matrix of the object function is solved by gradient method, Jacobi's system of linear equations is formed, the line is solved Property equation group calculate the direction of steepest descent of the object function;Wherein, the target is calculated according to following 4th relational expression The gradient of function:J (x) is r (x) Jacobian matrix;
Using golden section approach, determine to change step-length on the direction of steepest descent;
Direction is changed according to the model of determination and changes step-length and determines model knots modification, and mould is changed according to the model knots modification Type, completes an iteration inverting;
According to default stopping criterion for iteration, recursive call iterative inversion, until the stopping criterion for iteration is met, output knot Really.
13. the method for the identification attitude of stratum in horizontal well environment as claimed in claim 8, it is characterised in that the basis The inversion result, stratigraphic sequence relation, the pendage situation and measurement point draw well rail to the distance on layer border Mark-position of stratum graph of a relation includes:
By AUTOMATIC ZONING to stratum division result, definitely order relation layer by layer;
The measurement point obtained by inverting is to layer frontier distance, stratigraphic dip inversion result, inverting formation true dip angle, thickness, well Eye track passes through bed boundary position and well track and bed boundary relative position relation;
According to the vertical continuity of bed boundary and stratum true dip angle situation of change, the ground stratum boundary outside instrument investigative range is estimated Face is moved towards and thickness change, and non-inverting subterranean formation zone outside instrument investigative range is continued according to existing inversion result, to protect Demonstrate,prove the continuity and integrality of attitude of stratum figure;
According to formation resistivity inversion result and Anisotropic inversion result, resistivity value is converted into chromatic value and is filled in correspondence Stratum in, and pass through color light and shade change qualitatively judge resistivity value change;
The attitude of stratum two dimensional image ultimately formed shows that in the horizontal Different Strata has different resistivity value, and correspondence is different Color, is changed by change in resistance, bed boundary change in location and thickness in color change marker;In the vertical, in layer Fill Color change represents the change occurred in layer with horizontal displacement changes in resistance rate, is used to refer to subterranean formation zone successional Quality.
14. the device of the identification attitude of stratum in a kind of horizontal well environment, it is characterised in that described device includes:Extraction module, Processing module, modeling module, computing module and display module;
Extraction module, orientation Natural Gamma-ray Logging Curves or Orientation density for being obtained according to horizontal well well logging are logged well bent Line drawing stratigraphic dip and formation azimuth angle;
Processing module, for obtain multiple sources away from the horizontal well well logging data under multiple frequencies;Using described in extraction Stratigraphic dip and formation azimuth angle are pre-processed to the horizontal well well logging data, the log data standardized; Log data to the standardization carries out curve relevant comparative, and carries out stratum AUTOMATIC ZONING using activity method, obtains stratum Interface;Wherein, the horizontal well well logging data are write music including instrumental azimuth resistivity log response value, density and hole Line;
Modeling module, for according to the horizontal well well logging data, offset well log data, the stratigraphic dip, automatic point The bed boundary of stratification point builds stratum initial model, and the stratum initial model is initialized and parameter setting;
Computing module, for calculating being rung with brill orientation instrument forward modeling for three-dimensional horizontal well stratum initial model using finite element method Should, contrast measured data and the forward response of stratum initial model set up object function, multi-parameter are carried out with least square method The model knots modification that iterative inversion solves the object function optimal solution obtains inversion result;
Display module, for according to the inversion result, stratigraphic sequence relation, the pendage situation and measurement point to layer side The distance on boundary, draws well track-position of stratum graph of a relation, and is coloured according to every layer of true formation resistivity inversion result, Display attitude of stratum information directly perceived, and export inversion result.
CN201710069496.1A 2017-02-08 2017-02-08 A kind of method and apparatus of identification attitude of stratum in horizontal well environment Pending CN107045154A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710069496.1A CN107045154A (en) 2017-02-08 2017-02-08 A kind of method and apparatus of identification attitude of stratum in horizontal well environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710069496.1A CN107045154A (en) 2017-02-08 2017-02-08 A kind of method and apparatus of identification attitude of stratum in horizontal well environment

Publications (1)

Publication Number Publication Date
CN107045154A true CN107045154A (en) 2017-08-15

Family

ID=59544534

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710069496.1A Pending CN107045154A (en) 2017-02-08 2017-02-08 A kind of method and apparatus of identification attitude of stratum in horizontal well environment

Country Status (1)

Country Link
CN (1) CN107045154A (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107765318A (en) * 2017-10-10 2018-03-06 中国石油天然气集团公司 A kind of method and device for determining geologic horizon
CN107861917A (en) * 2017-11-29 2018-03-30 中国石油集团长城钻探工程有限公司 Method of the multivariate data combined calculation well to stratigraphic boundary distance in horizontal well
CN107942400A (en) * 2017-11-10 2018-04-20 中国石油天然气股份有限公司 A kind of method and device for predicting drilling risk
CN108073765A (en) * 2017-12-11 2018-05-25 中海石油(中国)有限公司湛江分公司 A kind of horizontal well routine well logging bed boundary identification is with border away from inversion method
CN108345719A (en) * 2018-01-17 2018-07-31 中国石油大学(华东) A kind of high angle hole/horizontal well is with boring inversion method after electromagnetic wave logging data is bored
CN108896004A (en) * 2018-08-01 2018-11-27 刘敬寿 A kind of fracture surface roughness anisotropy characterizing method
CN108915676A (en) * 2018-07-20 2018-11-30 陕西延长石油(集团)有限责任公司研究院 A kind of compact reservoir hole movable fluid invasion profile imaging method
CN108979629A (en) * 2018-07-10 2018-12-11 中国海洋石油集团有限公司 It is a kind of based on density, neutron, nuclear magnetic resonance log gas-bearing formation density calculation method
CN109901239A (en) * 2019-03-08 2019-06-18 中国石油天然气集团有限公司 The measurement method and system of density of earth formations
CN110552689A (en) * 2018-05-15 2019-12-10 中国石油化工股份有限公司 method for determining distance from while-drilling instrument to stratum boundary

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345002A (en) * 2013-07-09 2013-10-09 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 Method for obtaining sandstone brine layer potassium ion content by means of logging information
CN106324689A (en) * 2016-06-24 2017-01-11 杭州迅美科技有限公司 Resistivity anisotropy recognition method in horizontal well stratum environment

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103345002A (en) * 2013-07-09 2013-10-09 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 Method for obtaining sandstone brine layer potassium ion content by means of logging information
CN106324689A (en) * 2016-06-24 2017-01-11 杭州迅美科技有限公司 Resistivity anisotropy recognition method in horizontal well stratum environment

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
汪中浩 等: "水平井地层电阻率各向异性研究及应用", 《石油物探》 *
胡松 等: "水平井随钻电磁波响应特征分析及应用实例", 《地质科技情报》 *
许泽瑞 等: "随钻电阻率测井响应分析与应用", 《岩性油气藏》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107765318A (en) * 2017-10-10 2018-03-06 中国石油天然气集团公司 A kind of method and device for determining geologic horizon
CN107765318B (en) * 2017-10-10 2019-06-11 中国石油天然气集团公司 A kind of method and device of determining geologic horizon
CN107942400A (en) * 2017-11-10 2018-04-20 中国石油天然气股份有限公司 A kind of method and device for predicting drilling risk
CN107942400B (en) * 2017-11-10 2019-06-11 中国石油天然气股份有限公司 A kind of method and device for predicting drilling risk
CN107861917B (en) * 2017-11-29 2018-11-02 中国石油集团长城钻探工程有限公司 Method of the multivariate data combined calculation wellbore to stratigraphic boundary distance in horizontal well
CN107861917A (en) * 2017-11-29 2018-03-30 中国石油集团长城钻探工程有限公司 Method of the multivariate data combined calculation well to stratigraphic boundary distance in horizontal well
CN108073765A (en) * 2017-12-11 2018-05-25 中海石油(中国)有限公司湛江分公司 A kind of horizontal well routine well logging bed boundary identification is with border away from inversion method
CN108345719A (en) * 2018-01-17 2018-07-31 中国石油大学(华东) A kind of high angle hole/horizontal well is with boring inversion method after electromagnetic wave logging data is bored
CN108345719B (en) * 2018-01-17 2021-09-17 中国石油大学(华东) Post-drilling inversion method for electromagnetic wave logging data while drilling of highly deviated well/horizontal well
CN110552689A (en) * 2018-05-15 2019-12-10 中国石油化工股份有限公司 method for determining distance from while-drilling instrument to stratum boundary
CN108979629A (en) * 2018-07-10 2018-12-11 中国海洋石油集团有限公司 It is a kind of based on density, neutron, nuclear magnetic resonance log gas-bearing formation density calculation method
CN108915676B (en) * 2018-07-20 2021-09-14 陕西延长石油(集团)有限责任公司研究院 Movable fluid invasion profile imaging method for tight reservoir pores
CN108915676A (en) * 2018-07-20 2018-11-30 陕西延长石油(集团)有限责任公司研究院 A kind of compact reservoir hole movable fluid invasion profile imaging method
CN108896004A (en) * 2018-08-01 2018-11-27 刘敬寿 A kind of fracture surface roughness anisotropy characterizing method
CN108896004B (en) * 2018-08-01 2020-03-20 中国石油大学(华东) Crack surface roughness anisotropy characterization method
CN109901239A (en) * 2019-03-08 2019-06-18 中国石油天然气集团有限公司 The measurement method and system of density of earth formations

Similar Documents

Publication Publication Date Title
CN107045154A (en) A kind of method and apparatus of identification attitude of stratum in horizontal well environment
CN106324689B (en) Resistivity anisotropy recognition methods under a kind of horizontal well ground environment
CN105938503B (en) A kind of direction signal multilayer interface identification method
Lai et al. A review on the applications of image logs in structural analysis and sedimentary characterization
Ahr et al. Confronting the carbonate conundrum
CN101438176A (en) Time lapse analysis with electromagnetic data
CN105093313B (en) A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device
CN104516018B (en) Porosity inversion method under lithological constraint in geophysical exploration
Poppelreiter et al. Borehole image log technology: application across the exploration and production life cycle
CN105005077A (en) Thin layer thickness prediction method with combination of real drilling wells and virtual wells under rare well condition
CN106951660A (en) A kind of marine clastics horizontal well reservoir log interpretation method and device
MX2014012042A (en) Formation volumetric evaluation using normalized differential data.
CN105074505A (en) Determination of true formation resistivity
Hurley Recognition of faults, unconformities, and sequence boundaries using cumulative dip plots
CN108073765A (en) A kind of horizontal well routine well logging bed boundary identification is with border away from inversion method
Mazaheri et al. Developing fracture measure as an index of fracture impact on well-logs
CN110685600B (en) Drill bit adjustment prediction method for geosteering
US20180238148A1 (en) Method For Computing Lithofacies Probability Using Lithology Proximity Models
CN107861917B (en) Method of the multivariate data combined calculation wellbore to stratigraphic boundary distance in horizontal well
Hurley AAPG Methods in Exploration, No. 16, Chapter 9: Borehole Images
CN107589469B (en) The determination method and apparatus of oil-water interfaces
CN108804728A (en) Horizontal well stratum reservoir grading analysis method and computer readable storage medium
CN105164729A (en) Image symmetry for dip determination
Nasseri et al. Evaluating Bangestan reservoirs and targeting productive zones in Dezful embayment of Iran
Poppelreiter et al. Applications of dipmeter and borehole image data in static models

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing

Applicant after: China Offshore Oil Group Co., Ltd.

Applicant after: China Oilfield Services Limited

Address before: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing

Applicant before: China National Offshore Oil Corporation

Applicant before: China Oilfield Services Limited

Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing

Applicant after: China Offshore Oil Group Co., Ltd.

Applicant after: China Oilfield Services Limited

Address before: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing

Applicant before: China National Offshore Oil Corporation

Applicant before: China Oilfield Services Limited

CB02 Change of applicant information
RJ01 Rejection of invention patent application after publication

Application publication date: 20170815

RJ01 Rejection of invention patent application after publication