CN107561574A - Shot position, the detection method of receiving point position and device - Google Patents

Shot position, the detection method of receiving point position and device Download PDF

Info

Publication number
CN107561574A
CN107561574A CN201710594250.6A CN201710594250A CN107561574A CN 107561574 A CN107561574 A CN 107561574A CN 201710594250 A CN201710594250 A CN 201710594250A CN 107561574 A CN107561574 A CN 107561574A
Authority
CN
China
Prior art keywords
grid chart
shot
receiving point
detection
detection grid
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
Application number
CN201710594250.6A
Other languages
Chinese (zh)
Inventor
杨敬磊
朱旭江
李伟波
柳兴刚
吴蔚
李永革
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201710594250.6A priority Critical patent/CN107561574A/en
Publication of CN107561574A publication Critical patent/CN107561574A/en
Pending legal-status Critical Current

Links

Abstract

The application embodiment provides a kind of shot position, the detection method of receiving point position and device, wherein, this method includes:Obtain shot position, the receiving point position in target area;According to shot position, receiving point position, detection grid chart is made;Shot point number and receiving point number in statistical unit grid;Shot point number and receiving point number in unit grids, it is determined that the quadrat method of detection grid chart;It is determined that the quadrat method of detection grid chart is located at shot position, receiving point position corresponding to the detection grid chart in preset range, meet preset requirement.Because the program is by detecting grid chart corresponding to making, and quadrat method is determined according to detection grid chart, and then judge whether target area shot position, receiving point position meet preset requirement using quadrat method, therefore solve whether can not meet that randomness requires to carry out the technical problem of accurate judgement present in existing method to shot position, receiving point position.

Description

Shot position, the detection method of receiving point position and device
Technical field
The application is related to technical field of geophysical exploration, the detection of more particularly to a kind of shot position, receiving point position Method and apparatus.
Background technology
During geophysical exploration is carried out, it is often necessary to first pass through big gun inspection method and target area progress data are adopted Collection, and then carry out specific data reconstruction according to the data of collection.This method is when it is implemented, in order to reach preferably reconstruct effect Fruit, generally require to lay substantial amounts of shot point, receiving point, while can also produce larger data processing amount.
In order to take into account production cost, i.e.,:Shot point and receiving point are arranged as few as possible;Meanwhile also for raising data The effect of reconstruct.In recent years, the development based on compressive sensing theory, it is proposed that a kind of new data acquisition, the side of data reconstruction Method.That is, data acquisition is carried out by irregular sampling (stochastical sampling), passes through the data recovering algorithms based on compressive sensing theory Carry out corresponding data reconstruction.By the above method, irregular sampling can be utilized on the premise of shot point and receiving point is reduced The characteristics of alias is converted into the low amplitude value noise for being easy to filter out, strategy is promoted to enter the lack sampling data gathered using sparse The accurate data reconstruction of row.But the above method is when it is implemented, higher for shot position, receiving point status requirement, Shot position, receiving point position need to meet preset requirement, i.e., spatially shot point, receiving point are random distributions.
But a kind of specific side effectively accordingly evaluated shot position, receiving point position is not yet proposed at present Method.That is, whether presence can not be put in existing method meets that randomness requires to carry out accurate judgement to shot position, receiving point position Technical problem.
In view of the above-mentioned problems, not yet propose effective solution at present.
The content of the invention
The application embodiment provides a kind of shot position, the detection method of receiving point position and device, existing to solve Whether have can not meet that the technology that randomness requires to carry out accurate judgement is asked to shot position, receiving point position present in method Topic.
The application embodiment provides a kind of shot position, the detection method of receiving point position, including:
Obtain shot position, the receiving point position in target area;
According to the shot position, the receiving point position, detection grid chart is made, wherein, the detection grid chart bag Include multiple unit grids;
Count shot point number and receiving point number in unit grids in the detection grid chart;
Shot point number and receiving point number in the unit grids, it is determined that the quadrat method of detection grid chart;
According to the quadrat method of the detection grid chart, confirm shot position corresponding to the detection grid chart, receive Whether point position meets preset requirement.
It is in one embodiment, described that grid chart is detected according to the shot position, the receiving point position, making, Including:
According to the shot position, the receiving point position, shooting distance, perpendicular offset of shot line are determined;
According to the shooting distance, the perpendicular offset of shot line, the detection grid chart is made.
In one embodiment, the length value of the unit grids is twice of the shooting distance, the unit grids The width value of point is twice of the perpendicular offset of shot line.
In one embodiment, the shot point number and receiving point number in the unit grids, it is determined that inspection The quadrat method of survey grid trrellis diagram, including:
Shot point number and receiving point number in the unit grids, it is determined that the points average of detection grid chart;
According to the shot point number and receiving point number in the detection points average of grid chart, the unit grids, really The points variance of regular inspection survey grid trrellis diagram;
According to the points average of the detection grid chart and the points variance of the detection grid chart, the detection net is determined The quadrat method of trrellis diagram.
In one embodiment, the shot point number and receiving point number in the unit grids, it is determined that inspection The points average of survey grid trrellis diagram, including:
According to below equation, the points average for detecting grid chart is determined:
In above formula, Mean be detection grid chart points average, XiFor numbering be i unit grids in shot point number and Receiving point number sum, M be detection grid chart in unit grids number, i be unit grid numbering, span be 1 to M。
In one embodiment, the big gun in the points average for detecting grid chart, the unit grids Point number and receiving point number, it is determined that the points variance of detection grid chart, including:
The points variance of the detection grid chart is determined according to below equation:
In above formula, Mean is the points average of detection grid chart, and Var is the points variance of detection grid chart.
In one embodiment, according to the points average of the detection grid chart and the points side of the detection grid chart Difference, the quadrat method of the detection grid chart is determined, including:
The quadrat method of the detection grid chart is determined according to below equation:
In above formula, Mean is the points average of detection grid chart, and Var is the points variance of detection grid chart, and VMR is detection The quadrat method of grid chart.
In one embodiment, according to the quadrat method of the detection grid chart, the detection grid chart pair is confirmed Whether shot position, the receiving point position answered meet preset requirement, including:
In the case where two-dimensional observational system is laid, the quadrat method for detecting grid chart is more than 0.5 and the inspection less than 1 Shot position, receiving point position corresponding to survey grid trrellis diagram, are defined as meeting the preset requirement;
In the case where stereo observing system is laid, the quadrat method for detecting grid chart is more than 0.7 and the inspection less than 1 Shot position, receiving point position corresponding to survey grid trrellis diagram, are defined as meeting the preset requirement.
In one embodiment, confirming that shot position, receiving point position are unsatisfactory for corresponding to the detection grid chart After the preset requirement, methods described also includes:
Distance is lacked by controlling, to the shot position, the receiving point position, carries out randomness adjustment.
In one embodiment, confirming that shot position, receiving point position meet institute corresponding to the detection grid chart After stating preset requirement, methods described also includes:
Excitaton source is laid in the shot position, receiver is laid in the receiving point position;
Stochastical sampling is carried out by the excitaton source, the receiver, obtains lack sampling data;
According to the lack sampling data, data reconstruction is carried out by compression sensing method.
The application embodiment further provides the detection means of a kind of shot position, receiving point position, including:
Acquisition module, for obtaining the shot position in target area, receiving point position;
Drawing module, for detecting grid chart according to the shot position, the receiving point position, making, wherein, it is described Detection grid chart includes multiple unit grids;
Statistical module, for counting shot point number and receiving point number in the detection grid chart in unit grids;
First determining module, for the shot point number and receiving point number in the unit grids, it is determined that detection net The quadrat method of trrellis diagram;
Second determining module, for the quadrat method according to the detection grid chart, confirm the detection grid chart pair Whether shot position, the receiving point position answered meet preset requirement.
In one embodiment, first determining module includes:
First determining unit, for the shot point number and receiving point number in the unit grids, it is determined that detection net The points average of trrellis diagram;
Second determining unit, for the shot point in the points average according to the detection grid chart, the unit grids Number and receiving point number, it is determined that the points variance of detection grid chart;
3rd determining unit, for the points average according to the detection grid chart and the points side of the detection grid chart Difference, determine the quadrat method of the detection grid chart.
In the application embodiment, variance is determined by detecting grid chart corresponding to making, and according to detection grid chart Average ratio, judge whether target area shot position, receiving point position meet preset requirement using quadrat method.Therefore, solve Determine not putting present in existing method and whether met that randomness requirement is accurately sentenced to shot position, receiving point position Disconnected technical problem, technique effect that is cost-effective, improving efficiency of construction is reached.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this Some embodiments described in application, for those of ordinary skill in the art, do not paying the premise of creative labor Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 be according to the application embodiment provide shot position, receiving point position detection method handling process Figure;
Fig. 2 be according to the application embodiment provide shot position, receiving point position detection means composition structure Figure;
Fig. 3 is obtained not using detection method/device of the application embodiment offer shot position, receiving point position With the evaluation result schematic diagram of layout scheme.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation Example only some embodiments of the present application, rather than whole embodiments.It is common based on the embodiment in the application, this area The every other embodiment that technical staff is obtained under the premise of creative work is not made, it should all belong to the application protection Scope.
Effectively shot position can not be directed in view of existing method, whether receiving point position meets applied compression sense Spatial distribution randomness required by knowing theory specifically detect, evaluated.And the specific laying of shot point, receiving point is not for advising Then sample (stochastical sampling), deficiency of data recovers have considerable influence.That is, often presence can not be put to shot point in existing method Whether position, receiving point position meet that randomness requires to carry out the technical problem of accurate judgement.For above-mentioned technical problem, with reference to The Specific construction situation of field of geophysical exploration, the application consideration can first obtain related to shot position, receiving point position Quadrat method, judge whether target area shot position, receiving point position meet preset requirement using the quadrat method. So as to, solve can not be put present in existing method whether shot position, receiving point position are met randomness require into The technical problem of row accurate judgement, technique effect that is cost-effective, improving efficiency of construction is reached.
Based on above-mentioned thinking thinking, the application embodiment provides a kind of shot position, the detection side of receiving point position Method.Refer to Fig. 1 according to the application embodiment provide shot position, receiving point position detection method handling process Figure.The shot position of the application embodiment offer, the detection method of receiving point position, can specifically include following step Suddenly.
Step S101:Obtain shot position, the receiving point position in target area.
In the present embodiment, in some target area before progress data acquisition, it is thus necessary to determine that specific to lay Scheme (for example, blown-out shot scheme, or distribution scheme), specific shot position, receiving point position are determined according to layout scheme.And then Excitaton source can be arranged on corresponding shot position, for example, collection big gun etc.;Arrange and receive on corresponding receiving point position Device, for example, wave detector etc..It is thus possible to the excitaton source and receiver by being arranged in target area carry out meeting compressed sensing The irregular sampling of theory calls, i.e. stochastical sampling.
In the present embodiment, above-mentioned stochastical sampling will can turn due to alias caused by lack part shot point, receiving point The characteristics of turning to the low amplitude value noise for being easy to filter out, then promote strategy to realize the recovery to lack sampling data by sparse, so as to Carry out data reconstruction.Less shot point, receiving point are laid so as to reach, obtains preferable data reconstruction effect.In this reality Apply in mode, it is necessary to illustrate, the alias (Alias) refers to obscuring in frequency caused by data from the sample survey.For example, certain When the sampling number in input signal each cycle of one frequency is less than two, the output end of system will be counted as another frequency letter Number sampling, in turn result in interference.The lack sampling refers to the sampling less than Nyquist sample frequencys, wherein, Nyquist Sample frequency (Nyquist frequency, nyquist frequency) is the frequency, its size etc. named with Harry Nyquist In sampling rate (sampling rate) half.
In the present embodiment, in order that must have the property that can preferably carry out data reconstruction by the data of stochastical sampling Can, it is desirable to which shot position, receiving point position have certain randomness.If shot point, receiving point are to be uniformly distributed or assemble If distribution, by the data that this layout scheme gathers when carrying out data reconstruction, precision can be fewer, and error can compare Greatly, i.e. the effect of data reconstruction is undesirable.Therefore, in the present embodiment, first to obtain in target area in layout scheme Shot position, receiving point position, to be made a concrete analysis of to shot position, receiving point position, according to analysis result, judge big gun Whether the laying of point, receiving point meets certain randomness, and then layout scheme is specifically evaluated.
Step S102:According to the shot position, the receiving point position, detection grid chart is made, wherein, the detection Grid chart includes multiple unit grids.
In one embodiment, in order to be put to emplacement, receiving point position whether meet randomness requirement specifically sentenced Break, it is necessary to according to shot position, receiving point position, first make corresponding detection grid chart.When it is implemented, can be according to following Mode performs.
S102-1:According to the shot position, the receiving point position, shooting distance, perpendicular offset of shot line are determined.
S102-2:According to the shooting distance, the perpendicular offset of shot line, the detection grid chart is made.
In the present embodiment, the detection grid chart includes multiple unit grids, is distributed again in each unit grids There are multiple shot points, receiving point.It should be noted that above-mentioned unit grids are specific in the follow-up processing for carrying out randomness judgement Handle unit.
In one embodiment, in order that obtaining unit grids as specific processing unit, there is preferable precision, simultaneously Treatment effeciency is taken into account, when it is implemented, when making detection grid chart, the unit grids can be specifically a kind of rectangle, Wherein, the length value of unit grids is twice of the shooting distance, and the width value of unit grids point is twice of the perpendicular offset of shot line. So, it is ensured that shot point, receiving point is distributed with each unit grids, while subsequent treatment efficiency can also be improved.Its In, above-mentioned shooting distance, perpendicular offset of shot line can specifically determine according to the shot position, receiving point position.
Step S103:Count shot point number and receiving point number in unit grids in the detection grid chart.
In the present embodiment, when it is implemented, can be using a unit grids as a sample space, respectively to inspection Shot point number in each grid and receiving point number in survey grid trrellis diagram in multiple unit grids are accordingly counted.
Step S104:Shot point number and receiving point number in the unit grids, it is determined that the side of detection grid chart Poor average ratio.
In one embodiment, in order to the layout scheme corresponding to detection grid chart, i.e., specific shot position, reception Whether point position meets that randomness requirement is judged, each unit grids that can be first in detection grid calculate detection The quadrat method of grid chart, judge index is used for the mean variance, shot position, receiving point position carried out specific Judge.In the present embodiment, in order to determine the quadrat method of detection grid chart, when it is implemented, can be according to following step It is rapid to perform.
Step S104-1:Shot point number and receiving point number in the unit grids, it is determined that detection grid chart Points average.
In one embodiment, in order to determine the points average of detection grid chart, when it is implemented, can be according to following Formula, calculate the points average of detection grid chart:
In above formula, Mean be detection grid chart points average, XiFor numbering be i unit grids in shot point number and Receiving point number sum, M be detection grid chart in unit grids number, i be unit grid numbering, span be 1 to M。
Step S104-2:Shot point number in the detection points average of grid chart, the unit grids and connect Sink number, it is determined that the points variance of detection grid chart.
In one embodiment, in order to determine the points variance of detection grid chart, when it is implemented, can be according to following Formula, calculate the points variance of detection grid chart:
In above formula, Mean is the points average of detection grid chart, and Var is the points variance of detection grid chart.
Step S104-3:According to it is described detection grid chart points average and it is described detect grid chart points variance, really The quadrat method of the fixed detection grid chart.
In one embodiment, in order to determine detect grid chart quadrat method, when it is implemented, can according to Lower formula, calculate the quadrat method of detection grid chart:
In above formula, Mean is the points average of detection grid chart, and Var is the points variance of detection grid chart, and VMR is detection The quadrat method of grid chart.
Step S105:According to the quadrat method of the detection grid chart, shot point corresponding to the detection grid chart is confirmed Whether position, receiving point position meet preset requirement.
In the present embodiment, in order to judge layout scheme, i.e. whether shot position, receiving point position meets that randomness will Ask, when it is implemented, mean variance can be used for into reference index.By judging the numerical value of quadrat method whether positioned at pre- If scope, to judge to detect whether shot position, receiving point position corresponding to grid chart meet preset requirement, i.e., whether meet Randomness requirement.Specifically, it may be determined that the quadrat method of the detection grid chart is located at the detection grid in preset range The corresponding shot position of figure, receiving point position meet preset requirement;Determine the quadrat method for detecting grid chart not position Preset requirement is unsatisfactory in the shot position corresponding to the detection grid chart in preset range, receiving point position.
In one embodiment, when it is determined that whether the quadrat method of detection grid chart is located at preset range, consider To specific different construction requirement and situation is laid, corresponding preset range can be set with construction requirement as the case may be, When it is implemented, it can perform in such a way.
S105-1:It is in the case where two-dimensional observational system is laid, the quadrat method for detecting grid chart is more than 0.5 and small The shot position corresponding to detection grid chart, receiving point position in 1, are defined as meeting the preset requirement.That is, for two dimension Preset range, according to construction requirement, can be arranged to by the situation that observation system is laid:0.5 to 1.
S105-2:It is in the case where stereo observing system is laid, the quadrat method for detecting grid chart is more than 0.7 and small The shot position corresponding to detection grid chart, receiving point position in 1, are defined as meeting the preset requirement.That is, for three-dimensional Preset range, according to construction requirement, can be arranged to by the situation that observation system is laid:0.7 to 1.
Certainly, it is necessary to explanation, when it is implemented, can as the case may be and specific requirement, setting removes above-mentioned institute Number range beyond the number range enumerated is as the preset range.On the other hand, the application does not repeat.
In the present embodiment, the two dimension in the two-dimensional observational system can specifically include:Longitudinal dimension spatially And time dimension;Three-dimensional in the stereo observing system can specifically include:Longitudinal dimension spatially, transverse direction spatially Dimension and time dimension.
In one embodiment, when it is implemented, it is considered that the quadrat method of detection grid chart is right closer to 1 institute Shot position, the receiving point position answered more meet preset requirement, i.e., more meet randomness requirement.In this way, it can preferably detect net The layout scheme that the mean variance of trrellis diagram is closer to 1 carries out stochastical sampling.
In the embodiment of the present application, by detecting grid chart corresponding to making, and determine that variance is equal according to detection grid chart It is worth ratio, judges whether target area shot position, receiving point position meet preset requirement using quadrat method.Therefore, solve It can not be put present in existing method and whether meet that to shot position, receiving point position randomness requires to carry out accurate judgement Technical problem, reached it is cost-effective, improve efficiency of construction technique effect.
In one embodiment, after it is determined that shot position, receiving point position be unsatisfactory for the preset requirement, in order that Shot position, receiving point position meet randomness requirement, when it is implemented, methods described can also include:
Distance is lacked by controlling, to the shot position, the receiving point position, carries out randomness adjustment.
In one embodiment, after it is determined that shot position, receiving point position meet the preset requirement, can utilize The layout scheme carries out specific excitation device, receiving device is laid, and then gathered data, and enters line number according to the data gathered According to reconstruct.When it is implemented, methods described can also comprise the following steps.
106:Excitaton source is laid in the shot position, receiver is laid in the receiving point position.
In the present embodiment, the excitaton source can be specifically collection big gun equal excitation equipment, and the reception implement body can To be the receiving devices such as wave detector.
107:Stochastical sampling is carried out by the excitaton source, the receiver, obtains lack sampling data.
In the present embodiment, can be swashed in order to obtain above-mentioned lack sampling data, during specific implementation by the excitaton source Hair, recycles the receiving device to be received, so as to obtain the lack sampling data in the space from subsurface reflective.
108:According to the lack sampling data, data reconstruction is carried out by compression sensing method.
In the present embodiment, above-mentioned data reconstruction can be specifically by carrying out being based on compressed sensing to lack sampling data Theoretical data reconstruction, and then the fully sampled geological data after being reconstructed.
In one embodiment, after the lack sampling data are obtained, first alias can also be converted into and is easy to filter out Low amplitude value noise, carry out denoising, obtain the higher lack sampling data of precision;The lack sampling data after processing are passed through again It is sparse to promote strategy to carry out data recovery;The complete geological data after data reconstruction can finally be obtained.
In the embodiment of the present application, by detecting grid chart corresponding to making, and determine that variance is equal according to detection grid chart It is worth ratio, judges whether target area shot position, receiving point position meet preset requirement using quadrat method.Therefore, solve It can not be put present in existing method and whether meet that to shot position, receiving point position randomness requires to carry out accurate judgement Technical problem, reached it is cost-effective, improve efficiency of construction technique effect;Also pass through the estimate of situation pair according to randomness Shot position, receiving point position carry out randomness adjustment, to cause shot position, receiving point position to possess certain randomness, And then improve the degree of accuracy of data reconstruction.
As can be seen from the above description, the detection of shot position, receiving point position that the application embodiment provides Method, quadrat method is determined by detecting grid chart corresponding to making, and according to detection grid chart, is sentenced using quadrat method Whether disconnected target area shot position, receiving point position meet preset requirement.Therefore, solves nothing present in existing method Whether method point meets that randomness requires to carry out the technical problem of accurate judgement to shot position, receiving point position, has reached saving Cost, the technique effect for improving efficiency of construction;Also by carrying out randomness to shot position, receiving point position according to judged result Adjustment, improves the degree of accuracy of data reconstruction.
Based on same inventive concept, the inspection of a kind of shot position, receiving point position is additionally provided in embodiment of the present invention Device is surveyed, as described in following embodiment.Solve the principle of problem and the detection of shot position, receiving point position due to device Method is similar, therefore the implementation of device may refer to the implementation of method, repeats part and repeats no more.It is used below, term " unit " or " module " can realize the combination of the software and/or hardware of predetermined function.Although described by following examples Device is preferably realized with software, but hardware, or software and hardware combination realization and may and be contemplated. Referring to Fig. 2, be the shot position of the application embodiment, receiving point position detection means a kind of composition structure chart, should Device can include:Acquisition module 201, drawing module 202, statistical module 203, the first determining module 204, the second determining module 205, the structure is specifically described below.
Acquisition module 201, it specifically can be used for obtaining shot position, receiving point position in target area.
It drawing module 202, specifically can be used for according to the shot position, the receiving point position, make detection grid Figure, wherein, the detection grid chart includes multiple unit grids.
Statistical module 203, it specifically can be used for counting the shot point number detected in grid chart in unit grids and connect Sink number.
First determining module 204, the shot point number and receiving point number that specifically can be used in the unit grids, It is determined that the quadrat method of detection grid chart.
Second determining module 205, it specifically can be used for the quadrat method according to the detection grid chart, confirm the inspection Whether shot position, receiving point position meet preset requirement corresponding to survey grid trrellis diagram, i.e., whether meet randomness requirement.
In one embodiment, in order to obtain corresponding detection grid chart, when it is implemented, above-mentioned drawing module 202 It can specifically include:
Range cells, it specifically can be used for, according to the shot position, the receiving point position, determining shooting distance, big gun line Away from.
Mapping unit, specifically the detection grid chart can be made according to the shooting distance, the perpendicular offset of shot line.
In the present embodiment, it is necessary to which explanation is, it is contemplated that processing accuracy takes into account treatment effeciency, the unit grids Length value is twice of the shooting distance, and the width value of the unit grids point is twice of the perpendicular offset of shot line.
In one embodiment for the shot point number and receiving point number in the unit grids, it is determined that detection The quadrat method of grid chart, first determining module can specifically include:
First determining unit, the shot point number and receiving point number that specifically can be used in the unit grids, really The points average of regular inspection survey grid trrellis diagram.
Second determining unit, it specifically can be used for the points average according to the detection grid chart, in the unit grids Shot point number and receiving point number, it is determined that detection grid chart points variance.
3rd determining unit, it specifically can be used for the points average according to the detection grid chart and the detection grid chart Points variance, determine it is described detection grid chart quadrat method.
In one embodiment, in order to determine the shot position under different situations, whether receiving point position meets presets It is required that whether meet randomness requirement, when it is implemented, second determining module can specifically include:
It 4th determining unit, specifically can be used in the case where two-dimensional observational system is laid, the side of grid chart will be detected Poor average ratio is more than 0.5 and shot position, the receiving point position detected corresponding to grid chart less than 1, is defined as described in satisfaction Preset requirement.That is, preset range, according to construction requirement, can be arranged to by situation about being laid for two-dimensional observational system:0.5 To 1.
It 5th determining unit, specifically can be used in the case where stereo observing system is laid, the side of grid chart will be detected Poor average ratio is more than 0.7 and shot position, the receiving point position detected corresponding to grid chart less than 1, is defined as described in satisfaction Preset requirement.That is, preset range, according to construction requirement, can be arranged to by situation about being laid for stereo observing system:0.7 To 1.
In one embodiment, it is corresponding in order to be carried out to the shot position for being unsatisfactory for preset requirement, receiving point position Adjustment, said apparatus can also include adjusting module.
Above-mentioned adjusting module, specifically it can be used for it is determined that shot position, receiving point position are unsatisfactory for the preset requirement In the case of, distance is lacked by controlling, to the shot position, the receiving point position, carries out randomness adjustment.So that The laying of shot point, receiving point has certain randomness, with meet stochastical sampling, imperfect recovery laying requirement.
In one embodiment, it is determined that layout scheme, i.e. shot position, receiving point position meet preset requirement, i.e., After meeting randomness requirement, data sampling and processing, data reconstruction can be carried out by module of constructing.Wherein, above-mentioned construction module It can specifically include:
Unit is laid, specifically can be used for laying excitaton source in the shot position, lay and connect in the receiving point position Receive device.
Sampling unit, it specifically can be used for carrying out stochastical sampling by the excitaton source, the receiver, obtain lack sampling Data.
Reconfiguration unit, it specifically can be used for, according to the lack sampling data, data reconstruction being carried out by compression sensing method.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment Divide mutually referring to what each embodiment stressed is the difference with other embodiment.It is real especially for system For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method Part explanation.
It should be noted that system, device, module or unit that above-mentioned embodiment illustrates, specifically can be by computer Chip or entity are realized, or are realized by the product with certain function.For convenience of description, in this manual, retouch It is divided into various units when stating apparatus above with function to describe respectively.Certainly, can be the function of each unit when implementing the application Realized in same or multiple softwares and/or hardware.
In addition, in this manual, adjective can be only used for an element or dynamic such as first and second Make to make a distinction with another element or action, without requiring or implying any actual this relation or order.Permit in environment Perhaps in the case of, be should not be interpreted as limited to reference to element or part or step (s) in only element, part or step one It is individual, and can be one or more of element, part or step etc..
As can be seen from the above description, the detection of shot position, receiving point position that the application embodiment provides Device, grid chart is detected corresponding to module making by charting, variance is determined according to detection grid chart by the first determining module Average ratio, judge whether target area shot position, receiving point position meet using quadrat method by the second determining module Preset requirement.Therefore, solving can not put present in existing method to whether shot position, receiving point position meet at random Property require carry out accurate judgement technical problem, reached it is cost-effective, improve efficiency of construction technique effect;Also pass through adjustment Module carries out randomness adjustment according to judged result to shot position, receiving point position, improves the degree of accuracy of data reconstruction.
Scene, the detection method/dress for shot position, receiving point position carried using the application is embodied at one Put, to shot position in three kinds of layout schemes not of the same race, receiving point position detected, judges whether to meet random distribution requirement (i.e. preset requirement or randomness requirement).Specific implementation process can perform according to following steps.
S1:According to regular blown-out shot scheme, random blown-out shot scheme 1, random blown-out shot scheme 2, it is corresponding that layout scheme is obtained respectively Shot position and receiving point point position.
S2:According to the shot position of different schemes and receiving point position, corresponding detection grid chart is generated:Detect grid chart 0th, grid chart 1, detection grid chart 2 are detected, wherein, each unit grids detected in grid chart are it is expected spatial sampling interval two (length value of i.e. described unit grids is twice of the shooting distance, and the width value of the unit grids point is the big gun line again Away from twice).
S3:Shot point number and receiving point number in each unit grids in each detection grid chart are counted respectively.
S4:Shot point number and receiving point number in each unit grids in each detection grid chart calculate each inspection The quadrat method of survey grid trrellis diagram.
S5:Stereo observing system is primarily directed in view of the construction and lays situation, and accordingly, preset range is set to 0.7 to 1.By the quadrat method of each detection grid chart compared with preset range, according to comparative result to each laying Scheme is accordingly evaluated.Can be refering to Fig. 3 detection that shot position, receiving point position are provided using the application embodiment Evaluation result schematic diagram under the different layout schemes that method, device obtains.Specifically, for regular blown-out shot scheme:Shot point position Put, receiving point position is unsatisfactory for preset requirement, that is, be unsatisfactory for randomness requirement, when it is implemented, follow-up data reconstruct effect Can be poor.Be all located at the mean variance of blown-out shot scheme 1, random blown-out shot scheme 2 in preset range, i.e., with blown-out shot scheme 1, with Shot position, receiving point position corresponding to machine blown-out shot scheme 2 all meet preset requirement, that is, meet randomness requirement, subsequently all Preferable data reconstruction effect can be obtained.Specific mean variance and the side of random blown-out shot scheme 2 relatively with blown-out shot scheme 1 again Poor average, it is known that:The mean variance of random blown-out shot scheme 2 relative to the mean variance of random blown-out shot scheme 1 closer to 1, because This, it is believed that the shot position, receiving point position corresponding to random blown-out shot scheme 2 are relative to corresponding to random blown-out shot scheme 1 Shot position, receiving point Randomness of position it is more preferable, i.e., using the data of random blown-out shot scheme 2 collection relative to using random empty The effect for the Data Data reconstruct that big gun scheme 1 gathers is more preferable.In summary, can be with excellent from above-mentioned three kinds of different layout schemes Random blown-out shot scheme 2 is selected to carry out random data sampling.
By above-mentioned example, the inspection for shot position, receiving point position provided using the application embodiment is demonstrated Surveying method, device can solve put present in existing method to whether shot position, receiving point position meet really Randomness requires the technical problem of progress accurate judgement, has reached technique effect that is cost-effective, improving efficiency of construction.
Although mentioning different embodiments in teachings herein, the application is not limited to be capable Situation described by industry standard or embodiment etc., some professional standards or the implementation described using self-defined mode or embodiment On the basis of embodiment amended slightly can also realize above-described embodiment it is identical, it is equivalent or it is close or deformation after it is anticipated that Implementation result.Using the embodiment of the data acquisition after these modifications or deformation, processing, output, judgment mode etc., still may be used To belong within the scope of the optional embodiment of the application.
Although this application provides the method operating procedure as described in embodiment or flow chart, based on conventional or noninvasive The means for the property made can include more or less operating procedures.The step of being enumerated in embodiment order is only numerous steps A kind of mode in execution sequence, does not represent unique execution sequence., can when device or client production in practice performs So that according to embodiment, either method shown in the drawings order is performed or parallel performed (such as at parallel processor or multithreading The environment of reason, even distributed data processing environment).Term " comprising ", "comprising" or its any other variant are intended to contain Lid nonexcludability includes, so that process, method, product or equipment including a series of elements not only will including those Element, but also the other element including being not expressly set out, or it is this process, method, product or equipment also to include Intrinsic key element.In the absence of more restrictions, be not precluded from the process including the key element, method, product or Other identical or equivalent elements in person's equipment also be present.
Device that above-described embodiment illustrates or module etc., it can specifically be realized by computer chip or entity, or by having There is the product of certain function to realize.For convenience of description, it is divided into various modules during description apparatus above with function to retouch respectively State.Certainly, the function of each module can be realized in same or multiple softwares and/or hardware when implementing the application, The module for realizing same function can be realized by the combination of multiple submodule etc..Device embodiment described above is only Schematically, for example, the division of the module, only a kind of division of logic function, can there is other draw when actually realizing The mode of dividing, such as multiple module or components can combine or be desirably integrated into another system, or some features can be ignored, Or do not perform.
It is also known in the art that in addition to realizing controller in a manner of pure computer readable program code, it is complete Entirely can by by method and step carry out programming in logic come controller with gate, switch, application specific integrated circuit, may be programmed The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind Hardware component, and what its inside was included is used to realize that the device of various functions can also to be considered as the structure in hardware component.Or Person even, not only can be able to will be the software module of implementation method but also can be hardware for realizing that the device of various functions is considered as Structure in part.
The application can be described in the general context of computer executable instructions, such as program Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type Part, data structure, class etc..The application can also be put into practice in a distributed computing environment, in these DCEs, By performing task by communication network and connected remote processing devices.In a distributed computing environment, program module can With in the local and remote computer-readable storage medium including storage device.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application On the part that is contributed in other words to prior art can be embodied in the form of software product, the computer software product It can be stored in storage medium, such as ROM/RAM, magnetic disc, CD, including some instructions are causing a computer equipment (can be personal computer, mobile terminal, server, either network equipment etc.) perform each embodiment of the application or implementation Method described in some parts of example.
Each embodiment in this specification is described by the way of progressive, same or analogous portion between each embodiment Divide mutually referring to what each embodiment stressed is the difference with other embodiment.The application can be used for crowd In more general or special purpose computing system environments or configuration.Such as:Personal computer, server computer, handheld device or Portable set, laptop device, multicomputer system, the system based on microprocessor, set top box, programmable electronics are set Standby, network PC, minicom, the DCE etc. of mainframe computer including any of the above system or equipment.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and Change is without departing from spirit herein, it is desirable to which appended embodiment includes these deformations and changed without departing from the application.

Claims (12)

1. a kind of shot position, the detection method of receiving point position, it is characterised in that including:
Obtain shot position, the receiving point position in target area;
According to the shot position, the receiving point position, detection grid chart is made, the detection grid chart includes multiple units Grid;
Count shot point number and receiving point number in unit grids described in the detection grid chart;
Shot point number and receiving point number in the unit grids, determine the quadrat method of the detection grid chart;
According to the quadrat method of the detection grid chart, shot position, receiving point position corresponding to the detection grid chart are confirmed Put and whether meet preset requirement.
2. according to the method for claim 1, it is characterised in that it is described according to the shot position, the receiving point position, Detection grid chart is made, including:
According to the shot position, the receiving point position, shooting distance, perpendicular offset of shot line are determined;
According to the shooting distance, the perpendicular offset of shot line, the detection grid chart is made.
3. according to the method for claim 2, it is characterised in that the length value of the unit grids is the two of the shooting distance Times, the width value of the unit grids point is twice of the perpendicular offset of shot line.
4. according to the method for claim 1, it is characterised in that the shot point number in the unit grids and connect Sink number, it is determined that the quadrat method of detection grid chart, including:
Shot point number and receiving point number in the unit grids, it is determined that the points average of detection grid chart;
According to the shot point number and receiving point number in the points average of the detection grid chart, the unit grids, it is determined that inspection The points variance of survey grid trrellis diagram;
According to the points average of the detection grid chart and the points variance of the detection grid chart, the detection grid chart is determined Quadrat method.
5. according to the method for claim 4, it is characterised in that the shot point number in the unit grids and connect Sink number, it is determined that the points average of detection grid chart, including:
According to below equation, the points average for detecting grid chart is determined:
<mrow> <mi>M</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </msubsup> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mi>M</mi> </mfrac> </mrow>
In above formula, Mean be detection grid chart points average, XiTo number shot point number and the reception in the unit grids for being i Point number sum, M are the number of unit grids in detection grid chart, and i is the numbering of unit grid, and span is 1 to M.
6. according to the method for claim 5, it is characterised in that described according to the points average for detecting grid chart, institute The shot point number and receiving point number in unit grids are stated, it is determined that the points variance of detection grid chart, including:
The points variance of the detection grid chart is determined according to below equation:
<mrow> <mi>V</mi> <mi>a</mi> <mi>r</mi> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </msubsup> <msup> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>M</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <mi>M</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </mrow>
In above formula, Mean is the points average of detection grid chart, and Var is the points variance of detection grid chart.
7. according to the method for claim 6, it is characterised in that according to the points average of the detection grid chart and the inspection The points variance of survey grid trrellis diagram, the quadrat method of the detection grid chart is determined, including:
The quadrat method of the detection grid chart is determined according to below equation:
<mrow> <mi>V</mi> <mi>M</mi> <mi>R</mi> <mo>=</mo> <mfrac> <mrow> <mi>V</mi> <mi>a</mi> <mi>r</mi> </mrow> <mrow> <mi>M</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> </mrow> </mfrac> </mrow>
In above formula, Mean is the points average of detection grid chart, and Var is the points variance of detection grid chart, and VMR is detection grid The quadrat method of figure.
8. according to the method for claim 1, it is characterised in that according to the quadrat method of the detection grid chart, confirm Whether shot position, receiving point position meet preset requirement corresponding to the detection grid chart, including:
In the case where two-dimensional observational system is laid, the quadrat method for detecting grid chart is more than 0.5 and the detection net less than 1 Shot position, receiving point position corresponding to trrellis diagram, are defined as meeting the preset requirement;
In the case where stereo observing system is laid, the quadrat method for detecting grid chart is more than 0.7 and the detection net less than 1 Shot position, receiving point position corresponding to trrellis diagram, are defined as meeting the preset requirement.
9. according to the method for claim 1, it is characterised in that confirm it is described detection grid chart corresponding to shot position, After receiving point position is unsatisfactory for the preset requirement, methods described also includes:
Distance is lacked by controlling, to the shot position, the receiving point position, carries out randomness adjustment.
10. according to the method for claim 1, it is characterised in that confirm it is described detection grid chart corresponding to shot position, After receiving point position meets the preset requirement, methods described also includes:
Excitaton source is laid in the shot position, receiver is laid in the receiving point position;
Stochastical sampling is carried out by the excitaton source, the receiver, obtains lack sampling data;
According to the lack sampling data, data reconstruction is carried out by compression sensing method.
11. a kind of shot position, the detection means of receiving point position, it is characterised in that including:
Acquisition module, for obtaining the shot position in target area, receiving point position;
Drawing module, for according to the shot position, the receiving point position, making detection grid chart, the detection grid Figure includes multiple unit grids;
Statistical module, for counting shot point number and receiving point number in the detection grid chart in unit grids;
First determining module, for the shot point number and receiving point number in the unit grids, determine the detection net The quadrat method of trrellis diagram;
Second determining module, for the quadrat method according to the detection grid chart, confirm corresponding to the detection grid chart Whether shot position, receiving point position meet preset requirement.
12. device according to claim 11, it is characterised in that first determining module includes:
First determining unit, for the shot point number and receiving point number in the unit grids, it is determined that detection grid chart Points average;
Second determining unit, for the shot point number in the detection points average of grid chart, the unit grids and Receiving point number, it is determined that the points variance of detection grid chart;
3rd determining unit, for according to it is described detection grid chart points average and it is described detect grid chart points variance, Determine the quadrat method of the detection grid chart.
CN201710594250.6A 2017-07-20 2017-07-20 Shot position, the detection method of receiving point position and device Pending CN107561574A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710594250.6A CN107561574A (en) 2017-07-20 2017-07-20 Shot position, the detection method of receiving point position and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710594250.6A CN107561574A (en) 2017-07-20 2017-07-20 Shot position, the detection method of receiving point position and device

Publications (1)

Publication Number Publication Date
CN107561574A true CN107561574A (en) 2018-01-09

Family

ID=60974627

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710594250.6A Pending CN107561574A (en) 2017-07-20 2017-07-20 Shot position, the detection method of receiving point position and device

Country Status (1)

Country Link
CN (1) CN107561574A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108732619A (en) * 2018-04-27 2018-11-02 中国科学院地质与地球物理研究所 A kind of submarine geophysics collecting method
CN111474574A (en) * 2019-01-23 2020-07-31 中国石油天然气集团有限公司 Compressed sensing-based generation method and device for seismic acquisition observation system
CN111965694A (en) * 2020-08-03 2020-11-20 中国石油天然气集团有限公司 Method and device for determining position of seismic physical point and seismic observation system
CN112099091A (en) * 2019-06-17 2020-12-18 中国石油天然气集团有限公司 Seismic exploration point data static correction method and device
CN112711066A (en) * 2019-10-25 2021-04-27 中国石油化工股份有限公司 Three-dimensional seismic exploration shot point arrangement uniformity evaluation method and device
CN112782781A (en) * 2019-11-07 2021-05-11 中国石油天然气集团有限公司 Evaluation method and device for earthquake observation system
CN112946736A (en) * 2019-11-26 2021-06-11 中国石油天然气集团有限公司 Three-dimensional observation system reconstruction method and system

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104422954A (en) * 2013-08-20 2015-03-18 中国石油化工股份有限公司 Sparse sampling method in seismic data regularization

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104422954A (en) * 2013-08-20 2015-03-18 中国石油化工股份有限公司 Sparse sampling method in seismic data regularization

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
戴小华 等: ""GIS支持下的种群分布格局分析"", 《中山大学学报(自然科学版)》 *
杨敬磊 等: ""基于数据重构的观测系统评价方法"", 《2017年物探技术研讨会》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108732619A (en) * 2018-04-27 2018-11-02 中国科学院地质与地球物理研究所 A kind of submarine geophysics collecting method
CN108732619B (en) * 2018-04-27 2021-02-26 中国科学院地质与地球物理研究所 Seabed geophysical data acquisition method
CN111474574A (en) * 2019-01-23 2020-07-31 中国石油天然气集团有限公司 Compressed sensing-based generation method and device for seismic acquisition observation system
CN111474574B (en) * 2019-01-23 2023-04-25 中国石油天然气集团有限公司 Method and device for generating seismic acquisition observation system based on compressed sensing
CN112099091A (en) * 2019-06-17 2020-12-18 中国石油天然气集团有限公司 Seismic exploration point data static correction method and device
CN112099091B (en) * 2019-06-17 2024-03-01 中国石油天然气集团有限公司 Static correction method and device for seismic exploration point data
CN112711066A (en) * 2019-10-25 2021-04-27 中国石油化工股份有限公司 Three-dimensional seismic exploration shot point arrangement uniformity evaluation method and device
CN112711066B (en) * 2019-10-25 2024-02-20 中国石油化工股份有限公司 Three-dimensional seismic exploration shot point layout uniformity evaluation method and device
CN112782781A (en) * 2019-11-07 2021-05-11 中国石油天然气集团有限公司 Evaluation method and device for earthquake observation system
CN112946736A (en) * 2019-11-26 2021-06-11 中国石油天然气集团有限公司 Three-dimensional observation system reconstruction method and system
CN111965694A (en) * 2020-08-03 2020-11-20 中国石油天然气集团有限公司 Method and device for determining position of seismic physical point and seismic observation system
CN111965694B (en) * 2020-08-03 2024-01-30 中国石油天然气集团有限公司 Method and device for determining position of physical point of earthquake and earthquake observation system

Similar Documents

Publication Publication Date Title
CN107561574A (en) Shot position, the detection method of receiving point position and device
CN106855640B (en) A kind of submarine cable seismic data first arrival calculation method and device
CN106199703B (en) A kind of automatic positioning of microseism focus and Reliability Synthesis evaluation method
CN103605161B (en) Analytical method and the device of the CMP bin attribute of earthquake three-dimensional observation system
CN105549083B (en) A kind of method and apparatus for determining seismic observation system repeatability and integrally measuring
CN105404726A (en) Method and device for reverting inter-well connectivity through Gaussian distribution-based capacitance model
CN104462695A (en) Weak signal detection method based on double-coupling Duffing vibrators and scale varying
CN106501841B (en) A kind of optimum design method and device of two dimension wode line seismic observation system
CN103777248A (en) TEM one-dimensional forward modeling method applicable to irregular transmitting loop
CN107576985B (en) A kind of method and apparatus of seismic inversion
CN103592684B (en) A kind of massive seismic data compression method keeping space attribute information and device
CN104391325B (en) Discontinuous heterogeneous geologic body detection method and device
CN107870361B (en) A kind of earthquake diving Wave chromatography imaging method, device and terminal device
CN106324671A (en) Method and device of examining shot point offset
CN106054252A (en) Pre-stack time migration method and device
CN105204067B (en) Seismic horizon method for tracing and device
CN106094020A (en) A kind of seismic inversion method and device
CN107807409B (en) The determination method and apparatus of density of earth formations and resistivity relation
CN103901475B (en) One attribute isogram method for drafting and device
CN107991705B (en) Log curve correction method and apparatus based on Two-dimensional Statistical feature
CN107703541A (en) The determination method and apparatus of stratigraphic dip
CN104865601B (en) The method for determining the type of the earthquake record gathered in seismic prospecting
Emery et al. Models for support and information effects: a comparative study
CN107085236B (en) The determination method and apparatus of maximum offset
CN107065008B (en) A kind of method and device of determining observation system degree of covering

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: 20180109