CN107765318B - A kind of method and device of determining geologic horizon - Google Patents
A kind of method and device of determining geologic horizon Download PDFInfo
- Publication number
- CN107765318B CN107765318B CN201710932748.9A CN201710932748A CN107765318B CN 107765318 B CN107765318 B CN 107765318B CN 201710932748 A CN201710932748 A CN 201710932748A CN 107765318 B CN107765318 B CN 107765318B
- Authority
- CN
- China
- Prior art keywords
- color
- geologic horizon
- color component
- geologic
- electrical resistivity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- 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/38—Processing data, e.g. for analysis, for interpretation, for correction
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the present application discloses a kind of method and device of determining geologic horizon.The method is provided with the resistivity inversion cross-sectional data of interval of interest, wherein the interval of interest includes at least two geologic horizons；The described method includes: obtaining the corresponding electrical resistivity range of geologic horizon in the interval of interest；Based on the corresponding electrical resistivity range of the geologic horizon, color mapping processing is carried out to the resistivity inversion cross-sectional data, obtains color mapping treated resistivity inversion cross-sectional data；According to the color mapping treated resistivity inversion cross-sectional data, the position of the geologic horizon is determined.The accuracy of the determined geologic horizon on electrical method inverting section can be improved in technical solution provided by the embodiments of the present application.
Description
Technical field
This application involves geophysical information interpretation technique field, in particular to a kind of the method and dress of determining geologic horizon
It sets.
Background technique
Geophysical information explanation is the important stage of geophysical exploration engineering.This stage can do for oil-gas exploration
Clear subsurface structure feature and indicate Favorable Reservoir.
In geophysical information explanation, the differentiation and tracking of geologic horizon are very important work, simultaneously because real
Border geologic horizon often has corresponding resistivity range of variation, so that carrying out chasing after for geologic horizon according to electrical method inverting section
Track is possibly realized.Geologic horizon tracking is one of the key technology in Interpretation of Electrical Anomalies, and good geologic horizon method for tracing can
To increase substantially the efficiency and precision of Interpretation of Electrical Anomalies.
Existing geologic horizon method for tracing is usually to be manually completed by Interpretation of Electrical Anomalies personnel.Specific method can be with
It is the corresponding resistivity range of variation of each geologic horizon of control, the resistivity data shown from electrical method inverting section identifies
Each geologic horizon.This method is not only time-consuming and laborious, but also to make by the knowledge and experience of people at least partly main
The property seen judgement, thus the tracking of geologic horizon may be influenced by artificial priori knowledge, the accuracy of geologic horizon tracking
It is difficult to ensure.
Summary of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device of determining geologic horizon, to improve in electrical method inverting
The accuracy of determined geologic horizon on section.
In order to solve the above technical problems, the method and device that the embodiment of the present application provides a kind of determining geologic horizon is in this way
It realizes:
A kind of method of determining geologic horizon, is provided with the resistivity inversion cross-sectional data of interval of interest, wherein the mesh
Interval include at least two geologic horizons；The described method includes:
Obtain the corresponding electrical resistivity range of geologic horizon in the interval of interest；
Based on the corresponding electrical resistivity range of the geologic horizon, color mapping is carried out to the resistivity inversion cross-sectional data
Processing obtains color mapping treated resistivity inversion cross-sectional data；
According to the color mapping treated resistivity inversion cross-sectional data, the position of the geologic horizon is determined.
It is described to be based on the corresponding electrical resistivity range of the geologic horizon in preferred embodiment, to the resistivity inversion section
Data carry out color mapping processing, obtain color mapping treated resistivity inversion cross-sectional data, comprising:
Based on the corresponding electrical resistivity range of the geologic horizon, designated color is set for the resistivity inversion cross-sectional data
Component value combination；
The new resistivity inversion cross-sectional data for being provided with designated color component value combination is reflected as the color
Resistivity inversion cross-sectional data of penetrating that treated.
It is described to be based on the corresponding electrical resistivity range of the geologic horizon in preferred embodiment, it is the resistivity inversion section
Data are arranged designated color component value and combine, comprising:
The numerical intervals of the corresponding electrical resistivity range of the geologic horizon are split point according to the type equal proportion of designated color；
Wherein, the designated color includes: one of R color in RGB color mark, G color and B color color；
According to the numerical intervals after described split point, by the corresponding electrical resistivity range of the geologic horizon resistivity value with
The designated color component value combination is mapped one by one；Wherein, the designated color component value combination includes default first face
The designated color component value of colouring component value, default second color component value and the designated color type；Default first face
Colouring component value indicates the pre-set color point of one of two kinds of colors in the RGB color mark in addition to designated color color
Magnitude；Default second color component value indicates in two kinds of colors in the RGB color mark in addition to the designated color
Another pre-set color component value of color.
In preferred embodiment, when the corresponding electrical resistivity range of the geologic horizon is greater than 1 ohm meter and less than 10 Europe
When momme, default first color component value is the pre-set color component value of the R color, and numerical value is 0, described default
Second color component value is the pre-set color component value of the G color, and numerical value is 202, and the designated color is the B face
Color.
In preferred embodiment, when the corresponding electrical resistivity range of the geologic horizon is more than or equal to 10 ohm meters and small
When 30 ohm meter, default first color component value is the pre-set color component value of the R color, and numerical value is 9, institute
The pre-set color component value that default second color component value is the B color is stated, and numerical value is 49, the designated color is described
G color.
In preferred embodiment, when the corresponding electrical resistivity range of the geologic horizon is more than or equal to 30 ohm meters and small
When 100 ohm meter, default first color component value is the pre-set color component value of the R color, and numerical value is
154, default second color component value is the pre-set color component value of the B color, and numerical value is 0, the designated color
For the G color.
In preferred embodiment, when the corresponding electrical resistivity range of the geologic horizon is more than or equal to 100 ohm meters and small
When 200 ohm meter, default first color component value is the pre-set color component value of the G color, and numerical value is
227, default second color component value is the pre-set color component value of the B color, and numerical value is 0, the designated color
For the R color.
It is described according to the color mapping treated resistivity inversion cross-sectional data in preferred embodiment, it determines describedly
The position of matter layer position, comprising:
Obtain electric logging data of the interval of interest at drilling well position；
According to the electric logging data, the relative position of geologic horizon in the interval of interest is determined；Wherein, the geology
The relative position of layer position indicates phase of the geologic horizon on the direction perpendicular to horizontal plane from shallow to deep in the interval of interest
To position；
According to the relative position of the geologic horizon and the color mapping treated resistivity inversion section number
According to determining the position of the geologic horizon.
In preferred embodiment, the color mapping treated resistivity inversion cross-sectional data includes color component combination and ground
The corresponding relationship of layer position and the corresponding relationship of color component combination and resistivity value；Correspondingly, described according to the geology
The relative position of layer position and the color mapping treated resistivity inversion cross-sectional data, determine the geologic horizon
Position, comprising:
When the corresponding electrical resistivity range of two non-conterminous geologic horizons is identical at least two geologic horizon
When, according to the correspondence of color component combination and corresponding relationship and the color component combination and resistivity value of position of stratum
Relationship determines the corresponding position of stratum of the identical electrical resistivity range, and according to the relative position of the geologic horizon, from institute
State the position for determining described two non-conterminous geologic horizons in the corresponding position of stratum of identical electrical resistivity range respectively；
When the corresponding electrical resistivity range of geologic horizon of any two at least two geologic horizon is not identical
When, according to the correspondence of color component combination and corresponding relationship and the color component combination and resistivity value of position of stratum
Relationship, determines the specified corresponding position of stratum of electrical resistivity range, and will the corresponding position of stratum work of the specified electrical resistivity range
For the position for specifying geologic horizon at least two geologic horizon；Wherein, the specified electrical resistivity range indicates the finger
Determine the corresponding electrical resistivity range of geologic horizon.
A kind of device of determining geologic horizon, described device provide the resistivity inversion cross-sectional data of interval of interest, wherein
The interval of interest includes at least two geologic horizons；Described device includes: that electrical resistivity range obtains module, color mapping processing
Module and geologic horizon position determination module；Wherein,
The electrical resistivity range obtains module, for obtaining the corresponding resistivity model of geologic horizon in the interval of interest
It encloses；
The color mapping processing module, for being based on the corresponding electrical resistivity range of the geologic horizon, to the resistance
Rate inverting cross-sectional data carries out color mapping processing, obtains color mapping treated resistivity inversion cross-sectional data；
The geologic horizon position determination module, for according to the color mapping treated resistivity inversion section number
According to determining the position of the geologic horizon.
The embodiment of the present application provides a kind of method and device of determining geologic horizon, in the available interval of interest
The corresponding electrical resistivity range of geologic horizon；It can be based on the corresponding electrical resistivity range of the geologic horizon, it is anti-to the resistivity
It drills cross-sectional data and carries out color mapping processing, available color mapping treated resistivity inversion cross-sectional data；It can root
According to the color mapping treated resistivity inversion cross-sectional data, determines the position of the geologic horizon, may be implemented in electricity
The position of automatic tracing geologic horizon on method inverting section, can be improved the standard of the determined geologic horizon on electrical method inverting section
Exactness.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is a kind of flow chart of the embodiment of the method for determining geologic horizon of the application；
Fig. 2 is to determine geologic horizon in the embodiment of the present application from color mapping treated resistivity inversion cross-sectional data
Position diagrammatic cross-section；
Fig. 3 is the composite structural diagram for the Installation practice that the application determines geologic horizon.
Specific embodiment
The embodiment of the present application provides a kind of method and device of determining geologic horizon.
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without making creative work
Range.
A kind of method of determining geologic horizon provided by the embodiments of the present application.The method of the determining geologic horizon is provided with
The resistivity inversion cross-sectional data of interval of interest.
In the present embodiment, the interval of interest refers to the still undetermined interval in the boundary line of geologic horizon.The purpose
Interval may include at least two geologic horizons.
In the present embodiment, the resistivity inversion cross-sectional data includes resistivity value and position of stratum.Wherein, described
Position of stratum and the resistivity value correspond.
Fig. 1 is a kind of flow chart of the embodiment of the method for determining geologic horizon of the application.As shown in Figure 1, it is described definitely
The method of matter layer position, includes the following steps.
Step S101: the corresponding electrical resistivity range of geologic horizon in the interval of interest is obtained.
It in the present embodiment, can be by the electric logging data and petrophysical data of the statistics interval of interest, really
The classification and the corresponding electrical resistivity range of each geologic horizon for the geologic horizon for including in the fixed interval of interest.For example,
According to the sequence of depth of stratum from shallow to deep, the interval of interest may include: the 4th system (Q), Pliocene (series) Zhang Jiapo group
(N2Z), Miocene series-Pliocene (series) Lantian Ba He group (N2lb-N1G), Paleogene System (E), ancient times circle (Pz2) He Xiagu circle (Pz1) etc. ground
Matter layer position.Table 1 is the corresponding electrical resistivity range of geologic horizon in the embodiment of the present application in interval of interest.As shown in table 1, institute is removed
State the corresponding electrical resistivity range of geologic horizon Q and geologic horizon N in interval of interest2Other than the corresponding electrical resistivity range of z is identical,
The corresponding electrical resistivity range of his geologic horizon is all different.
The corresponding electrical resistivity range of 1 geologic horizon of table
| Geologic horizon | Q | N2z | N2lb-N1g | E | Pz2 | Pz1 |
| Electrical resistivity range/ohm meter | 10-30 | 1-10 | 10-30 | 30-100 | 100-200 | >200 |
Step S102: be based on the corresponding electrical resistivity range of the geologic horizon, to the resistivity inversion cross-sectional data into
The processing of row color mapping obtains color mapping treated resistivity inversion cross-sectional data.
In the present embodiment, it is based on the corresponding electrical resistivity range of the geologic horizon, to the resistivity inversion section
Data carry out color mapping processing, obtain color mapping treated resistivity inversion cross-sectional data, can specifically include, can be with
Based on the corresponding electrical resistivity range of the geologic horizon, designated color component value group is set for the resistivity inversion cross-sectional data
It closes.The new resistivity inversion cross-sectional data of designated color component value combination can will be provided with as the color mapping
Resistivity inversion cross-sectional data that treated.
In the present embodiment, it is based on the corresponding electrical resistivity range of the geologic horizon, is the resistivity inversion section
Data are arranged designated color component value and combine, and can specifically include, can be by the corresponding electrical resistivity range of the geologic horizon
Numerical intervals are split point according to the type equal proportion of designated color.Wherein, the designated color includes: R color in RGB color mark, G face
One of color and B color color.The designated color component value can be determined according to the type of the designated color.Example
Such as, the type of the designated color is 256 kinds, then the designated color component value can be 0~255.It can be according to described
Numerical intervals after splitting point, by the resistivity value and the designated color component in the corresponding electrical resistivity range of the geologic horizon
Value combination is mapped one by one.Wherein, the designated color component value combination includes default first color component value, default second
The designated color component value of color component value and the designated color type；Described in the default first color component value expression
The pre-set color component value of one of two kinds of colors in RGB color mark in addition to designated color color；Described default the
Second colors component value indicates the default of another color in two kinds of colors in the RGB color mark in addition to the designated color
Color component value.
In the present embodiment, when the corresponding electrical resistivity range of the geologic horizon is greater than 1 ohm meter and less than 10
When ohm meter, default first color component value can be the pre-set color component value of the R color, and numerical value is 0, institute
The pre-set color component value that default second color component value can be the G color is stated, and numerical value is 202, the designated color
It can be the B color.
In the present embodiment, when the corresponding electrical resistivity range of the geologic horizon be more than or equal to 10 ohm meters,
And when less than 30 ohm meter, default first color component value can be the pre-set color component value of the R color, and count
Value is 9, and default second color component value can be the pre-set color component value of the B color, and numerical value is 49, the finger
Determining color can be the G color.
In the present embodiment, when the corresponding electrical resistivity range of the geologic horizon be more than or equal to 30 ohm meters,
And when less than 100 ohm meter, default first color component value can be the pre-set color component value of the R color, and
Numerical value is 154, and default second color component value can be the pre-set color component value of the B color, and numerical value is 0, institute
Stating designated color can be the G color.
In the present embodiment, when the corresponding electrical resistivity range of the geologic horizon be more than or equal to 100 ohm meters,
And when less than 200 ohm meter, default first color component value can be the pre-set color component value of the G color, and
Numerical value is 227, and default second color component value can be the pre-set color component value of the B color, and numerical value is 0, institute
Stating designated color can be the R color.
In the present embodiment, described when the corresponding electrical resistivity range of the geologic horizon is greater than 200 ohm meter
Default first color component value can be the pre-set color component value of the G color, and numerical value is 144, default second face
Colouring component value can be the pre-set color component value of the B color, and numerical value is 0, and the designated color can be the R face
Color.
Step S103: according to the color mapping treated resistivity inversion cross-sectional data, the geologic horizon is determined
Position.
In the present embodiment, it according to the color mapping treated resistivity inversion cross-sectional data, determines describedly
The position of matter layer position, can specifically include, electric logging data of the available interval of interest at drilling well position.According to institute
Electric logging data are stated, can determine the relative position of geologic horizon in the interval of interest.Wherein, the geologic horizon is opposite
Position indicates relative position of the geologic horizon on the direction perpendicular to horizontal plane from shallow to deep in the interval of interest.According to
The relative position of the geologic horizon and the color mapping treated resistivity inversion cross-sectional data, can determine institute
State the position of geologic horizon.
In the present embodiment, treated that resistivity inversion cross-sectional data may include color component for the color mapping
The corresponding relationship of combination and corresponding relationship and the color component combination and resistivity value of position of stratum.
In the present embodiment, according to the relative position of the geologic horizon and the color mapping treated electricity
Resistance rate inverting cross-sectional data, determines the position of the geologic horizon, can specifically include, when at least two geologic horizon
When the corresponding electrical resistivity range of two non-conterminous geologic horizons is identical, according to color component combination and position of stratum
Corresponding relationship and color component combination and the corresponding relationship of resistivity value, can determine the identical electrical resistivity range
Corresponding position of stratum, and according to the relative position of the geologic horizon, it can be corresponding from the identical electrical resistivity range
The position of described two non-conterminous geologic horizons is determined in position of stratum respectively.Alternatively, working as at least two geologic horizon
When the corresponding electrical resistivity range of the geologic horizon of middle any two is not identical, according to color component combination and stratum position
The corresponding relationship of corresponding relationship and the color component combination and resistivity value set can determine that specified electrical resistivity range is corresponding
Position of stratum, and can be using the corresponding position of stratum of the specified electrical resistivity range as at least two geologic horizon
The position of specified geologic horizon.Wherein, the specified electrical resistivity range indicates the corresponding resistivity model of the specified geologic horizon
It encloses.
For example, Fig. 2 be in the embodiment of the present application from color mapping treated resistivity inversion cross-sectional data definitely
The diagrammatic cross-section of the position of matter layer position.Abscissa and ordinate in Fig. 2 respectively indicate interval of interest position in the horizontal direction
Set the position with depth of stratum direction.The color actually filled in Fig. 2 is colour, is replaced in figure using gray level image.In figure
Black and white curved line indicate the line of demarcation between each geologic horizon, black skew lines in figure indicates in interval of interest
Tomography.As shown in Fig. 2, the sequence of depth of stratum from shallow to deep can be from color mapping according to the direction with horizontal plane
Geologic horizon Q in the interval of interest, geologic horizon N are determined in resistivity inversion cross-sectional data that treated respectively2Z, geology
Layer position N2lb-N1G, geologic horizon E, geologic horizon Pz2With geologic horizon Pz1Position.
The embodiment of the method for the determining geologic horizon, to obtain the corresponding resistivity of geologic horizon in the interval of interest
Range；It can be based on the corresponding electrical resistivity range of the geologic horizon, color is carried out to the resistivity inversion cross-sectional data and is reflected
Penetrate processing, available color mapping treated resistivity inversion cross-sectional data；After being handled according to the color mapping
Resistivity inversion cross-sectional data, determine the position of the geologic horizon, the automatic tracing on electrical method inverting section may be implemented
The accuracy of the determined geologic horizon on electrical method inverting section can be improved in the position of geologic horizon.
Fig. 3 is the composite structural diagram for the Installation practice that the application determines geologic horizon.The dress of the determining geologic horizon
It sets and the resistivity inversion cross-sectional data of interval of interest is provided, wherein the interval of interest includes at least two geologic horizons.Such as figure
Shown in 3, the device of the determining geologic horizon includes: that electrical resistivity range obtains module 100,200 and of color mapping processing module
Geologic horizon position determination module 300.
The electrical resistivity range obtains module 100, can be used for obtaining the corresponding electricity of geologic horizon in the interval of interest
Resistance rate range.
The color mapping processing module 200 can be used for based on the corresponding electrical resistivity range of the geologic horizon, to institute
It states resistivity inversion cross-sectional data and carries out color mapping processing, obtain color mapping treated resistivity inversion cross-sectional data.
The geologic horizon position determination module 300, can be used for that treated that resistivity is anti-according to the color mapping
Cross-sectional data is drilled, determines the position of the geologic horizon.
The Installation practice of the determining geologic horizon is corresponding with the embodiment of the method for the determining geologic horizon, can be with
Realize the technical solution for determining the embodiment of the method for geologic horizon, and the technical effect of adquisitiones embodiment.
In the 1990s, the improvement of a technology can be distinguished clearly be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So
And with the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow to be programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, programmable logic device
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, logic function determines device programming by user.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, designs and makes without asking chip maker
Dedicated IC chip.Moreover, nowadays, substitution manually makes IC chip, this programming is also used instead mostly " is patrolled
Volume compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development,
And the source code before compiling also write by handy specific programming language, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but there are many kind, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed is most generally used at present
Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art
It will be apparent to the skilled artisan that only needing method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, so that it may it is readily available the hardware circuit for realizing the logical method process.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable
Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions for including in it can also be considered as in hardware component.Or
Even, can will be considered as realizing the device of various functions either the software module of implementation method can be Hardware Subdivision again
Structure in part.
Device that above-described embodiment illustrates, module can specifically realize by computer chip or entity, or by having certain
The product of function is planted to realize.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module can be realized in the same or multiple software and or hardware when application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
It realizes by means of software and necessary general hardware platform.Based on this understanding, the technical solution essence of the application
On in other words the part that contributes to existing technology can be embodied in the form of software products, in a typical configuration
In, calculating equipment includes one or more processors (CPU), input/output interface, network interface and memory.The computer is soft
Part product may include that some instructions are used so that a computer equipment (can be personal computer, server or network
Equipment etc.) execute method described in certain parts of each embodiment of the application or embodiment.The computer software product can
To be stored in memory, memory may include the non-volatile memory in computer-readable medium, random access memory
(RAM) and/or the forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer
The example of readable medium.Computer-readable medium includes that permanent and non-permanent, removable and non-removable media can be by
Any method or technique come realize information store.Information can be computer readable instructions, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only memory
(ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory techniques, CD-ROM are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or
Other magnetic storage devices or any other non-transmission medium, can be used for storage can be accessed by a computing device information.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation
Data-signal and carrier wave.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for device reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as: personal computer, clothes
Business device computer, handheld device or portable device, laptop device, multicomputer system, microprocessor-based system, set
Top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer including any of the above system or equipment
Distributed computing environment etc..
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application there are many deformation and
Variation is without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application's
Spirit.
Claims (9)
1. a kind of method of determining geologic horizon, which is characterized in that the resistivity inversion cross-sectional data of interval of interest is provided,
In, the interval of interest includes at least two geologic horizons；The described method includes:
Obtain the corresponding electrical resistivity range of geologic horizon in the interval of interest；
Based on the corresponding electrical resistivity range of the geologic horizon, the resistivity inversion cross-sectional data is carried out at color mapping
Reason obtains color mapping treated resistivity inversion cross-sectional data；Wherein, the color mapping treated resistivity inversion
Cross-sectional data includes the correspondence of color component combination with corresponding relationship and the color component combination and resistivity value of position of stratum
Relationship；
According to the color mapping treated resistivity inversion cross-sectional data, the position of the geologic horizon is determined, comprising: when
When the corresponding electrical resistivity range of two non-conterminous geologic horizons is identical at least two geologic horizon, according to described
The corresponding relationship of color component combination and corresponding relationship and the color component combination and resistivity value of position of stratum, determines phase
The same corresponding position of stratum of electrical resistivity range, and according to the relative position of the geologic horizon, from the identical resistivity
The position of described two non-conterminous geologic horizons is determined in the corresponding position of stratum of range respectively；When at least two geology
When the corresponding electrical resistivity range of geologic horizon of any two is not identical in layer position, according to color component combination and ground
The corresponding relationship of layer position and the corresponding relationship of color component combination and resistivity value, determine that specified electrical resistivity range is corresponding
Position of stratum, and specified using the corresponding position of stratum of the specified electrical resistivity range as at least two geologic horizon
The position of geologic horizon；Wherein, the specified electrical resistivity range indicates the corresponding electrical resistivity range of the specified geologic horizon.
2. a kind of method of determining geologic horizon according to claim 1, which is characterized in that described to be based on the geological stratification
The corresponding electrical resistivity range in position carries out color mapping processing to the resistivity inversion cross-sectional data, obtains color mapping processing
Resistivity inversion cross-sectional data afterwards, comprising:
Based on the corresponding electrical resistivity range of the geologic horizon, designated color component is set for the resistivity inversion cross-sectional data
Value combination；
Using be provided with the designated color component value combination new resistivity inversion cross-sectional data as the color mapping at
Resistivity inversion cross-sectional data after reason.
3. a kind of method of determining geologic horizon according to claim 2, which is characterized in that described to be based on the geological stratification
The corresponding electrical resistivity range in position is arranged designated color component value for the resistivity inversion cross-sectional data and combines, comprising:
The numerical intervals of the corresponding electrical resistivity range of the geologic horizon are split point according to the type equal proportion of designated color；Its
In, the designated color includes: one of R color in RGB color mark, G color and B color color；
According to the numerical intervals after splitting point, by the corresponding electrical resistivity range of the geologic horizon resistivity value with it is described
The combination of designated color component value is mapped one by one；Wherein, the designated color component value combination includes default first color point
The designated color component value of magnitude, default second color component value and the designated color type；Default first color point
Magnitude indicates the pre-set color component of one of two kinds of colors in the RGB color mark in addition to designated color color
Value；Default second color component value indicates another in two kinds of colors in the RGB color mark in addition to the designated color
A kind of pre-set color component value of color.
4. a kind of method of determining geologic horizon according to claim 3, which is characterized in that when the geologic horizon is corresponding
Electrical resistivity range when being greater than 1 ohm meter and less than 10 ohm meter, default first color component value is the R face
The pre-set color component value of color, and numerical value is 0, default second color component value is the pre-set color component of the G color
Value, and numerical value is 202, the designated color is the B color.
5. a kind of method of determining geologic horizon according to claim 3, which is characterized in that when the geologic horizon is corresponding
Electrical resistivity range when being more than or equal to 10 ohm meters and less than 30 ohm meter, default first color component value
For the pre-set color component value of the R color, and numerical value is 9, and default second color component value is the default of the B color
Color component value, and numerical value is 49, the designated color is the G color.
6. a kind of method of determining geologic horizon according to claim 3, which is characterized in that when the geologic horizon is corresponding
Electrical resistivity range when being more than or equal to 30 ohm meters and less than 100 ohm meter, default first color component value
For the pre-set color component value of the R color, and numerical value is 154, and default second color component value is the pre- of the B color
If color component value, and numerical value is 0, the designated color is the G color.
7. a kind of method of determining geologic horizon according to claim 3, which is characterized in that when the geologic horizon is corresponding
Electrical resistivity range when being more than or equal to 100 ohm meters and less than 200 ohm meter, default first color component
Value is the pre-set color component value of the G color, and numerical value is 227, and default second color component value is the B color
Pre-set color component value, and numerical value is 0, the designated color is the R color.
8. a kind of method of determining geologic horizon according to claim 1, which is characterized in that described to be reflected according to the color
Resistivity inversion cross-sectional data of penetrating that treated, determines the position of the geologic horizon, comprising:
Obtain electric logging data of the interval of interest at drilling well position；
According to the electric logging data, the relative position of geologic horizon in the interval of interest is determined；Wherein, the geologic horizon
Relative position indicate the opposite position of the geologic horizon in the interval of interest on the direction perpendicular to horizontal plane from shallow to deep
It sets；
According to the relative position of the geologic horizon and the color mapping treated resistivity inversion cross-sectional data, really
The position of the fixed geologic horizon.
9. a kind of device of determining geologic horizon, which is characterized in that the resistivity inversion section of described device offer interval of interest
Data, wherein the interval of interest includes at least two geologic horizons；Described device includes: that electrical resistivity range obtains module, face
Color maps processing module and geologic horizon position determination module；Wherein,
The electrical resistivity range obtains module, for obtaining the corresponding electrical resistivity range of geologic horizon in the interval of interest；
The color mapping processing module, it is anti-to the resistivity for being based on the corresponding electrical resistivity range of the geologic horizon
It drills cross-sectional data and carries out color mapping processing, obtain color mapping treated resistivity inversion cross-sectional data；Wherein, the face
The corresponding relationship and color of color mapping treated resistivity inversion cross-sectional data includes color component combination and position of stratum
The corresponding relationship of component combination and resistivity value；
The geologic horizon position determination module, for according to the color mapping treated resistivity inversion cross-sectional data,
Determine the position of the geologic horizon；Wherein, the geologic horizon position determination module is used to work as at least two geological stratification
When the corresponding electrical resistivity range of two non-conterminous geologic horizons is identical in position, according to color component combination and stratum
The corresponding relationship of corresponding relationship and the color component combination and resistivity value of position, determines that identical electrical resistivity range is corresponding
Position of stratum, and according to the relative position of the geologic horizon, from the corresponding position of stratum of the identical electrical resistivity range
The middle position for determining described two non-conterminous geologic horizons respectively；When the ground of any two at least two geologic horizon
When the corresponding electrical resistivity range in matter layer position is not identical, the corresponding relationship with position of stratum is combined according to the color component,
And color component combination and the corresponding relationship of resistivity value, determine the specified corresponding position of stratum of electrical resistivity range, and by institute
The corresponding position of stratum of specified electrical resistivity range is stated as the position for specifying geologic horizon at least two geologic horizon；Its
In, the specified electrical resistivity range indicates the corresponding electrical resistivity range of the specified geologic horizon.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710932748.9A CN107765318B (en) | 2017-10-10 | 2017-10-10 | A kind of method and device of determining geologic horizon |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710932748.9A CN107765318B (en) | 2017-10-10 | 2017-10-10 | A kind of method and device of determining geologic horizon |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107765318A CN107765318A (en) | 2018-03-06 |
| CN107765318B true CN107765318B (en) | 2019-06-11 |
Family
ID=61267235
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710932748.9A Active CN107765318B (en) | 2017-10-10 | 2017-10-10 | A kind of method and device of determining geologic horizon |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107765318B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112763567B (en) * | 2020-12-30 | 2023-11-14 | 核工业北京地质研究院 | Element geochemistry method suitable for variegated stratum division |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102042010A (en) * | 2010-09-07 | 2011-05-04 | 中国石油天然气股份有限公司 | Method for determining development position of fractured-vuggy carbonate reservoir |
| CN103809216A (en) * | 2012-11-12 | 2014-05-21 | 中国石油天然气集团公司 | Resistivity data and earthquake data combined speed field construction method |
| CN103954991A (en) * | 2014-03-13 | 2014-07-30 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Multi-attribute earthquake data inversion method and device |
| CN104854479A (en) * | 2013-02-05 | 2015-08-19 | 哈利伯顿能源服务公司 | Apparatus and methods to visualize formation related features |
| CN106802432A (en) * | 2017-01-17 | 2017-06-06 | 安徽理工大学 | The method for surveying and device of soil layer construction |
| CN107045154A (en) * | 2017-02-08 | 2017-08-15 | 中国海洋石油总公司 | A kind of method and apparatus of identification attitude of stratum in horizontal well environment |
| CN107144878A (en) * | 2017-04-14 | 2017-09-08 | 中国石油天然气股份有限公司 | A kind of fault recognizing method and device |
| CN107656309A (en) * | 2017-08-15 | 2018-02-02 | 中国石油天然气股份有限公司 | A kind of method and device for determining fault throw |
-
2017
- 2017-10-10 CN CN201710932748.9A patent/CN107765318B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102042010A (en) * | 2010-09-07 | 2011-05-04 | 中国石油天然气股份有限公司 | Method for determining development position of fractured-vuggy carbonate reservoir |
| CN103809216A (en) * | 2012-11-12 | 2014-05-21 | 中国石油天然气集团公司 | Resistivity data and earthquake data combined speed field construction method |
| CN104854479A (en) * | 2013-02-05 | 2015-08-19 | 哈利伯顿能源服务公司 | Apparatus and methods to visualize formation related features |
| CN103954991A (en) * | 2014-03-13 | 2014-07-30 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Multi-attribute earthquake data inversion method and device |
| CN106802432A (en) * | 2017-01-17 | 2017-06-06 | 安徽理工大学 | The method for surveying and device of soil layer construction |
| CN107045154A (en) * | 2017-02-08 | 2017-08-15 | 中国海洋石油总公司 | A kind of method and apparatus of identification attitude of stratum in horizontal well environment |
| CN107144878A (en) * | 2017-04-14 | 2017-09-08 | 中国石油天然气股份有限公司 | A kind of fault recognizing method and device |
| CN107656309A (en) * | 2017-08-15 | 2018-02-02 | 中国石油天然气股份有限公司 | A kind of method and device for determining fault throw |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107765318A (en) | 2018-03-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107765298B (en) | A kind of method and device of determining near-surface velocity model | |
| CN107367762B (en) | A kind of method and device of determining reservoir parameter | |
| CN107942404B (en) | A kind of method and device in determining crack | |
| US20200341160A1 (en) | Identifying potential hydrocarbon traps in a subterranean region using recursive anisotropic erosion of seismic data | |
| CN104975850A (en) | Carbonate rock microfacies type identification method and sedimentary facies description method thereof | |
| CN106761672B (en) | Every interlayer inter-wells predicting method and device | |
| Kanaujia et al. | 1D velocity structure and characteristics of contemporary local seismicity around the Tehri region, Garhwal Himalaya | |
| CN106970423B (en) | A kind of method and device of determining shale gas stratigraphic structure model | |
| CN104637075A (en) | Automatic rapid mapping method for sand body | |
| CN107765318B (en) | A kind of method and device of determining geologic horizon | |
| CN109143375B (en) | Method and device for determining drilling direction | |
| CN108090656A (en) | A kind of method and device of definite sand body connectedness | |
| CN107144878A (en) | A kind of fault recognizing method and device | |
| CN107246900A (en) | A kind of method and device for determining oil-water interface position | |
| CN106772599B (en) | A kind of method and device calculating formation shear speed | |
| CN108254786A (en) | A kind of reservoir sequence recognition methods and device | |
| CN106950601A (en) | Static correcting method and device | |
| CN107656309B (en) | A kind of method and device of determining fault throw | |
| CN107703537B (en) | The methods of exhibiting and device of a kind of inspection point in three-dimensional earth's surface | |
| CN106908842B (en) | A kind of method and device of determining individual well sedimentary micro curve | |
| CN106772596B (en) | A kind of method and device of determining pre-stack time migration velocity field | |
| CN107939371B (en) | A kind of method and device of determining well pattern thickening feasibility | |
| Hu et al. | Automatic fracture identification using ant tracking in Tahe oilfield | |
| CN108492014A (en) | A kind of data processing method and device of determining geological resources | |
| CN108107470B (en) | A kind of seismic data processing technique and device |
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 |