CN108873059A - Construct the method and device of loess thickness model - Google Patents
Construct the method and device of loess thickness model Download PDFInfo
- Publication number
- CN108873059A CN108873059A CN201810592993.4A CN201810592993A CN108873059A CN 108873059 A CN108873059 A CN 108873059A CN 201810592993 A CN201810592993 A CN 201810592993A CN 108873059 A CN108873059 A CN 108873059A
- Authority
- CN
- China
- Prior art keywords
- loess
- elevation
- shot point
- rock
- point
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/104—Generating seismic energy using explosive charges
- G01V1/13—Arrangements or disposition of charges to produce a desired pattern in space or time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/129—Source location
- G01V2210/1299—Subsurface, e.g. in borehole or below weathering layer or mud line
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Radar, Positioning & Navigation (AREA)
- Processing Or Creating Images (AREA)
Abstract
The embodiment of the present application provides a kind of method and device for constructing loess thickness model, and this method includes：The remote sensing image comprising altitude data in loess tableland work area is obtained, and determines that the surface layer loess-rock interface in the loess tableland work area is distributed according to the remote sensing image；The quantity and distribution at control point are determined according to the surface layer interface distributions；Determine loess-rock interface elevation at the control point；It is carried out at interpolation according to the loess at the control point-rock interface elevation or continuation is managed, obtain loess-rock initial interface elevation in the loess tableland work area at each shot point；Loess at each shot point-rock initial interface elevation is smoothed, the corresponding loess-rock interface elevation obtained at each shot point；It determines the earth's surface elevation at each shot point, according to the loess at each shot point-rock interface elevation and the corresponding loess thickness determined at each shot point of earth's surface elevation, and constructs loess thickness model accordingly.The modeling efficiency of loess thickness model can be improved in the embodiment of the present application.
Description
Technical field
This application involves seismic exploration technique fields, more particularly, to a kind of method and dress for constructing loess thickness model
It sets.
Background technique
Loess tableland seismic data quality is mainly restricted by loose loess formation attenuation by absorption and strong surface wave interference.To improve money
Item matter generally requires by the way of shooting on group in this area, according to the loess thickness of shot point design point by point combination well number,
Well depth and dose equal excitation parameter.And before loess tableland seismic prospecting wellbore construction, need to complete full work area loess thickness mode
Type building is to get the loess thickness for arriving each excitation point, for instructing the point-by-point design of Huangtuyuan area shot point excitation condition.
It is mainly at present lithology investigation method in the method for Huangtuyuan area building loess thickness model, one is designed in work area
The points for investigation for determining density is dug a well using Luoyang Spade and investigates the loess thickness for obtaining points for investigation, is similar to the point as control point
A certain range of loess thickness, but by loess thickness change is larger and cost is influenced, the depth of near-surface investigation point and close
Degree is all unable to satisfy the exact requirements for implementing loess thickness.
Summary of the invention
The embodiment of the present application is designed to provide a kind of method and device for constructing loess thickness model, to improve building
The efficiency and precision of loess thickness model.
In order to achieve the above objectives, on the one hand, the embodiment of the present application provides a kind of method for constructing loess thickness model, packet
It includes：
The remote sensing image comprising altitude data in loess tableland work area is obtained, and the loess is determined according to the remote sensing image
The surface layer loess in plateau work area-rock interface distribution；
The quantity and distribution for determining control point are distributed according to the surface layer loess-rock interface；
Determine loess-rock interface elevation at the control point；
Interpolation or continuation processing are carried out to the loess at the control point-rock interface elevation, determine the loess tableland work
Loess-rock initial interface elevation in area at each shot point；
Loess at each shot point-rock initial interface elevation is smoothed, it is corresponding to obtain at each shot point
Loess-rock interface elevation；
The earth's surface elevation at each shot point is determined, according to the loess at each shot point-rock interface elevation and earth's surface elevation
The corresponding loess thickness determined at each shot point, and loess thickness model is constructed accordingly.
In the method for the building loess thickness model of the embodiment of the present application, described determined according to the surface layer interface distributions is controlled
The quantity and distribution of point are made, including：
According in the surface layer interface distributions, loess-rock interface exposure region inner boundary fluctuating situation determines control point
Quantity and distribution.
In the method for the building loess thickness model of the embodiment of the present application, loess-rock at the determination control point
Stone interface elevation, including：
From the loess-rock interface elevation extracted in remote sensing image from control point.
In the method for the building loess thickness model of the embodiment of the present application, the loess-rock at the control point
Interface elevation carries out interpolation or continuation processing, determines loess-rock initial interface in the loess tableland work area at each shot point
Elevation, including：
Determine that a shot point and its periphery meet the control point of preset condition；
Using specified interpolation algorithm or continuation algorithm, default item is met with the position coordinates of the shot point, the shot point periphery
For the position coordinates and loess at the control point of part-rock interface elevation as input, the loess-rock interface obtained at the shot point is high
Journey；
And so on, obtain loess-rock interface elevation at each shot point.
In the method for the building loess thickness model of the embodiment of the present application, the earth's surface elevation at each shot point of determination,
According to the loess at each shot point-rock interface elevation and the corresponding loess thickness determined at each shot point of earth's surface elevation, including：
From the earth's surface elevation extracted in remote sensing image from each shot point；
With the corresponding loess-rock interface elevation subtracted at each shot point of earth's surface elevation at each shot point, obtain each
Loess thickness at shot point.
On the other hand, the embodiment of the present application also provides a kind of devices for constructing loess thickness model, including：
Surface layer interface obtains module, for obtaining the remote sensing image comprising altitude data in loess tableland work area, and according to institute
State the surface layer interface distributions that remote sensing image determines the loess tableland work area；
Control point determining module, for determining the quantity and distribution at control point according to the surface layer interface distributions；
First interface elevation determining module, for determining loess-rock interface elevation at the control point；
Second contact surface elevation determining module, for the loess at the control point-rock interface elevation carry out interpolation or
Continuation processing, determines loess-rock initial interface elevation in the loess tableland work area at each shot point；
Interface elevation Leveling Block, for being smoothed to the loess at each shot point-rock initial interface elevation,
Corresponding loess-rock interface the elevation obtained at each shot point；
Thickness model constructs module, for determining the earth's surface elevation at each shot point, according to loess-rock at each shot point
Stone interface elevation and the corresponding loess thickness determined at each shot point of earth's surface elevation, and loess thickness model is constructed accordingly.
In the device of the building loess thickness model of the embodiment of the present application, described determined according to the surface layer interface distributions is controlled
The quantity and distribution of point are made, including：
According in the surface layer interface distributions, loess-rock interface exposure region inner boundary fluctuating situation determines control point
Quantity and distribution.
In the device of the building loess thickness model of the embodiment of the present application, loess-rock at the determination control point
Stone interface elevation, including：
From the loess-rock interface elevation extracted in remote sensing image from control point.
In the device of the building loess thickness model of the embodiment of the present application, the loess-rock at the control point
Interface elevation carries out interpolation or continuation processing, determines loess-rock initial interface in the loess tableland work area at each shot point
Elevation, including：
Determine that a shot point and its periphery meet the control point of preset condition；
Using specified interpolation algorithm or continuation algorithm, default item is met with the position coordinates of the shot point, the shot point periphery
For the position coordinates and loess at the control point of part-rock interface elevation as input, the loess-rock interface obtained at the shot point is high
Journey；
And so on, obtain loess-rock interface elevation at each shot point.
In the device of the building loess thickness model of the embodiment of the present application, the earth's surface elevation at each shot point of determination,
According to the loess at each shot point-rock interface elevation and the corresponding loess thickness determined at each shot point of earth's surface elevation, including：
From the earth's surface elevation extracted in remote sensing image from each shot point；
With the corresponding loess-rock interface elevation subtracted at each shot point of earth's surface elevation at each shot point, obtain each
Loess thickness at shot point.
As can be seen from the technical scheme provided by the above embodiments of the present application, the embodiment of the present application can utilize software and loess tableland work
The remote sensing image comprising altitude data in area determines the loess thickness at each shot point, and constructs loess thickness model accordingly, with
The loess thickness and modeling scheme on the basis of investigation of existing manual research point are compared, and building loess thickness is substantially increased
The efficiency of model.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The some embodiments recorded in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, it is also possible to obtain other drawings based on these drawings.In the accompanying drawings：
Fig. 1 is the flow chart that the method for loess thickness model is constructed in one embodiment of the application；
Fig. 2 is the schematic diagram for determining control point in one embodiment of the application in aerophotograph；
Fig. 3 is the space interpolation schematic diagram in one embodiment of the application；
Fig. 4 is that the control point of one embodiment Air China on piece of the application is distributed；
Fig. 5 is that the shot point of one embodiment Air China on piece of the application is distributed；
Fig. 6 is the structural block diagram that the device of loess thickness model is constructed in one embodiment of the application；
Fig. 7 is the structural block diagram that the device of loess thickness model is constructed in another embodiment of the application.
Specific embodiment
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in example is applied, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.Based on the embodiment in the application, this field is common
The application protection all should belong in technical staff's every other embodiment obtained without creative efforts
Range.
Refering to what is shown in Fig. 1, the method for the building loess thickness model of the embodiment of the present application may comprise steps of：
S101, the remote sensing image comprising altitude data for obtaining loess tableland work area, and institute is determined according to the remote sensing image
State the surface layer interface distributions in loess tableland work area.
In recent years, with the reduction of remote sensing image cost, the remote sensing images such as high-precision aerophotograph using increasingly extensive.
The plane precision of high-precision aerophotograph can achieve 0.25 × 0.25m at present, and height accuracy reaches 2m, and therefore, aerophotograph precision can be with
Meet seismic prospecting production requirement well.During realizing the application, the inventor of the present application discovered that：Due to from aerophotograph
On can intuitively identify the interface of loess and lower overlying strata stone, therefore by picking up loess-rock interface elevation information, it can structure
Build loess thickness model.Wherein, interface of the loess-rock interface between Huangtuyuan area upper layer loess and underlying rock.
Therefore, in some embodiments of the application, the remote sensing image comprising altitude data in loess tableland work area can be obtained first, and according to
Remote sensing image can be clearly seen the surface layer loess in the loess tableland work area-rock interface distribution, for example, interface undulation form, walk
To, quantity, size, whether there is or not loess-rock interface exposures etc..
S102, the quantity and distribution that control point is determined according to the surface layer interface distributions.
In the embodiment of the present application, the control point is the reference point selected from remote sensing image, which is from
(such as shown in Fig. 2) selected on loess-rock interface of exposing, for the Huang at the subsequent calculating each shot point in loess tableland work area
Soil-rock interface elevation.
It is described to determine that the quantity at control point and distribution according to the surface layer interface distributions in some embodiments of the application
To be according in the surface layer interface distributions, loess-rock interface exposure region interface fluctuating situation determines the quantity at control point
And distribution.Such as in loess-rock interface exposure region, each protrusion can at least select a control point；When protrusion is certain
Partial fluctuations than it is more gentle when, control point can be suitably reduced on the part, to reduce while guaranteeing modeling accuracy
Calculation amount and cost of implementation；When the fluctuations of the certain parts of protrusion are bigger, it can be fitted on the part and increase control point, to mention
High modeling accuracy, to may make the model of foundation more true and reliable.In an illustrative embodiments, according to surface layer interface point
Cloth can determine that the quantity at control point and distribution are as shown in Figure 4；In Fig. 4, grey dot indicates control point.
S103, loess-rock interface elevation at the control point is determined.
In some embodiments of the application, since remote sensing image itself includes the information such as altitude data, work as control point
After place is selected, correspondingly, can directly extract the loess-rock interface elevation and position coordinates letter at control point by remote sensing image
Breath.
S104, interpolation or continuation processing are carried out to the loess at the control point-rock interface elevation, determine the loess
Loess-rock initial interface elevation in plateau work area at each shot point.
In some embodiments of the application, the quantity of shot point and distribution are the needs according to seismic prospecting in loess tableland work area
It has been predefined that, therefore, the position of each shot point, example can be marked out in remote sensing image according to the position coordinates of these shot points
As shown in figure 5, black dot indicates shot point in Fig. 5.Under normal circumstances, partial region belongs to Huang in loess tableland work area
Soil-rock interface exposure region.Therefore, those are located at the shot point in the non-exposure region of loess-rock interface, although and position
In loess-rock interface exposure region, but the position not shot point on loess-rock interface, it can not be straight by remote sensing image
Connect the loess-rock interface elevation extracted at corresponding shot point.This just needs to predict at its position so that control point is as a reference point
Loess-rock interface elevation.
In some embodiments of the application, using the loess at control point-rock interface elevation, and pass through interpolation algorithm
Or continuation algorithm predicts the loess at shot point-rock interface elevation, specifically：
Firstly, determining that a shot point and its periphery meet the control point of preset condition；In some exemplary embodiments, institute
It states and meets preset condition for example and can be control point and shot point is located in the same protrusion, and/or, the spacing at control point and shot point
In default spacing range etc..
Then, using specified interpolation algorithm or continuation algorithm, met with the position coordinates of the shot point, the shot point periphery pre-
If the position coordinates and loess at the control point of condition-rock interface elevation obtain loess-rock circle at the shot point as input
Face elevation；
And so on, repeating above step finally can get the loess at each shot point-rock interface elevation.
In one exemplary embodiment, as shown in figure 3, control point 1., 2. 3. control point belongs to shot point with control point 4. all
While meet the control point of preset condition, extracted from aerophotograph obtain control point 1., control point 2. with control point 3. from loess-
After rock interface elevation, the loess-rock interface elevation of shot point 4. can be calculated by interpolation algorithm.In an exemplary implementation
In example, there are 15332 shot points (as shown in Figure 5) in certain loess tableland work area, lays 646, control point altogether in the work area and extract and obtain
The interface elevation (as shown in Figure 4) for obtaining this 646 points, the loess-rock of 15332 shot points of the whole district can be obtained by interpolation algorithm
Stone interface elevation.
In some embodiments of the application, preferably, for 3-d seismic exploration, may be selected spatial interpolation algorithm (such as
Kriging regression method, nearest neighbor point interpolation method etc.) etc. primary Calculations obtain loess-rock interface elevation of each shot point of the whole district；
And for two-dimensional seismic survey, the loess-of each shot point of the whole district can be obtained by primary Calculations such as linear interpolation or continuation algorithms
Rock interface elevation.
S105, the loess at each shot point-rock initial interface elevation is smoothed, it is corresponding to obtain each shot point
The loess at place-rock interface elevation.
It, can be to the Huang at each shot point of previous step acquisition in order to improve modeling accuracy in some embodiments of the application
Soil-rock interface elevation is smoothed.Mean value smoothing method can be used in some exemplary embodiments to carry out smoothly, such as
Different smooth radiuses (200m, 300m, 400m etc.) can be selected high to the initial loess at the shot point of part-rock interface respectively
Cheng Jinhang is smooth, and the control point interface elevation that Statistical Comparison picks up after smooth takes error in the error size of smooth front and back
Reckling is as suitable smooth radius, and then the suitable smooth radius of selection is to initial loess-rock circle at all shot points
Face elevation is smoothed, to obtain loess final at all shot points-rock interface elevation.
S106, earth's surface elevation at each shot point is determined, according to the loess at each shot point-rock interface elevation and earth's surface
The corresponding loess thickness determined at each shot point of elevation, and loess thickness model is constructed accordingly.
In some embodiments of the application, earth's surface elevation at each shot point of determination, according to the Huang at each shot point
Soil-rock interface elevation and the corresponding loess thickness determined at each shot point of earth's surface elevation for example can be and first mentions from remote sensing image
The earth's surface elevation at each shot point is taken, then with the corresponding loess-rock subtracted at each shot point of earth's surface elevation at each shot point
Stone interface elevation, obtains the loess thickness at each shot point, and specific formula for calculation is as follows：
Tij=Hij1-Hij2
Wherein, TijLoess thickness i.e. at point (i, j), Hij1Earth's surface elevation i.e. at point (i, j), Hij2I.e. at point (i, j)
Loess-rock interface elevation, (i, j) is geographical coordinate, and i is northern coordinate, and j is eastern coordinate.
And on the basis of determining the loess thickness at each shot point, it can easily construct corresponding loess thickness mode
Type.
In other embodiments of the application, after establishing loess thickness model, loess thickness model can also be verified
With precision analysis.Such as can be by a certain number of shot points near remote sensing image pickup loess-rock interface as check post, it will
The earth's surface elevation of these shot points subtracts each other with the loess near it-rock interface elevation (can directly extract on remote sensing image), from
And accurate loess thickness at these shot points is obtained, and remove in verifying model loess at these check posts in this, as reference thickness
The accuracy of thickness.The loess of corresponding position is thick in loess thickness at these check posts, with above-mentioned loess thickness model
Error between degree within the allowable range when, then can be using above-mentioned loess thickness model as final loess thickness model.And work as
Error in loess thickness at these check posts, with above-mentioned loess thickness model between the loess thickness of corresponding position exceeds
, can be on the basis of changing number of control points and position when allowed band (such as error is no more than 1% etc.), and/or changing
On the basis of smoothing parameter, loess Thickness modeling is re-started, and carries out verifying and precision analysis again, is wanted until obtaining and meeting
Until the loess thickness model asked.
Refering to what is shown in Fig. 7, a kind of device of building loess thickness model of the embodiment of the present application may include memory, place
Reason device and the computer program being stored on the memory, the computer program execute when being run by the processor
Following steps：
The remote sensing image comprising altitude data in loess tableland work area is obtained, and the loess is determined according to the remote sensing image
The surface layer interface distributions in plateau work area；
The quantity and distribution at control point are determined according to the surface layer interface distributions；
Determine loess-rock interface elevation at the control point；
Interpolation or continuation processing are carried out to the loess at the control point-rock interface elevation, determine the loess tableland work
Loess-rock initial interface elevation in area at each shot point；
Loess at each shot point-rock initial interface elevation is smoothed, it is corresponding to obtain at each shot point
Loess-rock interface elevation；
The earth's surface elevation at each shot point is determined, according to the loess at each shot point-rock interface elevation and earth's surface elevation
The corresponding loess thickness determined at each shot point, and loess thickness model is constructed accordingly.
Although procedures described above process includes the multiple operations occurred with particular order, it should however be appreciated that understand,
These processes may include more or fewer operations, these operations can be executed sequentially or be executed parallel (such as using parallel
Processor or multi-thread environment).
Refering to what is shown in Fig. 6, the device of another building loess thickness model of the embodiment of the present application may include：
Surface layer interface obtains module 61, can be used for obtaining the remote sensing image comprising altitude data in loess tableland work area, and
The surface layer interface distributions in the loess tableland work area are determined according to the remote sensing image；
Control point determining module 62 can be used for determining the quantity and distribution at control point according to the surface layer interface distributions；
First interface elevation determining module 63, the loess-rock interface elevation being determined at the control point；
Second contact surface elevation determining module 64 can be used for carrying out the loess at the control point-rock interface elevation
Interpolation or continuation processing, determine loess-rock initial interface elevation in the loess tableland work area at each shot point；
Interface elevation Leveling Block 65 can be used for carrying out the loess at each shot point-rock initial interface elevation flat
Sliding processing, the corresponding loess-rock interface elevation obtained at each shot point；
Thickness model constructs module 66, the earth's surface elevation being determined at each shot point, at each shot point
Loess-rock interface elevation and the corresponding loess thickness determined at each shot point of earth's surface elevation, and loess thickness mode is constructed accordingly
Type.
For convenience of description, it is divided into various units when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each unit can be realized in the same or multiple software and or hardware when application.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM) and/or
The forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer-readable medium
Example.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data.
The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves
State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable
Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM),
Digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetic storage devices
Or any other non-transmission medium, can be used for storage can be accessed by a computing device information.As defined in this article, it calculates
Machine readable medium does not include temporary computer readable media (transitory media), such as the data-signal and carrier wave of modulation.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
Include, so that process, method or equipment including a series of elements not only include those elements, but also including not having
There is the other element being expressly recited, or further includes for this process, method or the intrinsic element of equipment.Do not having more
In the case where more limitations, the element that is limited by sentence "including a ...", it is not excluded that including process, the side of the element
There is also other identical elements in method or equipment.
It will be understood by those skilled in the art that embodiments herein can provide as method, system or computer program product.
Therefore, complete hardware embodiment, complete software embodiment or embodiment combining software and hardware aspects can be used in the application
Form.It is deposited moreover, the application can be used to can be used in the computer that one or more wherein includes computer usable program code
The shape for the computer program product implemented on storage media (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
Formula.
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The above description is only an example of the present application, is not intended to limit this application.For those skilled in the art
For, various changes and changes are possible in this application.All any modifications made within the spirit and principles of the present application are equal
Replacement, improvement etc., should be included within the scope of the claims of this application.
Claims (10)
1. a kind of method for constructing loess thickness model, which is characterized in that including：
The remote sensing image comprising altitude data in loess tableland work area is obtained, and the loess tableland work is determined according to the remote sensing image
The surface layer loess in area-rock interface distribution；
The quantity and distribution for determining control point are distributed according to the surface layer loess-rock interface；
Determine loess-rock interface elevation at the control point；
Interpolation or continuation processing are carried out to the loess at the control point-rock interface elevation, determined in the loess tableland work area
Loess-rock initial interface elevation at each shot point；
Loess at each shot point-rock initial interface elevation is smoothed, the corresponding loess-obtained at each shot point
Rock interface elevation；
Determine the earth's surface elevation at each shot point, it is corresponding according to the loess at each shot point-rock interface elevation and earth's surface elevation
It determines the loess thickness at each shot point, and constructs loess thickness model accordingly.
2. the method for building loess thickness model as described in claim 1, which is characterized in that described according to the surface layer interface
It is distributed the quantity and distribution for determining control point, including：
According in the surface layer interface distributions, loess-rock interface exposure region inner boundary fluctuating situation determines the number at control point
Amount and distribution.
3. the method for building loess thickness model as described in claim 1, which is characterized in that at the determination control point
Loess-rock interface elevation, including：
From the loess-rock interface elevation extracted in remote sensing image from control point.
4. as described in claim 1 building loess thickness model method, which is characterized in that it is described to the control point at
Loess-rock interface elevation carries out interpolation or continuation processing, determines loess-rock in the loess tableland work area at each shot point
Initial interface elevation, including：
Determine that a shot point and its periphery meet the control point of preset condition；
Using specified interpolation algorithm or continuation algorithm, preset condition is met with the position coordinates of the shot point, the shot point periphery
The position coordinates and loess at control point-rock interface elevation obtain loess-rock interface elevation at the shot point as input；
And so on, obtain loess-rock interface elevation at each shot point.
5. the method for building loess thickness model as described in claim 1, which is characterized in that at each shot point of determination
Earth's surface elevation, according to the loess at each shot point-rock interface elevation and the corresponding loess determined at each shot point of earth's surface elevation
Thickness, including：
From the earth's surface elevation extracted in remote sensing image from each shot point；
With the corresponding loess-rock interface elevation subtracted at each shot point of earth's surface elevation at each shot point, each shot point is obtained
The loess thickness at place.
6. a kind of device for constructing loess thickness model, which is characterized in that including：
Surface layer interface obtains module, for obtaining the remote sensing image comprising altitude data in loess tableland work area, and according to described distant
Sense image determines the surface layer interface distributions in the loess tableland work area；
Control point determining module, for determining the quantity and distribution at control point according to the surface layer interface distributions；
First interface elevation determining module, for determining loess-rock interface elevation at the control point；
Second contact surface elevation determining module, for carrying out interpolation or continuation to the loess at the control point-rock interface elevation
Processing, determines loess-rock initial interface elevation in the loess tableland work area at each shot point；
Interface elevation Leveling Block, it is corresponding for being smoothed to the loess at each shot point-rock initial interface elevation
Obtain loess-rock interface elevation at each shot point；
Thickness model constructs module, for determining the earth's surface elevation at each shot point, according to loess-rock circle at each shot point
Face elevation and the corresponding loess thickness determined at each shot point of earth's surface elevation, and loess thickness model is constructed accordingly.
7. the device of building loess thickness model as claimed in claim 6, which is characterized in that described according to the surface layer interface
It is distributed the quantity and distribution for determining control point, including：
According in the surface layer interface distributions, loess-rock interface exposure region inner boundary fluctuating situation determines the number at control point
Amount and distribution.
8. the device of building loess thickness model as claimed in claim 6, which is characterized in that at the determination control point
Loess-rock interface elevation, including：
From the loess-rock interface elevation extracted in remote sensing image from control point.
9. as claimed in claim 6 building loess thickness model device, which is characterized in that it is described to the control point at
Loess-rock interface elevation carries out interpolation or continuation processing, determines loess-rock in the loess tableland work area at each shot point
Initial interface elevation, including：
Determine that a shot point and its periphery meet the control point of preset condition；
Using specified interpolation algorithm or continuation algorithm, preset condition is met with the position coordinates of the shot point, the shot point periphery
The position coordinates and loess at control point-rock interface elevation obtain loess-rock interface elevation at the shot point as input；
And so on, obtain loess-rock interface elevation at each shot point.
10. the device of building loess thickness model as claimed in claim 6, which is characterized in that at each shot point of determination
Earth's surface elevation, according to the corresponding Huang determined at each shot point of the loess at each shot point-rock interface elevation and earth's surface elevation
Soil thickness, including：
From the earth's surface elevation extracted in remote sensing image from each shot point；
With the corresponding loess-rock interface elevation subtracted at each shot point of earth's surface elevation at each shot point, each shot point is obtained
The loess thickness at place.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810592993.4A CN108873059A (en) | 2018-06-11 | 2018-06-11 | Construct the method and device of loess thickness model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810592993.4A CN108873059A (en) | 2018-06-11 | 2018-06-11 | Construct the method and device of loess thickness model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108873059A true CN108873059A (en) | 2018-11-23 |
Family
ID=64337746
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810592993.4A Pending CN108873059A (en) | 2018-06-11 | 2018-06-11 | Construct the method and device of loess thickness model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108873059A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114429728A (en) * | 2020-10-15 | 2022-05-03 | 中国石油化工股份有限公司 | Loess tableland stratum physical model material, loess tableland stratum physical model and preparation method |
-
2018
- 2018-06-11 CN CN201810592993.4A patent/CN108873059A/en active Pending
Non-Patent Citations (3)
| Title |
|---|
| 刘红星 等: "综合高精度航片的黄土厚度模型构建及应用", 《2016年煤层气学术研讨会论文集》 * |
| 肖梦雄 等: "应用高精度航片构建黄土厚度模型研究", 《2016中国地球科学联合学术年会论文集（十九）——专题40：油气田与煤田地球物理勘探》 * |
| 邓域才: "《铁路规划与机助设计》", 30 April 1996, 中国铁道出版社 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114429728A (en) * | 2020-10-15 | 2022-05-03 | 中国石油化工股份有限公司 | Loess tableland stratum physical model material, loess tableland stratum physical model and preparation method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106990438B (en) | Complicated earth surface region micro logging point and well depth determine method and device | |
| AU2010308501B2 (en) | Method for optimization with gradient information | |
| AU2011339017B2 (en) | Constructing geologic models from geologic concepts | |
| CN106896406B (en) | Method and apparatus based on impedance Domain Properties prediction High-quality Reservoir | |
| US9915742B2 (en) | Method and system for geophysical modeling of subsurface volumes based on label propagation | |
| CN107783186B (en) | A kind of near-surface velocity modeling method and device | |
| CN103077558A (en) | Method for establishing large-sized karst cave reservoir body distribution model of fractured-vuggy carbonate rock reservoir | |
| CN107817525B (en) | A kind of near surface constant speed fill method and device | |
| US10387583B2 (en) | Rotations from gradient directions | |
| WO2014099204A1 (en) | Method and system for geophysical modeling of subsurface volumes based on computed vectors | |
| CN108254784A (en) | A kind of fault recognizing method, apparatus and system based on two-dimension earthquake data | |
| CN107144877B (en) | Anomalous structure recognition methods and device | |
| CN102057368A (en) | Distribution of properties in a 3D volumetric model using a maximum continuity field | |
| CN108445532B (en) | A kind of Depth Domain inverse migration method and device | |
| CN109444955A (en) | Interpolation method is disturbed when the bilinearity of three dimensional seismic raytracing is walked | |
| CN107956465A (en) | The more borehole logging tool log normalization method and devices of the whole district based on association well | |
| CN105549084A (en) | Method and system for building three-dimensional high-precision velocity model | |
| CN108845350A (en) | The method and device of inverting two-dimension speed model | |
| CN109975872B (en) | Method for directly compiling balance profile by using Corel DRAW software | |
| CN109598068B (en) | Ancient structure constraint modeling method, device and equipment | |
| CN108873059A (en) | Construct the method and device of loess thickness model | |
| CN106842314A (en) | The determination method of formation thickness | |
| CN106443829A (en) | Method and apparatus for constructing near-surface model | |
| CN109633781A (en) | Geological property acquisition methods and device, electronic equipment, storage medium | |
| CN107797148B (en) | A kind of aeromagnetic anomaly field separation method and system based on three-dimensional geological modeling |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181123 |