Background technology
The near surface structure investigation is an important step of seismic prospecting, all plays an important role in instructing the open-air earthquake stimulation reception parameter designing of producing of seismic prospecting and indoor earthquake data static corrections processing.The data processing method of the near surface structure investigation of adopting at present is mainly by 0.25~1 near surface structure points for investigation of every square kilometre of laying in enumeration district, gather the near surface structure data, and when processing, adopt the pure mathematics data interpolating method, the zone that does not gather the near surface structure data is generated the Similarity interpolation data automatically, finally produce and generate the near surface structure figure in the whole enumeration district, wherein, described near surface structure data owner will comprise speed data and the thickness data of each layer of near surface, simultaneously also can comprise other data, such as the depth data of water table, the Lithology data, topographic condition data etc.Said method changes the little north and the little Continental Facies Stratigraphy district of rock property variation in the earth's surface landform, and the precision of the near surface structure data that calculate is relatively high.
Yet marine strata in south China has following features: the 1. old stratum of the main exposure in marine deposit stratum, mostly be carbonatite (that is, limestone) and sandstone formation crop out, and near-surface velocity and variation in thickness are large; 2. southern landform is mainly the mountain region, and surface relief is large, and the structural change of near surface is fast.Therefore, the near surface structure data that adopt the data processing method of the near surface structure investigation of above-mentioned routine to calculate can not truly reflect the structural change situation of near surface, and error is larger, thus the needs of seismic exploration fully 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: the near surface structure points for investigation in the work area is gathered the near surface structure data; From near surface structure data sample storehouse, seek the data interval of the corresponding near surface structure data of at least one feature of the not points for investigation in the work area, and not the near surface structure data constraint of points for investigation in described data interval; Judge that 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, the feature that then increases points for investigation not is so that the near surface structure data of points for investigation are not further retrained, until the data interval after the constraint is less than predetermined value; And with the near surface structure data of points for investigation and not the data interval of the near surface structure data of points for investigation obtain near surface structure data profile in the work area.
The near surface structure data can comprise at least one in speed data, thickness data and the water table depth data of each layer of near surface.
The predetermined value of 100m/s can be set speed data, thickness data be arranged the predetermined value of 1m, and the water table degree of depth be arranged the predetermined value of 0.5m.
Can be in advance gather and obtain near surface structure data sample storehouse by various near surface structures being carried out the near surface structure investigation, perhaps with the near surface structure data of the near surface structure points for investigation in the work area as near surface structure data sample storehouse.
The feature of points for investigation can not comprise rock property, topographic condition and the water table degree of depth.
Topographic condition can comprise earth's surface elevation, the gradient and from the distance of cliff.
Near surface structure data profile in the work area can comprise at least one in speed data distribution plan, thickness data distribution plan and the water table depth data distribution plan of each layer of near surface.
But the situation of change of the needs processed of points for investigation based on data or near surface structure and being distributed not.
Data interval can be the set of at least one value data.
Can be by obtaining value data to getting intermediate value less than two endpoint values of the data interval after the constraint of predetermined value, and obtain near surface structure data profile in the work area with the near surface structure data of points for investigation and described value data.
Will be in ensuing description part set forth the present invention other aspect and/or advantage, some will be clearly by describing, and perhaps 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 illustrates the according to an exemplary embodiment of the present invention process flow diagram of marine strata in south China near surface structure data processing method.
With reference to Fig. 1, at step S110, the near surface structure points for investigation in the work area is gathered the near surface structure data.Here, described near surface structure data can comprise at least one in speed data, thickness data and the water table depth data of each layer of near surface, but should understand and the invention is not restricted to this, described near surface structure data also can comprise other data, such as Lithology data, topographic condition data etc.; Can be in enumeration district every square kilometre lay 0.25~1 near surface structure points for investigation, but the invention is not restricted to this, those skilled in the art can in the light of actual conditions adjust points for investigation with needs fully accordingly.In addition, can come near surface structure investigation is gathered by existing methods such as little refraction, micro logging, will repeat no more here.
At step S120, from near surface structure data sample storehouse, seek the data interval of the corresponding near surface structure data of at least one feature of the not points for investigation in the not image data district in the work area, and not the near surface structure data constraint of points for investigation in described data interval.Here, the feature of described not 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 from the distance of cliff, and these features of 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 a plurality of value data, but and the situation of change of the needs processed of those skilled in the art's based on data or the near surface structure described not points for investigation that distributes.In addition, can come the near surface structure data are retrained according to the order of rock property, topographic condition and the water table degree of depth, but this only is example, those skilled in the art can adjust the order of constraint and change fully as required; Only as example, if with the rock property of a certain not points for investigation as feature, and described rock property is limestone, then near surface structure data sample storehouse, obtain the corresponding near surface structure data of limestone data interval (such as, speed data is interval, thickness data is interval and the depth data interval of water table etc.), and near surface structure data constraint that will described not points for investigation is in the data interval of acquisition.
At step S130, whether the variation range of the described data interval among the 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 the predetermined value of 1m (rice), and the water table degree of depth is arranged the predetermined value of 0.5m (rice).
If the variation range of described data interval is less than predetermined value, then at step S140, with the near surface structure data of points for investigation and not the data interval of the near surface structure data of points for investigation obtain near surface structure data profile in the work area.Here, the near surface structure data profile in the work area can include but not limited at least one in the depth data distribution plan of speed data distribution plan, thickness data distribution plan and water table of each layer of near surface.In addition, can by obtaining the near surface structure data profile such as existing methods such as data Cheng Tu, will repeat no more here.
If the variation range of described data interval is not less than predetermined value, then increase the not feature of points for investigation at step S150, and step S120 is returned in operation, namely, from near surface structure data sample storehouse, seek the described feature of points for investigation not and described increase feature the data interval of common corresponding near surface structure data, so that the near surface structure data of points for investigation are not further retrained, until the variation range of the data interval after the constraint is less than predetermined value.
Only as example, if above-mentioned a certain not points for investigation pass through rock property constraint after, obtain the data interval of the corresponding near surface structure data of limestone and the variation range of this data interval and be not less than predetermined value, then can proceed constraint to described not points for investigation, such as, if the elevation of described not points for investigation is 800m, then can near surface structure data sample storehouse, obtain rock property and be limestone and elevation and be the corresponding near surface structure data of 800m data interval (such as, speed data is interval, the depth data interval of thickness data interval and water table etc.), if the variation range of described data interval still is not less than predetermined value, then proceed constraint, such as, if the water table degree of depth of described not points for investigation is 10m, then can obtain rock property near surface structure data sample storehouse is limestone, elevation is that 800m and the water table degree of depth are the data interval of the corresponding near surface structure data of 10m, etc., until the variation range of data interval after the constraint less than predetermined value less than predetermined value.
Should understand, the data interval of the near surface structure data of not points for investigation recited above can be the set of at least one value data, and only as example be: if the data interval of the near surface structure data of points for investigation is not the set of value data, then can be by such as obtaining value data to getting intermediate value less than two endpoint values of the data interval after the constraint of predetermined value, and obtain near surface structure data profile in the work area with the near surface structure data of points for investigation and described value data.
According to exemplary embodiment of the present invention, by the interpolation method replacement traditional pure mathematics numerical interpolation method of employing with geological condition constraints, introduced simultaneously statistical principle, thereby the geology correlativity is strengthened greatly, improved marine strata in south China near surface structure enquiry data processing accuracy, solved the inaccurate problem of marine strata in south China near surface structure enquiry data deal with data, also indirectly improve seismic prospecting later stage field construction and excited the precision that receives parameter designing precision and the indoor processing static correction of seismic data, have workable, the characteristics that precision is high; In addition, the present invention is not only applicable to the formation of marine facies in southern China area, can also promote in other similar area.
Although illustrated and described the present invention with reference to certain exemplary embodiments, but it will be apparent to one skilled in the art that the various changes that in the situation that does not break away from the spirit and scope of the present invention that scope limits by claim and equivalent thereof, can make on form and the details.