CN105868482A - Method and device for calculating depositional stage micro-amplitude palaeohigh - Google Patents
Method and device for calculating depositional stage micro-amplitude palaeohigh Download PDFInfo
- Publication number
- CN105868482A CN105868482A CN201610213582.0A CN201610213582A CN105868482A CN 105868482 A CN105868482 A CN 105868482A CN 201610213582 A CN201610213582 A CN 201610213582A CN 105868482 A CN105868482 A CN 105868482A
- Authority
- CN
- China
- Prior art keywords
- structural
- data
- point
- deposition phase
- recent
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
The invention provides a method and device for calculating depositional stage micro-amplitude palaeohigh and relates to the technical field of geological exploration. The method comprises the steps of determining current structural trend-surface map data according to current structural map data, and determining a to-be-constructed amplitude of each inter-well grid point, wherein the to-be-constructed amplitude is vertical direction amplitude of a point on a current structural map relative to the current structural trend-surface; determining depositional stage structural trend-surface map data of a to-be-calculated regional stratum according to a ni-out line of the current structural map; determining relatively high and low data of the depositional stage structure according to the depositional stage structural trend-surface map data and the to-be-constructed amplitude, generating depositional stage micro-amplitude palaeohigh map data according to the relatively high and low data of the depositional stage structure, and determining map data of sedimentary microfacies. According to the method and device for calculating the depositional stage micro-amplitude palaeohigh, the problems that in the prior art, only palaeohigh with the uplift amplitude of about 50 m and above can be calculated and the precision of calculating depositional stage micro-amplitude palaeohigh with the uplift amplitude of only 20 m and even below is poor can be solved.
Description
Technical field
The present invention relates to technical field of geological exploration, particularly relate to a kind of deposit the phase spend a little palaeohigh projectional technique and
Device.
Background technology
Currently, worldwide, directly related with ancient landform oil and gas reserves is more than 5 × 108The big oil of bbl oil equivalent
Gas field is numerous.Ancient landform refers to during earth history interim by factors such as structural deformation, deposit filling, differential compaction and weathering and erosions
Form and remain so far or be embedded in some geomorphic type of underground.The hydrocarbon trap relevant with ancient landform is divided into two big classes: a class
For with the deposition relevant hydrocarbon trap of phase ancient landform, such as: reefs, granule beach, Sha Ba, sand-flood, erosion river valley etc.；Another kind of be with
Degrade the hydrocarbon trap that phase ancient landform is relevant, such as: the carbonate rock produced because of ancient landform or the raw eluviation of table, hill reservoir
Deng.And current ancient landform inversion method is also broadly divided into and degrades ancient landform and deposition ancient landform inverting two broad aspect.Wherein, heavy
Long-pending ancient landform inversion method mainly has impression method, sedimentology method, sequence stratigraphy method and flattening method.Impression method advantage is simple
Easily operation, and shortcoming to be superstratum reference layer be difficult to determine, compaction correction is difficult to determine, can only sxemiquantitative to characterize the deposition phase ancient
Landforms.Sedimentology method advantage is that theoretical basis is enriched, and shortcoming is to be only capable of palaeogeomorphic spatial and temporal pattern of qualitative displaying deposition phase.Sequence
Stratigraphy and flattening method advantage are made by Based On Sequence Stratigraphic Theory and geophysical prospecting technology realizes palaeogeomorphic portraying, and shortcoming is
The quality requiring seismic data wants height, and physical prospecting longitudinal resolution simultaneously is low, deep layer inferior for Use of Geophysical Data relatively thin strate
In deposition phase micro-amplitude scintillation portray unclear.Basin modelling law technology is that Consideration is complete for the advantage of paleogeographic reconstruction,
Data type is many, it is assumed that condition is relatively reasonable, but its shortcoming to be parameter too many, by many factors control, hardly result in reliably
Parameter.
It is the key factor controlling the distribution of tidal-flat facies carbonate rock granule beach that subaqueous deposit phase deposition phase spends palaeohigh a little.
Granule beach is again later stage white cloudsization and grows for the best microfacies of High-quality Reservoir.Therefore, portraying the palaeohigh under water of deposition phase is
Find the prerequisite of tidal-flat facies high-quality dolostone reservoirs.Prior art one uses the side that sequence stratigraphy combines with impression method
Method, chooses destination layer according to Based On Sequence Stratigraphic Theory, is asked for by residual thickness, true thickness correction, compacting recovers, buried hill degrades
The steps such as district determines, paleao-water depth recovery, recover the ancient landform of dam, s_4 formation SQ2 main beach, Shubei area deposition period.According to extensive
Multiple palaeogeomorphic feature, is divided into littoral zone, buial hill zone, north slope zone, south slope belt and 5 landforms of Depression Belt by study area
Unit.Ach-bar sandbody is mainly distributed on littoral zone and south, north slope zone, grows the most with south slope belt.Ancient landform passes through
Deposition landform shape and the water body distribution that controlled thereof and relative water depth, hydrodynamic force, thing source transport passage for transporting and unloading place etc. because of
Element controls the growth of ach-bar sandbody.Dam, chiltern beach, study area major developmental, in SQ2 low stand system tract deposition period, invades system to lake
Area deposition water body in period is deepened, and grows carbonate bank dam in buial hill zone.Prior art two utilizes basin modelling inverting ancient landform
The residual thickness figure on the different times stratum that technology and foundation 2D and 3D seismic data interpretation and drilling well individual-layer data obtain is with existing
Modern buried depth figure, on the basis of to main structural unconformable surface erosion restoration and decompaction correction, carries out paleostructure, ancient landform
Recover；Calculating the ancient buried depth in different configuration period of Different Strata and ancient thickness, system have extensively studied Bachu-Markit
Area paleostructure history of evolution, illustrate Bachu uplift with and field palaeohigh structural evolution relation and Transport, combine Austria and make pottery
It is karst reservoir development characteristics, has inquired into relation between palaeohigh evolution and karst reservoir growth, establish karst reservoir and grow
Structural-lithological reservoir, calculated the Ordovician system favourable karst reservoir Distribution Area.
Inventor studying the present invention and find during above-mentioned prior art, prior art one and prior art two be all for
Calculate that protuberance amplitude is about 50 meters and above palaeohigh.Protuberance amplitude is only had to the deposition phase of 20 meters the least
Spending palaeohigh a little, the projection accuracy of prior art is poor.
Summary of the invention
Embodiments of the invention provide a kind of and deposit the phase and spend projectional technique and the device of palaeohigh a little, to solve existing skill
Art is only capable of calculating that protuberance amplitude, about 50 meters and above palaeohigh, only has 20 meters the least heavy for protuberance amplitude
The long-pending phase spends the problem that the projection accuracy of palaeohigh prior art is poor a little.
For reaching above-mentioned purpose, the present invention adopts the following technical scheme that
A kind of deposit the phase and spend the projectional technique of palaeohigh a little, including:
Obtain the Recent Structural diagram data waiting to calculate stratum, area；
Recent Structural isogram data are determined according to described Recent Structural diagram data；
The plan of mesh point between each well is determined according to described Recent Structural diagram data and described Recent Structural isogram data
Structural range；Described plan Structural range is the vertical direction amplitude relative to Recent Structural trend surface of the point on Recent Structural figure；
Obtain Recent Structural figure degrades line data, and determines wait to calculate according to the line data of degrading of described Recent Structural figure
The deposition phase structural trend surface diagram data on stratum, area；
Determine that the deposition phase constructs according to described deposition phase structural trend surface diagram data and described plan Structural range relatively high
Low data, and spend palaeohigh diagram data a little according to the data genaration deposition phase of height relatively of structure of described deposition phase；
Spend palaeohigh diagram data a little according to the described deposition phase and determine sedimentary micro diagram data.
Concrete, determine Recent Structural isogram data according to described Recent Structural diagram data, including:
The first data point and the second data point is obtained in the Recent Structural figure that Recent Structural diagram data is corresponding；Wherein, institute
The coordinate stating the first data point is (x0, y0, z0), the coordinate of described second data point is (xn, yn, zn)；Described first data point
For the minimum point in Recent Structural figure；Described second data point be in Recent Structural figure centered by peak, with preset away from
From for the minimum point in the range of radius；
The straight line that first data point and the second data point connect presets discrepancy in elevation Δ Z every one and determines one the 3rd data
Point；The coordinate of described 3rd data point is (xk, yk, zk)；
According to formula:
In determining space, each point coordinates (z), and inserted by x, y according to kriging analysis algorithm by point coordinates each in space
Value, generates Recent Structural isogram data；At straight lineOn the height value of point equal
For zk。
Concrete, determine grid between each well according to described Recent Structural diagram data and described Recent Structural isogram data
The plan Structural range of point, including:
By the coordinate of mesh point between each well in Recent Structural diagram data and each point in Recent Structural isogram data
Coordinate in the vertical direction one_to_one corresponding；
Determine in the Recent Structural isogram data that between each well in Recent Structural diagram data, mesh point is corresponding
The amplitude difference of some in the vertical direction, and the plan Structural range of mesh point between each well is determined according to each amplitude difference.
Concrete, obtain Recent Structural figure degrades line data, and according to described Recent Structural figure to degrade line data true
Surely wait to calculate the deposition phase structural trend surface diagram data on stratum, area, including:
Any two points is chosen in line data corresponding degrading of degrading at Recent Structural figure on line, and with wherein as coordinate
Initial point sets up azimuthal coordinates system；The direct north of described azimuthal coordinates system is 0 °, and direction, due east is 90 °；Wherein, in described orientation
In coordinate system origin be (x '0, y '0), another point coordinates in described any two points be (x '1, y '1)；
Will from (x '0, y '0) point to (x '1, y '1) vector direction as degrade line trend and deposition phase structural trend surface
Trend；Will from (x '0, y '0) point to (x '1, y '1) vector vertical direction as deposition phase structural trend surface tendency；
According to formula:
Determine each point coordinates on deposition phase structural trend surface (x ', y ', z ')；Wherein, (x'i,y'i,z'i) it was the deposition phase
Certain point coordinates on structural trend surface；(x'j,y'j,z'j) for along certain point that deposition phase structural trend surface tendency predeterminable range is l '
Coordinate；α is the deposition phase structural trend surface inclination angle pre-set；θ is the direct north degrading line trend with azimuthal coordinates system
Angle；For crossing (x'j,y'j,z'j) point, and move towards direction along deposition phase structural trend surface
Linear equation；zx”,y”=z'jFor crossing (x'j,y'j,z'j) put and along depositing the elevation that phase structural trend surface moves towards the straight line in direction
Value；
Deposition phase structural trend surface diagram data is constituted according to each point coordinates (x ', y ', z ') on deposition phase structural trend surface.
Concrete, determine what the deposition phase constructed according to described deposition phase structural trend surface diagram data and described plan Structural range
Height data relatively, including:
With it, deposition phase structural trend surface diagram data of each mesh point is intended Structural range be added, determine what the deposition phase constructed
Height data relatively.
A kind of deposit the phase and spend the estimation device of palaeohigh a little, including:
Recent Structural diagram data acquiring unit, for obtaining the Recent Structural diagram data waiting to calculate stratum, area；
According to described Recent Structural diagram data, Recent Structural isogram data determination unit, for determining that Recent Structural becomes
Gesture face diagram data；
Intend Structural range and determine unit, for according to described Recent Structural diagram data and described Recent Structural isogram number
According to the plan Structural range determining mesh point between each well；Described plan Structural range is that the point on Recent Structural figure is relative to Recent Structural
The vertical direction amplitude of trend surface；
Deposition phase structural trend surface diagram data determines unit, degrades line data for obtain Recent Structural figure, and according to
The line data of degrading of described Recent Structural figure determine the deposition phase structural trend surface diagram data waiting to calculate stratum, area；
The deposition phase spends palaeohigh diagram data signal generating unit a little, for according to described deposition phase structural trend surface diagram data and
Described plan Structural range determines the data of height relatively that the deposition phase constructs, and according to the data of height relatively of structure of described deposition phase
The generation deposition phase spends palaeohigh diagram data a little；
For spending palaeohigh diagram data a little according to the described deposition phase, sedimentary microfacies map data determination unit, determines that deposition is micro-
Phase diagram data.
Additionally, described Recent Structural isogram data determination unit, specifically for:
The first data point and the second data point is obtained in the Recent Structural figure that Recent Structural diagram data is corresponding；Wherein, institute
The coordinate stating the first data point is (x0, y0, z0), the coordinate of described second data point is (xn, yn, zn)；Described first data point
For the minimum point in Recent Structural figure；Described second data point be in Recent Structural figure centered by peak, with preset away from
From for the minimum point in the range of radius；
The straight line that first data point and the second data point connect presets discrepancy in elevation Δ Z every one and determines one the 3rd data
Point；The coordinate of described 3rd data point is (xk, yk, zk)；
According to formula:
In determining space, each point coordinates (z), and inserted by x, y according to kriging analysis algorithm by point coordinates each in space
Value, generates Recent Structural isogram data；At straight lineOn the height value of point equal
For zk。
Additionally, described plan Structural range determines unit, specifically for:
By the coordinate of mesh point between each well in Recent Structural diagram data and each point in Recent Structural isogram data
Coordinate in the vertical direction one_to_one corresponding；
Determine in the Recent Structural isogram data that between each well in Recent Structural diagram data, mesh point is corresponding
The amplitude difference of some in the vertical direction, and the plan Structural range of mesh point between each well is determined according to each amplitude difference.
Additionally, deposition phase structural trend surface diagram data determines unit, specifically for:
Any two points is chosen in line data corresponding degrading of degrading at Recent Structural figure on line, and with wherein as coordinate
Initial point sets up azimuthal coordinates system；The direct north of described azimuthal coordinates system is 0 °, and direction, due east is 90 °；Wherein, in described orientation
In coordinate system origin be (x '0, y '0), another point coordinates in described any two points be (x '1, y '1)；
Will from (x '0, y '0) point to (x '1, y '1) vector direction as degrade line trend and deposition phase structural trend surface
Trend；Will from (x '0, y '0) point to (x '1, y '1) vector vertical direction as deposition phase structural trend surface tendency；
According to formula:
Determine each point coordinates on deposition phase structural trend surface (x ', y ', z ')；Wherein, (x'i,y'i,z'i) it was the deposition phase
Certain point coordinates on structural trend surface；(x'j,y'j,z'j) for along certain point that deposition phase structural trend surface tendency predeterminable range is l '
Coordinate；α is the deposition phase structural trend surface inclination angle pre-set；θ is the direct north degrading line trend with azimuthal coordinates system
Angle；For crossing (x'j,y'j,z'j) point, and move towards direction along deposition phase structural trend surface
Linear equation；zx”,y”=z'jFor crossing (x'j,y'j,z'j) put and along depositing the elevation that phase structural trend surface moves towards the straight line in direction
Value；
Deposition phase structural trend surface diagram data is constituted according to each point coordinates (x ', y ', z ') on deposition phase structural trend surface.
Additionally, the described deposition phase spend palaeohigh diagram data signal generating unit a little, specifically for:
With it, deposition phase structural trend surface diagram data of each mesh point is intended Structural range be added, determine what the deposition phase constructed
Height data relatively.
What the embodiment of the present invention provided a kind of deposit the phase and spends projectional technique and the device of palaeohigh, its targeted ring a little
Border is the deposition phase of the buried Tibetan tidal-flat facies carbonate rock in weak structural deformation district.The embodiment of the present invention is according to Recent Structural figure number
According to determining Recent Structural isogram data；And then determine respectively according to Recent Structural diagram data and Recent Structural isogram data
The plan Structural range of mesh point between well；Wherein, this plan Structural range is that the point on Recent Structural figure is relative to Recent Structural trend
The vertical direction amplitude in face；Obtain Recent Structural figure afterwards degrades line data, and degrades line data according to Recent Structural figure
Determine the deposition phase structural trend surface diagram data waiting to calculate stratum, area；According to deposition phase structural trend surface diagram data and plan structure
Amplitude determines the data of height relatively that the deposition phase constructs, and according to data genaration deposition phase of height relatively of deposition phase structure a little
Degree palaeohigh diagram data；Spend palaeohigh diagram data a little according to the deposition phase and determine sedimentary micro diagram data.Permissible by the present invention
Protuberance amplitude only deposition phase of 20 meters the least being spent palaeohigh a little calculate, the deposition phase formed spends Gu a little
Protuberance figure precision is higher.Can solve prior art be only capable of calculating protuberance amplitude about 50 meters and above palaeohigh, right
The deposition phase only having 20 meters the least in protuberance amplitude spends palaeohigh a little, the problem that the projection accuracy of prior art is poor.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, also may be used
To obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 provides a kind of for the embodiment of the present invention and deposits the flow chart one that the phase spends the projectional technique of palaeohigh a little；
Fig. 2 provides a kind of for the embodiment of the present invention and deposits the flowchart 2 that the phase spends the projectional technique of palaeohigh a little；
Fig. 3 is Soviet Union and Eastern Europe's block Gu in Ordovician period deposition phase physical features general layout and structure Triassic period the in of Soviet Union and Eastern Europe's block in the embodiment of the present invention
The schematic diagram of the physical features general layout after reversion；
Fig. 4 is the schematic diagram of Marine Carbonate Rocks depositional environment in the embodiment of the present invention；
Fig. 5 is Marine Carbonate Rocks depositional environment rough schematic view in the embodiment of the present invention；
Fig. 6 be in the embodiment of the present invention Recent Structural figure degrade line schematic diagram；
Fig. 7 is the azimuthal coordinates system in the embodiment of the present invention and degrades line schematic diagram；
Fig. 8 is the schematic diagram in an application example region in the embodiment of the present invention；
Fig. 9 is the open-flow capacity in the embodiment of the present invention and the relation schematic diagram of granule beach reservoir thickness；
Figure 10 a kind of deposits the structural representation that the phase spends the estimation device of palaeohigh a little for what the embodiment of the present invention provided.
Detailed description of the invention
Below in conjunction with 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
Describe, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments wholely.Based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under not making creative work premise
Embodiment, broadly falls into the scope of protection of the invention.
When realizing the embodiment of the present invention, inventor finds, existing paleotopography restoration methods is primarily directed to the most several
The scope secondary structural element of ten thousand square kilometres or the biggest third-level tectonic element of sub terrain differential magnitude.For area only
Have hundreds of square kilometre, only in monoclinal structure, growth amplitude 20 meters even more small nose like structure, well spacing have reached 1000 meters very
To less gas field, it is still desirable to understand the anisotropism of High-quality Reservoir distribution, to improving the success rate of Encryption Well.And
The above-mentioned inapplicable area of ancient landform inversion method and the least gas field of Structural range.Its reason is, work area target zone buries
Deep at about 3000 meters, seimic wave propagation distance is remote, central reflector layer is many, causes reflecting that the signal of target zone is weak, and resolution drops
Low；Earth's surface is loose desert region, causes high-frequency signal deep fades during earthquake-capturing, causes seismic resolution to reduce；
Seismic wave propagation speed is low, form strong amplitude, high continuous print seismic reflection, cause its underpart to have again coal seam to cause on target zone
Formation signal is weak, resolution reduces.Earth's surface loosens, bury deep, coal seam interference causes target zone seismic data quality very poor, it is impossible to
Effectively reflection thickness only has the stratigraphic sequence relation of about 50 meters.Therefore, ancient landform inversion method based on seismic data cannot
Application.
Deposit the phase spend the projectional technique of palaeohigh a little as it is shown in figure 1, the embodiment of the present invention provides a kind of, including:
Step 101, obtain the Recent Structural diagram data waiting to calculate stratum, area.
Step 102, determine Recent Structural isogram data according to Recent Structural diagram data.
Step 103, determine the plan of mesh point between each well according to Recent Structural diagram data and Recent Structural isogram data
Structural range.
Wherein, intending Structural range is the vertical direction amplitude relative to Recent Structural trend surface of the point on Recent Structural figure.
Step 104, obtain Recent Structural figure degrade line data, and determine treat according to the line data of degrading of Recent Structural figure
Calculate the deposition phase structural trend surface diagram data on stratum, area.
Step 105, according to deposition phase structural trend surface diagram data with intend Structural range to determine that the deposition phase constructs relatively high
Low data, and spend palaeohigh diagram data a little according to the data genaration deposition phase of height relatively of deposition phase structure.
Step 106, spend palaeohigh diagram data according to the deposition phase a little and determine sedimentary micro diagram data.
What the embodiment of the present invention provided a kind of deposit the phase and spends the projectional technique of palaeohigh a little, and its targeted environment is weak
The deposition phase of the buried Tibetan tidal-flat facies carbonate rock in structural deformation district.The embodiment of the present invention determines according to Recent Structural diagram data
Recent Structural isogram data；And then determine net between each well according to Recent Structural diagram data and Recent Structural isogram data
The plan Structural range of lattice point；Wherein, this plan Structural range is point the erecting relative to Recent Structural trend surface on Recent Structural figure
Nogata is to amplitude；Obtain Recent Structural figure afterwards degrades line data, and determines treat according to the line data of degrading of Recent Structural figure
Calculate the deposition phase structural trend surface diagram data on stratum, area；True according to deposition phase structural trend surface diagram data and plan Structural range
Surely the data of height relatively of deposition phase structure, and according to the data genaration deposition phase Du Gulong a little of height relatively of deposition phase structure
Play diagram data；Spend palaeohigh diagram data a little according to the deposition phase and determine sedimentary micro diagram data.Can be to protuberance by the present invention
Amplitude only deposition phase of 20 meters the least is spent palaeohigh a little and calculates, the deposition phase formed spends palaeohigh figure a little
Precision is higher.Prior art can be solved and be only capable of reckoning protuberance amplitude about 50 meters and above palaeohigh, for protuberance
Amplitude only has the deposition phase of 20 meters the least and spends palaeohigh a little, the problem that the projection accuracy of prior art is poor.
In order to make those skilled in the art be better understood by the present invention, a more detailed embodiment is set forth below,
Depositing the phase spend the projectional technique of palaeohigh a little as in figure 2 it is shown, the embodiment of the present invention provides a kind of, its targeted environment is weak
The deposition phase of the buried Tibetan tidal-flat facies carbonate rock in structural deformation district, the most weak structural deformation district refers to that stratum does not produces substantially
The region of bending, fold or fracture.In order to the deposition phase realized as shown in Figure 2 spends the projectional technique of palaeohigh a little, herein
There is hypothesis below condition in the present invention: as it is shown on figure 3, in weak structural deformation district, the variation of structural trend surface does not has change office
The plan Structural range of portion's point, as Soviet Union and Eastern Europe's block Gu in Ordovician period deposition phase physical features general layout is becoming Soviet Union and Eastern Europe's block after tectonic in-version
Triassic period tectonic in-version after physical features general layout.Or the variation of structural trend surface changes journey for the plan Structural range of well point, local
Spend identical.It is to say, each point now intends the plan Structural range that Structural range is equal to or is proportional to deposition phase each point.
The projectional technique that this deposition phase spends palaeohigh a little includes:
Step 201, obtain the Recent Structural diagram data waiting to calculate stratum, area.
The formation technology of Recent Structural figure herein is the most ripe, such as: utilize reference lamina method, cycle of sedimentation method,
The multiple methods such as uniform thickness method, characteristic curve method, first select the reference lamina on stratum, resettle Comparison of standards section and set up stratum pair
Ratio network, in the range of finally determining gas field, the step such as the degree of depth on top circle, all drilling wells rational stratum and end circle carries out stratum pair
Ratio, and it is converted into the height above sea level degree of depth.Utilize kriging analysis algorithm to calculate the height above sea level degree of depth of grid between well, identical height above sea level line is obtained
To subsurface contour, thus form Recent Structural figure.
Step 202, in the Recent Structural figure that Recent Structural diagram data is corresponding, obtain the first data point and the second data point.
Wherein, the coordinate of the first data point is (x0, y0, z0), the coordinate of the second data point is (xn, yn, zn)。
First data point is the minimum point in Recent Structural figure；Described second data point be in Recent Structural figure with the highest
Centered by Dian, with predeterminable range (such as 5000m) as radius in the range of minimum point；In order to make the plan structure finally given
Amplitude is on the occasion of, z0The minimum value of contour being typically less than in Recent Structural figure, znTypically it is less than in Recent Structural figure
The minima of high point regional area.
Step 203, the straight line that the first data point and the second data point connect presets discrepancy in elevation Δ Z every one determine one the
Three data points.
Wherein, the coordinate of the 3rd data point is (xk, yk, zk)。
Represent the region number divided between the first data point and the second data point.
Step 204, according to formula:
In determining space, each point coordinates (z), and inserted by x, y according to kriging analysis algorithm by point coordinates each in space
Value, generates Recent Structural isogram data；At straight lineOn the height value of point equal
For zk。
Wherein, k ∈ [1, n-1].
Represented the first data point and the linear equation of the second data point.
Represent through the 3rd data point, and be perpendicular to straight line
Linear equation.
Step 205, by the coordinate of mesh point between each well in Recent Structural diagram data and Recent Structural isogram data
In each point coordinates in the vertical direction one_to_one corresponding.
The Recent Structural isogram that between step 206, each well determined in Recent Structural diagram data, mesh point is corresponding
The amplitude difference of the some in the vertical direction in data, and the plan Structural range of mesh point between each well is determined according to each amplitude difference.
What deserves to be explained is, as shown in Figure 4, Marine Carbonate Rocks depositional environment is generally by the slope on land to ocean,
Growing the protuberance of local on this slope, these swell or are submerged under water or surface.The end face of these protuberances is relative
Trend surface in slope has certain forward amplitude.
This deposition landforms of Marine Carbonate Rocks can simplify as it is shown in figure 5, the heave line such as the curve in Fig. 5 on ground
Shown in, the trend surface on slope is as shown in the straight line in Fig. 5.In embodiments of the present invention, intending Structural range is on Recent Structural figure
Point relative to the vertical direction amplitude of Recent Structural trend surface, i.e. point on curve in Fig. 5 is relative to the point on straight line
Vertical direction distance 301.This be different from routine Structural range of the prior art 302, i.e. with the point on the curve in Fig. 5 to
Straight line does vertical line section, and the length of this vertical line section is Structural range 302.
Step 207, choose any two points on line in corresponding the degrading of line data of degrading of Recent Structural figure, and with wherein one
Point sets up azimuthal coordinates system for zero.
The direct north of described azimuthal coordinates system is 0 °, and direction, due east is 90 °；Wherein, in described azimuthal coordinates system Central Plains
Point coordinates be (x '0, y '0), another point coordinates in described any two points be (x '1, y '1)。
What deserves to be explained is, as shown in Figure 6, the line 401 that degrades in Recent Structural figure refers on the ground, side degrading line 401
Layer thickness is more than 0, and is 0 at its opposite side formation thickness.
It addition, this orientation coordinate system can be as it is shown in fig. 7, its direct north be 0 °, direction, due east is 90 °, then (x '0, y
’0) with (x '1, y '1The line angle of strike that degrades represented by) is
Step 208, will from (x '0, y '0) point to (x '1, y '1) vector direction as degrade line trend and deposition phase structure
Make trend surface trend；Will from (x '0, y '0) point to (x '1, y '1) vector vertical direction as deposition phase structural trend surface incline
To.
Step 209, according to formula:
Determine each point coordinates on deposition phase structural trend surface (x ', y ', z '), and according to each point on deposition phase structural trend surface
Coordinate (x ', y ', z ') constitute deposition phase structural trend surface diagram data.
Wherein, (x'i,y'i,z'i) it was certain point coordinates on deposition phase structural trend surface；(x'j,y'j,z'j) it is along heavy
Long-pending phase structural trend surface tendency predeterminable range is the coordinate of certain point of l '；α is the deposition phase structural trend surface inclination angle pre-set；
θ is the angle degrading line trend with the direct north of azimuthal coordinates system；For crossing (x'j,y'j,z
'j) point, and the linear equation in direction is moved towards along deposition phase structural trend surface；zx”,y”=z'jFor crossing (x'j,y'j,z'j) put and edge
Deposition phase structural trend surface moves towards the height value of the straight line in direction.
Step 210, the deposition phase structural trend surface diagram data of each mesh point is intended Structural range with it it is added, determine deposition
The data of height relatively of phase structure, and spend palaeohigh figure number a little according to the data genaration deposition phase of height relatively of deposition phase structure
According to.
Step 211, spend palaeohigh diagram data according to the deposition phase a little and determine sedimentary micro diagram data.
It is that direction, northwest is high that the such as deposition phase spends palaeohigh figure corresponding to palaeohigh diagram data a little, low total of southeastern direction
Physique office.If the rock sample taken out when drilling well, also referred to as the rock color of core show as grey black, then it is ring under water
Border, can incorporate subtidal zone into；If rock type is grainstone or the dolomite with residual particles structure, then its sedimentary micro
For granule beach；If dolomitic white clouds type is evaporation pump, then it is offshore environment, can incorporate supratidal zone into.By this kind
Mode can obtain sedimentary microfacies map accurately.
In order to make those skilled in the art be better understood by the technique effect of the present invention, it is set forth below a present invention's
Application example:
It is theoretical for instructing, utilize logging geology method to follow the trail of each well horse 55 end face and bottom surface depth data with sedimentology,
And combine well head reinforcing the heart height above sea level drilling depth is converted into the height above sea level degree of depth, in the case of considering regional structure pattern, in utilizing gram
Lattice method calculates the height above sea level degree of depth of each point, and forms Recent Structural figure (part A of Fig. 8).The well spacing of study area is to reach 800 meters
Left and right, utilizes the precision of the structural map that well data obtains can meet to portray and spend nose structure a little.From Recent Structural region feature
Set out, calculate mean obliquity and the tendency on stratum, utilize mean obliquity and tendency, calculate Recent Structural trend surface (the B portion of Fig. 8
Point).Utilize Recent Structural to deduct Recent Structural trend, i.e. obtain plan Structural range (C portion of Fig. 8) now.Study area exists
The most overall uplift erosion after horse five sections deposition, the trend of horse 55 strata denudation line is paleotectonic average trend, then paleostructure
It is inclined to just vertical with degrading line average trend.After determining paleotectonic tendency and inclination angle, calculate deposition phase structure trend
Face (the D part of Fig. 8).Calculating deposition phase structural trend surface, the terrain slope less than 3 ° of tidal-flat facies deposition landforms, uses now
Certain angle less than 3 ° is as deposition phase structural trend surface inclination angle.
Based on the present invention it is assumed that changed to Recent Structural trend surface by deposition phase structural trend surface do not change local
The plan Structural range of point, or the variation of structural trend surface is identical for the plan Structural range change degree of partial points.Namely
Saying, each point now intends the plan Structural range that Structural range is equal to or is proportional to deposition phase each point.Recent Structural is utilized to be intended
Structural range adds deposition phase structural trend surface, i.e. can obtain deposition phase paleotopography (the E part of Fig. 8).Gu is obtained in conjunction with inverting
Landform, coring well petrology, the data such as white clouds type, can further inverting weak deformation structure district sedimentary microfacies map (Fig. 8's
F part).Study area has 266 mouthfuls of wells and bores chance horse 55 stratum, and Microfacies Types has supratidal zone cloud level ground, granule beach, Hui Ping, ash balance model
Four kinds, in this work area physical property and reservoir development probability show as table 1 below:
Table 1:
From table 1 it follows that granule beach physical property is better than other microfacies physical property；Study area has 266 mouthfuls of wells and bores chance horse
55 stratum, wherein 76 mouthfuls of wells bore and meet granule beach microfacies, and granule beach Drilling ratio is only 28.6%, but meets 76 mouthfuls of granule beach at brill
In well, having the granule beach of 68 mouthfuls of wells to be evolved into reservoir, the probability finding reservoir after boring chance granule beach is 89.5%.
Surveyed with horse 55 monolayer by the dolostone reservoirs thickness (being called for short granule beach reservoir thickness) developed by granule beach
From the point of view of the open-flow capacity graph of a relation (Fig. 9) of examination high yield well, granule beach reservoir thickness and open-flow capacity are in positive strongly
Pass relation, its correlation coefficient can reach: R=0.95.Therefore, in horse 55 stratum find high yield gas well it is crucial that find
Granule beach dolostone reservoirs, granule beach reservoir is composed again and is stored in the microfacies of granule beach.Known by the practice of oil-gas exploration and development, granule
Beach microfacies is primarily present in the position that the lowest protuberance top water power capacity is stronger.So by the method for the present invention, finding
To deposition phase the lowest protuberance distribution, then can search out the distribution on granule beach, and then search out high yield gas well
Distribution site.As a example by the present embodiment result, find granule beach and to have the probability of high yield gas well be 89.5% × 95%=
85%.
What the embodiment of the present invention provided a kind of deposit the phase and spends the projectional technique of palaeohigh a little, and its targeted environment is weak
The deposition phase of the buried Tibetan tidal-flat facies carbonate rock in structural deformation district.The embodiment of the present invention determines according to Recent Structural diagram data
Recent Structural isogram data；And then determine net between each well according to Recent Structural diagram data and Recent Structural isogram data
The plan Structural range of lattice point；Wherein, this plan Structural range is point the erecting relative to Recent Structural trend surface on Recent Structural figure
Nogata is to amplitude；Obtain Recent Structural figure afterwards degrades line data, and determines treat according to the line data of degrading of Recent Structural figure
Calculate the deposition phase structural trend surface diagram data on stratum, area；True according to deposition phase structural trend surface diagram data and plan Structural range
Surely the data of height relatively of deposition phase structure, and according to the data genaration deposition phase Du Gulong a little of height relatively of deposition phase structure
Play diagram data；Spend palaeohigh diagram data a little according to the deposition phase and determine sedimentary micro diagram data.Can be to protuberance by the present invention
Amplitude only deposition phase of 20 meters the least is spent palaeohigh a little and calculates, the deposition phase formed spends palaeohigh figure a little
Precision is higher.Prior art can be solved and be only capable of reckoning protuberance amplitude about 50 meters and above palaeohigh, for protuberance
Amplitude only has the deposition phase of 20 meters the least and spends palaeohigh a little, the problem that the projection accuracy of prior art is poor.
Corresponding to the embodiment of the method described in above-mentioned Fig. 1 and Fig. 2, as shown in Figure 10, the embodiment of the present invention also provides for one
The deposition phase spends the estimation device of palaeohigh a little, including:
Recent Structural diagram data acquiring unit 51, for obtaining the Recent Structural diagram data waiting to calculate stratum, area.
Recent Structural isogram data determination unit 52, for determining Recent Structural according to described Recent Structural diagram data
Isogram data.
Intend Structural range and determine unit 53, for according to described Recent Structural diagram data and described Recent Structural isogram
Data determine the plan Structural range of mesh point between each well；Described plan Structural range is that the point on Recent Structural figure is relative to structure now
Make the vertical direction amplitude of trend surface.
Deposition phase structural trend surface diagram data determines unit 54, degrades line data, and root for obtain Recent Structural figure
The deposition phase structural trend surface diagram data waiting to calculate stratum, area is determined according to the line data of degrading of described Recent Structural figure.
The deposition phase spends palaeohigh diagram data signal generating unit 55 a little, for according to described deposition phase structural trend surface diagram data
The relative height data that the deposition phase constructs are determined with described plan Structural range, and according to the number of height relatively of structure of described deposition phase
Palaeohigh diagram data is spent a little according to generating the deposition phase.
Sedimentary microfacies map data determination unit 56, determines deposition for spending palaeohigh diagram data a little according to the described deposition phase
Microfacies diagram data.
Additionally, described Recent Structural isogram data determination unit 52, specifically can:
The first data point and the second data point is obtained in the Recent Structural figure that Recent Structural diagram data is corresponding；Wherein, institute
The coordinate stating the first data point is (x0, y0, z0), the coordinate of described second data point is (xn, yn, zn)；Described first data point
For the minimum point in Recent Structural figure；Described second data point be in Recent Structural figure centered by peak, with preset away from
From for the minimum point in the range of radius.
The straight line that first data point and the second data point connect presets discrepancy in elevation Δ Z every one and determines one the 3rd data
Point；The coordinate of described 3rd data point is (xk, yk, zk)；
According to formula:
In determining space, each point coordinates (z), and inserted by x, y according to kriging analysis algorithm by point coordinates each in space
Value, generates Recent Structural isogram data；At straight lineOn the height value of point equal
For zk。
Additionally, this plan Structural range determines unit 53, specifically can:
By the coordinate of mesh point between each well in Recent Structural diagram data and each point in Recent Structural isogram data
Coordinate in the vertical direction one_to_one corresponding.
Determine in the Recent Structural isogram data that between each well in Recent Structural diagram data, mesh point is corresponding
The amplitude difference of some in the vertical direction, and the plan Structural range of mesh point between each well is determined according to each amplitude difference.
Additionally, this deposition phase structural trend surface diagram data determines unit 54, specifically can:
Any two points is chosen in line data corresponding degrading of degrading at Recent Structural figure on line, and with wherein as coordinate
Initial point sets up azimuthal coordinates system；The direct north of described azimuthal coordinates system is 0 °, and direction, due east is 90 °；Wherein, in described orientation
In coordinate system origin be (x '0, y '0), another point coordinates in described any two points be (x '1, y '1)；
Will from (x '0, y '0) point to (x '1, y '1) vector direction as degrade line trend and deposition phase structural trend surface
Trend；Will from (x '0, y '0) point to (x '1, y '1) vector vertical direction as deposition phase structural trend surface tendency；
According to formula:
Determine each point coordinates on deposition phase structural trend surface (x ', y ', z ')；Wherein, (x'i,y'i,z'i) it was the deposition phase
Certain point coordinates on structural trend surface；(x'j,y'j,z'j) for along certain point that deposition phase structural trend surface tendency predeterminable range is l '
Coordinate；α is the deposition phase structural trend surface inclination angle pre-set；θ is the direct north degrading line trend with azimuthal coordinates system
Angle；For crossing (x'j,y'j,z'j) point, and move towards direction along deposition phase structural trend surface
Linear equation；zx”,y”=z'jFor crossing (x'j,y'j,z'j) put and along depositing the elevation that phase structural trend surface moves towards the straight line in direction
Value；
Deposition phase structural trend surface diagram data is constituted according to each point coordinates (x ', y ', z ') on deposition phase structural trend surface.
Additionally, this deposition phase spend palaeohigh diagram data signal generating unit 55 a little, specifically can be by the deposition phase of each mesh point
Structural trend surface diagram data is intended Structural range with it and is added, and determines the data of height relatively that the deposition phase constructs.
What deserves to be explained is, what the embodiment of the present invention provided a kind of deposit the phase and spends estimation device concrete of palaeohigh a little
Implementation may refer to above-mentioned embodiment of the method, and here is omitted.
What the embodiment of the present invention provided a kind of deposit the phase and spends the estimation device of palaeohigh a little, and its targeted environment is weak
The deposition phase of the buried Tibetan tidal-flat facies carbonate rock in structural deformation district.The embodiment of the present invention determines according to Recent Structural diagram data
Recent Structural isogram data；And then determine net between each well according to Recent Structural diagram data and Recent Structural isogram data
The plan Structural range of lattice point；Wherein, this plan Structural range is point the erecting relative to Recent Structural trend surface on Recent Structural figure
Nogata is to amplitude；Obtain Recent Structural figure afterwards degrades line data, and determines treat according to the line data of degrading of Recent Structural figure
Calculate the deposition phase structural trend surface diagram data on stratum, area；True according to deposition phase structural trend surface diagram data and plan Structural range
Surely the data of height relatively of deposition phase structure, and according to the data genaration deposition phase Du Gulong a little of height relatively of deposition phase structure
Play diagram data；Spend palaeohigh diagram data a little according to the deposition phase and determine sedimentary micro diagram data.Can be to protuberance by the present invention
Amplitude only deposition phase of 20 meters the least is spent palaeohigh a little and calculates, the deposition phase formed spends palaeohigh figure a little
Precision is higher.Prior art can be solved and be only capable of reckoning protuberance amplitude about 50 meters and above palaeohigh, for protuberance
Amplitude only has the deposition phase of 20 meters the least and spends palaeohigh a little, the problem that the projection accuracy of prior art is poor.
Those skilled in the art are it should be appreciated that embodiments of the invention can be provided as method, system or computer program
Product.Therefore, the reality in terms of the present invention can use complete hardware embodiment, complete software implementation or combine software and hardware
Execute the form of example.And, the present invention can use at one or more computers wherein including computer usable program code
The upper computer program product implemented of usable storage medium (including but not limited to disk memory, CD-ROM, optical memory etc.)
The form of product.
The present invention is with reference to method, equipment (system) and the flow process of computer program according to embodiments of the present invention
Figure and/or block diagram describe.It should be understood that can the most first-class by computer program instructions flowchart and/or block diagram
Flow process in journey and/or square frame and flow chart and/or block diagram and/or the combination of square frame.These computer programs can be provided
Instruction arrives the processor of general purpose computer, special-purpose computer, Embedded Processor or other programmable data processing device to produce
A raw machine so that the instruction performed by the processor of computer or other programmable data processing device is produced for real
The device of the function specified in one flow process of flow chart or multiple flow process and/or one square frame of block diagram or multiple square frame now.
These computer program instructions may be alternatively stored in and computer or other programmable data processing device can be guided with spy
Determine in the computer-readable memory that mode works so that the instruction being stored in this computer-readable memory produces and includes referring to
Make the manufacture of device, this command device realize at one flow process of flow chart or multiple flow process and/or one square frame of block diagram or
The function specified in multiple square frames.
These computer program instructions also can be loaded in computer or other programmable data processing device so that at meter
Perform sequence of operations step on calculation machine or other programmable devices to produce computer implemented process, thus at computer or
The instruction performed on other programmable devices provides for realizing at one flow process of flow chart or multiple flow process and/or block diagram one
The step of the function specified in individual square frame or multiple square frame.
The present invention applies specific embodiment principle and the embodiment of the present invention are set forth, above example
Explanation be only intended to help to understand method and the core concept thereof of the present invention；Simultaneously for one of ordinary skill in the art,
According to the thought of the present invention, the most all will change, in sum, in this specification
Hold and should not be construed as limitation of the present invention.
Claims (10)
1. one kind deposits the phase and spends the projectional technique of palaeohigh a little, it is characterised in that including:
Obtain the Recent Structural diagram data waiting to calculate stratum, area；
Recent Structural isogram data are determined according to described Recent Structural diagram data；
The plan structure of mesh point between each well is determined according to described Recent Structural diagram data and described Recent Structural isogram data
Amplitude；Described plan Structural range is the vertical direction amplitude relative to Recent Structural trend surface of the point on Recent Structural figure；
Obtain Recent Structural figure degrades line data, and determines wait to calculate area according to the line data of degrading of described Recent Structural figure
The deposition phase structural trend surface diagram data on stratum；
The relative height number that the deposition phase constructs is determined with described plan Structural range according to described deposition phase structural trend surface diagram data
According to, and spend palaeohigh diagram data a little according to the data genaration deposition phase of height relatively of structure of described deposition phase；
Spend palaeohigh diagram data a little according to the described deposition phase and determine sedimentary micro diagram data.
The deposition phase the most according to claim 1 spends the projectional technique of palaeohigh a little, it is characterised in that according to described now
Structural map data determine Recent Structural isogram data, including:
The first data point and the second data point is obtained in the Recent Structural figure that Recent Structural diagram data is corresponding；Wherein, described
The coordinate of one data point is (x0, y0, z0), the coordinate of described second data point is (xn, yn, zn)；Described first data point is existing
Minimum point in modern structural map；Described second data point be in Recent Structural figure centered by peak, with predeterminable range be
Minimum point in the range of radius；
The straight line that first data point and the second data point connect presets discrepancy in elevation Δ Z every one and determines one the 3rd data point；Institute
The coordinate stating the 3rd data point is (xk, yk, zk)；
According to formula:
In determining space, (x, y, z), and carry out interpolation according to kriging analysis algorithm to point coordinates each in space to each point coordinates, raw
Become Recent Structural isogram data；At straight lineOn the height value of point be zk。
The deposition phase the most according to claim 2 spends the projectional technique of palaeohigh a little, it is characterised in that according to described now
Structural map data and described Recent Structural isogram data determine the plan Structural range of mesh point between each well, including:
By the coordinate of mesh point between each well in Recent Structural diagram data and each point coordinates in Recent Structural isogram data
In the vertical direction one_to_one corresponding；
Determine that the point in the Recent Structural isogram data that between each well in Recent Structural diagram data, mesh point is corresponding exists
Amplitude difference on vertical direction, and the plan Structural range of mesh point between each well is determined according to each amplitude difference.
The deposition phase the most according to claim 3 spends the projectional technique of palaeohigh a little, it is characterised in that obtain Recent Structural
Figure degrade line data, and determine wait to calculate the deposition phase structure on stratum, area according to the line data of degrading of described Recent Structural figure
Isogram data, including:
Any two points is chosen in line data corresponding degrading of degrading at Recent Structural figure on line, and with wherein as zero
Set up azimuthal coordinates system；The direct north of described azimuthal coordinates system is 0 °, and direction, due east is 90 °；Wherein, at described azimuthal coordinates
In system origin be (x '0, y '0), another point coordinates in described any two points be (x '1, y '1)；
Will from (x '0, y '0) point to (x '1, y '1) vector direction as degrade line trend and deposition phase structural trend surface trend；
Will from (x '0, y '0) point to (x '1, y '1) vector vertical direction as deposition phase structural trend surface tendency；
According to formula:
Determine each point coordinates on deposition phase structural trend surface (x ', y ', z ')；Wherein, (x'i,y'i,z'i) it was deposition phase structure
Certain point coordinates on trend surface；(x'j,y'j,z'j) for along certain seat put that deposition phase structural trend surface tendency predeterminable range is l '
Mark；α is the deposition phase structural trend surface inclination angle pre-set；θ is the folder degrading line trend with the direct north of azimuthal coordinates system
Angle；For crossing (x'j,y'j,z'j) point, and the straight line in direction is moved towards along deposition phase structural trend surface
Equation；zx”,y”=z'jFor crossing (x'j,y'j,z'j) put and along depositing the height value that phase structural trend surface moves towards the straight line in direction；
Deposition phase structural trend surface diagram data is constituted according to each point coordinates (x ', y ', z ') on deposition phase structural trend surface.
The deposition phase the most according to claim 4 spends the projectional technique of palaeohigh a little, it is characterised in that according to described deposition
Phase structural trend surface diagram data determines the relative height data that the deposition phase constructs with described plan Structural range, including: by each grid
The deposition phase structural trend surface diagram data of point is intended Structural range with it and is added, determines the relative just data that the deposition phase constructs.
6. one kind deposits the phase and spends the estimation device of palaeohigh a little, it is characterised in that including:
Recent Structural diagram data acquiring unit, for obtaining the Recent Structural diagram data waiting to calculate stratum, area；
Recent Structural isogram data determination unit, for determining Recent Structural trend surface according to described Recent Structural diagram data
Diagram data；
Intend Structural range and determine unit, for true according to described Recent Structural diagram data and described Recent Structural isogram data
The plan Structural range of mesh point between fixed each well；Described plan Structural range is that the point on Recent Structural figure is relative to Recent Structural trend
The vertical direction amplitude in face；
Deposition phase structural trend surface diagram data determines unit, degrades line data for obtain Recent Structural figure, and according to described
The line data of degrading of Recent Structural figure determine the deposition phase structural trend surface diagram data waiting to calculate stratum, area；
The deposition phase spends palaeohigh diagram data signal generating unit a little, for according to described deposition phase structural trend surface diagram data and described
Intend the data of height relatively that Structural range determines that the deposition phase construct, and according to the most just data genaration of structure of described deposition phase
The deposition phase spends palaeohigh diagram data a little；
Sedimentary microfacies map data determination unit, determines sedimentary microfacies map for spending palaeohigh diagram data a little according to the described deposition phase
Data.
The deposition phase the most according to claim 6 spends the estimation device of palaeohigh a little, it is characterised in that described Recent Structural
Isogram data determination unit, specifically for: in the Recent Structural figure that Recent Structural diagram data is corresponding, obtain the first data
Point and the second data point；Wherein, the coordinate of described first data point is (x0, y0, z0), the coordinate of described second data point is
(xn, yn, zn)；Described first data point is the minimum point in Recent Structural figure；Described second data point is in Recent Structural figure
Centered by peak, with predeterminable range as radius in the range of minimum point；
The straight line that first data point and the second data point connect presets discrepancy in elevation Δ Z every one and determines one the 3rd data point；Institute
The coordinate stating the 3rd data point is (xk, yk, zk)；
According to formula:
In determining space, (x, y, z), and carry out interpolation according to kriging analysis algorithm to point coordinates each in space to each point coordinates, raw
Become Recent Structural isogram data；At straight lineOn the height value of point be zk。
The deposition phase the most according to claim 7 spends the estimation device of palaeohigh a little, it is characterised in that described plan structure width
Degree determines unit, specifically for:
By the coordinate of mesh point between each well in Recent Structural diagram data and each point coordinates in Recent Structural isogram data
In the vertical direction one_to_one corresponding；
Determine that the point in the Recent Structural isogram data that between each well in Recent Structural diagram data, mesh point is corresponding exists
Amplitude difference on vertical direction, and the plan Structural range of mesh point between each well is determined according to each amplitude difference.
The deposition phase the most according to claim 8 spends the estimation device of palaeohigh a little, it is characterised in that deposition phase structure becomes
Gesture face diagram data determines unit, specifically for:
Any two points is chosen in line data corresponding degrading of degrading at Recent Structural figure on line, and with wherein as zero
Set up azimuthal coordinates system；The direct north of described azimuthal coordinates system is 0 °, and direction, due east is 90 °；Wherein, at described azimuthal coordinates
In system origin be (x '0, y '0), another point coordinates in described any two points be (x '1, y '1)；
Will from (x '0, y '0) point to (x '1, y '1) vector direction as degrade line trend and deposition phase structural trend surface trend；
Will from (x '0, y '0) point to (x '1, y '1) vector vertical direction as deposition phase structural trend surface tendency；
According to formula:
Determine each point coordinates on deposition phase structural trend surface (x ', y ', z ')；Wherein, (x'i,y'i,z'i) it was deposition phase structure
Certain point coordinates on trend surface；(x'j,y'j,z'j) for along certain seat put that deposition phase structural trend surface tendency predeterminable range is l '
Mark；α is the deposition phase structural trend surface inclination angle pre-set；θ is the folder degrading line trend with the direct north of azimuthal coordinates system
Angle；For crossing (x'j,y'j,z'j) point, and the straight line in direction is moved towards along deposition phase structural trend surface
Equation；zx”,y”=z'jFor crossing (x'j,y'j,z'j) put and along depositing the height value that phase structural trend surface moves towards the straight line in direction；
Deposition phase structural trend surface diagram data is constituted according to each point coordinates (x ', y ', z ') on deposition phase structural trend surface.
The deposition phase the most according to claim 9 spends the estimation device of palaeohigh a little, it is characterised in that the described deposition phase
Degree palaeohigh diagram data signal generating unit a little, specifically for: the deposition phase structural trend surface diagram data of each mesh point is intended with it
Structural range is added, and determines the most just data that the deposition phase construct.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610213582.0A CN105868482B (en) | 2016-04-07 | 2016-04-07 | A kind of deposition phase spends the projectional technique and device of palaeohigh a little |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610213582.0A CN105868482B (en) | 2016-04-07 | 2016-04-07 | A kind of deposition phase spends the projectional technique and device of palaeohigh a little |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105868482A true CN105868482A (en) | 2016-08-17 |
| CN105868482B CN105868482B (en) | 2018-09-07 |
Family
ID=56636986
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610213582.0A Active CN105868482B (en) | 2016-04-07 | 2016-04-07 | A kind of deposition phase spends the projectional technique and device of palaeohigh a little |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105868482B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107315195A (en) * | 2017-05-31 | 2017-11-03 | 中国石油天然气股份有限公司 | The determination method and apparatus of microrelief |
| CN109870719A (en) * | 2019-01-25 | 2019-06-11 | 中国石油天然气集团有限公司 | A kind of carbonate rock thin tight reservoir distribution determination method, apparatus and system |
| CN110119578A (en) * | 2019-05-16 | 2019-08-13 | 中国人民解放军海军工程大学 | Seabed trend surface building method, seabed Trend Prediction method and device |
| CN111260100A (en) * | 2018-11-30 | 2020-06-09 | 中国石油天然气股份有限公司 | Reservoir production capacity prediction method and device |
| CN111411947A (en) * | 2020-05-14 | 2020-07-14 | 中国石油天然气集团有限公司 | Method for predicting formation pressure based on HT L parameter |
| CN113970796A (en) * | 2020-07-23 | 2022-01-25 | 中国石油化工股份有限公司 | Method for accurately recovering ancient water depth of sedimentary basin |
| US11542816B2 (en) | 2021-01-28 | 2023-01-03 | Saudi Arabian Oil Company | System and method for hydrocarbon exploration based on imaging tunnel valleys |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2183332C1 (en) * | 2000-11-08 | 2002-06-10 | Зубков Михаил Юрьевич | Method for prediction of development zones of secondary fractured reservoirs in sedimentary mantle |
| CN105137482A (en) * | 2015-09-09 | 2015-12-09 | 中国海洋石油总公司 | Method for calculating ancient gradient of sedimentary body |
| CN105425315A (en) * | 2015-12-21 | 2016-03-23 | 中国石油大学(北京) | Inversion method of microscale deposition stage palaeogeomorphology in weak construction deformation range |
| CN105426316A (en) * | 2015-11-09 | 2016-03-23 | 北京大学 | Quota control temperature based racetrack memory chip and control method therefor |
-
2016
- 2016-04-07 CN CN201610213582.0A patent/CN105868482B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2183332C1 (en) * | 2000-11-08 | 2002-06-10 | Зубков Михаил Юрьевич | Method for prediction of development zones of secondary fractured reservoirs in sedimentary mantle |
| CN105137482A (en) * | 2015-09-09 | 2015-12-09 | 中国海洋石油总公司 | Method for calculating ancient gradient of sedimentary body |
| CN105426316A (en) * | 2015-11-09 | 2016-03-23 | 北京大学 | Quota control temperature based racetrack memory chip and control method therefor |
| CN105425315A (en) * | 2015-12-21 | 2016-03-23 | 中国石油大学(北京) | Inversion method of microscale deposition stage palaeogeomorphology in weak construction deformation range |
Non-Patent Citations (1)
| Title |
|---|
| 林畅松 等: "塔里木盆地古生代中央隆起带古构造地貌及其对沉积相发育分布的制约", 《中国科学 D 辑:地球科学》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107315195A (en) * | 2017-05-31 | 2017-11-03 | 中国石油天然气股份有限公司 | The determination method and apparatus of microrelief |
| CN107315195B (en) * | 2017-05-31 | 2019-05-07 | 中国石油天然气股份有限公司 | The determination method and apparatus of microrelief |
| CN111260100A (en) * | 2018-11-30 | 2020-06-09 | 中国石油天然气股份有限公司 | Reservoir production capacity prediction method and device |
| CN109870719A (en) * | 2019-01-25 | 2019-06-11 | 中国石油天然气集团有限公司 | A kind of carbonate rock thin tight reservoir distribution determination method, apparatus and system |
| CN109870719B (en) * | 2019-01-25 | 2021-01-29 | 中国石油天然气集团有限公司 | Well position layout method, device and system for carbonate rock tight thin reservoir |
| CN110119578A (en) * | 2019-05-16 | 2019-08-13 | 中国人民解放军海军工程大学 | Seabed trend surface building method, seabed Trend Prediction method and device |
| CN111411947A (en) * | 2020-05-14 | 2020-07-14 | 中国石油天然气集团有限公司 | Method for predicting formation pressure based on HT L parameter |
| CN111411947B (en) * | 2020-05-14 | 2023-09-15 | 中国石油天然气集团有限公司 | Method for predicting formation pressure based on HTL parameters |
| CN113970796A (en) * | 2020-07-23 | 2022-01-25 | 中国石油化工股份有限公司 | Method for accurately recovering ancient water depth of sedimentary basin |
| US11542816B2 (en) | 2021-01-28 | 2023-01-03 | Saudi Arabian Oil Company | System and method for hydrocarbon exploration based on imaging tunnel valleys |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105868482B (en) | 2018-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105868482B (en) | A kind of deposition phase spends the projectional technique and device of palaeohigh a little | |
| CN107526101B (en) | A kind of acquisition and processing method obtaining earthquake reflected wave | |
| CN104765065B (en) | Method and device for determining sedimentary stratum ancient landform | |
| CN105425315B (en) | A kind of small scale depositional phase ancient landform inversion method in weak structural deformation area | |
| CN102243678A (en) | Method for analyzing sand bodies in reservoirs based on inversion technique of sedimentary dynamics | |
| CN103293553B (en) | Upper and lower cable seismic acquisition data boundary element continuation bearing calibration at the bottom of a kind of Complex Sea | |
| US11163093B2 (en) | Constructing stratigraphic images of sediment distribution in a subsurface of the earth | |
| CN105093307A (en) | Lower palaeozoic tilted stratum true thickness calculation method | |
| Playton et al. | Architecture and genesis of prograding deep boundstone margins and debris-dominated carbonate slopes: Examples from the Permian Capitan Formation, Southern Guadalupe Mountains, West Texas | |
| Sacchi et al. | The use and beauty of ultra-high-resolution seismic reflection imaging in Late Quaternary marine volcaniclastic settings, Napoli Bay, Italy | |
| Yarus et al. | Facies simulation in practice: Lithotype proportion mapping and Plurigaussian simulation, a powerful combination | |
| Rush et al. | Geostatistical facies modeling trends for oolitic tidal sand shoals | |
| Vittori et al. | The Quaternary Congo deep-sea fan: Preliminary results on reservoir complexity in turbiditic systems using 2D high resolution seismic and multibeam data | |
| CN109085644A (en) | True earth's surface imaging method when being walked based on dual-beam | |
| Di Fiore et al. | Joining up land and marine seismic data: Case study from Procida and Ventotene Islands (Tyrrhenian Sea, Italy) | |
| Leslie et al. | Internal structure of Puna Ridge: evolution of the submarine East Rift Zone of Kilauea Volcano, Hawai ̀i | |
| Punzo et al. | Ensuring very shallow-water sediment properties: case study from Capo Granitola harbour, Sicily (Italy) | |
| Hou et al. | Analysis of Sandbody Development Model in Lishui Depression | |
| Alotaibi | Seismic Structural Investigation and Reservoir Characterization of the Moki Formation in Maari Field, Taranaki Basin, New Zealand | |
| Huang et al. | Seismic facies in the Pearl River Submarine Canyon, northern South China Sea | |
| Marini et al. | Multi-scale characterization of the Pleistocene-Holocene Tiber delta deposits as a depositional analogue for hydrocarbon reservoirs | |
| Javadova | EARTH SCIENCES | |
| CN104932013A (en) | Delta sub-phase identification method based on earthquake high-frequency cycle thickness | |
| Karlsen | 3D Subsurface Modeling based on Ground Penetrating Radar Survey of Glaciomarine Deposits–Hagadrag Aquifer, Norway. A sedimentological, hydrogeological and geophysical study | |
| Yu et al. | Theory and Methods for Studying Clastic Sequence Stratigraphy |
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 |