CN105607145B - The compaction correction method and device of Clastic Stratum of Country Rocks - Google Patents
The compaction correction method and device of Clastic Stratum of Country Rocks Download PDFInfo
- Publication number
- CN105607145B CN105607145B CN201510564658.XA CN201510564658A CN105607145B CN 105607145 B CN105607145 B CN 105607145B CN 201510564658 A CN201510564658 A CN 201510564658A CN 105607145 B CN105607145 B CN 105607145B
- Authority
- CN
- China
- Prior art keywords
- compaction
- individual well
- current
- lithology
- depth
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of compaction correction method and device of Clastic Stratum of Country Rocks, the compaction correction method of the Clastic Stratum of Country Rocks includes：For each individual well in current Clastic Stratum of Country Rocks, according to the log data of the individual well, the Clastic Stratum of Country Rocks residing for the individual well is divided into different lithology and is compacted unit, log data includes：Porosity, depth and shale content；Unit is compacted according to lithology compaction correction is carried out to current Clastic Stratum of Country Rocks, and calculate the compaction correction amount of each individual well respectively；Planar interpolation is carried out to the compaction correction amount of each individual well, obtains amount of compaction planar interpolation distribution map；According to amount of compaction planar interpolation distribution map and Depth Domain thickness chart, it is determined that the Depth Domain compaction correction thickness chart of current Clastic Stratum of Country Rocks.Compacting dividing elements are carried out based on shale content, lithology can be segmented, improve individual well compaction correction precision；Three-dimensional interpolation compaction correction based on well amount of compaction, it can realize that the compacting of three-dimensional Clastic Stratum of Country Rocks recovers, improve clastic rock paleogeographic reconstruction precision.
Description
Technical field
The present invention relates to clastic rock paleogeographic reconstruction technical field, more particularly to a kind of compaction correction side of Clastic Stratum of Country Rocks
Method and device.
Background technology
Ancient landform controls the spatial of sedimentary system, also controls the source-reservoir-seal assemblage of oil reservoir to a certain extent,
It is the important process for carrying out sedimentology research that ancient landform, which recover,.For clastic deposited sediments stratum, except denudation, structure
Motion, diagenesis and abnormal pressure are made, Compaction is also to influence an important factor for ancient landform is formed and developed.
At present, for the research of Clastic Stratum of Country Rocks compaction correction, it is mostly based on empirical index number of the porosity with change in depth
Model or other porosity-depth models.Different lithology in Clastic Stratum of Country Rocks, its compacted property are also different.At present by chip
Rock stratum is roughly divided into two kinds of lithology of sandstone and mud stone, summarizes Clastic Stratum of Country Rocks using the compacted property of both lithology
The compacted property of all lithology, such as the compacted property of silty, argillaceous siltstoue, clean sandstone and pure shale, do not account for
Influence of the lithology for porosity evolution, have impact on the precision of compaction correction.
The content of the invention
The invention provides a kind of compaction correction method and device of Clastic Stratum of Country Rocks, at least to solve existing clastic rock
In formation compaction trimming process, Clastic Stratum of Country Rocks is roughly divided into two kinds of lithology of sandstone and mud stone, compaction correction precision is not
The problem of high.
According to an aspect of the invention, there is provided a kind of compaction correction method of Clastic Stratum of Country Rocks, including：For current
Each individual well in Clastic Stratum of Country Rocks, according to the log data of the individual well, the Clastic Stratum of Country Rocks residing for the individual well is divided into not
Same lithology compacting unit, wherein, the log data includes：Porosity, depth and shale content；It is compacted according to the lithology
Unit carries out compaction correction to the current Clastic Stratum of Country Rocks, and calculates the compaction correction amount of each individual well respectively；To each list
The compaction correction amount of well carries out planar interpolation, obtains the amount of compaction planar interpolation distribution map of the current Clastic Stratum of Country Rocks；According to
The Depth Domain thickness chart of the amount of compaction planar interpolation distribution map and the current Clastic Stratum of Country Rocks, determines the current clastic rock
The Depth Domain compaction correction thickness chart on stratum.
In one embodiment,, will according to the log data of the individual well for each individual well in current Clastic Stratum of Country Rocks
Clastic Stratum of Country Rocks residing for the individual well is divided into different lithology compacting units, including：, will according to shale content and default step-length
The log data is divided at least two groups, and at least two group is ranked up according to preset order；To institute
State at least two groups and be combined processing, obtain at least two combinations, wherein, the number of the number of the combination and the group
Mesh is identical, and the numbering of the numbering and at least two group of at least two combination corresponds, and present combination I includes：
Group 1 is to group i data, and i is the integer more than or equal to 2, and I is the integer more than or equal to 2；For the hole in each combination
Degrees of data and depth data, according to porosity depth compaction model, the initial porosity and compacting factor of the combination are fitted, and
Calculate the residual error between porosity real data and fitting data；The residual error of more each combination, is determined corresponding to least residual
Combination, and be compacted unit using the lithology in shale content section corresponding to the combination as lithology；Judge the Minimum Residual
Whether last group at least two group is included in combination corresponding to difference；If it is not, then from the minimum
Next group of combination corresponding to residual error starts, and combined treatment is re-started to group, calculate combination residual error and really
Determine lithology compacting unit；If it is, stop calculating.
In one embodiment, the porosity depth compaction model is：φ (z)=φ0e-cz, wherein, φ0Represent initial
Porosity, c represent compacting factor, and z represents depth, and φ (z) represents the porosity size under different depth.
In one embodiment, each individual well in for current Clastic Stratum of Country Rocks, according to the log data of the individual well,
Clastic Stratum of Country Rocks residing for the individual well is divided into before different lithology compacting units, methods described also includes：By individual well mud
Matter content log is demarcated with well log interpretation shale content data, with reduce the individual well shale content log with
Difference between the well log interpretation shale content data.
In one embodiment, it is compacted unit according to the lithology and compaction correction is carried out to the current Clastic Stratum of Country Rocks,
And the compaction correction amount of each individual well is calculated respectively, including：For each individual well, following steps are performed both by determine the individual well
Compaction correction amount：The Clastic Stratum of Country Rocks residing for the current individual well is compacted according to the lithology compacting unit at current individual well
Correction, obtain the compaction correction thickness of each lithology compacting unit；By all lithology compacting unit at the current individual well
Compaction correction thickness is added, and obtains the compaction correction thickness of the current individual well；Compaction correction using the current individual well is thick
Degree divided by original stratigraphic thickness, the compaction correction amount of the current individual well is calculated.
In one embodiment, planar interpolation is carried out to the compaction correction amount of each individual well, including：Utilize Kriging regression
Method carries out planar interpolation to the compaction correction amount of each individual well.
In one embodiment, according to the depth of the amount of compaction planar interpolation distribution map and the current Clastic Stratum of Country Rocks
Domain thickness chart, the Depth Domain compaction correction thickness chart of the current Clastic Stratum of Country Rocks is determined, including：Institute is obtained using seismic interpretation
State the time-domain thickness chart of current Clastic Stratum of Country Rocks；By settling time domain depth domain model, the time-domain thickness chart is turned
It is changed to Depth Domain thickness chart；According to the amount of compaction planar interpolation distribution map and the Depth Domain thickness chart, the depth is determined
Domain compaction correction thickness chart.
According to another aspect of the present invention, there is provided a kind of compaction correction device of Clastic Stratum of Country Rocks, including：Division is single
Member, for for each individual well in current Clastic Stratum of Country Rocks, according to the log data of the individual well, by the chip residing for the individual well
Rock STRATIGRAPHIC DIVISION is compacted unit for different lithology, wherein, the log data includes：Porosity, depth and shale content；Pressure
Real correction and computing unit, compaction correction is carried out to the current Clastic Stratum of Country Rocks for being compacted unit according to the lithology, and
The compaction correction amount of each individual well is calculated respectively；Interpolating unit, for carrying out planar interpolation to the compaction correction amount of each individual well,
Obtain the amount of compaction planar interpolation distribution map of the current Clastic Stratum of Country Rocks；Determining unit, for according to the amount of compaction plane
The Depth Domain thickness chart of interpolation distribution and the current Clastic Stratum of Country Rocks, determine the Depth Domain pressure of the current Clastic Stratum of Country Rocks
Real correction thickness chart.
In one embodiment, the division unit includes：Division module, for according to shale content and default step-length,
The log data is divided at least two groups, and at least two group is ranked up according to preset order；Group
Processing module is closed, for being combined processing at least two group, obtains at least two combinations, wherein, the combination
Number it is identical with the number of the group, the numbering of at least two combination with the numbering of at least two group one by one
Corresponding, present combination I includes：Group 1 is to group i data, and i is the integer more than or equal to 2, and I is the integer more than or equal to 2；
The Fitting Calculation module, for for the porosity data in each combination and depth data, according to porosity depth compaction model,
The initial porosity and compacting factor of the combination are fitted, and calculates the residual error between porosity real data and fitting data；
Comparison module, for the residual error of more each combination, determine the combination corresponding to least residual, and by shale corresponding to the combination
The lithology of content interval is compacted unit as a lithology；Judge module, for judging the combination corresponding to the least residual
In whether include at least two group in last group；If it is not, then from the group corresponding to the least residual
The next group closed starts, and combined treatment is re-started to group, calculates the residual error of combination and determines lithology compacting unit；
If it is, stop calculating.
In one embodiment, the porosity depth compaction model is：φ (z)=φ0e-cz, wherein, φ0Represent initial
Porosity, c represent compacting factor, and z represents depth, and φ (z) represents the porosity size under different depth.
In one embodiment, described device also includes：Demarcate unit, for by individual well shale content log with survey
Well explains that shale content data are demarcated, and is contained with reducing the individual well shale content log with the well log interpretation shale
Measure the difference between data.
In one embodiment, the compaction correction and computing unit, specifically for for each individual well, being performed both by following
Step is to determine the compaction correction amount of the individual well：According to the lithology compacting unit at current individual well to residing for the current individual well
Clastic Stratum of Country Rocks carries out compaction correction, obtains the compaction correction thickness of each lithology compacting unit；At the current individual well
The compaction correction thickness of all lithology compacting units is added, and obtains the compaction correction thickness of the current individual well；Worked as using described
The compaction correction thickness divided by original stratigraphic thickness of preceding individual well, the compaction correction amount of the current individual well is calculated.
In one embodiment, the interpolating unit is specifically used for utilizing compacting of the Kriging regression method to each individual well
Correcting value carries out planar interpolation.
In one embodiment, the determining unit includes：Acquisition module, it is described current for being obtained using seismic interpretation
The time-domain thickness chart of Clastic Stratum of Country Rocks；Modular converter is thick by the time-domain for by settling time domain depth domain model
Degree figure is converted to Depth Domain thickness chart；Determining module, for according to the amount of compaction planar interpolation distribution map and the Depth Domain
Thickness chart, determine the Depth Domain compaction correction thickness chart.
By the compaction correction method and device of the Clastic Stratum of Country Rocks of the present invention, clastic rock compacting is carried out based on shale content
Unit divides automatically, can segment lithology, establishes different lithology compacting units, effectively improves individual well compaction correction and (or is
Compacting recover) precision；Unit is compacted according to the lithology of division and calculates the compaction correction amount of individual well, and then is entered based on well amount of compaction
Row three-dimensional interpolation compaction correction, it can realize that the compacting of three-dimensional Clastic Stratum of Country Rocks recovers, effectively improve clastic rock paleogeographic reconstruction
Precision.
Brief description of the drawings
Accompanying drawing described herein is used for providing a further understanding of the present invention, forms the part of the application, this hair
Bright schematic description and description is used to explain the present invention, does not form limitation of the invention.In the accompanying drawings：
Fig. 1 is the flow chart of the compaction correction method of the Clastic Stratum of Country Rocks of the embodiment of the present invention；
Fig. 2 is the schematic flow sheet of the compacting dividing elements method of the lithology based on shale content of the embodiment of the present invention；
Fig. 3 is the lithology compacting dividing elements schematic diagram of the embodiment of the present invention；
Fig. 4 is the amount of compaction planar interpolation distribution map of the embodiment of the present invention；
Fig. 5 is the structured flowchart of the compaction correction device of the Clastic Stratum of Country Rocks of the embodiment of the present invention；
Fig. 6 is the compacting dividing elements schematic diagram of the individual well lithology based on shale content of the embodiment of the present invention；
Fig. 7 A are the ancient landform three-dimensional displays before the compaction correction of the embodiment of the present invention；
Fig. 7 B are the ancient landform three-dimensional displays after the compaction correction of the embodiment of the present invention.
Embodiment
With reference to the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Ground describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.Based on this
The embodiment of invention, the every other implementation that those of ordinary skill in the art are obtained under the premise of creative work is not made
Example, belongs to protection scope of the present invention.
The embodiments of the invention provide a kind of compaction correction method of Clastic Stratum of Country Rocks, Fig. 1 is the broken of the embodiment of the present invention
Consider the flow chart of the compaction correction method on rock stratum to be worth doing.As shown in figure 1, this method comprises the following steps：
Step S101, for each individual well in current Clastic Stratum of Country Rocks, according to the log data of the individual well, by the individual well
Residing Clastic Stratum of Country Rocks is divided into different lithology compacting units, wherein, log data includes：Porosity, depth and shale
Content.
Step S102, it is compacted unit according to lithology and compaction correction is carried out to current Clastic Stratum of Country Rocks, and calculates respectively each
The compaction correction amount of individual well.
Step S103, planar interpolation is carried out to the compaction correction amount of each individual well, obtains the compacting of current Clastic Stratum of Country Rocks
Measure planar interpolation distribution map.Preferably, plane can be carried out to the compaction correction amount of each individual well using Kriging regression method
Interpolation.
Step S104, according to the Depth Domain thickness chart of amount of compaction planar interpolation distribution map and current Clastic Stratum of Country Rocks (i.e.
Layer original thickness chart), it is determined that the Depth Domain compaction correction thickness chart of current Clastic Stratum of Country Rocks.
By above-described embodiment, clastic rock compacting unit is carried out based on shale content and divided automatically, lithology can be segmented, built
Vertical different lithology compacting unit, effectively improve the precision of individual well compaction correction (or recovering for compacting)；According to the rock of division
Property compacting unit calculate individual well compaction correction amount, and then based on well amount of compaction carry out three-dimensional interpolation compaction correction, three can be realized
The compacting for tieing up Clastic Stratum of Country Rocks recovers, and effectively improves the precision of clastic rock paleogeographic reconstruction.
In one embodiment, step S101 divides lithology compacting unit and can be realized by the following method：
(1) according to shale content and default step-length, log data is divided at least two groups, and according to preset order
At least two group is ranked up.Wherein, default step-length can be configured as needed, for example, 10%, preset order
Can be the size according to shale content, for example, from small to large or from big to small.
(2) processing is combined at least two group, obtains at least two combinations, wherein, the number and group of combination
Other number is identical, and the numbering of at least two combination corresponds with the numbering of at least two group, present combination I bags
Include：Group 1 is to group i data, and i is the integer more than or equal to 2, and I is the integer more than or equal to 2.
For example, being combined processing to i group, obtain I and combine, i=I, wherein, group 1 is as combination 1, group
1st, 2 as combination 2, and group 1,2,3 is as combination 3, and by that analogy, group 1 to i is as combination I.
(3) for the porosity data in each combination and depth data, according to porosity depth compaction model (i.e. pressure
Real mould is based on porosity and depth), the initial porosity and compacting factor of the combination are fitted, and calculate porosity reality
Residual error between border data and fitting data；Wherein, porosity real data is the porosity inputted, and porosity depth is compacted
Model can be：φ (z)=φ0e-cz, φ0Initial porosity is represented, c represents compacting factor, and z represents depth, and φ (z) is represented not
With the porosity size under depth.
(4) residual error of more each combination, the combination corresponding to least residual is determined, and shale corresponding to the combination is contained
The lithology for measuring section is compacted unit as a lithology.
(5) judge in the combination corresponding to least residual whether to include in all groups for being sorted according to preset order
Last group.If it is not, then since next group of the combination corresponding to least residual, group is re-started to group
Conjunction is handled, and is calculated the residual error of combination and is determined lithology compacting unit；If it is, stop calculating.
In actual applications, for a user, porosity, shale content, depth data can be inputted on interface, is led to
Compacting unit can be divided automatically according to the data that user inputs by crossing preset program.
In order to clearly illustrate that the lithology based on shale content is compacted dividing elements method, enter with reference to an example
Row explanation.
Fig. 2 is the schematic flow sheet of the compacting dividing elements method of the lithology based on shale content of the embodiment of the present invention, such as
Shown in Fig. 2, comprise the following steps：
Step S201：Data are classified
By porosity-depth-shale content data according to shale content size, according to certain step-length, ascending point
For different groups, for example, X (1), X (2), X (3), X (4) ... X (N).Wherein it is possible to step-length is selected according to actual conditions,
Too small step-length make it that packet is excessive, and every group of porosity data is very few and influences the fitting of porosity depth compaction curve, mistake
Big step-length does not reach the effect of subdivision compacting unit.Preferably, step-length can be 10%.
For example, minimum shale content is 5%, it is step-length with 10%, data is divided into following group：5%-15%, 15%-
25%, 25%-35%, 35%-45%, 45%-55% etc..
Step S202：Data combine
Processing is combined to the obtained groups of step S201, obtains following combination：(X (1)), (X (1), X (2)), (X
(1), X (2), X (3)), (X (1), X (2), X (3), X (4)) ... ..., until last combination contains the number of all groups
According to.
For example, 5%-15% is combined into, and 5%-25%, 5%-35%, the combination such as 5%-45%, 5%-55%.
Step S203：Data are fitted
Porosity data and depth data in the different combinations that step S202 is obtained, it is compacted according to porosity depth
Model, such as φ (z)=φ0e-cz, fit initial porosity φ0With compacting factor c, and porosity real data and fitting are returned to
Residual error between data.
Step S204：Discriminating data
The residual error for each combination that comparison step S203 is returned, to be combined as breakpoint corresponding to least residual, it is used as the
One lithology is compacted unit.If combination corresponding to least residual include last group (according to the sequence of above-mentioned group,
Last group is shale content highest group), then return to step S201 next group at the breakpoint, again
It is combined calculating；Otherwise step S205 is performed, calculating terminates.
For example, combination (X (1), X (2), X (3)), i.e. the data residual error of 5%-35% is minimum, then by mud corresponding to the combination
The lithology of matter content interval is compacted unit as first lithology.Return to step S201 continues to operate, and remaining group is X (4),
X (5), X (6), X (7) ... X (N), obtain following combination：(X (4)), (X (4), X (5)), (X (4), X (5), X (6)), (X
(4), X (5), X (6), X (7)) ... ..., and carry out step S203 the Fitting Calculation.
Clastic Stratum of Country Rocks is subdivided into different lithology, it is that raising clastic rock compacting is extensive to establish different lithology compacting units
The needs of multiple precision.Lithology be compacted dividing elements schematic diagram as shown in figure 3, leftmost Clastic Stratum of Country Rocks be divided into the right it is several
Kind different lithology compacting unit (lithology 1 as shown in Figure 3, lithology 2, lithology 3, lithology 4).Each of which compacting unit
Longitudinal direction is expressed as depth, from top to bottom deepens；It is laterally porosity, increases from left to right.Circle in each lithology compacting unit
Point reflects under the lithology porosity with the evolution of depth.Formula in Fig. 3 is φ (z)=φ for each lithology0e-cz, such as the φ of lithology 11(z)=φ01e-c1z。
Because inherent mechanism be present between the change of depth in the shale content of clastic rock, initial porosity size, porosity
On contact, and can be proved in terms of rock physicses theoretical model, petrophysics experiment and real data three.Meanwhile
Shale content does not have tendency in depth, does not change with depth.Therefore, using shale content as according to the lithology list segmented out
Member, there can be each independent compacted property.Clastic rock compacting unit is carried out based on shale content to divide automatically, at utmost
Reduce the residual error between measured data and fitting data, the precision that compacting recovers can be effectively improved.
In one embodiment, before step S101, the above method can also include：Individual well shale content is logged well bent
Line is demarcated with well log interpretation shale content data, to reduce individual well shale content log and well log interpretation shale content
Difference between data.The above-mentioned control for the quality of data, ensure that the data differences of two different data sources are little, make to be based on
The porosity of well log interpretation-shale content data can reasonably be applied to individual well.If two data source difference are excessive, abandon
The use of the well data.
In one embodiment, step S102 can be directed to each individual well, be performed both by following steps to determine the individual well
Compaction correction amount：Compacting school is carried out to the Clastic Stratum of Country Rocks residing for current individual well according to the lithology compacting unit at current individual well
Just, the compaction correction thickness of each lithology compacting unit is obtained；By the compacting school of all lithology compacting unit at current individual well
Positive thickness is added, and obtains the compaction correction thickness of current individual well；Utilize the compaction correction thickness divided by prime stratum of current individual well
Thickness, the compaction correction amount of current individual well is calculated.
In one embodiment, step S104 can include：The time of current Clastic Stratum of Country Rocks is obtained using seismic interpretation
Domain thickness chart；By settling time domain depth domain model, time-domain thickness chart is converted into Depth Domain thickness chart；According to amount of compaction
Planar interpolation distribution map and Depth Domain thickness chart, determine Depth Domain compaction correction thickness chart.
By above-mentioned introduction, the three-dimensional interpolation compaction correction flow based on well amount of compaction is as follows：
Step 1：Data quality control
Individual well shale content log and well log interpretation shale content data are demarcated, ensure two different pieces of informations
The data differences in source are little, the porosity based on well log interpretation-shale content data is reasonably applied to individual well.
Step 2：Individual well lithology is compacted unit subdivision
The porosity through quality control, shale content data are inputted, carrying out clastic rock compacting unit based on shale content draws
Point, subdivision lithology compacting unit.
Step 3：Individual well amount of compaction calculates
Compacting recovery is carried out to individual well clastic rock target zone, calculates amount of compaction (the compaction correction thickness/primitively of each individual well
Thickness degree).
Step 4：Kriging method planar interpolation
Planar interpolation is carried out to the compaction correction amount of each individual well using Kriging regression method, obtains amount of compaction planar interpolation
Distribution map, as shown in figure 4, the amount of compaction planar interpolation that the amount of compaction planar interpolation of W1 wellblocks is 1.6, W2 wellblocks is 1.37, W3
The amount of compaction planar interpolation of wellblock is that the amount of compaction planar interpolation of 1.31, W4 wellblocks is that the amount of compaction planes of 1.51, W5 wellblocks is inserted
The amount of compaction planar interpolation being worth for 1.63, W6 wellblocks is 1.55.
Step 5：Compacting recovers thickness chart and calculated
The time-domain thickness chart of target zone can be obtained through seismic interpretation, passes through depth model when establishing, can be by time-domain thickness
Figure is converted to Depth Domain thickness chart.Depth Domain thickness chart with reference to amount of compaction planar interpolation distribution map and after time and depth transfer,
Obtain final Depth Domain compaction correction thickness chart.
Three-dimensional interpolation compaction correction technology based on well amount of compaction, which is fully extracted and applied, is accurately compacted information on well,
And earthquake section information is make use of, give full play to the advantage that seismic data is expanded in the plane.Basic skills is to pass through input
The porosity of each well, shale content data, different lithology is segmented out using the compacting dividing elements method based on shale content
It is compacted unit.Compaction correction is carried out based on this, calculates the compaction correction amount of each individual well.Utilize Kriging regression method pair
The compaction correction amount of each well carries out planar interpolation, obtains amount of compaction planar interpolation distribution map.With reference to the interpolation distribution and process
Depth Domain slice thickness information after time and depth transfer, obtains final compaction correction thickness chart.
The method of the embodiment of the present invention be related to calculate link it is relatively fewer, at utmost remain well amount of compaction information, drop
Low well shake combine caused by error.Meanwhile in Kriging regression link, rationally inserted using amount of compaction distribution trend
Value, interpolation result have higher confidence level.And this method has the characteristics of small, fast, clever, possesses industrialized popularization value.
Based on same inventive concept, the embodiment of the present invention additionally provides a kind of compaction correction device of Clastic Stratum of Country Rocks, can
For realizing the method described by above-described embodiment, as described in the following examples.Due to the compaction correction of Clastic Stratum of Country Rocks
Device solve problem principle it is similar to the compaction correction method of Clastic Stratum of Country Rocks, therefore the implementation of the device may refer to it is above-mentioned
The implementation of method, repeat part and repeat no more.Used below, term " unit " or " module " can realize predetermined function
Software and/or hardware combination.Although system described by following examples is preferably realized with software, hardware,
Or the realization of the combination of software and hardware is also what may and be contemplated.If for example, being realized with hardware, this area can be used
Any one of known following technology or their combination are realized：With for realizing that logic function is patrolled to data-signal
The discrete logic of gate circuit is collected, the application specific integrated circuit with suitable combinational logic gate circuit, programmable gate array
(PGA), field programmable gate array (FPGA) etc..
Fig. 5 is the structured flowchart of the compaction correction device of the Clastic Stratum of Country Rocks of the embodiment of the present invention, as shown in figure 5, this is broken
The compaction correction device on bits rock stratum includes：Division unit 51, compaction correction and computing unit 52, interpolating unit 53 and determination are single
Member 54, the structure is specifically described below.
Division unit 51, for for each individual well in current Clastic Stratum of Country Rocks, according to the log data of the individual well, inciting somebody to action
Clastic Stratum of Country Rocks residing for the individual well is divided into different lithology compacting units, wherein, log data includes：Porosity, depth
And shale content；
Compaction correction and computing unit 52, for carrying out compacting school to current Clastic Stratum of Country Rocks according to lithology compacting unit
Just, the compaction correction amount of each individual well and is calculated respectively；
Interpolating unit 53, for carrying out planar interpolation to the compaction correction amount of each individual well, obtain current Clastic Stratum of Country Rocks
Amount of compaction planar interpolation distribution map；
Determining unit 54, for the Depth Domain thickness according to amount of compaction planar interpolation distribution map and current Clastic Stratum of Country Rocks
Figure, it is determined that the Depth Domain compaction correction thickness chart of current Clastic Stratum of Country Rocks.
By above-described embodiment, clastic rock compacting unit is carried out based on shale content and divided automatically, lithology can be segmented, built
Vertical different lithology compacting unit, effectively improve the precision of individual well compaction correction；Unit is compacted according to the lithology of division and calculates list
The compaction correction amount of well, and then three-dimensional interpolation compaction correction is carried out based on well amount of compaction, the pressure of three-dimensional Clastic Stratum of Country Rocks can be realized
It is real to recover, improve the precision of its paleogeographic reconstruction.
In one embodiment, division unit 51 can include：Division module, combined treatment module, the Fitting Calculation module,
Comparison module and judge module.
Division module, for according to shale content and default step-length, log data being divided into at least two groups, and press
At least two group is ranked up according to preset order；
Combined treatment module, for being combined processing at least two group, at least two combinations are obtained, wherein,
The number of combination is identical with the number of group, a pair of number and the numbering 1 of at least two group of at least two combination
Should, present combination I includes：Group 1 is to group i data, and i is the integer more than or equal to 2, and I is the integer more than or equal to 2；
The Fitting Calculation module, for for the porosity data in each combination and depth data, according to porosity depth
Compaction model, fits the initial porosity and compacting factor of the combination, and calculate porosity real data and fitting data it
Between residual error；
Comparison module, for the residual error of more each combination, determine the combination corresponding to least residual, and by the combination pair
The lithology in the shale content section answered is compacted unit as a lithology；
Judge module, for judging whether include the institute to be sorted according to preset order in the combination corresponding to least residual
There is last group in group；If it is not, then since next group of the combination corresponding to least residual, to group
Combined treatment is re-started, the residual error of combination is calculated and determines lithology compacting unit；If it is, stop calculating.
In one embodiment, porosity depth compaction model can be：φ (z)=φ0e-cz, wherein, φ0Represent initial
Porosity, c represent compacting factor, and z represents depth, and φ (z) represents the porosity size under different depth.
In one embodiment, said apparatus can also include：Unit is demarcated, for by individual well shale content log
Demarcated with well log interpretation shale content data, to reduce individual well shale content log and well log interpretation shale content number
Difference between.
In one embodiment, compaction correction and computing unit 52 are specifically used for being performed both by following steps for each individual well
To determine the compaction correction amount of the individual well：According to the lithology compacting unit at current individual well to the clastic rock residing for current individual well
Layer carries out compaction correction, obtains the compaction correction thickness of each lithology compacting unit；All lithology at current individual well are compacted
The compaction correction thickness of unit is added, and obtains the compaction correction thickness of current individual well；Utilize the compaction correction thickness of current individual well
Divided by original stratigraphic thickness, the compaction correction amount of current individual well is calculated.
In one embodiment, interpolating unit 53 is specifically used for utilizing compacting school of the Kriging regression method to each individual well
Positive quantity carries out planar interpolation.
In one embodiment, determining unit 54 includes：Acquisition module, for obtaining current clastic rock using seismic interpretation
The time-domain thickness chart on stratum；Modular converter, for by settling time domain depth domain model, time-domain thickness chart to be converted to
Depth Domain thickness chart；Determining module, for according to amount of compaction planar interpolation distribution map and Depth Domain thickness chart, determining Depth Domain pressure
Real correction thickness chart.
Certainly, the division of said units or module is a kind of signal division, and the invention is not limited in this.As long as can be real
The Module Division of the existing purpose of the present invention, all should belong to protection scope of the present invention.
In order to carry out apparent explanation to the compaction correction method and device of above-mentioned Clastic Stratum of Country Rocks, with reference to tool
The embodiment of body illustrates, however, it should be noted that the embodiment is merely to be better described the present invention, not structure
The present invention improperly limits in pairs.
By taking Dagang Oilfield individual well Clastic Stratum of Country Rocks as an example, clastic rock compacting dividing elements are carried out based on shale content.Tool
Body, by the porosity of the well, depth, shale content data entry program, step-length 5%, automatic division result such as Fig. 6 institutes
Show, be divided into three groups of different lithology (lithology 1, lithology 2 and lithology 3) automatically according to shale content.In Fig. 6, lines represent to intend
The porosity of conjunction-depth compaction curve；Point outside fitting lines represents the input porosity data of different shale content components
(i.e. real data)；Point on fitting lines represents to be fitted porosity corresponding to porosity-depth compaction curve of fitting.From
Numerical value corresponding to dynamic division is as shown in table 1.
The clastic rock lithology of table 1 compacting unit divides table automatically
As can be seen from Table 1, when not dividing lithology, porosity real data and fitting data residual error are 47.3.Using broken
Lithology is subdivided into three groups by the bits automatic splitting scheme of rock pressure reality unit, and regression criterion is much smaller than the residual error not being grouped.Therefore, base
In clastic rock shale content, different compacting units is divided by segmenting lithology, compaction correction precision can be effectively increased.
By taking the clastic rock target zone of Dagang Oilfield work area as an example, using the three-dimensional interpolation compaction correction skill based on well amount of compaction
Art.Fig. 7 A be the embodiment of the present invention compaction correction before ancient landform three-dimensional display, as shown in Figure 7 A, W1 wellblocks is flatly
Shape, W1 wellblocks and W2 wellblocks physical features height are close.Fig. 7 B be the embodiment of the present invention compaction correction after ancient landform three-dimensional display,
As shown in Figure 7 B, W1 wellblocks show as indentation section, and physical features is significantly lower than W2 wellblocks, and thus W1 wellblocks can determine whether as lake basin center,
Mutually confirmed with the correspondingly deposited phasor of target zone.Numerical value in Fig. 7 A and Fig. 7 B represents thickness, and left side shows that different colors is deep
The shallow different-thickness scope corresponded in 0~900.The well shake for being compacted dividing elements based on shale content by application combines chip
Rock pressure reality bearing calibration, the precision of compaction correction is improved, has really recovered the ancient landform of clastic rock work area target zone.
In summary, the well shake provided in an embodiment of the present invention based on shale content compacting dividing elements combines chip rock pressure
Real bearing calibration, on the basis of traditional porosity-depth empirical index number model, by introducing shale content this parameter, carry
Go out to build porosity-shale content-depth compacting function, devise the compacting unit of the clastic rock based on shale content and divide automatically
Method, to reach subdivision lithology, the purpose of different lithology compacting unit is established, effectively improve the precision of individual well compaction correction, and
Using convenient；The three-dimensional interpolation compaction correction side based on well amount of compaction is applied in the division for being then based on above-mentioned lithology compacting unit
Method, realize that the compacting of three-dimensional Clastic Stratum of Country Rocks recovers, should three-dimensional interpolation compaction correction method based on well amount of compaction have it is easy,
Fast, the characteristics of with a high credibility, possesses oil-gas exploration and development industrial application value.The above method can be applied to oil gas and survey
Visit the research of paleogeographic reconstruction in exploitation.High-precision paleogeographic reconstruction research is the important of raising oil-gas exploration and development success rate
Factor.The present invention can effectively improve the precision of clastic rock paleogeographic reconstruction, finally can be for oil-gas exploration and development target preferably and well
Position deployment provides important references.
Any process or method described otherwise above description in flow chart or herein is construed as, and represents to include
Module, fragment or the portion of the code of the executable instruction of one or more the step of being used to realize specific logical function or process
Point, and the scope of the preferred embodiment of the present invention includes other realization, wherein can not press shown or discuss suitable
Sequence, including according to involved function by it is basic simultaneously in the way of or in the opposite order, carry out perform function, this should be of the invention
Embodiment person of ordinary skill in the field understood.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or the spy for combining the embodiment or example description
Point is contained at least one embodiment or example of the present invention.In this manual, to the schematic representation of above-mentioned term not
Necessarily refer to identical embodiment or example.Moreover, specific features, structure, material or the feature of description can be any
One or more embodiments or example in combine in an appropriate manner.
Particular embodiments described above, the purpose of the present invention, technical scheme and beneficial effect are carried out further in detail
Describe in detail it is bright, should be understood that the foregoing is only the present invention specific embodiment, the guarantor being not intended to limit the present invention
Scope is protected, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc., should be included in this
Within the protection domain of invention.
Claims (12)
1. a kind of compaction correction method of Clastic Stratum of Country Rocks, it is characterised in that including：
For each individual well in current Clastic Stratum of Country Rocks, according to the log data of the individual well, by the clastic rock residing for the individual well
STRATIGRAPHIC DIVISION is compacted unit for different lithology, wherein, the log data includes：Porosity, depth and shale content；
Unit is compacted according to the lithology compaction correction is carried out to the current Clastic Stratum of Country Rocks, and calculate each individual well respectively
Compaction correction amount；
Planar interpolation is carried out to the compaction correction amount of each individual well, obtains the amount of compaction planar interpolation of the current Clastic Stratum of Country Rocks
Distribution map；
According to the amount of compaction planar interpolation distribution map and the Depth Domain thickness chart of the current Clastic Stratum of Country Rocks, it is determined that described work as
The Depth Domain compaction correction thickness chart of preceding Clastic Stratum of Country Rocks；
For each individual well in current Clastic Stratum of Country Rocks, according to the log data of the individual well, by the clastic rock residing for the individual well
STRATIGRAPHIC DIVISION is compacted unit for different lithology, including：
According to shale content and default step-length, the log data is divided at least two groups, and according to preset order pair
At least two group is ranked up；
Processing is combined at least two group, obtains at least two combinations, wherein, the number of the combination with it is described
The number of group is identical, and the numbering of the numbering and at least two group of at least two combination corresponds, current group
Closing I includes：Group 1 is to group i data, and i is the integer more than or equal to 2, and I is the integer more than or equal to 2；
For the porosity data in each combination and depth data, according to porosity depth compaction model, the combination is fitted
Initial porosity and compacting factor, and calculate the residual error between porosity real data and fitting data；
The residual error of more each combination, determines the combination corresponding to least residual, and by shale content section corresponding to the combination
Lithology be compacted unit as lithology；
Judge whether include last group at least two group in the combination corresponding to the least residual；
If it is not, then since next group of the combination corresponding to the least residual, group is re-started at combination
Reason, calculate the residual error of combination and determine lithology compacting unit；
If it is, stop calculating.
2. according to the method for claim 1, it is characterised in that the porosity depth compaction model is：φ (z)=φ0e-cz, wherein, φ0Initial porosity is represented, c represents compacting factor, and z represents depth, and φ (z) represents the porosity under different depth
Size.
3. according to the method for claim 1, it is characterised in that each individual well in for current Clastic Stratum of Country Rocks, root
According to the log data of the individual well, the Clastic Stratum of Country Rocks residing for the individual well is divided into before different lithology compacting units, it is described
Method also includes：
Individual well shale content log and well log interpretation shale content data are demarcated, contained with reducing the individual well shale
Measure the difference between log and the well log interpretation shale content data.
4. according to the method for claim 1, it is characterised in that be compacted unit to the current clastic rock according to the lithology
Stratum carries out compaction correction, and calculates the compaction correction amount of each individual well respectively, including：
For each individual well, following steps are performed both by determine the compaction correction amount of the individual well：
Compaction correction is carried out to the Clastic Stratum of Country Rocks residing for the current individual well according to the lithology compacting unit at current individual well, obtained
To the compaction correction thickness of each lithology compacting unit；
The compaction correction thickness of all lithology compacting unit at the current individual well is added, obtains the pressure of the current individual well
Real correction thickness；
Using the compaction correction thickness divided by original stratigraphic thickness of the current individual well, the compacting of the current individual well is calculated
Correcting value.
5. according to the method for claim 1, it is characterised in that planar interpolation is carried out to the compaction correction amount of each individual well,
Including：Planar interpolation is carried out to the compaction correction amount of each individual well using Kriging regression method.
6. method according to any one of claim 1 to 5, it is characterised in that according to the amount of compaction planar interpolation point
The Depth Domain thickness chart of Butut and the current Clastic Stratum of Country Rocks, determine the Depth Domain compaction correction of the current Clastic Stratum of Country Rocks
Thickness chart, including：
The time-domain thickness chart of the current Clastic Stratum of Country Rocks is obtained using seismic interpretation；
By settling time domain depth domain model, the time-domain thickness chart is converted into Depth Domain thickness chart；
According to the amount of compaction planar interpolation distribution map and the Depth Domain thickness chart, the Depth Domain compaction correction thickness is determined
Figure.
A kind of 7. compaction correction device of Clastic Stratum of Country Rocks, it is characterised in that including：
Division unit, for for each individual well in current Clastic Stratum of Country Rocks, according to the log data of the individual well, by the individual well
Residing Clastic Stratum of Country Rocks is divided into different lithology compacting units, wherein, the log data includes：Porosity, depth and
Shale content；
Compaction correction and computing unit, compacting school is carried out to the current Clastic Stratum of Country Rocks for being compacted unit according to the lithology
Just, the compaction correction amount of each individual well and is calculated respectively；
Interpolating unit, for carrying out planar interpolation to the compaction correction amount of each individual well, obtain the current Clastic Stratum of Country Rocks
Amount of compaction planar interpolation distribution map；
Determining unit, for the Depth Domain thickness according to the amount of compaction planar interpolation distribution map and the current Clastic Stratum of Country Rocks
Figure, determine the Depth Domain compaction correction thickness chart of the current Clastic Stratum of Country Rocks；
The division unit includes：
Division module, for according to shale content and default step-length, the log data being divided into at least two groups, and press
At least two group is ranked up according to preset order；
Combined treatment module, for being combined processing at least two group, at least two combinations are obtained, wherein, institute
State that the number of combination is identical with the number of the group, the volumes of the numbering and at least two group of at least two combination
Number correspond, present combination I includes：Group 1 is to group i data, and i is the integer more than or equal to 2, and I is more than or equal to 2
Integer；
The Fitting Calculation module, for for the porosity data in each combination and depth data, being compacted according to porosity depth
Model, fits the initial porosity and compacting factor of the combination, and calculates between porosity real data and fitting data
Residual error；
Comparison module, for the residual error of more each combination, the combination corresponding to least residual is determined, and by corresponding to the combination
The lithology in shale content section is compacted unit as a lithology；
Judge module, for judging whether include in the combination corresponding to the least residual at least two group most
The latter group；If it is not, then since next group of the combination corresponding to the least residual, group is re-started
Combined treatment, calculate the residual error of combination and determine lithology compacting unit；If it is, stop calculating.
8. device according to claim 7, it is characterised in that the porosity depth compaction model is：φ (z)=φ0e-cz, wherein, φ0Initial porosity is represented, c represents compacting factor, and z represents depth, and φ (z) represents the porosity under different depth
Size.
9. device according to claim 7, it is characterised in that described device also includes：
Unit is demarcated, for individual well shale content log and well log interpretation shale content data to be demarcated, to reduce
Difference between the individual well shale content log and the well log interpretation shale content data.
10. device according to claim 7, it is characterised in that the compaction correction and computing unit, specifically for for
Each individual well, following steps are performed both by determine the compaction correction amount of the individual well：
Compaction correction is carried out to the Clastic Stratum of Country Rocks residing for the current individual well according to the lithology compacting unit at current individual well, obtained
To the compaction correction thickness of each lithology compacting unit；
The compaction correction thickness of all lithology compacting unit at the current individual well is added, obtains the pressure of the current individual well
Real correction thickness；
Using the compaction correction thickness divided by original stratigraphic thickness of the current individual well, the compacting of the current individual well is calculated
Correcting value.
11. device according to claim 7, it is characterised in that the interpolating unit is specifically used for utilizing Kriging regression
Method carries out planar interpolation to the compaction correction amount of each individual well.
12. the device according to any one of claim 7 to 11, it is characterised in that the determining unit includes：
Acquisition module, for obtaining the time-domain thickness chart of the current Clastic Stratum of Country Rocks using seismic interpretation；
Modular converter, for by settling time domain depth domain model, the time-domain thickness chart to be converted into Depth Domain thickness
Figure；
Determining module, for according to the amount of compaction planar interpolation distribution map and the Depth Domain thickness chart, determining the depth
Domain compaction correction thickness chart.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510564658.XA CN105607145B (en) | 2015-09-07 | 2015-09-07 | The compaction correction method and device of Clastic Stratum of Country Rocks |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510564658.XA CN105607145B (en) | 2015-09-07 | 2015-09-07 | The compaction correction method and device of Clastic Stratum of Country Rocks |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105607145A CN105607145A (en) | 2016-05-25 |
| CN105607145B true CN105607145B (en) | 2018-03-13 |
Family
ID=55987220
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510564658.XA Active CN105607145B (en) | 2015-09-07 | 2015-09-07 | The compaction correction method and device of Clastic Stratum of Country Rocks |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105607145B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106405642B (en) * | 2016-07-19 | 2018-06-01 | 西安石油大学 | Seismic inversion method for predicting reservoir based on de- compacting acoustic speed |
| CN109100281A (en) * | 2018-09-10 | 2018-12-28 | 中国石油天然气股份有限公司 | Oil saturation determines method and apparatus |
| CN112083483A (en) * | 2020-09-17 | 2020-12-15 | 北京中恒利华石油技术研究所 | Sedimentary ancient landform restoration method based on rock physical model de-compaction correction |
| CN113777657B (en) * | 2021-09-07 | 2023-10-17 | 中国石油大学(华东) | Slope zone well-connected profile seismic phase characterization industrialization flow based on differential compaction recovery |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007131356A1 (en) * | 2006-05-12 | 2007-11-22 | Stra Terra Inc. | Information characterization system and methods |
| CN103206207A (en) * | 2013-03-01 | 2013-07-17 | 陕西延长石油(集团)有限责任公司研究院 | Favorable oil reservoir region identification method based on hydrocarbon generation period ancient structure |
| CN104122602A (en) * | 2014-07-31 | 2014-10-29 | 中国石油大学(北京) | Method and device for acquiring erosion thickness of sedimentary basin |
| CN104181603A (en) * | 2014-07-24 | 2014-12-03 | 中国石油大学(华东) | Identification method of deposition and diagenetic integrated phase of clastic rocks |
| CN104765966A (en) * | 2015-04-17 | 2015-07-08 | 中国地质大学（北京） | Method and device for measuring original thickness of longshore bar in ancient lake |
-
2015
- 2015-09-07 CN CN201510564658.XA patent/CN105607145B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007131356A1 (en) * | 2006-05-12 | 2007-11-22 | Stra Terra Inc. | Information characterization system and methods |
| CN103206207A (en) * | 2013-03-01 | 2013-07-17 | 陕西延长石油(集团)有限责任公司研究院 | Favorable oil reservoir region identification method based on hydrocarbon generation period ancient structure |
| CN104181603A (en) * | 2014-07-24 | 2014-12-03 | 中国石油大学(华东) | Identification method of deposition and diagenetic integrated phase of clastic rocks |
| CN104122602A (en) * | 2014-07-31 | 2014-10-29 | 中国石油大学(北京) | Method and device for acquiring erosion thickness of sedimentary basin |
| CN104765966A (en) * | 2015-04-17 | 2015-07-08 | 中国地质大学（北京） | Method and device for measuring original thickness of longshore bar in ancient lake |
Non-Patent Citations (2)
| Title |
|---|
| 关于碎屑岩层的去压实校正方法的讨论;漆家福 等;《石油实验地质》;20010930;第23卷(第3期);第351-356页 * |
| 碎屑岩地层的压实校正方法;杨桥 等;《石油实验地质》;20030430;第25卷(第2期);第206-208页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105607145A (en) | 2016-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105607145B (en) | The compaction correction method and device of Clastic Stratum of Country Rocks | |
| Mehrabi et al. | Paleoclimate and tectonic controls on the depositional and diagenetic history of the Cenomanian–early Turonian carbonate reservoirs, Dezful Embayment, SW Iran | |
| Wen et al. | Use of border regions for improved permeability upscaling | |
| CN105913494B (en) | The modeling of multi-scale facture fine geology and method for numerical simulation and device | |
| US20160123121A1 (en) | Method of searching for oil-gas reservoir based on trap-3d software | |
| CN104199124B (en) | Mixed phase stratigraphic analysis method and apparatus | |
| CN110056346B (en) | Oil reservoir three-dimensional original water saturation simulation method based on trend change function | |
| CN106407503B (en) | Forecast Means of Reservoir Fractures and device | |
| CN109783944B (en) | Three-dimensional flow unit model building method, device and equipment | |
| Amorosi et al. | Sand composition changes across key boundaries of siliciclastic and hybrid depositional sequences | |
| Ielpi | Anatomy of major coal successions: Facies analysis and sequence architecture of a brown coal-bearing valley fill to lacustrine tract (Upper Valdarno Basin, Northern Apennines, Italy) | |
| CN104504754A (en) | Multipoint statistic modeling method and device | |
| Palermo et al. | Reservoir properties and petrophysical modelling of carbonate sand bodies: outcrop analogue study in an epicontinental basin (Triassic, Germany) | |
| Daraei et al. | Factors shaping reservoir architecture in the Jurassic Arab carbonates: a case from the Persian Gulf | |
| CN103377196B (en) | A kind of lithofacies palaeogeography data processing method based on GIS-Geographic Information System | |
| Riegel et al. | Petrophysical properties and microstructural analysis of faulted heterolithic packages: a case study from miocene turbidite successions, Italy | |
| Tyler et al. | Faster history matching and uncertainty in predicted production profiles with stochastic modeling | |
| Hornung et al. | Reservoir architecture in a terminal alluvial plain: an outcrop analogue study (Upper Triassic, Southern Germany) Part II: cyclicity, controls and models | |
| CN107219564B (en) | A kind of method and device handling vertical even well profile | |
| Whitaker et al. | Simulating flow through forward sediment model stratigraphies: insights into climatic control of reservoir quality in isolated carbonate platforms | |
| Yarus et al. | Facies simulation in practice: Lithotype proportion mapping and Plurigaussian simulation, a powerful combination | |
| CN107832482A (en) | Compact reservoir multi-scale facture network modelling and analogy method | |
| CN111208558B (en) | Method and device for establishing ultra-deep low-amplitude three-dimensional geological structure | |
| Poland et al. | Sedimentology of an erg to an erg-margin depositional system, the Rush Springs Sandstone of western Oklahoma, USA: Implications for paleowinds across northwestern Pangea during the Guadalupian (Middle Permian) | |
| Burs et al. | Developing a structural and conceptual model of a tectonically limited karst aquifer: a hydrogeological study of the Hastenrather Graben near Aachen, Germany |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |