CN108107469A - Distributary channel method for determining position and device - Google Patents

Distributary channel method for determining position and device Download PDF

Info

Publication number
CN108107469A
CN108107469A CN201711145503.8A CN201711145503A CN108107469A CN 108107469 A CN108107469 A CN 108107469A CN 201711145503 A CN201711145503 A CN 201711145503A CN 108107469 A CN108107469 A CN 108107469A
Authority
CN
China
Prior art keywords
seismic
well
distributary channel
strata
distributary
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
Application number
CN201711145503.8A
Other languages
Chinese (zh)
Other versions
CN108107469B (en
Inventor
黄军平
刘震华
杨占龙
洪忠
向倪娇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201711145503.8A priority Critical patent/CN108107469B/en
Publication of CN108107469A publication Critical patent/CN108107469A/en
Application granted granted Critical
Publication of CN108107469B publication Critical patent/CN108107469B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6161Seismic or acoustic, e.g. land or sea measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/40Protecting water resources
    • Y02A20/402River restoration

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)

Abstract

The present invention provides a kind of distributary channel method for determining position and device, wherein, this method includes:Obtain the drilling well log of purpose series of strata;According to the drilling well log, the well logging phase character of distributary channel is determined;The position of distributary channel is calibrated on seismic profile according to the well logging phase character, obtains the characteristics of seismic of distributary channel;Seismic waveform inversion is carried out based on the log, obtains the position of distributary channel in purpose series of strata;Tectonic interpretation is carried out to the purpose series of strata according to the characteristics of seismic, to recover to obtain micro- ancient landform;The position of the distributary channel obtained based on the characteristics of seismic, inverting and micro- ancient landform, obtain the spread of distributary channel in the purpose series of strata.The present invention solves the existing distribution characteristics that can not be simple and efficient and identify distributary channel in shallow-circulation flow, and reservoir is caused to teach difficult the technical issues of predicting, has reached the technique effect for being simple and efficient predicting reservoir.

Description

Distributary channel method for determining position and device
Technical field
The present invention relates to technical field of geological exploration, more particularly to a kind of distributary channel method for determining position and device.
Background technology
With the continuous progress of oil-gas exploration technology, the discovery of structural deposit is more and more difficult, and the reserves of discovery are more next It is smaller.Under this trend, the exploration of lithologic deposit is more and more important.In recent years, in each oil field, it is found that in succession very More lithologic deposits, especially some oil and gas productions can reach 50,000,000 tons or more of annual output.In the past for breaking fold belt development Fault depressed lacustrine basin, can be by identifying Slope break belt, and then predicts the Distribution Characteristics of Reservoir Body.However, for breaking fold belt agensis Shallow-water lacustrine basin, can be by predicting delta distribution, and then predict the Distribution Characteristics of Reservoir Body.
However, for the place of developement of the delta, how by high-precision three-dimensional whether Reservoir Body centainly develop and Reservoir Body distribution characteristics in shake data meticulous depiction shallow-circulation flow is required for enough prospect pit quantity and foot at present Enough seismic datas could be realized.However, in a practical situation, often degree of prospecting than it is relatively low, prospect pit is less, and earthquake provides Material is also inadequate.
Based on existing existing research and exploration data, shallow water triangle is had found in Paleozoic group, the Mesozoic group and Cenozoic Continent deposits, and more with terrestrial facies lake shallow-circulation flow.Wherein, shallow-circulation flow is (1m-6m) shallow in water body, landform it is gentle (< 0.1 °) position formed the deltaic reservoirs based on distributary channel.The Reservoir Body formed under this sedimentation setting The difficult point of Geophysical recognition essentially consists in the following:
1) ancient lake landform is relatively flat, the more difficult feature for identifying ancient stream channel on seismic profile;
2) sandstone development, single sand body is longitudinally stacked on top of each other and thickness change is larger, and lateral migration is fast, and intersection is compound, stablizes Property is poor;
3) net sandstone thinner thickness, vertical upper distribution is more scattered, and cross directional variations are big;
4) reservoir change is big, and anisotropism is strong, and hydrocarbon distribution is strictly by sand body morphology and physical characteristics control.
For such paleotopography feature, stored up using three dimensional seismic data it is difficult to accurately depict in shallow-circulation flow The Distribution Characteristics of collective.
In view of the above-mentioned problems, currently no effective solution has been proposed.
The content of the invention
An embodiment of the present invention provides a kind of distributary channel method for determining position, and Reservoir Body is depicted with what is be simple and efficient Distribution Characteristics, this method includes:
Obtain the drilling well log of purpose series of strata;
According to the drilling well log, the well logging phase character of distributary channel is determined;
The position of distributary channel is calibrated on seismic profile according to the well logging phase character, obtains the earthquake of distributary channel Phase character;
Seismic waveform inversion is carried out based on the log, obtains the position of distributary channel in purpose series of strata;
Tectonic interpretation is carried out to the purpose series of strata according to the characteristics of seismic, to recover to obtain micro- ancient landform;
The position of the distributary channel obtained based on the characteristics of seismic, inverting and micro- ancient landform, obtain the mesh Series of strata in distributary channel spread.
In one embodiment, the position of distributary channel is being demarcated on seismic profile according to the well logging phase character, The characteristics of seismic of distributary channel is obtained, including:
According to the well logging phase character of distributary channel, the well logging phase character type of distributary channel is determined;
According to the well logging phase character, distributary channel is demarcated onto corresponding seismic profile, to identify sympathetic earthquakes The characteristics of seismic of distributary channel on section.
In one embodiment, after recovering to obtain micro- ancient landform, further include:
Ask for the compacting factor of every mouthful of well in micro- ancient landform;
Planar interpolation is carried out to the compacting factor of every mouthful of well, obtains compacting factor plan view;
Micro- ancient landform with the compacting factor plan view is multiplied, obtains micro- ancient landform after compaction correction.
In one embodiment, the compacting factor of every mouthful of well in micro- ancient landform is asked for, including:
By porosity depth compaction model, ask for being compacted the compacting factor of every mouthful of well purpose series of strata:
φ (z)=φ0e-cz
Wherein, φ0Represent initial porosity, c represents compacting factor, and z represents depth.
In one embodiment, seismic waveform inversion is carried out based on the log, obtains shunting in purpose series of strata The position in river, including:
Load the seismic data and well-log information of the purpose series of strata;
The seismic data and the well-log information are demarcated;
The related well sample sheet of selection, sample for reference spatial distribution distance and curve distribution feature establish initial model;
Selected in the preset range of well point with the most similar seismic waveform of indicatrix, using global optimization approach, The seismic waveform of selection is generalized to no wellblock, the company's of obtaining well seismic profile;
By the even well seismic profile compared with characterizing the log value of distributary channel and carrying out repeatedly, until being met The geological model of preset requirement;
The geological model of foundation is applied to entire seismic data cube, to realize seismic waveform inversion, obtains purpose The position of distributary channel in series of strata.
The embodiment of the present invention additionally provides a kind of determining device of distributary channel position, is preserved with depicting of being simple and efficient The Distribution Characteristics of body, the device include:
Acquisition module, for obtaining the drilling well log of purpose series of strata;
Second determining module, for according to the drilling well log, determining the well logging phase character of distributary channel;
Demarcating module for calibrating the position of distributary channel on seismic profile according to the well logging phase character, obtains The characteristics of seismic of distributary channel;
Inverting module carries out seismic waveform inversion for being based on the log, obtains distributary channel in purpose series of strata Position;
Recovery module, for carrying out tectonic interpretation to the purpose series of strata according to the characteristics of seismic, to recover Obtain micro- ancient landform;
Second determining module, for the position of distributary channel that is obtained based on the characteristics of seismic, inverting and described micro- Ancient landform obtains the spread of distributary channel in the purpose series of strata.
In one embodiment, the demarcating module includes:
Determination unit for the well logging phase character according to distributary channel, determines the well logging phase character type of distributary channel;
Unit is demarcated, for according to the well logging phase character, distributary channel being demarcated onto corresponding seismic profile, to know Do not go out the characteristics of seismic of distributary channel on sympathetic earthquakes section.
In one embodiment, above device further includes:
Module is asked for, for after recovering to obtain micro- ancient landform, asking for the compacting system of every mouthful of well in micro- ancient landform Number;
Interpolating module for carrying out planar interpolation to the compacting factor of every mouthful of well, obtains compacting factor plan view;
Processing module, for micro- ancient landform to be multiplied with the compacting factor plan view, after obtaining compaction correction Micro- ancient landform.
In one embodiment, the module of asking for is specifically used for, by porosity depth compaction model, asking for being compacted The compacting factor of every mouthful of well purpose series of strata:
φ (z)=φ0e-cz
Wherein, φ0Represent initial porosity, c represents compacting factor, and z represents depth.
In one embodiment, stating inverting module includes:
Loading unit, for loading the seismic data and well-log information of the purpose series of strata;
Unit is demarcated, for being demarcated to the seismic data and the well-log information;
Unit is established, for selecting related well sample sheet, sample for reference spatial distribution distance and curve distribution feature are established just Beginning model;
Unit is promoted, for being selected in the preset range of well point and the most similar seismic waveform of indicatrix, use The seismic waveform of selection is generalized to no wellblock, the company's of obtaining well seismic profile by global optimization approach;
Comparing unit, for the even well seismic profile repeatedly to be compared with characterizing the log value of distributary channel Compared with until obtaining the geological model for meeting preset requirement;
Applying unit, the geological model for that will establish is applied to entire seismic data cube, to realize seismic waveform Inverting obtains the position of distributary channel in purpose series of strata.
In embodiments of the present invention, organically combined by the phase that will log well, seismic facies, seismic waveform inversion three, it is accurate Distribution of the distributary channel on seismic profile really is had identified, further, is had identified point using micro- ancient landform analytical technology The position of fluid stream channel and distribution characteristics identify distributary channel in shallow-circulation flow so as to solve existing can not be simple and efficient Distribution characteristics, and the technical issues of cause reservoir religion difficult prediction, reached the technique effect for being simple and efficient predicting reservoir.
Description of the drawings
Attached drawing described herein is used for providing a further understanding of the present invention, forms the part of the application, not Form limitation of the invention.In the accompanying drawings:
Fig. 1 is the method flow diagram of distributary channel method for determining position according to embodiments of the present invention;
Fig. 2 is the another method flow chart of distributary channel method for determining position according to embodiments of the present invention;
Fig. 3 is the well logging facies type schematic diagram of distributary channel according to embodiments of the present invention;
Fig. 4 is the seismic facies identification schematic diagram for the L1 well distributary channels that top circle according to embodiments of the present invention evens up;
Fig. 5 is the seismic facies identification schematic diagram for the L1 well distributary channels that top circle according to embodiments of the present invention does not even up;
Fig. 6 is the characteristics of seismic schematic diagram of distributary channel according to embodiments of the present invention;
Fig. 7 is L2 well seismic waveform inversion sections excessively according to embodiments of the present invention;
Fig. 8 is the ancient landform schematic diagram before compaction correction according to embodiments of the present invention;
Fig. 9 is the ancient landform schematic diagram after compaction correction according to embodiments of the present invention;
Figure 10 is ancient landform and distributary channel three-dimensional overlay display schematic diagram after compaction correction according to embodiments of the present invention;
Figure 11 is another structure diagram of the determining device of distributary channel position according to embodiments of the present invention.
Specific embodiment
Understand to make the object, technical solutions and advantages of the present invention clearer, it is right with reference to embodiment and attached drawing The present invention is described in further details.Here, the exemplary embodiment and its explanation of the present invention be for explaining the present invention, but simultaneously It is not as a limitation of the invention.
In this example, it is contemplated that can distributary channel be identified, with reference to the distributary river of well-log facies recognition by electrofacies analysis Road position on seismic profile, with reference to seismic facies and seismic waveform inversion technology identification distributary channel position and distribution characteristics, and It can recover ancient landform by compaction correction using three dimensional seismic data, it is genuine that front distributary river combined under ancient landform background Seismic phase and inverting profile features are accurately identified distributary channel position.
Based on this, in this example, a kind of distributary channel method for determining position is provided, as shown in Figure 1, can include such as Lower step:
Step 101:Obtain the drilling well log of purpose series of strata;
Step 102:According to the drilling well log, the well logging phase character of distributary channel is determined;
Specifically, it according to the well logging phase character of distributary channel, can determine the well logging phase character type of distributary channel;Root According to the well logging phase character, by distributary channel calibration to corresponding seismic profile, shunted with identifying on sympathetic earthquakes section The characteristics of seismic in river.
I.e., it is possible to according to the electrofacies analysis to distributary channel, cutd open by fine synthesis record calibration to corresponding earthquake On face, and then it can identify distributary channel characteristics of seismic on sympathetic earthquakes section.It demarcates several mouthfuls of prospect pits more, can will summarize The well logging facies type of distributary channel be converted to earthquake facies type.
Step 103:The position of distributary channel is calibrated on seismic profile according to the well logging phase character, obtains distributary river The characteristics of seismic in road;
Step 104:Seismic waveform inversion is carried out based on the log, obtains the position of distributary channel in purpose series of strata It puts;
Step 105:Tectonic interpretation is carried out to the purpose series of strata according to the characteristics of seismic, it is micro- to recover to obtain Ancient landform;
After recovering to obtain micro- ancient landform, amount of compaction correction can be carried out to it.For example, it can recover to obtain micro- Gu After landforms, the compacting factor of every mouthful of well in micro- ancient landform is asked for;Planar interpolation is carried out to the compacting factor of every mouthful of well, is obtained To compacting factor plan view;Micro- ancient landform with the compacting factor plan view is multiplied, obtains micro- Gu after compaction correction Landforms.
Wherein, compacting factor can be asked for by porosity depth compaction model:
φ (z)=φ0e-cz
Wherein, φ0Represent initial porosity, c represents compacting factor, and z represents depth.
Step 106:The position of the distributary channel obtained based on the characteristics of seismic, inverting and micro- ancient landform, are obtained The spread of distributary channel into the purpose series of strata.
In above-mentioned steps 104, seismic waveform inversion is carried out based on log, obtains distributary channel in purpose series of strata Position can include:
S1:Load the seismic data and well-log information of the purpose series of strata;
S2:The seismic data and the well-log information are demarcated;
S3:The related well sample sheet of selection, sample for reference spatial distribution distance and curve distribution feature establish initial model;
S4:Selected in the preset range of well point with the most similar seismic waveform of indicatrix, calculated using global optimization The seismic waveform of selection is generalized to no wellblock, the company's of obtaining well seismic profile by method;
S5:By the even well seismic profile compared with characterizing the log value of distributary channel and carrying out repeatedly, until obtaining Meet the geological model of preset requirement;
S6:The geological model of foundation is applied to entire seismic data cube, to realize seismic waveform inversion, obtains mesh Series of strata in distributary channel position.
The above method is illustrated with reference to a specific embodiment, it should be noted, however, that the specific implementation Example does not form the improper restriction to the application merely to the application is better described.
In this example, in a manner that micro- paleogeographic reconstruction, seismic waveform inversion and electrofacies analysis are combined, three-dimensional is utilized Seismic data can effectively identify the position that distributary channel is developed in shallow-circulation flow, so as to promote under the sedimentation setting Effective identification of Reservoir Body and the further exploration of lithologic deposit, solve and are difficult to effectively know by conventional reservoir Predicting Technique The problem of position that distributary channel is developed under not micro- ancient landform.
Specifically, it as shown in Fig. 2, may include steps of:
Step 201:The well-log facies recognition of distributary channel.
It is that distributary channel is the most directly recognized in view of drilling well, therefore the well logging by meeting purpose series of strata drilling well to boring Curve is analyzed, and the well logging phase character of the distributary channel in certain depth bounds at the well point can be determined on well, then, can To sum up the well logging phase character type of distributary channel based on the well logging phase character determined.
Step 202:Seismic facies and seismic waveform inversion synthesis identification distributary channel.
According to the calibration of fine synthesis record to the electrofacies analysis of distributary channel, can be passed through in above-mentioned steps 201 to phase On the seismic profile answered, and then it can identify the characteristics of seismic of distributary channel on corresponding seismic profile.
For example, several mouthfuls of prospect pits can be demarcated, then can be by the well logging facies type for the distributary channel summarized in step 201 more Be converted to earthquake facies type;Secondly, the more apparent seismic profile of several distributary channel characteristics of seismic can be selected, with ground Seismic wave shape inversion technique obtains distributary channel inverting section.
The purpose of inverting is primarily due to, using well-log information, can only obtain the distribution of longitudinal distributary channel, and Quantity is fewer.However, the position of more distributary channel can be obtained based on seismic inversion.Therefore, in this example, it is to combine The position for the distributary channel that inverting determines is as an accurately reference, to cause the spread of definite distributary channel more Accurately.
Wherein, seismic waveform inversion in the horizontal, can make full use of the seismic waveform information that spatial distribution is intensive. On longitudinal direction, high-resolution drilling information can be made full use of.Seismic waveform is densely distributed space structure data, reflection The spatial variations of depositional environment and lithology combination.Seismic waveform inversion is usually preferably to go out correlation using seismic waveform similitude Well sample sheet, sample for reference spatial distribution distance and curve distribution feature establish initial model, instead of variogram analysis space Variant structure carries out radio-frequency component the optimal estimation of unbiased.
Specifically, the inversion method can be divided into following steps:
S1:Loading earthquake, well-log information are simultaneously demarcated.
S2:It is preferred that related well sample sheet, sample for reference spatial distribution distance and curve distribution feature establish initial model.
S3:Preferably go out near well point with the most similar seismic waveform of indicatrix, using global optimization approach, promote To no wellblock, to enhance inverting certainty, realize from completely random to progressively definite transformation.By company's well seismic profile of inverting Comparison is repeated in log value with characterizing distributary channel, suitable geological model is obtained until establishing, by the ground of foundation Matter model is applied to entire seismic data cube, to achieve the purpose that inverting distributary channel.
Step 203:Micro- paleogeographic reconstruction technology identification distributary channel position and distribution characteristics.
It, can be on three dimensional seismic data to shunting after seismic facies identification by distributary channel in above-mentioned steps 202 The top bottom of series of strata residing for river carries out Fine structural interpretation, then, then does micro- paleogeographic reconstruction.
Palaeotectonic setting in view of shallow-circulation flow development be landform it is gentle (<0.1 °) and the gradient it is smaller, therefore, can be with Compaction correction is done to ready-made ancient landform so that the smaller distributary channel of the gradient can also be embodied as out.For through overcompaction Micro- ancient landform after amount correction, in conjunction with the characteristics of seismic of distributary channel, can preferably identify distributary channel position.In order to To preferably be compared before compaction correction with micro- ancient landform after compaction correction, can will before compaction correction and compaction correction after Two micro- ancient landforms be standardized, so as to unified dimension.
Wherein, the principle of amount of compaction correction is as follows:
Research for Clastic Stratum of Country Rocks compaction correction, be mostly based on porosity with the empirical index number model of change in depth or Other porosity-depth models.Different lithology, compacted property are also different in Clastic Stratum of Country Rocks.It specifically, can be 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.For example, it is special to characterize compacting that purpose series of strata compacting factor can be asked for by using log Sign.
Specifically, compacting factor can be asked for by the following steps:
S1:First, the compacting factor of every mouthful of well is asked for.Such as:According to the log data of individual well, by the chip residing for individual well Rock STRATIGRAPHIC DIVISION is compacted unit at least two different lithology, passes through porosity depth compaction model:φ (z)=φ0e-czIt asks The compacting factor of pressure every mouthful of well purpose series of strata in fact, wherein, φ0Represent initial porosity, c represents compacting factor, and z represents deep Degree.
S2:Planar interpolation is carried out to the compacting factor of every mouthful of individual well, specifically, Kriging regression method can be utilized to each The compacting factor of a individual well carries out planar interpolation, and after interpolation completion, the compacting factor after interpolation and research area can sink Product background compares and analyzes, to determine whether the plane distribution of compacting factor meets research area's geologic setting, otherwise, then Think to modify to compacting factor.
S3:The micro- ancient landform recovered before compaction correction with the compacting factor plan view asked for is multiplied, obtains compaction correction Micro- ancient landform afterwards.
It is described as follows with reference to specific experiment scene:
S1:Well logging recognition distributary channel.
The well logging facies type of distributary channel is illustrated in figure 3, the well logging that distributary channel can be summed up by the figure is mutually special Sign, i.e.,:Bell, box-shaped, tooth box-shaped and complex.Pass through the analysis for phase of logging well to a large amount of drilling wells of research area's purpose series of strata, hair In present shallow-circulation flow, distributary channel is developed the most, is secondly deposited between development distributary channel, river mouth bar is less, gravity Stream deposition is less.
S2:Seismic facies and seismic waveform inversion synthesis identification distributary channel.
Electrofacies analysis through step S1, is recorded by fine synthesis, and correspondence can be preferably demarcated on seismic profile The location of distributary channel.The earthquake of distributary channel on the seismic profile for illustrating L1 wells shown in Fig. 4 and Fig. 5 preferably Phase character evens up purpose series of strata top circle, from Fig. 4 in, can preferably identify the seismic facies of distributary channel, i.e. earthquake A series of micro- ancient landforms on section, similar to Downcutting valley Seismic reflection character;Then, if purpose series of strata top circle not drawn It is flat, it is from Fig. 5, then more difficult to identify on seismic profile that ancient landform is exactly distributary channel slightly for this.According to these contents, in Fig. 6 It is middle purpose series of strata top circle is evened up after, distributary channel can be preferably identified, according to shape of the distributary channel on seismic profile State can be divided into 3 type of U-shaped, V-type and L-type.
The distributary channel gone out according to well-log facies recognition,, can be fine using seismic waveform inversion technology with reference to seismic data Ground depicts feature of the distributary channel on inverting section, as can be seen from Figure 7 the local sandstone of these distributary channels development Reservoir Body is more developed.
S3:Micro- paleogeographic reconstruction identification distributary channel position and distribution characteristics.
It, can be on three dimensional seismic data to distributary channel institute after seismic facies identification by the distributary channel of step S2 The top bottom for locating series of strata carries out Fine structural interpretation, and variation is spent in particularly purpose series of strata bottom a little, and it is special can to combine seismic facies well Sign identifies, then does micro- paleogeographic reconstruction.
Due to the palaeotectonic setting of shallow-circulation flow development be landform it is gentle (<0.1 °), the gradient is smaller.Compacting factor it is big The small feature that can reflect rich sand and landforms under rich mud background well, therefore, it should which compaction correction is done to ancient landform.For warp The micro- ancient landform recovered after over-pressed actual quantities correction, with reference to the characteristics of seismic of distributary channel, can preferably identify distributary river The position in road and distribution characteristics.It, can be in order to will preferably be compared before compaction correction with micro- ancient landform after compaction correction Micro- ancient landform is standardized, unified dimension.As shown in figure 8, it is the ancient landform before compaction correction, as shown in figure 9, being Ancient landform after compaction correction.As shown in Figure 10, it is ancient landform after compaction correction and distributary channel three-dimensional overlay display schematic diagram.
The method provided by upper example solves what distributary channel in shallow-circulation flow was difficult on seismic profile Problem.Specifically, being combined by phase of logging well, seismic facies, seismic waveform inversion, distributary channel has been accurately identified in seismic profile On distribution, position and the distribution characteristics of distributary channel are identified using micro- ancient landform analytical technology, so as to be recognized accurately The distribution characteristics of distributary channel (also crying ancient stream channel) in shallow-circulation flow, and then solve the problems such as reservoir prediction is difficult.
The technology application has been obtained into good application effect, has played preferable production actual effect.By to Old Well Reexamination, It is found that the more active drilling well of oil/gas show is all got on distributary channel substantially, the poor drilling well of oil/gas show has all been got into point Fluid stream channel edge.
In upper example, known using electrofacies analysis technology, seismic facies technology, seismic inversion and micro- paleogeographic reconstruction technology The position of distributary channel is not gone out.Specifically, having identified distributary channel by electrofacies analysis technology, combine seismic facies technology And seismic inversion has identified distributary channel position;It is special that distributary channel distribution has been accurately identified by micro- paleogeographic reconstruction It seeks peace development degree.It is influenced by seismic data quality, the more difficult identification distributary channel position on seismic profile, comprehensive logging phase, Seismic facies and seismic waveform inversion technology identification distributary channel position, after compaction correction, with micro- paleogeographic reconstruction technology Distributary channel position and distribution characteristics can be accurately identified.
Based on same inventive concept, a kind of determining device of distributary channel position is additionally provided in the embodiment of the present invention, such as Described in the following examples.The principle solved the problems, such as due to the determining device of distributary channel position is determined with distributary channel position Method is similar, therefore the implementation of the determining device of distributary channel position may refer to the reality of distributary channel method for determining position It applies, overlaps will not be repeated.Used below, term " unit " or " module " can realize the software of predetermined function And/or the combination of hardware.Although following embodiment described device is preferably realized with software, hardware or soft The realization of the combination of part and hardware is also what may and be contemplated.Figure 11 is determining for the distributary channel position of the embodiment of the present invention A kind of structure diagram of device as shown in figure 11, can include:Acquisition module 1101, the first determining module 1102, demarcating module 1103rd, inverting module 1104,1105 and second determining module 1106 of recovery module, below illustrate the structure.
Acquisition module 1101, for obtaining the drilling well log of purpose series of strata;
First determining module 1102, for according to the drilling well log, determining the well logging phase character of distributary channel;
Demarcating module 1103, for the position of distributary channel to be calibrated on seismic profile according to the well logging phase character, Obtain the characteristics of seismic of distributary channel;
Inverting module 1104 carries out seismic waveform inversion for being based on the log, obtains shunting in purpose series of strata The position in river;
Recovery module 1105, for carrying out tectonic interpretation to the purpose series of strata according to the characteristics of seismic, with Recovery obtains micro- ancient landform;
Second determining module 1106, for the position of distributary channel obtained based on the characteristics of seismic, inverting and institute Micro- ancient landform is stated, obtains the spread of distributary channel in the purpose series of strata.
In one embodiment, demarcating module 1103 can include:Determination unit, for the well logging according to distributary channel Phase character determines the well logging phase character type of distributary channel;Unit is demarcated, for according to the well logging phase character, will shunt River is demarcated onto corresponding seismic profile, to identify the characteristics of seismic of distributary channel on sympathetic earthquakes section.
In one embodiment, above device can also include:Ask for module, for recover obtain micro- ancient landform it Afterwards, the compacting factor of every mouthful of well in micro- ancient landform is asked for;Interpolating module, for carrying out plane to the compacting factor of every mouthful of well Interpolation obtains compacting factor plan view;Processing module, for micro- ancient landform to be multiplied with the compacting factor plan view, Obtain micro- ancient landform after compaction correction.
In one embodiment, module is asked for specifically to can be used for, by porosity depth compaction model, asking for being compacted The compacting factor of every mouthful of well purpose series of strata:
φ (z)=φ0e-cz
Wherein, φ0Represent initial porosity, c represents compacting factor, and z represents depth.
In one embodiment, inverting module 1104 can include:Loading unit, for loading the purpose series of strata Seismic data and well-log information;Unit is demarcated, for being demarcated to the seismic data and the well-log information;It establishes single Member, for selecting related well sample sheet, sample for reference spatial distribution distance and curve distribution feature establish initial model;It promotes single Member, for selected in the preset range of well point with the most similar seismic waveform of indicatrix, will using global optimization approach The seismic waveform of selection is generalized to no wellblock, the company's of obtaining well seismic profile;Comparing unit, for by it is described even well seismic profile with The log value of characterization distributary channel is repeatedly compared, until obtaining the geological model for meeting preset requirement;Applying unit, The geological model for that will establish is applied to entire seismic data cube, to realize seismic waveform inversion, obtains purpose series of strata The position of middle distributary channel.
In another embodiment, a kind of software is additionally provided, which is used to perform above-described embodiment and preferred reality Apply the technical solution described in mode.
In another embodiment, a kind of storage medium is additionally provided, above-mentioned software is stored in the storage medium, it should Storage medium includes but not limited to:CD, floppy disk, hard disk, scratch pad memory etc..
It can be seen from the above description that the embodiment of the present invention realizes following technique effect:By the phase that will log well, Seismic facies, seismic waveform inversion three are organically combined, and distribution of the distributary channel on seismic profile has been recognized accurately, into One step, position and the distribution characteristics of distributary channel are had identified using micro- ancient landform analytical technology, it is existing so as to solve The distribution characteristics for identifying distributary channel in shallow-circulation flow can not be simple and efficient, and the technology of the difficult prediction of reservoir religion is caused to be asked Topic, has reached the technique effect for being simple and efficient predicting reservoir.
Obviously, those skilled in the art should be understood that each module of the above-mentioned embodiment of the present invention or each step can be with It is realized with general computing device, they can concentrate on single computing device or be distributed in multiple computing devices On the network formed, optionally, they can be realized with the program code that computing device can perform, it is thus possible to by it Store and performed in the storage device by computing device, and in some cases, can be held with the order being different from herein They are either fabricated to each integrated circuit modules or will be multiple in them by the shown or described step of row respectively Module or step are fabricated to single integrated circuit module to realize.In this way, the embodiment of the present invention be not restricted to it is any specific hard Part and software combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to limit the invention, for the skill of this field For art personnel, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles of the invention, made Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.

Claims (10)

1. a kind of distributary channel method for determining position, which is characterized in that including:
Obtain the drilling well log of purpose series of strata;
According to the drilling well log, the well logging phase character of distributary channel is determined;
The position of distributary channel is calibrated on seismic profile according to the well logging phase character, the seismic facies for obtaining distributary channel is special Sign;
Seismic waveform inversion is carried out based on the log, obtains the position of distributary channel in purpose series of strata;
Tectonic interpretation is carried out to the purpose series of strata according to the characteristics of seismic, to recover to obtain micro- ancient landform;
The position of the distributary channel obtained based on the characteristics of seismic, inverting and micro- ancient landform, obtain the target zone The spread of distributary channel in system.
2. the method as described in claim 1, which is characterized in that according to the well logging phase character in the calibration point on seismic profile The position of fluid stream channel obtains the characteristics of seismic of distributary channel, including:
According to the well logging phase character of distributary channel, the well logging phase character type of distributary channel is determined;
According to the well logging phase character, distributary channel is demarcated onto corresponding seismic profile, to identify sympathetic earthquakes section The characteristics of seismic of upper distributary channel.
3. the method as described in claim 1, which is characterized in that after recovering to obtain micro- ancient landform, further include:
Ask for the compacting factor of every mouthful of well in micro- ancient landform;
Planar interpolation is carried out to the compacting factor of every mouthful of well, obtains compacting factor plan view;
Micro- ancient landform with the compacting factor plan view is multiplied, obtains micro- ancient landform after compaction correction.
4. method as claimed in claim 3, which is characterized in that ask for the compacting factor of every mouthful of well in micro- ancient landform, wrap It includes:
By porosity depth compaction model, ask for being compacted the compacting factor of every mouthful of well purpose series of strata:
φ (z)=φ0e-cz
Wherein, φ0Represent initial porosity, c represents compacting factor, and z represents depth.
5. the method as described in claim 1, which is characterized in that seismic waveform inversion is carried out based on the log, is obtained The position of distributary channel in purpose series of strata, including:
Load the seismic data and well-log information of the purpose series of strata;
The seismic data and the well-log information are demarcated;
The related well sample sheet of selection, sample for reference spatial distribution distance and curve distribution feature establish initial model;
Selected in the preset range of well point with the most similar seismic waveform of indicatrix, using global optimization approach, will select The seismic waveform taken is generalized to no wellblock, the company's of obtaining well seismic profile;
By the even well seismic profile compared with characterizing the log value of distributary channel and carrying out repeatedly, until obtaining meeting default It is required that geological model;
The geological model of foundation is applied to entire seismic data cube, to realize seismic waveform inversion, obtains purpose series of strata The position of middle distributary channel.
6. a kind of determining device of distributary channel position, which is characterized in that including:
Acquisition module, for obtaining the drilling well log of purpose series of strata;
First determining module, for according to the drilling well log, determining the well logging phase character of distributary channel;
Demarcating module for calibrating the position of distributary channel on seismic profile according to the well logging phase character, is shunted The characteristics of seismic in river;
Inverting module carries out seismic waveform inversion for being based on the log, obtains the position of distributary channel in purpose series of strata It puts;
Recovery module, for carrying out tectonic interpretation to the purpose series of strata according to the characteristics of seismic, to recover to obtain Micro- ancient landform;
Second determining module, for the position of distributary channel obtained based on the characteristics of seismic, inverting and micro- ancient ground Looks obtain the spread of distributary channel in the purpose series of strata.
7. device as claimed in claim 6, which is characterized in that the demarcating module includes:
Determination unit for the well logging phase character according to distributary channel, determines the well logging phase character type of distributary channel;
Unit is demarcated, for according to the well logging phase character, distributary channel being demarcated onto corresponding seismic profile, to identify The characteristics of seismic of distributary channel on sympathetic earthquakes section.
8. device as claimed in claim 6, which is characterized in that further include:
Module is asked for, for after recovering to obtain micro- ancient landform, asking for the compacting factor of every mouthful of well in micro- ancient landform;
Interpolating module for carrying out planar interpolation to the compacting factor of every mouthful of well, obtains compacting factor plan view;
Processing module for micro- ancient landform to be multiplied with the compacting factor plan view, obtains micro- Gu after compaction correction Landforms.
9. device as claimed in claim 8, which is characterized in that the module of asking for is specifically used for being compacted by porosity depth Model asks for being compacted the compacting factor of every mouthful of well purpose series of strata:
φ (z)=φ0e-cz
Wherein, φ0Represent initial porosity, c represents compacting factor, and z represents depth.
10. device as claimed in claim 6, which is characterized in that the inverting module includes:
Loading unit, for loading the seismic data and well-log information of the purpose series of strata;
Unit is demarcated, for being demarcated to the seismic data and the well-log information;
Unit is established, for selecting related well sample sheet, sample for reference spatial distribution distance and curve distribution feature establish introductory die Type;
Promote unit, for selected in the preset range of well point with the most similar seismic waveform of indicatrix, using the overall situation The seismic waveform of selection is generalized to no wellblock, the company's of obtaining well seismic profile by optimization algorithm;
Comparing unit, compared with the even well seismic profile is carried out repeatedly with the log value of characterization distributary channel, directly To obtaining the geological model for meeting preset requirement;
Applying unit, the geological model for that will establish is applied to entire seismic data cube, to realize seismic waveform inversion, Obtain the position of distributary channel in purpose series of strata.
CN201711145503.8A 2017-11-17 2017-11-17 Distributary channel method for determining position and device Active CN108107469B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711145503.8A CN108107469B (en) 2017-11-17 2017-11-17 Distributary channel method for determining position and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711145503.8A CN108107469B (en) 2017-11-17 2017-11-17 Distributary channel method for determining position and device

Publications (2)

Publication Number Publication Date
CN108107469A true CN108107469A (en) 2018-06-01
CN108107469B CN108107469B (en) 2019-07-09

Family

ID=62207620

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711145503.8A Active CN108107469B (en) 2017-11-17 2017-11-17 Distributary channel method for determining position and device

Country Status (1)

Country Link
CN (1) CN108107469B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111624651A (en) * 2019-02-28 2020-09-04 中国石油天然气股份有限公司 Reservoir prediction method and device based on ancient landform constraint
CN111766630A (en) * 2020-07-02 2020-10-13 中国地质大学(北京) Restoration method for ancient landform of basin
CN111983671A (en) * 2019-05-23 2020-11-24 中国石油天然气股份有限公司 Shallow water lake basin reservoir prediction method and device based on micro-ancient landform restoration
CN111997604A (en) * 2020-10-16 2020-11-27 陕西延长石油(集团)有限责任公司 Method for identifying single-stage river channel boundary
CN113703044A (en) * 2020-05-20 2021-11-26 中国石油化工股份有限公司 Correction method and device for width of ancient river channel, electronic equipment and storage medium
CN113721298A (en) * 2021-05-25 2021-11-30 中国石油化工股份有限公司 Reservoir prediction description method for sediment end underwater diversion river channel
US11542816B2 (en) 2021-01-28 2023-01-03 Saudi Arabian Oil Company System and method for hydrocarbon exploration based on imaging tunnel valleys

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101408624A (en) * 2007-10-08 2009-04-15 陶庆学 Forecasting and evaluating technologies of three-dimensional earthquake optimum time window river course sand body storage layer
CN102650702A (en) * 2012-05-03 2012-08-29 中国石油天然气股份有限公司 Seismic waveform analysis and reservoir prediction method and device
CN106772680A (en) * 2016-12-30 2017-05-31 中国石油天然气股份有限公司 Characterized based on two-dimentional poststack seismic data micro crack and determine method and apparatus

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101408624A (en) * 2007-10-08 2009-04-15 陶庆学 Forecasting and evaluating technologies of three-dimensional earthquake optimum time window river course sand body storage layer
CN102650702A (en) * 2012-05-03 2012-08-29 中国石油天然气股份有限公司 Seismic waveform analysis and reservoir prediction method and device
CN106772680A (en) * 2016-12-30 2017-05-31 中国石油天然气股份有限公司 Characterized based on two-dimentional poststack seismic data micro crack and determine method and apparatus

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
HUANG JUNPING ET AL.: "The origin and accumulation model of crude oils from oil reservoirs Chang 9 and Chang 10 in the Yanchang Formation of the Ordos Basin", 《CHIN.J.GEOCHEM》 *
HUANG JUNPING ET AL.: "The Palaeoclimate and Microscopic Reservoir Characteristics of Shanxi Formation Shales from Lower Permian in Qinshui Basin", 《ACTA GEOLOGICA SINICA》 *
吕乐: "苏里格气田盒8段河道砂体储层地震预测分析", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
姜正龙等: "古地貌恢复方法及应用—以济阳坳陷桩西地区沙二段为例", 《现代地质》 *
张荣军等: "《低渗透油藏有利区预测及工程设计——以西峰油田延安组油藏为例》", 31 August 2013, 北京:石油工业出版社 *
杨占龙等: ""三相"联合解释技术在岩性油气藏勘探中的应用—以吐哈盆地胜北地区为例", 《天然气地球科学》 *
许杰等: "鄂尔多斯盆地大牛地气田盒3段辫状河道沉积相预测", 《岩性油气藏》 *
邹新宁等: "鄂尔多斯盆地苏里格气田河道砂体识别研究", 《中国优秀博硕士学位论文全文数据库(博士) 基础科学辑》 *
郭川等: "塔河油田12区奥陶系油藏东部区域岩溶古河道识别及其意义", 《石油地质与工程》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111624651A (en) * 2019-02-28 2020-09-04 中国石油天然气股份有限公司 Reservoir prediction method and device based on ancient landform constraint
CN111624651B (en) * 2019-02-28 2023-08-22 中国石油天然气股份有限公司 Reservoir prediction method and device based on paleo-topography constraint
CN111983671A (en) * 2019-05-23 2020-11-24 中国石油天然气股份有限公司 Shallow water lake basin reservoir prediction method and device based on micro-ancient landform restoration
CN113703044A (en) * 2020-05-20 2021-11-26 中国石油化工股份有限公司 Correction method and device for width of ancient river channel, electronic equipment and storage medium
CN113703044B (en) * 2020-05-20 2024-04-09 中国石油化工股份有限公司 Correction method and device for ancient river channel width, electronic equipment and storage medium
CN111766630A (en) * 2020-07-02 2020-10-13 中国地质大学(北京) Restoration method for ancient landform of basin
CN111766630B (en) * 2020-07-02 2022-04-22 中国地质大学(北京) Restoration method for ancient landform of basin
CN111997604A (en) * 2020-10-16 2020-11-27 陕西延长石油(集团)有限责任公司 Method for identifying single-stage river channel boundary
CN111997604B (en) * 2020-10-16 2023-12-08 陕西延长石油(集团)有限责任公司 Single-stage river boundary identification method
US11542816B2 (en) 2021-01-28 2023-01-03 Saudi Arabian Oil Company System and method for hydrocarbon exploration based on imaging tunnel valleys
CN113721298A (en) * 2021-05-25 2021-11-30 中国石油化工股份有限公司 Reservoir prediction description method for sediment end underwater diversion river channel
CN113721298B (en) * 2021-05-25 2023-11-03 中国石油化工股份有限公司 Reservoir prediction description method for underwater diversion river channel at deposition end

Also Published As

Publication number Publication date
CN108107469B (en) 2019-07-09

Similar Documents

Publication Publication Date Title
CN108107469B (en) Distributary channel method for determining position and device
Grammer et al. Integration of outcrop and modern analogs in reservoir modeling: Overview with examples from the Bahamas
Regone et al. Geologic model building in SEAM Phase II—Land seismic challenges
CN107817535A (en) The determination method and apparatus of short lap
US20110231164A1 (en) Generating facies probablity cubes
CN108802812A (en) A kind of formation lithology inversion method of well shake fusion
CA2674089A1 (en) Method for generating reservoir models utilizing synthetic stratigraphic columns
CN106405637B (en) A kind of method and system for the prediction of thin sand-mud interbed sand thickness
CN105116449B (en) A kind of recognition methods of weak reflection reservoir
Hart Whither seismic stratigraphy?
Grijalba-Cuenca et al. Geostatistical inversion of 3D seismic data to extrapolate wireline petrophysical variables laterally away from the well
Mayall et al. Seismic interpretation workflows for deep-water systems: A practical guide for the subsurface
Embry et al. Seismic facies, stratigraphy, geomorphology, and seismic modelling of a Lower Cretaceous carbonate platform
Nwaezeapu et al. Sequence stratigraphic approach to hydrocarbon exploration: A case study of Chiadu field at eastern onshore Niger Delta basin, Nigeria
CN105093308B (en) A kind of well trajectory design method and system
Sverdrup et al. Sealing properties of faults and their influence on water-alternating-gas injection efficiency in the Snorre field, northern North Sea
Saeid et al. A depositional model for the Carbonera Formation, Llanos Foothills, Colombia, from workflow of a sequence stratigraphic framework and interpretation from well-log stacking patterns, well cuttings, and three-dimensional seismic spectral decomposition
Pyrcz et al. Representative input parameters for geostatistical simulation
Abbaszadeh et al. Development of an integrated reservoir model for a naturally fractured volcanic reservoir in China
Dim et al. Hydrocarbon leads and prospects opportunities across a cluster of fields in parts of Onshore Niger Delta Basin
Méndez et al. Predicting and 3D modeling of karst zones using seismic facies analysis in Ordovician carbonates of the Tahe oilfield, China
Ma et al. Coupling spatial and frequency uncertainty analyses in reservoir modeling: example of Judy Creek Reef complex in Swan Hills, Alberta, Canada
Akingbade et al. An integrated approach to identifying and reducing reservoir uncertainties in a turbidite deep water environment
Pennington et al. Calibration of seismic attributes for reservoir characterization
Latimer et al. Integrated seismic reservoir characterization and modeling: A Gulf of Mexico 3D case history

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant