CN104808259A - Detection method and device of gypsum rock distribution - Google Patents
Detection method and device of gypsum rock distribution Download PDFInfo
- Publication number
- CN104808259A CN104808259A CN201510151976.3A CN201510151976A CN104808259A CN 104808259 A CN104808259 A CN 104808259A CN 201510151976 A CN201510151976 A CN 201510151976A CN 104808259 A CN104808259 A CN 104808259A
- Authority
- CN
- China
- Prior art keywords
- data
- cream rock
- distribution
- point
- rock
- 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.)
- Pending
Links
Abstract
The invention provides a detection method and device of gypsum rock distribution and relates to the technical field of geological prospecting. The method includes that equal interval sampling is performed on a well logging density curve, a natural gamma curve and an interval transit time curve of each well in a preset area so that a plurality of sampling points are obtained; the sampling points of each well are joined to generate a plurality of cross plots under a orthogonal coordinate system; data points for representing gypsum rock features are determined according to density values and natural gamma values of all points in the cross plots; the data points for representing the gypsum rock features are projected on an area in a depth range that each data point corresponds to, and gypsum rock longitudinal distribution feature data of each well are generated; seismic amplitude attribute data of the preset area are obtained, and gypsum rock horizontal distribution feature data are generated according to the seismic amplitude attribute data; gypsum rock space distribution feature data are generated according to the gypsum rock longitudinal distribution feature data and the gypsum rock horizontal distribution feature data of each well. By means of the detection method and device of gypsum rock distribution, space distribution of gypsum rocks can be accurately detected.
Description
Technical field
The present invention relates to technical field of geological exploration, particularly relate to detection method and the device of the distribution of a kind of cream rock.
Background technology
Current, numerous with the large oil and gas fields that cream rock is cap rock, such as all there is this kind of oil gas field the area such as North America, the Middle East, Siberia and Chinese Tarim Basin, the Sichuan Basin etc.
Cream rock serves the effect of capping to oil gas, simultaneously or the mark of the Favorable Areas of source rock development.Because cream rock is formed and the environment of deposition often salinity is very high, the sulfate bacterium of anaerobism not easily survives, and forms weak oxide-reducing environment, and this environment is beneficial to organic preservation and to the hydrocarbon conversion.In addition, cream rock or the favourable development area of High-quality Reservoir.Because cream rock stratum thermal conductivity is high, kapillary diverts oneself from loneliness or boredom pressure greatly, effectively can slows down the compaction on stratum, impel the factor of porosity that stratum, its underpart keeps higher.Along with buried depth of strata increases, gypsum can be sloughed bulk crystallization water and change to muriacite, and these water of crystallization are rich in organic acid, and its dissolution enhances fluid-rock reaction, corrosion mineral and form secondary pores, effectively can improve reservoir properties.In addition, cream rock also has very strong plasticity.Simultaneously due to the existence of cream rock, also can cause in basin and complicated thermochemical sulfate reduction (being called for short TSR) occurs between hydrocarbon-water-rock three, can hydrogen sulfide gas be produced, have a series of impact to oil-gas exploration, Hydrocarbon Formation Reservoirs etc.
Visible, if the space distribution of cream rock accurately can be detected, will to oil-gas mining advantageously.
Summary of the invention
The detection method that embodiments of the invention provide a kind of cream rock to distribute and device, to realize the space distribution accurately detecting cream rock.
For achieving the above object, the present invention adopts following technical scheme:
A detection method for cream rock distribution, comprising:
The well logging densimetric curve of well each in one predeterminable area, gamma ray curve and interval transit time curve are equidistantly sampled, acquires multiple sampled point;
The sampled point of each well is crossed, the multiple figure that cross under generating orthogonal coordinate system; The horizontal ordinate of the described figure of crossing is interval transit time value, ordinate is density value, the parameter that crosses is natural gamma value;
According to density value and the natural gamma value of each point in the described figure that crosses, determine the data point representing cream rock feature;
By described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well;
Obtain the seismic amplitude attribute data of described predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data;
According to cream rock Vertical Distribution Characteristics data and the described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data of described each well.
Concrete, the density value of each point and natural gamma value in the figure that crosses described in described basis, determine the data point representing cream rock feature, comprising:
Obtain the density threshold and natural gamma threshold value that pre-set;
Choose density value from each point the described figure of crossing and be greater than described density threshold and natural gamma value is less than the number of targets strong point of described natural gamma threshold value, using described number of targets strong point as being used for the data point representing cream rock feature.
In addition, by described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, after generating the cream rock Vertical Distribution Characteristics data of each well, comprising:
Obtain the well lithological column data of actual measurement, described histogram data and described cream rock Vertical Distribution Characteristics data are compared, determines the coincidence rate of described cream rock Vertical Distribution Characteristics data and described histogram data.
Concrete, the seismic amplitude attribute data of the described predeterminable area of described acquisition, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data, comprising:
Pass through formula:
y＝5.022742+0.000374x
Described seismic amplitude attribute data is converted to cream rock planar characteristics of distribution data; Wherein x is the amplitude in seismic amplitude attribute, and y is cream rock thickness.
Described seismic amplitude attribute data is converted to cream rock planar characteristics of distribution data.
Concrete, described seismic amplitude attribute data comprises average absolute value amplitude, average trough amplitude, maximum value amplitude, accumulated absolute values amplitude, RMS amplitude, amplitude variations data.
A pick-up unit for cream rock distribution, comprising:
Sampling unit, equidistantly samples for the well logging densimetric curve to well each in a predeterminable area, gamma ray curve and interval transit time curve, acquires multiple sampled point;
Cross figure generation unit, for crossing to the sampled point of each well, and the multiple figure that cross under generating orthogonal coordinate system; The horizontal ordinate of the described figure of crossing is interval transit time value, ordinate is density value, the parameter that crosses is natural gamma value;
Data point determining unit, for density value and the natural gamma value of each point in the figure that crosses described in basis, determines the data point representing cream rock feature;
Vertical Distribution Characteristics data generating unit, for by described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well;
Planar characteristics of distribution data generating unit, for obtaining the seismic amplitude attribute data of described predeterminable area, and generates cream rock planar characteristics of distribution data according to described seismic amplitude attribute data;
Spatial distribution characteristic data generating unit, for according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.
Concrete, this data point determining unit, comprising:
Threshold value acquisition module, for obtaining the density threshold and natural gamma threshold value that pre-set;
Data point determination module, is greater than described density threshold and natural gamma value is less than the number of targets strong point of described natural gamma threshold value, using described number of targets strong point as being used for the data point representing cream rock feature for choosing density value from each point in the described figure of crossing.
In addition, the pick-up unit of this cream rock distribution, also comprises:
Coincidence rate authentication unit, for obtaining the well lithological column data of actual measurement, comparing described histogram data and described cream rock Vertical Distribution Characteristics data, determining the coincidence rate of described cream rock Vertical Distribution Characteristics data and described histogram data.
In addition, described planar characteristics of distribution data generating unit, specifically for:
Pass through formula:
y＝5.022742+0.000374x
Described seismic amplitude attribute data is converted to cream rock planar characteristics of distribution data; Wherein x is the amplitude in seismic amplitude attribute, and y is cream rock thickness.
The detection method of cream rock that the embodiment of the present invention provides distribution and device, equidistantly sample to the well logging densimetric curve of well each in a predeterminable area, gamma ray curve and interval transit time curve, acquire multiple sampled point; And the sampled point of each well is crossed, the multiple figure that cross under generating orthogonal coordinate system; Thus according to the density value of each point in the described figure that crosses and natural gamma value, determine the data point representing cream rock feature; By described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well; And obtain the seismic amplitude attribute data of described predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data; Thus according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.The cream rock spatial distribution characteristic data generated can represent the space distribution of cream rock, thus the present invention can accurately provide the space distribution of cream rock for oil-gas mining.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
The process flow diagram one of the detection method of a kind of cream rock distribution that Fig. 1 provides for the embodiment of the present invention;
The flowchart 2 of the detection method of a kind of cream rock distribution that Fig. 2 provides for the embodiment of the present invention;
Fig. 3 is that log response in the embodiment of the present invention crosses schematic diagram;
Fig. 4 is that the gypsum-salt rock in the embodiment of the present invention grows longitudinal profile feature schematic diagram;
Fig. 5 is the gypsum-salt rock flat distribution map of the different attribute in the embodiment of the present invention;
Fig. 6 is the flat distribution map of the different earthquake attribute in the embodiment of the present invention;
Fig. 7 is the imperial hillock district gypsum-salt rock thickness prediction figure in the embodiment of the present invention;
The structural representation one of the pick-up unit of a kind of cream rock distribution that Fig. 8 provides for the embodiment of the present invention;
The structural representation two of the pick-up unit of a kind of cream rock distribution that Fig. 9 provides for the embodiment of the present invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
The invention provides the detection method of a kind of cream rock distribution, as shown in Figure 1, comprising:
Step 101, well logging densimetric curve, gamma ray curve and interval transit time curve to well each in a predeterminable area are equidistantly sampled, and acquire multiple sampled point.
Step 102, the sampled point of each well to be crossed, the multiple figure that cross under generating orthogonal coordinate system.
Wherein, the horizontal ordinate of figure of crossing described in is interval transit time value, ordinate is density value, the parameter that crosses is natural gamma value.
Step 103, according to described in cross the density value of each point in figure and natural gamma value, determine the data point representing cream rock feature.
Step 104, by described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well.
Step 105, obtain the seismic amplitude attribute data of described predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data.
Step 106, according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.
The detection method of cream rock that the embodiment of the present invention provides distribution, equidistantly samples to the well logging densimetric curve of well each in a predeterminable area, gamma ray curve and interval transit time curve, acquires multiple sampled point; And the sampled point of each well is crossed, the multiple figure that cross under generating orthogonal coordinate system; Thus according to the density value of each point in the described figure that crosses and natural gamma value, determine the data point representing cream rock feature; By described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well; And obtain the seismic amplitude attribute data of described predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data; Thus according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.The cream rock spatial distribution characteristic data generated can represent the space distribution of cream rock, thus the present invention can accurately provide the space distribution of cream rock for oil-gas mining.
In order to make those skilled in the art better understand the present invention, enumerate a more specifically embodiment below, as shown in Figure 2, the detection method of a kind of cream rock distribution that the embodiment of the present invention provides, comprising:
Step 201, well logging densimetric curve, gamma ray curve and interval transit time curve to well each in a predeterminable area are equidistantly sampled, and acquire multiple sampled point.
Step 202, the sampled point of each well to be crossed, the multiple figure that cross under generating orthogonal coordinate system.
Wherein, the horizontal ordinate of figure of crossing described in is interval transit time value, ordinate is density value, the parameter that crosses is natural gamma value.
The density threshold that step 203, acquisition pre-set and natural gamma threshold value.
Step 204, to choose density value from each point the described figure of crossing and be greater than described density threshold and natural gamma value is less than the number of targets strong point of described natural gamma threshold value, using described number of targets strong point as being used for the data point representing cream rock feature.
Step 205, by described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well.
The well lithological column data of step 206, acquisition actual measurement, compare described histogram data and described cream rock Vertical Distribution Characteristics data, determine the coincidence rate of described cream rock Vertical Distribution Characteristics data and described histogram data.
By determining that the coincidence rate of described cream rock Vertical Distribution Characteristics data and described histogram data can know that whether cream rock Vertical Distribution Characteristics data are accurate.
Step 207, obtain the seismic amplitude attribute data of described predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data.
Concrete, this step 207 can be pass through formula: y=5.022742+0.000374x, described seismic amplitude attribute data is converted to cream rock planar characteristics of distribution data; Wherein x is the amplitude in seismic amplitude attribute, and y is cream rock thickness.
In one embodiment, seismic amplitude attribute data comprises average absolute value amplitude, average trough amplitude, maximum value amplitude, accumulated absolute values amplitude, RMS amplitude, amplitude variations data.
Step 208, according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.
The detection method of cream rock that the embodiment of the present invention provides distribution, equidistantly samples to the well logging densimetric curve of well each in a predeterminable area, gamma ray curve and interval transit time curve, acquires multiple sampled point; And the sampled point of each well is crossed, the multiple figure that cross under generating orthogonal coordinate system; Thus according to the density value of each point in the described figure that crosses and natural gamma value, determine the data point representing cream rock feature; By described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well; And obtain the seismic amplitude attribute data of described predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data; Thus according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.The cream rock spatial distribution characteristic data generated can represent the space distribution of cream rock, thus the present invention can accurately provide the space distribution of cream rock for oil-gas mining.
Corresponding to above-mentioned embodiment of the method, enumerate the instantiation of the detection method of an above-mentioned cream rock distribution below:
First according to the well logging of 3 sections, COMPONENT IN SOUTH SICHUAN BASIN land mine densimetric curve, natural gamma curve and interval transit time curve, different well is equidistantly sampled, the sampled point of every mouthful of well is crossed, set up the multiple figure that cross under orthogonal coordinate system, as shown in Figure 3, cream rock is distinguished from carbonatite (limestone and pierite).Concrete, as shown in Figure 3 cross in figure, horizontal ordinate is interval transit time value, ordinate is density value, GR value is as the third parameter that crosses, characterize different natural gamma values by gray scale, gray scale high (inclined black) represent that GR value is lower, gray scale is low (partially white) to represent GR value higher.Cross on figure as can be seen from 12 mouthfuls of wells, within the scope of study area, the natural gamma value of thunder three sections and density value are roughly in negative correlativing relation.With the high point of gray scale generally within the higher upper half of density value, with the low point of gray scale generally within the lower bottom half of density value, and point different situation of horizontal ordinate direction sampled point is poor, and the data point that gray scale is high and low is mixed in together.Therefore, gray scale high and be positioned at the usual corresponding cream rock of point of upper half or gypsum content higher, gray scale is low and be positioned at the point then usual corresponding non-cream rock development area of bottom half.
Afterwards, according to the density in every mouthful of well corresponding to cream rock and natural gamma Distribution value scope, natural gamma threshold value and a density threshold is pre-set.Choose density value from each point the described figure of crossing and be greater than described density threshold and natural gamma value is less than the number of targets strong point of described natural gamma threshold value, using described number of targets strong point as being used for the data point representing cream rock feature.By for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well, such as by target data spot projection in the corresponding degree of depth, form a curve prediction lithological column (as shown in first row on the right side of scatter diagram in Fig. 3), this histogram is the individual well cream rock predicted by density and natural gamma and longitudinally grows histogram.In addition, as shown in Figure 4, this figure is that gypsum-salt rock grows longitudinal profile feature schematic diagram.
Afterwards, natural gamma and density value is utilized to predict the thunder three sections of cream rock genesis analysis histograms drawn, contrast with actual well logging lithological column, in selected 12 mouthfuls of wells, the cream rock utilizing logging trace to obtain predicts the outcome more satisfactory with actual lithology corresponding relation, prediction coincidence rate more than 90%.If corresponding relation is not good, may be then because lithology is (as near Lg18 well 4248m, near Lg6 well 3030m) that the cream matter pierite containing muriacite mineral causes, have plenty of (as near Lg12 well 4170m, near Lg20 well 4320m) because cream matter limestone causes, what have is then because logging trace interrupts caused (as dark in Lg173 well 3338m).
Afterwards, obtain the seismic amplitude attribute data of predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data.Specific implementation can be such as utilize imperial hillock district three dimensional seismic data, explain Ding Di circle of thunder three sections, then the seismic properties (Fig. 5) that thunder three sections amounts to 9 types is extracted, comprising 6 kinds of amplitude generic attributes, 2 kinds of energy generic attributes and a kind of reflection strength attribute.Amplitude generic attribute comprises Average AbsoluteAmplitude (average absolute value amplitude), Average Trough Amplitude (average trough amplitude), MaximumAbsolute Amplitude (maximum value amplitude), Total Absolute Amplitude (accumulated absolute values amplitude), Root-Mean-Square (RMS amplitude) and Variance in Amplitude (amplitude variations).Energy generic attribute comprises Average Energy (average energy) and Total Energy (gross energy).Reflection strength attribute is AverageReflection Strength (average reflection intensity).
In order to strengthen the corresponding relation of amplitude attribute and cream rock further, improve precision of prediction, select four kinds of amplitude association attributeses to carry out superposition, the comprehensive earthquake amplitude attribute (on the left of Fig. 6) obtained discloses cream rock thickness distribution trend (on the right side of Fig. 6) in planar characteristics of distribution with drilling well comparatively identical.
Although the response ratio of seismic amplitude attribute to cream rock is more responsive, precision of prediction is also higher, still only relies on earthquake prediction cream rock plane distribution can not obtain the thickness spread situation of cream rock.Therefore, utilize the linear fit relational expression (1) that related coefficient is the highest in the present embodiment, earthquake synthesized attribute figure is converted to the cream rock thickness distribution prognostic chart having more geological Significance, as shown in Figure 7.
y＝5.022742+0.000374x (1)
X-seismic properties (amplitude); Y-cream rock thickness (m).
As shown in Figure 7, in plane, cream rock variation in thickness is remarkable, and the thickest in hundred meters, the thinnest only has tens meters.Plane is seen, main exist three cream rock developmental centers (A, B, C), and area is the widest and that thickness is maximum one (A) is positioned at the imperial hilllock 37-dragon wellblock, hilllock 106 of southwest corner, work area, and in triangular shape, its distribution area is about 530km
2, cream rock thickness exceedes 100m more.Another cream rock developmental center (B) is positioned near the wellblock, imperial hilllock 10 of Middle-north Area, work area, and in strip, area is about 360km
2, thickness is between 90-100m.Last cream rock developmental center (C) is positioned at the southeast, work area, and its distributed areas comparatively disperse, and thickness is thin compared with the first two center.
Other position cream rocks in study area sporadicly distribute, thickness 50-80m, and cream rock is in " concentrating type " large area spread.The weak area that cream rock is grown is ribbon, and along northwest (NW)-east southeast trend, area is about 800km
2, the cream rock thickness in this band is less than 50m substantially.Only near the well of imperial hilllock 172,161,21 and 11, cream rock is grown among a small circle, and thickness maximum is close to hundred meters.
Even dragon hillock district thunder three sections does not find obvious Salt deposition in the area that cream rock thickness is maximum yet, only extraneous in work area, on the south imperial Gang37Jing and imperial hilllock 160 well about 16km imperial hilllock 20 well in grow the thick rock salt of nearly 30m.This illustrates thunder three sections of deposition period, sedimentary environment in study area belongs to the cream rock lake deposit environment of half limitation, deposited gypsolyte and pierite, the imperial Gang20Jing being in study area southwestward then belongs to the cream lopolith sedimentary environment of subtidal zone, mainly with carbonatite, cream rock and Salt deposition for feature.
Corresponding to the embodiment of the method for Fig. 1 and Fig. 2, the pick-up unit of a kind of cream rock distribution provided by the invention, as shown in Figure 8, comprising:
Sampling unit 31, equidistantly can sample to the well logging densimetric curve of well each in a predeterminable area, gamma ray curve and interval transit time curve, acquire multiple sampled point.
Cross figure generation unit 32, can cross to the sampled point of each well, the multiple figure that cross under generating orthogonal coordinate system; The horizontal ordinate of the described figure of crossing is interval transit time value, ordinate is density value, the parameter that crosses is natural gamma value.
Data point determining unit 33, according to the density value of each point in the described figure that crosses and natural gamma value, can determine the data point representing cream rock feature.
Vertical Distribution Characteristics data generating unit 34, by described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, can generate the cream rock Vertical Distribution Characteristics data of each well.
Planar characteristics of distribution data generating unit 35, can obtain the seismic amplitude attribute data of described predeterminable area, and generates cream rock planar characteristics of distribution data according to described seismic amplitude attribute data.
Spatial distribution characteristic data generating unit 36, can according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.
Concrete, as shown in Figure 9, this data point determining unit 33, can comprise:
Threshold value acquisition module 331, can obtain the density threshold and natural gamma threshold value that pre-set.
Data point determination module 332, can choose density value and be greater than described density threshold and natural gamma value is less than the number of targets strong point of described natural gamma threshold value, using described number of targets strong point as being used for the data point representing cream rock feature from each point the described figure of crossing.
In addition, as shown in Figure 9, the pick-up unit of this cream rock distribution, also comprises:
Coincidence rate authentication unit 37, can obtain the well lithological column data of actual measurement, described histogram data and described cream rock Vertical Distribution Characteristics data is compared, determines the coincidence rate of described cream rock Vertical Distribution Characteristics data and described histogram data.
In addition, described planar characteristics of distribution data generating unit 35, specifically can by formula: described seismic amplitude attribute data is converted to cream rock planar characteristics of distribution data by y=5.022742+0.000374x; Wherein x is the amplitude in seismic amplitude attribute, and y is cream rock thickness.
What deserves to be explained is, the specific implementation of the pick-up unit of the cream rock distribution that the embodiment of the present invention provides, see above-mentioned embodiment of the method, repeats no more herein.
The pick-up unit of cream rock that the embodiment of the present invention provides distribution, equidistantly samples to the well logging densimetric curve of well each in a predeterminable area, gamma ray curve and interval transit time curve, acquires multiple sampled point; And the sampled point of each well is crossed, the multiple figure that cross under generating orthogonal coordinate system; Thus according to the density value of each point in the described figure that crosses and natural gamma value, determine the data point representing cream rock feature; By described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well; And obtain the seismic amplitude attribute data of described predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data; Thus according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.The cream rock spatial distribution characteristic data generated can represent the space distribution of cream rock, thus the present invention can accurately provide the space distribution of cream rock for oil-gas mining.
Apply specific embodiment in the present invention to set forth principle of the present invention and embodiment, the explanation of above embodiment just understands method of the present invention and core concept thereof for helping; Meanwhile, for one of ordinary skill in the art, according to thought of the present invention, all will change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention.
Claims (9)
1. a detection method for cream rock distribution, is characterized in that, comprising:
The well logging densimetric curve of well each in one predeterminable area, gamma ray curve and interval transit time curve are equidistantly sampled, acquires multiple sampled point;
The sampled point of each well is crossed, the multiple figure that cross under generating orthogonal coordinate system; The horizontal ordinate of the described figure of crossing is interval transit time value, ordinate is density value, the parameter that crosses is natural gamma value;
According to density value and the natural gamma value of each point in the described figure that crosses, determine the data point representing cream rock feature;
By described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well;
Obtain the seismic amplitude attribute data of described predeterminable area, and generate cream rock planar characteristics of distribution data according to described seismic amplitude attribute data;
According to cream rock Vertical Distribution Characteristics data and the described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data of described each well.
2. the detection method of cream rock according to claim 1 distribution, is characterized in that, the density value of each point and natural gamma value in the figure that crosses described in described basis, determine the data point representing cream rock feature, comprising:
Obtain the density threshold and natural gamma threshold value that pre-set;
Choose density value from each point the described figure of crossing and be greater than described density threshold and natural gamma value is less than the number of targets strong point of described natural gamma threshold value, using described number of targets strong point as being used for the data point representing cream rock feature.
3. the detection method of cream rock according to claim 2 distribution, is characterized in that, by described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, after generating the cream rock Vertical Distribution Characteristics data of each well, comprising:
Obtain the well lithological column data of actual measurement, described histogram data and described cream rock Vertical Distribution Characteristics data are compared, determines the coincidence rate of described cream rock Vertical Distribution Characteristics data and described histogram data.
4. the detection method of cream rock distribution according to claim 3, is characterized in that, the seismic amplitude attribute data of the described predeterminable area of described acquisition, and generates cream rock planar characteristics of distribution data according to described seismic amplitude attribute data, comprising:
Pass through formula:
y＝5.022742+0.000374x
Described seismic amplitude attribute data is converted to cream rock planar characteristics of distribution data; Wherein x is the amplitude in seismic amplitude attribute, and y is cream rock thickness.
5. the detection method of cream rock distribution according to claim 4, it is characterized in that, described seismic amplitude attribute data comprises average absolute value amplitude, average trough amplitude, maximum value amplitude, accumulated absolute values amplitude, RMS amplitude, amplitude variations data.
6. a pick-up unit for cream rock distribution, is characterized in that, comprising:
Sampling unit, equidistantly samples for the well logging densimetric curve to well each in a predeterminable area, gamma ray curve and interval transit time curve, acquires multiple sampled point;
Cross figure generation unit, for crossing to the sampled point of each well, and the multiple figure that cross under generating orthogonal coordinate system; The horizontal ordinate of the described figure of crossing is interval transit time value, ordinate is density value, the parameter that crosses is natural gamma value;
Data point determining unit, for density value and the natural gamma value of each point in the figure that crosses described in basis, determines the data point representing cream rock feature;
Vertical Distribution Characteristics data generating unit, for by described for representing that the data point of cream rock feature projects to the depth range corresponding to each data point, generate the cream rock Vertical Distribution Characteristics data of each well;
Planar characteristics of distribution data generating unit, for obtaining the seismic amplitude attribute data of described predeterminable area, and generates cream rock planar characteristics of distribution data according to described seismic amplitude attribute data;
Spatial distribution characteristic data generating unit, for according to the cream rock Vertical Distribution Characteristics data of described each well and described cream rock planar characteristics of distribution data genaration cream rock spatial distribution characteristic data.
7. the pick-up unit of cream rock distribution according to claim 6, it is characterized in that, data point determining unit, comprising:
Threshold value acquisition module, for obtaining the density threshold and natural gamma threshold value that pre-set;
Data point determination module, is greater than described density threshold and natural gamma value is less than the number of targets strong point of described natural gamma threshold value, using described number of targets strong point as being used for the data point representing cream rock feature for choosing density value from each point in the described figure of crossing.
8. the pick-up unit of cream rock distribution according to claim 7, is characterized in that, also comprise:
Coincidence rate authentication unit, for obtaining the well lithological column data of actual measurement, comparing described histogram data and described cream rock Vertical Distribution Characteristics data, determining the coincidence rate of described cream rock Vertical Distribution Characteristics data and described histogram data.
9. the pick-up unit of cream rock distribution according to claim 8, is characterized in that, described planar characteristics of distribution data generating unit, specifically for:
Pass through formula:
y＝5.022742+0.000374x
Described seismic amplitude attribute data is converted to cream rock planar characteristics of distribution data; Wherein x is the amplitude in seismic amplitude attribute, and y is cream rock thickness.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510151976.3A CN104808259A (en) | 2015-04-01 | 2015-04-01 | Detection method and device of gypsum rock distribution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510151976.3A CN104808259A (en) | 2015-04-01 | 2015-04-01 | Detection method and device of gypsum rock distribution |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104808259A true CN104808259A (en) | 2015-07-29 |
Family
ID=53693244
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510151976.3A Pending CN104808259A (en) | 2015-04-01 | 2015-04-01 | Detection method and device of gypsum rock distribution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104808259A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110083725A (en) * | 2019-04-15 | 2019-08-02 | 青岛海洋科学与技术国家实验室发展中心 | Sample element data processing method and system |
| CN111665574A (en) * | 2019-03-05 | 2020-09-15 | 中国石油化工集团公司 | Well logging interpretation method and system for rock thermal parameters |
| CN111983189A (en) * | 2020-08-18 | 2020-11-24 | 成都理工大学 | Method for quantitatively analyzing ultra-deep gypsum salt rock and dolomite coupling mechanism |
| CN112485843A (en) * | 2019-09-12 | 2021-03-12 | 中国石油天然气股份有限公司 | Prediction method, device and equipment for overflow layer in gypsum rock |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101315428A (en) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | High resolution prospecting method of large-ventage clastic rock deposition stratum based on horizontal wave velocity |
| CN102312669A (en) * | 2010-07-01 | 2012-01-11 | 中国石油集团长城钻探工程有限公司 | Method (Ip) for evaluating properties of reservoir based on gas survey logging data |
| CN102455437A (en) * | 2010-10-14 | 2012-05-16 | 中国石油天然气集团公司 | Method for determining reservoir layer and fluid body |
| CN102749651A (en) * | 2012-07-12 | 2012-10-24 | 恒泰艾普石油天然气技术服务股份有限公司 | Calibration method for fracture logging of reservoir |
| CN102759748A (en) * | 2012-07-04 | 2012-10-31 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic reservoir prediction method by cross analysis based on target analysis |
| CN102914793A (en) * | 2012-10-12 | 2013-02-06 | 中国石油天然气股份有限公司 | Nonequivalent modeling method and device for heterogeneous reservoir stratum |
| CN102967883A (en) * | 2012-11-20 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for predicting rock brittleness probability through pre-stack elasticity parameter inversion of shale gas |
-
2015
- 2015-04-01 CN CN201510151976.3A patent/CN104808259A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101315428A (en) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | High resolution prospecting method of large-ventage clastic rock deposition stratum based on horizontal wave velocity |
| CN102312669A (en) * | 2010-07-01 | 2012-01-11 | 中国石油集团长城钻探工程有限公司 | Method (Ip) for evaluating properties of reservoir based on gas survey logging data |
| CN102455437A (en) * | 2010-10-14 | 2012-05-16 | 中国石油天然气集团公司 | Method for determining reservoir layer and fluid body |
| CN102759748A (en) * | 2012-07-04 | 2012-10-31 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Seismic reservoir prediction method by cross analysis based on target analysis |
| CN102749651A (en) * | 2012-07-12 | 2012-10-24 | 恒泰艾普石油天然气技术服务股份有限公司 | Calibration method for fracture logging of reservoir |
| CN102914793A (en) * | 2012-10-12 | 2013-02-06 | 中国石油天然气股份有限公司 | Nonequivalent modeling method and device for heterogeneous reservoir stratum |
| CN102967883A (en) * | 2012-11-20 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for predicting rock brittleness probability through pre-stack elasticity parameter inversion of shale gas |
Non-Patent Citations (3)
| Title |
|---|
| 张水山等: "潜江凹陷蚌湖向斜周缘岩性油藏识别技术及效果", 《岩性油气藏》 * |
| 徐美娥等: "四川盆地东南部中、下寒武统膏岩盖层分布特征及封盖有效性", 《石油天然气地质》 * |
| 魏水建等: "川东北通南巴地区三叠系膏盐岩盖层预测", 《石油实验地质》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111665574A (en) * | 2019-03-05 | 2020-09-15 | 中国石油化工集团公司 | Well logging interpretation method and system for rock thermal parameters |
| CN111665574B (en) * | 2019-03-05 | 2023-04-07 | 中国石油化工集团公司 | Well logging interpretation method and system for rock thermal parameters |
| CN110083725A (en) * | 2019-04-15 | 2019-08-02 | 青岛海洋科学与技术国家实验室发展中心 | Sample element data processing method and system |
| CN110083725B (en) * | 2019-04-15 | 2023-04-07 | 青岛海洋科技中心 | Sample element data processing method and system |
| CN112485843A (en) * | 2019-09-12 | 2021-03-12 | 中国石油天然气股份有限公司 | Prediction method, device and equipment for overflow layer in gypsum rock |
| CN111983189A (en) * | 2020-08-18 | 2020-11-24 | 成都理工大学 | Method for quantitatively analyzing ultra-deep gypsum salt rock and dolomite coupling mechanism |
| CN111983189B (en) * | 2020-08-18 | 2021-06-08 | 成都理工大学 | Method for quantitatively analyzing ultra-deep gypsum salt rock and dolomite coupling mechanism |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106869790B (en) | A kind of quick fine geology guidance method of shale gas horizontal well | |
| CN104612675B (en) | A kind of carbonate formation Lithology while drilling method for quickly identifying | |
| CN103670383B (en) | A kind of method and apparatus identifying shale oil deposit Effective Reservoirs | |
| CN102052070B (en) | Complex lithology clastic rock stratum constituent logging processing method | |
| CN109388817A (en) | A kind of Reservoir Fracture three-dimensional modeling method | |
| CN104655821B (en) | A kind of carbonate reservoir origin cause of formation sentences knowledge method | |
| CN101403301B (en) | Method for acquiring container layer coarseness area distribution rule | |
| CN105986815A (en) | Method for recognizing geological sweet spots in shale stratum | |
| CN105590012A (en) | Estimation method for favorable sand of sandstone-type uranium deposit adapted to interlayer oxidation zone | |
| Yu et al. | Fault system impact on paleokarst distribution in the Ordovician Yingshan Formation in the central Tarim basin, northwest China | |
| CN105676294B (en) | A kind of carbonate rock fractured cave type Reservoir Body quantitative description | |
| CN104808259A (en) | Detection method and device of gypsum rock distribution | |
| CN105242307B (en) | Carbonate complex seismic reservoir porosity acquisition methods and device | |
| CN105240005A (en) | Low-hole extra-low-permeability unconventional reservoir dissolution hole identification method | |
| CN107515957A (en) | Mud shale SEQUENCE STRATIGRAPHIC DIVISION method | |
| CN105425315A (en) | Inversion method of microscale deposition stage palaeogeomorphology in weak construction deformation range | |
| CN104632201A (en) | Determination method for lower limiting value of organic carbon content in shale oil and gas 'dessert area' | |
| CN103291287A (en) | Effectiveness level measuring method for hole-type reservoirs | |
| Barnett et al. | Seismic characterization and origin of clinoforms in lacustrine depositional environments: a case study from the Cretaceous of the South Atlantic | |
| CN111967631A (en) | Method for predicting prospect area of early mineralization of sandstone-type uranium deposit exploration | |
| CN108804728A (en) | Horizontal well stratum reservoir grading analysis method and computer readable storage medium | |
| CN105425309A (en) | Method for determining tufaceous content of high-gamma sandstone | |
| CN106199754B (en) | Oil gas drilling target integrates optimizing evaluation method | |
| Ping et al. | Distribution pattern and origin of the columnar pull-down anomalies in Gaoshiti Block of Central Sichuan Basin, SW China | |
| CN107092995A (en) | Geothermal resources assessment method and apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150729 |