Background technology
Near surface structure investigation is an important step of seismic prospecting, all plays an important role in the earthquake stimulation receiving parameter design instructing seismic prospecting field to produce and indoor earthquake data static corrections processing.The data processing method of the near surface structure investigation of current employing is mainly through every square kilometre of laying 0.25 ~ 1 near surface structure points for investigation in enumeration district, gather near surface structure data, and adopt pure mathematics data interpolating method when processing, Similarity interpolation data are generated automatically to the region not gathering near surface structure data, finally produce and generate the near surface structure figure in whole enumeration district, wherein, described near surface structure data mainly comprise speed data and the thickness data of each layer of near surface, also can comprise other data simultaneously, the depth data of such as water table, Lithology data, topographic condition data etc.Said method changes little Continental Facies Stratigraphy district in the little north of earth's surface topography variation and rock property, and the precision of the near surface structure data calculated is relatively high.
But marine strata in south China has following features: the 1. old stratum of the main exposure in marine deposit stratum, mostly is carbonatite (that is, limestone) and sandstone formation crop out, near-surface velocity and variation in thickness large; 2. southern landform is mainly mountain region, and surface relief is large, and the structure change of near surface is fast.Therefore, the near surface structure data that the data processing method adopting the near surface structure of above-mentioned routine to investigate calculates truly can not reflect the structure situation of change of near surface, and error is comparatively large, thus can not the needs of seismic exploration completely contentedly.
Summary of the invention
The object of the present invention is to provide a kind of marine strata in south China near surface structure data processing method, comprising: near surface structure data are gathered to the near surface structure points for investigation in work area; The data interval of the near surface structure data corresponding at least one feature of the non-points for investigation in work area is found near surface structure data sample storehouse, and by the near surface structure data constraint of non-points for investigation in described data interval; Judge whether the variation range of described data interval is less than predetermined value, if the variation range of described data interval is not less than predetermined value, then increase the feature of non-points for investigation to retrain further the near surface structure data of non-points for investigation, until the data interval after constraint is less than predetermined value; And the data interval of the near surface structure data of utility efficiency point and the near surface structure data of non-points for investigation obtains the near surface structure data profile in work area.
Near surface structure data can comprise at least one in the speed data of each layer of near surface, thickness data and water table depth data.
The predetermined value of 100m/s can be set speed data, thickness data be arranged to the predetermined value of 1m, and the predetermined value of 0.5m is arranged to the water table degree of depth.
Near surface structure data sample storehouse can be obtained in advance by carrying out near surface structure investigation collection to various near surface structure, or using the near surface structure data of the near surface structure points for investigation in work area as near surface structure data sample storehouse.
The feature of non-points for investigation can comprise rock property, topographic condition and the water table degree of depth.
Topographic condition can comprise earth's surface elevation, the gradient and the distance from cliff.
Near surface structure data profile in work area can comprise at least one in the speed data distribution plan of each layer of near surface, thickness data distribution plan and water table depth data distribution plan.
Non-points for investigation can be distributed based on the situation of change of the needs of data processing or near surface structure.
Data interval can be the set of at least one value data.
By getting intermediate value to obtain value data to two endpoint values of the data interval after the constraint being less than predetermined value, and the near surface structure data of utility efficiency point and described value data obtain the near surface structure data profile in work area.
Part in ensuing description is set forth the present invention other in and/or advantage, some will be clearly by describing, or can learn through enforcement of the present invention.
Embodiment
Now, describe exemplary embodiment of the present invention in detail, its example represents in the accompanying drawings, and wherein, identical label represents identical parts all the time.
Fig. 1 is the process flow diagram that marine strata in south China near surface structure data processing method is according to an exemplary embodiment of the present invention shown.
With reference to Fig. 1, in step S110, near surface structure data are gathered to the near surface structure points for investigation in work area.Here, described near surface structure data can comprise at least one in the speed data of each layer of near surface, thickness data and water table depth data, but should understand and the present invention is not limited thereto, described near surface structure data also can comprise other data, such as Lithology data, topographic condition data etc.; Can in enumeration district every square kilometre lay 0.25 ~ 1 near surface structure points for investigation, but the present invention is not limited thereto, those skilled in the art can in the light of actual conditions adjust points for investigation with needs completely accordingly.In addition, by existing method such as such as little refraction, micro logging etc., near surface structure investigation is gathered, will repeat no more here.
In step S120, the data interval of the near surface structure data corresponding at least one feature of the non-points for investigation in the non-image data district in work area is found near surface structure data sample storehouse, and by the near surface structure data constraint of non-points for investigation in described data interval.Here, the feature of described non-points for investigation can include but not limited to rock property, topographic condition and the water table degree of depth, described topographic condition can include but not limited to earth's surface elevation, the gradient and the distance from cliff, and these features of non-points for investigation can be determined by methods such as geologic examinations by those skilled in the art; Described data interval can be the set of multiple value data, and those skilled in the art can distribute based on the situation of change of the needs of data processing or near surface structure described in non-points for investigation.In addition, can retrain, but this is only example according to the order of rock property, topographic condition and the water table degree of depth near surface structure data, those skilled in the art completely can adjust the order of constraint and change as required; Only exemplarily, if using the rock property of a certain non-points for investigation as feature, and described rock property is limestone, then near surface structure data sample storehouse, obtain the data interval of the near surface structure data corresponding to limestone (such as, speed data is interval, thickness data is interval and the depth data of water table interval etc.), and by the near surface structure data constraint of described non-points for investigation in the data interval obtained.
In step S130, whether the variation range of the described data interval in determining step S120 is less than predetermined value.Here, those skilled in the art can arrange described predetermined value according to actual needs, such as, the predetermined value of 100m/s (meter per second) can be set speed data, thickness data is arranged to the predetermined value of 1m (rice), and the predetermined value of 0.5m (rice) is arranged to the water table degree of depth.
If the variation range of described data interval is less than predetermined value, then in step S140, the data interval of the near surface structure data of utility efficiency point and the near surface structure data of non-points for investigation obtains the near surface structure data profile in work area.Here, the near surface structure data profile in work area can include but not limited at least one in the depth data distribution plan of the speed data distribution plan of each layer of near surface, thickness data distribution plan and water table.In addition, obtain near surface structure data profile by existing methods such as such as data Cheng Tu, will repeat no more here.
If the variation range of described data interval is not less than predetermined value, the feature of non-points for investigation is then increased in step S150, and operation returns step S120, namely, find near surface structure data sample storehouse the feature increased described in the described characteristic sum of non-points for investigation the data interval of common corresponding near surface structure data, to retrain further the near surface structure data of non-points for investigation, until the variation range of data interval after constraint is less than predetermined value.
Only exemplarily, if after being retrained by rock property of above-mentioned a certain non-points for investigation, obtain the data interval of near surface structure data corresponding to limestone and the variation range of this data interval is not less than predetermined value, then can proceed constraint to described non-points for investigation, such as, if the elevation of described non-points for investigation is 800m, then can obtain near surface structure data sample storehouse rock property be limestone and elevation be the data interval of near surface structure data corresponding to 800m (such as, speed data is interval, thickness data depth data that is interval and water table is interval), if the variation range of described data interval is still not less than predetermined value, then proceed constraint, such as, if the water table degree of depth of described non-points for investigation is 10m, then can obtain rock property near surface structure data sample storehouse is limestone, elevation is 800m and the data interval of the near surface structure data of the water table degree of depth corresponding to 10m, etc., until the variation range of data interval after constraint is less than predetermined value be less than predetermined value.
Should understand, the data interval of the near surface structure data of non-points for investigation recited above can be the set of at least one value data, and only exemplarily: if the data interval of the near surface structure data of non-points for investigation is not the set of value data, then by such as getting intermediate value to obtain value data to two endpoint values of the data interval after the constraint being less than predetermined value, and the near surface structure data of utility efficiency point and described value data obtain the near surface structure data profile in work area.
According to exemplary embodiment of the present invention, traditional pure mathematics numerical interpolation method is replaced by adopting the interpolation method of band geological condition constraints, introduce statistical principle simultaneously, thus geology correlativity is strengthened greatly, improve marine strata in south China near surface structure enquiry data processing accuracy, solve the problem that marine strata in south China near surface structure enquiry data process data are inaccurate, also indirectly improve the precision that seismic prospecting later stage field construction excites receiving parameter design accuracy and the indoor process static correction of seismic data, have workable, the feature that precision is high, in addition, the present invention is not only applicable to formation of marine facies in southern China area, can also promote in other similar area.
Although illustrate and describe the present invention with reference to certain exemplary embodiments, but it will be apparent to one skilled in the art that the various changes can made when not departing from the spirit and scope of the present invention of scope by claim and equivalents thereof in form and details.