CN108693562A - A method of with electric imaging logging measure paleocurrent to - Google Patents
A method of with electric imaging logging measure paleocurrent to Download PDFInfo
- Publication number
- CN108693562A CN108693562A CN201810465781.XA CN201810465781A CN108693562A CN 108693562 A CN108693562 A CN 108693562A CN 201810465781 A CN201810465781 A CN 201810465781A CN 108693562 A CN108693562 A CN 108693562A
- Authority
- CN
- China
- Prior art keywords
- data
- paleocurrent
- sector
- pole plate
- conductivity
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention belongs to geologic structure determination techniques fields, be related to it is a kind of with electric imaging logging measure paleocurrent to method, using electric imaging logging analysis determine downhole in reservoir paleocurrent to method, first carry out the processing of electric imaging logging data, data are carried out again to divide sector count, then according to statistical disposition chart measure paleocurrent to;Its processing step include determine sandstone reservoir type, processing imaging data, generate Electrical imaging image, deduct shale layer, divide sector count, draw rose figure, determine paleocurrent to verify result of calculation;Paleocurrent direction is determined using electric imaging logging data, a large amount of geologic features are identified and picked up by manual method in Electrical imaging figure without taking a significant amount of time, the many restrictions of outcrop observation have been broken away from simultaneously, it is reliable that it measures technological principle science, measurement accuracy is good, and computer disposal intelligence degree is high, saves manpower financial capacity, processing speed is fast, and application environment is friendly.
Description
Technical field:
The invention belongs to geologic structure determination techniques field, be related to it is a kind of with electric imaging logging measure paleocurrent to method,
Especially a kind of use passes through pretreated electric imaging logging data, then by dividing the data statistics result of sector to measure
Type water course Sandbody Reservoirs paleocurrent to technical solution.
Background technology:
In petroleum geology exploration technical field paleocurrent direction abbreviation paleocurrent direction analysis and measurement be to research block into
One of the important means of row depositional environment and Analysis on depositional system, paleocurrent is to the main source side for determining a certain period deposit
To and dispersal direction, be related to the regularity of distribution of sand body and the spread of sedimentary system, for predicting favorably to store up the distribution model of body
Enclose, optimize injection-production program, improve recovery ratio etc. is of great significance.Paleocurrent to information can be in sedimentary formation with specific
Form remain, paleocurrent direction is may determine that and restored according to the mark of these forms.For a long time, paleocurrent direction analysis
Rely primarily on the observation to outcrop and measurement, by analyze the deposition directional structure such as a variety of stratifications, level feature and
Aligning for gravel, Skeletal etc. can therefrom extract paleocurrent to information.However, the application of such method is by stratum exposure
The limitation of many objective condition such as situation, when target zone without appear or without corresponding paleocurrent to mark when, field inspection can not be passed through
Determine paleocurrent to.
In the prior art, judge paleocurrent to being a kind of common technological means using well-log information.In the 1990s,
The be widely used in paleocurrent direction analysis and sand body sedimentary structure etc. of dipmeter log have played positive effect, and stratum is inclined
Angle, which is logged well, can measure the inclination angle and azimuth on stratum at wellbore, summarize the response modes of various geologic features, then judge accordingly
The orientation constructed in the layers such as the occurrence on stratum and stratification, then can determine paleocurrent direction on this basis.Dipmeter log
Continuous formation information is can get, has higher longitudinal resolution, but experiment curv is only 4~8, relevant comparative can only be used
Method is handled, preferable in the stable and continuous effect of distribution of strata, complicated and changeable or heterogeneous on ground layer vertical
Property it is strong when less reliable, multi-solution enhancing, explain and be frequently subjected to query using result.Chinese Patent Application No.
CN201510073967.7 discloses a kind of bad ground clay content measurement method based on electric imaging logging data, this method
Accurate clay content data can be provided, to lay the foundation for formation evaluation and fluid identification, but do not refer to measurement paleocurrent
To technical solution;So people are continually looking for new technical method, and using electric imaging logging data determine paleocurrent to
Technology not yet has been reported that.
With the progress of logging technique, the electric imaging logging technology developed by dipmeter log is because that can provide high score
Resolution borehole wall image is used widely, and compared with stratigraphic dip data, Electrical imaging figure can not only provide the inclination angle and side on stratum
Position, common sedimentary structure also have display in figure, therefore, by identifying sedimentary structure, are presented with dip vector result map
Blue pattern, green pattern draw rose figure to determine paleocurrent to being possibly realized, but this method is the problem is that treatment people needs
Have certain explanation experience, specify a variety of paleocurrents to response modes of the mark in image, by way of human-computer interaction
A large amount of sedimentary structures with correlated characteristic are identified and picked up in figure, the time is of high cost, analysis result have it is certain not really
It is qualitative.So to overcome outcrop observe limitation, improve by electric imaging logging data judge paleocurrent to efficiency and
Accuracy rate, find it is a kind of it is new from electric imaging logging data determine paleocurrent to method there is important practical value.
Invention content:
It is an object of the invention to overcome disadvantage of the existing technology, seek to provide a kind of utilization electric imaging logging data
Analysis determine downhole in reservoir paleocurrent to method, first carry out electric imaging logging data processing, then to data carry out layering point fan
Area count, then according to statistical disposition chart measure paleocurrent to.
To achieve the goals above, the technical solution adopted by the present invention comprises the following steps that:
(1) sandstone reservoir type is determined:
The variation of paleocurrent direction (abbreviation paleocurrent to) can lead to the anisotropy of sand body plane, show as the side of permeability
Tropism determines paleocurrent to the direction of permeability is then determined by electric imaging logging data by analyzing the direction of permeability;For
The sand body of fluvial sandbody and various type water courses, permeability direction and paleocurrent are measured by electric imaging logging data to consistent
The paleocurrent of such sand body to;According to this principle, depositional environment and sandstone reservoir class are first determined by existing areal geology data
Type;
(2) imaging data is handled:
First by electric imaging logging instrument (abbreviation instrument) using the stratum conductance in multiple pole plates difference production wells week each orientation
Rate information;Each pole plate is equipped with multiple button electrodes, and each button electrode measurement obtains a conductance profile, suitable according to pole plate
Sequence obtains the two-dimensional conductivity image in ring well week with the expansion of positive north orientation;The bearing data of each pole plate is by being mounted on electric imaging logging
Three orthogonal magnetometer measures in instrument obtain;Since geographical south poles and earth magnetism south poles are misaligned, to what is measured
Orientation is corrected, to exempt finally determining paleocurrent to deviation;Wherein directional correction formula is as follows:
φ'=φ+δ-γ (1)
γ=Δ λ sin ψs (2)
In formula, φ 'For the azimuth after correction;φ is the azimuth measured;δ is magnetic declination, and eastern magnetic declination is just western magnetic
Drift angle is negative;γ is the convergence of meridians of Gaussian plane, and eastern convergency is just that western convergency is negative;Δ λ is to calculate point in
Entreat the difference of longitude between meridian;ψ is to calculate latitude where point;
Again by the computer system equipped with electric imaging logging data processing software, added using integration method or filter method
Velocity correction;Since hole deviation, hole anomaly and wall cake are to the damping action of instrument, electric imaging logging instrument is made not in underground
Regular motion, or even there is the case where card stops, cause surface seismic records depth and instrument real depth to mismatch, shows logging diagram
Upper stretching and compression for image, raw measurement results cannot reflect the electrical variation in corresponding depth stratum, can not be used for really
Determine paleocurrent to, need to carry out acceleration correction to electric imaging logging data, eliminate the stretching in image and compression phenomena, while it is extensive
The real depth of multiple log data;Software wherein of the present invention is the business software of maturation, including the maximum oil field in the whole world
The Techlog softwares of service company Schlumberger Schlumberger, the CIFLog softwares of petrochina and the big soft skill of oil of Beijing stone
The Forward.NET softwares of art Co., Ltd;
Then pole plate data balancing is carried out, electric imaging logging instrument is in underground work, due to the power supply electricity of each pole plate
Pressure has differences, and button electrode surface can form mud film or oil film, affected by these factors, the conductivity value measured with it is true
For real value there are relatively large deviation, it is obviously brighter or dark than remaining pole plate image that the image that ultimately generates will appear some pole plate image
The case where, or occur black stripe or the illusion of vertical nicking in some pole plate image, cause instrument to same stratum
Response there is notable difference;Determine paleocurrent to handling result will be by some pole plate or certain electricity using uncorrected data
The influence of pole, it is difficult to obtain true paleocurrent direction, using data balancing algorithm, the measurement result of all electrodes is carried out
Statistical correction;
(3) Electrical imaging image is generated:
Then electric imaging logging data after treatment, then using coloration standardization scale are mapped at two-dimentional colour code data
It is different colors by colour code data scale, further according to No.1 pole plate orientation 2-D data along just to specific imaging chromatography
The north is unfolded to 360 °, that is, generates the intuitive Electrical imaging image of vision;Coloration standardization used be static color scale determine method with
Dynamic coloration standardization carries out scale, the Electrical imaging image of generation to the data being entirely tested in well section and certain window length respectively
Respectively static map and Dynamic Graph;Static map reflects that the overall variation on stratum, Dynamic Graph are then more clear to the description of local feature
It is clear, subsequent statistical analysis is carried out using static scale data, supplementary explanation is carried out using Dynamic Graph;By multiple pole plate Data Synthesis
Complete X-Y scheme realizes the intuitive display to data, is used for follow-up geologic interpretation;Its Electrical imaging image is according to pole plate orientation edge
Direct north is presented, and has been subjected to directional correction, realizes accurate data positioning, maximum conductivity side can be reflected after processing
To, and then determine paleocurrent direction;Full well section scale is carried out to data, realizes the normalization of data, facilitates subsequent statistical point
Analysis;
(4) shale layer is deducted:
The shale layer being created in Lentic environment cannot reflect the direction of paleocurrent, should be by mud stone before processing and measurement
Layer is deducted from the well section of processing, and conclusion is explained using conventional logging, and dynamic, static map in conjunction with step (3) generation identify
Shale layer therein, and it is deducted from tested interval;
(5) divide sector count:
It is divided into 0.1in (0.00254m) between electric imaging logging data sampling, first chooses 40 sampled points, the system of a length of 0.1m
Window length (i.e. lift height) is counted, using the depth indexing of intermediate point as the depth indexing of statistical result;To improve result of calculation
Longitudinal resolution, it is long that sliding step is chosen to be 1/2 window, i.e. 0.05m;
360 ° of Electrical imaging image is divided into 24 sectors again, 15 ° of each sector calculates separately current window, current
Sector internal conductance rate average value;There are problems that two when practical operation:First, since electric imaging logging instrument is using the survey of multiple pole plates
Amount, wellbore coverage rate is always less than 100%, and there are the blank tapes of no data between pole plate, and with the increase blank bandwidth of hole diameter
Degree gradually increases;Second is that since instrument can rotate during underground survey, the pole plate data and blank tape that measure
Place orientation changes at any time;The two factors cause the data amount check in the long different sectors of same window unequal;If current
Sector is covered by a certain pole plate completely, and the data amount check for statistics is 600 (40 × 15);If there are blank for current sector
Band, data amount check are then less than 600;The most extreme case that can occur is all blank tapes in current sector, and no valid data are used for
Statistics;In order to eliminate due to do not have data or number of samples it is very few caused by statistical error, be arranged that single window is long, single sector number
It is 100 according to number threshold value, is counted if effective data amount check is more than threshold value, calculate conductivity average value;Such as valid data
Number is less than threshold value, then abandons current sector data, do not count;
After obtaining each sector conductivity mean value, maximum value therein is selected, then using sector where maximum value as current
The output valve of window;According to the long moving window of 1/2 window, each sector conductivity mean value is calculated successively, determines maximum conductivity mean value
Terminate until currently processed well section place sector;It is 10m to be provided with single processing well segment length;
(6) rose figure is drawn:
Using the identical division methods with step (5), whole circumference is divided into 24 sectors, 15 ° of each sector will walk
Suddenly it includes maximum electricity that the maximum conductivity mean value position calculated in (5) adds up according to sector numbers and obtains each sector successively
The data amount check of conductance mean value, result is plotted in circumference, forms the rose for determining conductivity maximum value direction
Figure, wherein the radius of each sector is directly proportional with the numerical values recited of corresponding statistical result;In order to protrude the difference between different sectors, respectively
The fill color or gray scale of sector are related to corresponding numerical values recited, and as a result similar sector has same or similar fill color or ash
Degree, it is on the contrary then difference is larger;
(7) determine paleocurrent to:
Determine that the maximum sector direction of radius is conductivity maximum value direction in rose figure;According to water-base mud electricity at
As logging principle obtains that conductivity is bigger to indicate that stratum permeability is better at this, and according to permeability direction and paleocurrent to consistent
Principle, conductivity maximum value direction are paleocurrent to using the maximum sector centers line direction of radius as final true
Fixed paleocurrent to, realize paleocurrent to measurement;
(8) result of calculation is verified:
Pie rock sample and standard plunger piston sample is made in full diameter rock sample, is redirected by remanent magnetism and determines rock core orientation;Use cake
Shape rock sample measures the conductivity of different direction, measures core permeability with standard plunger piston sample, result is drawn in polar coordinates, really
Maximum conductivity direction and permeability principal direction are determined, with paleocurrent determining in step (7) to comparing, with both verifications
Consistency, realization with electric imaging logging measure paleocurrent to the effect of.
Compared with prior art, the present invention paleocurrent direction is determined using electric imaging logging data, when without spending a large amount of
Between in Electrical imaging figure by manual method identify and pick up a large amount of geologic features, while broken away from outcrop observation it is many
Limitation has the unrivaled advantage of other methods in terms of paleocurrent direction determination, has application and promotional value;It is measured
Technological principle science is reliable, and measurement accuracy is good, and computer disposal intelligence degree is high, saves manpower financial capacity, and processing speed is fast,
Application environment is friendly.
Description of the drawings:
Fig. 1 is the processing step schematic process flow diagram that the method for the present invention is related to.
Fig. 2 is static and dynamic electric image of the present invention.
Fig. 3 divides sector count principle schematic to be of the present invention.
Fig. 4 be it is of the present invention for determine paleocurrent to rose figure.
Fig. 5 is that circumferential conductivity and permeability survey result of the present invention correspond to schematic diagram.
Specific implementation mode:
Technical solution to further illustrate the present invention with reference to the accompanying drawings and examples.
Embodiment 1:
The present embodiment be related to using electric imaging logging data determine type water course Sandbody Reservoirs paleocurrent to method, realize
Specific process step it is as follows:
(1) sandstone reservoir type is determined:
For the present embodiment according to related geologic information, tested well region is located at Ordos Basin her Shan slope Zhong Bei
Portion, processing interval are 1 section of mountain, and comprehensive Petrographic Features and sedimentary structure feature determine that the interval is sedimentary environment of fluvial facies, symbol
Close the method for the invention application conditions, can be determined by Image Logging Data paleocurrent to;
(2) imaging data is handled:
The tested original electric imaging logging data with well is obtained, first the azimuthal curves measured are corrected;It should
Well point magnetic declination is -2.63 °, and it is 1.88 ° to calculate the difference of longitude between point and central meridian, and latitude where calculating point is north
38.968 ° of latitude, the convergence of meridians that Gaussian plane is calculated by formula (2) are:
γ=1.88 × sin38.968 ° -1.182 ° (3)
It is herein western convergency, negative value should be taken, therefore the convergence of meridians is -1.182 °, after calculating correction by formula (1)
Azimuth is:
φ'=φ -2.63 ° of -1.448 ° of+1.182 ° of-φ (4)
The orientation curve of instrument record is corrected according to formula (4), eliminates the error that is introduced by magnetic declination, guarantee at
As data are accurately positioned;
Acceleration is carried out to the data after directional correction in the computer equipped with electric imaging logging data processing software
Correction and data equalization processing are eliminated and meet card as instrument, compression of images caused by unfreezing, stretch phenomenon and non_uniform response
Situation provides the input data of high quality for subsequent step;Wherein involved electric imaging logging data processing software is maturation
Business software, including the Techlog softwares of the maximum oilfield service companies Schlumberger Schlumberger in the whole world, petrochina
CIFLog softwares and Beijing Shi great You soft techniques Co., Ltd Forward.NET softwares;
(3) Electrical imaging image is generated:
Conventional static state is respectively adopted to pretreated electric imaging logging data and dynamic coloration scaling method is carved
Degree is two-dimentional colour code data, and when static demarcating carries out scale to full well section, and dynamic calibration uses length for the window of 1m;After calibration
Original conductivity data is converted to the colour code data of 0-255, and colour code data are then mapped as corresponding color with Heat chromatographies,
Form color matrix data;Finally according to the pole plate orientation curve after stringent directional correction by data along direct north exhibition
It opens, that is, generates the intuitive Electrical imaging image of vision, as shown in Fig. 2, it is bent to be from left to right followed successively by depth track, natural gamma in figure
Line, quiescent imaging figure and dynamic imaging figure;
(4) shale layer is deducted:
Label 1 is coring position in Fig. 2, and depth 2975.11m respectively takes 5m as processing up and down using centered on the position
Well section, i.e. 2970~2980m;It is a continuous Reservoir Section, no mud that comprehensive Conventional Logs and Electrical imaging image, which obtain the well section,
Rock stratum takes the well section static imaging data to carry out next step processing so there is no need to deduct;
(5) divide sector count:
Layered shaping is carried out to processing well section static imaging data downwards since initial depth, each small thickness includes
40 sampled points, example well Electrical imaging sampling interval are 0.1in (0.00254m), a length of 0.1m of statistic window;Most from still image
Left side (direct north) starts to be divided into a region for every 15 ° to the right, is the sector region of an equal angular in polar coordinates;Such as
Fig. 3 show layering and divides sector count principle schematic, and black squares 2 indicate that valid data, blank square 3 indicate in figure
The position indicates current window and current sector without valid data, rectangle frame 4;It calculates by longitudinal 40 sampled points, laterally 15 ° of structures
At rectangle frame in include the numbers of valid data, valid data are less than 100, illustrate current region largely or entirely by blank
Band occupies, which is abandoned;Valid data are more than 100, then calculate the average value of data in the region;Window moves right
It is 15 ° dynamic, current region data mean value is calculated successively and is finished up to 24 sectors all calculate, and it is electric to obtain each sector in current window length
Conductance mean value;Maximum value therein is selected from 24 result of calculation, records sector number where it;It is grown according to 1/2 window
Step-length moves down statistical window until processing well section cut-off depth, is calculated successively further according to preceding method and step average
Sector number where conductivity maximum value;
(6) rose figure is drawn:
Based on the result of calculation of step (5), the number for falling the average conductivity maximum value in each sector is counted respectively, and
Rose figure is drawn with radius and area filling color (or gray scale size) that the value is each sector;
Average conductivity maximum value number in 1 each sector of table
| Sector number | Maximum value number | Sector number | Maximum value number |
| 0 | 6 | 12 | 11 |
| 1 | 16 | 13 | 11 |
| 2 | 6 | 14 | 8 |
| 3 | 5 | 15 | 0 |
| 4 | 0 | 16 | 3 |
| 5 | 4 | 17 | 15 |
| 6 | 8 | 18 | 10 |
| 7 | 19 | 19 | 42 |
| 8 | 4 | 20 | 8 |
| 9 | 3 | 21 | 12 |
| 10 | 0 | 22 | 2 |
| 11 | 1 | 23 | 3 |
The rose to be drawn according to 1 data of table is schemed as shown in Figure 4, and 0 ° is direct north, and each sector region occupies 15 °
Range, fan-shaped radius size is directly proportional to average conductivity maximum value number in the region, fan-shaped fill color (or gray scale size)
It is related to maximum value number;
(7) determine paleocurrent to:
Will processing well section Electrical imaging data statistics result be depicted as rose figure after, can determine from figure paleocurrent to;In Fig. 4
The signified fan-shaped radius of label 5 is maximum, indicates that most maximum conductivity average value is fallen in the area, label corresponding thereto
For 6 fan-shaped radius secondly, and the two differs 180 °, illustrate two sector directions for permeability principal direction, i.e., paleocurrent to;
The present embodiment is using angle where fan-shaped center line as final paleocurrent to therefore, from the paleocurrent that Fig. 4 is determined to being 292 °
(or 112 °);Realize with electric imaging logging data determination paleocurrent to function.
Microresistivity log anisotropy caused by anisotropy of the present embodiment using type water course Sandbody Reservoirs itself come
It determines permeability principal direction, and then determines paleocurrent direction, expanded electric imaging logging data answering in paleocurrent direction analysis field
With being of great significance to the research of research area's lithofacies palaeogeography, reservoir prediction and production development.
Embodiment 2:
The present embodiment is related to a kind of method of the measurement result of verification embodiment 1, and it includes coring data, Fig. 2 first to use well
Middle label 1 is coring position, makes pie rock sample and standard plunger piston sample respectively of full-hole core, then redirect really by remanent magnetism
Determine the orientation of rock sample;Pie rock sample is divided into 24 sectors, each sector side cuts out smooth section, measures and calculates every
The conductivity of one group of relative sector, result is drawn in polar coordinates, as shown in label 7 in Fig. 4;Rock is measured with standard plunger piston sample
Measurement result is drawn in polar coordinates, as shown in label 8 in Fig. 5 by heart permeability according to orientation where it;It is determined by Fig. 5
Maximum conductivity direction and permeability principal direction are 112.5 ° (or 292.5 °), compare the present embodiment acquired results and rock physics
Experimental result, the two have good consistency, demonstrate the validity and accuracy of 1 the method for embodiment, illustrate this Shen
The technical solution related to has practicability.
Claims (2)
1. it is a kind of with electric imaging logging measure paleocurrent to method, it is characterised in that the technical solution of use includes following technique step
Suddenly:
(1) sandstone reservoir type is determined:
Paleocurrent direction or paleocurrent to variation can lead to the anisotropy of sand body plane, show as the directionality of permeability,
Determine paleocurrent to and the direction of permeability is then determined by electric imaging logging data by analyzing the direction of permeability;For river
The sand body of quicksand body and various type water courses, permeability direction and paleocurrent can be surveyed by electric imaging logging data to consistent
The paleocurrent of such fixed sand body to;According to this principle, depositional environment and sandstone reservoir type are first determined by areal geology data;
(2) imaging data is handled:
The formation conductivity information in production wells week each orientation is first distinguished using multiple pole plates by electric imaging logging instrument or instrument;
Each pole plate is equipped with multiple button electrodes, and each button electrode measurement obtains a conductance profile, according to pole plate sequence with just
North orientation expansion obtains the two-dimensional conductivity image in ring well week;The bearing data of each pole plate is by electric imaging logging instrument
Three orthogonal magnetometer measures obtain;Since geographical south poles and earth magnetism south poles are misaligned, to the orientation that measures into
Row correction, to exempt finally determining paleocurrent to deviation;Wherein directional correction formula is as follows:
φ'=φ+δ-γ (1)
γ=Δ λ sin ψs (2)
In formula, φ 'For the azimuth after correction;φ is the azimuth measured;δ is magnetic declination, and eastern magnetic declination is just western magnetic declination
It is negative;γ is the convergence of meridians of Gaussian plane, and eastern convergency is just that western convergency is negative;Δ λ is to calculate point and center
Difference of longitude between noon line;ψ is to calculate latitude where point;
Again by the computer system equipped with electric imaging logging data processing software, acceleration is carried out using integration method or filter method
Correction;Since hole deviation, hole anomaly and wall cake are to the damping action of instrument, electric imaging logging instrument is made irregular in underground
Movement, or even occur card the case where stopping, surface seismic records depth and instrument real depth being caused to mismatch, shows in logging diagram and is
The stretching and compression of image, raw measurement results cannot reflect the electrical variation in corresponding depth stratum, can not be used to determine Gu
Flow direction needs to carry out acceleration correction to electric imaging logging data, eliminates the stretching in image and compression phenomena, while restoring to survey
The real depth of well data;
Then pole plate data balancing is carried out, electric imaging logging instrument is in underground work, since the supply voltage of each pole plate is deposited
In difference, button electrode surface can form mud film or oil film, affected by these factors, the conductivity value and actual value measured
There are relatively large deviation, the image that ultimately generates will appear some pole plate image obviously feelings brighter or dark than remaining pole plate image
Condition, or occur black stripe or the illusion of vertical nicking in some pole plate image, lead to sound of the instrument to same stratum
Notable difference should occur;Determine that paleocurrent will be by some pole plate or certain electrodes to, handling result using uncorrected data
It influences, it is difficult to obtain true paleocurrent direction, using data balancing algorithm, the measurement result of all electrodes is counted
Correction;
(3) Electrical imaging image is generated:
Then electric imaging logging data after treatment, then using coloration standardization scale are mapped to spy at two-dimentional colour code data
Colour code data scale is different colors by fixed imaging chromatography, further according to No.1 pole plate orientation 2-D data along the positive north
It is unfolded to 360 °, that is, generates the intuitive Electrical imaging image of vision;Coloration standardization used is that static color scale determines method and dynamic
Coloration standardization carries out scale, the Electrical imaging image difference of generation to the data being entirely tested in well section and certain window length respectively
For static map and Dynamic Graph;Static map reflects that the overall variation on stratum, Dynamic Graph are then more clear the description of local feature, adopt
Subsequent statistical analysis is carried out with static scale data, supplementary explanation is carried out using Dynamic Graph;Multiple pole plate Data Synthesis are complete
X-Y scheme, realize intuitive display to data, be used for follow-up geologic interpretation;Its Electrical imaging image is according to pole plate orientation along due north
Direction is presented, and has been subjected to directional correction, is realized accurate data positioning, maximum conductivity direction can be reflected after processing, in turn
Determine paleocurrent direction;Full well section scale is carried out to data, realizes the normalization of data, is convenient for subsequent statistical analysis;
(4) shale layer is deducted:
The shale layer being created in Lentic environment cannot reflect paleocurrent direction, should be by shale layer from before processing and measurement
It is deducted in the well section of reason, conclusion is explained using conventional logging, Dynamic Graph, the static map generated in conjunction with step (3) identifies wherein
Shale layer, and it is deducted from tested interval;
(5) divide sector count:
It is divided into 0.1in (0.00254m) between electric imaging logging data sampling, first chooses 40 sampled points, the statistic window of a length of 0.1m
Long (i.e. lift height), using the depth indexing of intermediate point as the depth indexing of statistical result;To improve the longitudinal direction of result of calculation
Resolution ratio, sliding step select 1/2 a length of 0.05m of window;
360 ° of Electrical imaging image is divided into 24 sectors again, 15 ° of each sector calculates separately current window, current sector
Internal conductance rate average value;There are problems that two when practical operation:First, since electric imaging logging instrument is using the measurement of multiple pole plates, well
For eye coverage rate always less than 100%, there are the blank tapes of no data between pole plate, and gradual with the increase blank bandwidth of hole diameter
Increase;Second is that the pole plate data measured since instrument can rotate during underground survey and blank tape place side
Position changes at any time;The two factors cause the data amount check in the long different sectors of same window unequal;Current sector quilt completely
A certain pole plate is covered, and the data amount check for statistics is 600;Or current sector, there are blank tape, data amount check is then less than
600;The most extreme case that can occur is all blank tapes in current sector, and no valid data are for counting;In order to eliminate due to
Do not have data or number of samples it is very few caused by statistical error, be arranged that single window is long, single sector data amount check threshold value is 100,
Valid data number is then counted more than threshold value, calculates conductivity average value;Or valid data number is less than threshold value, then abandons
Current sector data, do not count;
After obtaining each sector conductivity mean value, maximum value therein is selected, then using sector where maximum value as current window
Output valve;According to the long moving window of 1/2 window, each sector conductivity mean value is calculated successively, determines maximum conductivity mean value place
Terminate until currently processed well section sector;It is 10m to be provided with single processing well segment length;
(6) rose figure is drawn:
Using the identical division methods with step (5), whole circumference is divided into 24 sectors, 15 ° of each sector, by step
(5) it includes maximum conductance that the maximum conductivity mean value position calculated in adds up according to sector numbers and obtains each sector successively
The data amount check of rate mean value, result is plotted in circumference, is formed the rose for determining conductivity maximum value direction and is schemed,
The radius of wherein each sector is directly proportional with the numerical values recited of corresponding statistical result;It is each to fan in order to protrude the difference between different sectors
The fill color or gray scale in area are related to corresponding numerical values recited, and as a result similar sector has same or similar fill color or ash
Degree, it is on the contrary then difference is larger;
(7) determine paleocurrent to:
Determine that the maximum sector direction of radius is conductivity maximum value direction in rose figure;It is surveyed according to water-base mud Electrical imaging
Well principle obtains that conductivity is bigger to indicate that stratum permeability is better at this, and according to permeability direction and paleocurrent to consistent original
Then, conductivity maximum value direction is paleocurrent to using the maximum sector centers line direction of radius as final determining
Paleocurrent to, realize paleocurrent to measurement.
2. a kind of verification measures method of the paleocurrent to measurement result with electric imaging logging, it is characterised in that full diameter rock sample to be made
Pie rock sample and standard plunger piston sample, then redirected by remanent magnetism and determine rock core orientation;The electricity of different direction is measured with pie rock sample
Conductance, with standard plunger piston sample measure core permeability, result is drawn in polar coordinates, determine maximum conductivity direction and
Permeability principal direction, with the paleocurrent that is measured with electric imaging logging to comparing, the consistency both verified, realization is surveyed with Electrical imaging
Well measures the effect of paleocurrent is to result of calculation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810465781.XA CN108693562B (en) | 2018-05-16 | 2018-05-16 | A method of with electric imaging logging measure paleocurrent to |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810465781.XA CN108693562B (en) | 2018-05-16 | 2018-05-16 | A method of with electric imaging logging measure paleocurrent to |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108693562A true CN108693562A (en) | 2018-10-23 |
| CN108693562B CN108693562B (en) | 2019-08-06 |
Family
ID=63847477
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810465781.XA Active CN108693562B (en) | 2018-05-16 | 2018-05-16 | A method of with electric imaging logging measure paleocurrent to |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108693562B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110080754A (en) * | 2019-04-25 | 2019-08-02 | 杭州迅美科技有限公司 | A kind of electric imaging logging image class PERIODIC INTERFERENCE processing method |
| CN112213465A (en) * | 2020-09-29 | 2021-01-12 | 西安石油大学 | Method for recovering palygorskite flow direction and original form of volcanic mechanism based on electric imaging logging |
| CN113393400A (en) * | 2021-06-22 | 2021-09-14 | 中海油田服务股份有限公司 | Processing method and device for eliminating noise characteristics of electric imaging image |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101630016A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | Method for improving imaging quality of vertical seismic profile |
| CN104989360A (en) * | 2015-07-10 | 2015-10-21 | 中国石油天然气股份有限公司 | Characterization method for dynamic fractures of low-permeability oil reservoir |
| US20150316666A1 (en) * | 2014-05-05 | 2015-11-05 | The Board Of Trustees Of The Leland Stanford Junior University | Efficient Similarity Search of Seismic Waveforms |
| CN105334536A (en) * | 2015-12-01 | 2016-02-17 | 中国石油大学(华东) | Effectiveness evaluation method for compact sandstone reservoir map cracking system |
| CN205422693U (en) * | 2016-04-14 | 2016-08-03 | 刘健 | Spike agent degree of depth logging instrument |
| CN107942393A (en) * | 2017-11-02 | 2018-04-20 | 中国科学院地质与地球物理研究所 | One kind is with brill orientation acoustic logging collecting method |
-
2018
- 2018-05-16 CN CN201810465781.XA patent/CN108693562B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101630016A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | Method for improving imaging quality of vertical seismic profile |
| US20150316666A1 (en) * | 2014-05-05 | 2015-11-05 | The Board Of Trustees Of The Leland Stanford Junior University | Efficient Similarity Search of Seismic Waveforms |
| CN104989360A (en) * | 2015-07-10 | 2015-10-21 | 中国石油天然气股份有限公司 | Characterization method for dynamic fractures of low-permeability oil reservoir |
| CN105334536A (en) * | 2015-12-01 | 2016-02-17 | 中国石油大学(华东) | Effectiveness evaluation method for compact sandstone reservoir map cracking system |
| CN205422693U (en) * | 2016-04-14 | 2016-08-03 | 刘健 | Spike agent degree of depth logging instrument |
| CN107942393A (en) * | 2017-11-02 | 2018-04-20 | 中国科学院地质与地球物理研究所 | One kind is with brill orientation acoustic logging collecting method |
Non-Patent Citations (1)
| Title |
|---|
| 李潮流 等: "利用测井信息判识古流向的方法探讨", 《测井技术》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110080754A (en) * | 2019-04-25 | 2019-08-02 | 杭州迅美科技有限公司 | A kind of electric imaging logging image class PERIODIC INTERFERENCE processing method |
| CN112213465A (en) * | 2020-09-29 | 2021-01-12 | 西安石油大学 | Method for recovering palygorskite flow direction and original form of volcanic mechanism based on electric imaging logging |
| CN113393400A (en) * | 2021-06-22 | 2021-09-14 | 中海油田服务股份有限公司 | Processing method and device for eliminating noise characteristics of electric imaging image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108693562B (en) | 2019-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108957549B (en) | Braided river sediment heterogeneous compact sandstone gas reservoir geological modeling method | |
| CN109061765B (en) | Trap evaluation method for heterogeneous thin sandstone interbed reservoir | |
| CN104977618B (en) | A kind of method evaluated shale gas reservoir and find dessert area | |
| Painter et al. | Fractional Lévy motion as a model for spatial variability in sedimentary rock | |
| CN104516018B (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
| CN105044770A (en) | Compact glutenite gas reservoir quantificational prediction method | |
| Williams-Stroud et al. | Microseismicity-constrained discrete fracture network models for stimulated reservoir simulation | |
| CN110412649B (en) | Identification method of single-stage underwater diversion river channel | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| CN108802812A (en) | A kind of formation lithology inversion method of well shake fusion | |
| CN109541685B (en) | River channel sand body identification method | |
| CN108693562B (en) | A method of with electric imaging logging measure paleocurrent to | |
| CN110056346B (en) | Oil reservoir three-dimensional original water saturation simulation method based on trend change function | |
| CN109725348A (en) | A method of sedimentary facies is identified based on seismic data | |
| CN110208859A (en) | Oil-base mud well crack quantitative parameter intelligence computation method based on ultrasonic imaging | |
| CN107589469B (en) | The determination method and apparatus of oil-water interfaces | |
| CN111352172B (en) | Method for acquiring spatial distribution position of uranium anomaly in sand body by well-seismic combination method | |
| CN107515957A (en) | Mud shale SEQUENCE STRATIGRAPHIC DIVISION method | |
| Xu et al. | 3D geostatistical modeling of Lascaux hill from ERT data | |
| CN108663714A (en) | A kind of sedimentary micro depicting method | |
| Bottari et al. | Coseismic damage at an archaeological site in sicily, italy: evidence of roman age earthquake surface faulting | |
| CN103278852B (en) | Utilize the method for seismic data volume waveform configuration characteristic model predicting oil/gas | |
| CN109975876A (en) | A kind of modeling method of the well shake fusion rate pattern based on tectonic level | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| Guderian et al. | Draugen field–successful reservoir management using 4d seismic |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |