CN107728212B - Transient State Rayleigh Wave detects bearing calibration and device - Google Patents

Transient State Rayleigh Wave detects bearing calibration and device Download PDF

Info

Publication number
CN107728212B
CN107728212B CN201710897962.5A CN201710897962A CN107728212B CN 107728212 B CN107728212 B CN 107728212B CN 201710897962 A CN201710897962 A CN 201710897962A CN 107728212 B CN107728212 B CN 107728212B
Authority
CN
China
Prior art keywords
measuring point
coordinate
abscissa
moving process
geophone group
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.)
Active
Application number
CN201710897962.5A
Other languages
Chinese (zh)
Other versions
CN107728212A (en
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.)
FUZHOU HUAHONG INTELLIGENT TECHNOLOGY Co.,Ltd.
Original Assignee
Anhui University of Science and Technology
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 Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN201710897962.5A priority Critical patent/CN107728212B/en
Publication of CN107728212A publication Critical patent/CN107728212A/en
Application granted granted Critical
Publication of CN107728212B publication Critical patent/CN107728212B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/52Move-out correction

Abstract

The present invention provides a kind of detection bearing calibration of Transient State Rayleigh Wave and device.This method comprises: moving geophone group according to pre-determined distance in the abscissa direction of datum level, the base position coordinate of multiple wave detectors in geophone group in each moving process is obtained;According in each moving process in geophone group each wave detector base position coordinate, determine that the base position coordinate of multiple measuring points respectively corresponds to the center position coordinates of geophone group in each moving process respectively;Carry out Transient State Rayleigh Wave detection processing respectively at the base position coordinate of each measuring point, determine each measuring point to bedrock surface depth measurement reference coordinate;According in each moving process in geophone group each wave detector actual position coordinate, obtain the corresponding matched curve of each measuring point respectively;According to the corresponding matched curve of each measuring point, processing is corrected to the depth measurement reference coordinate of each measuring point respectively, obtains the correction of soundings coordinate of each measuring point.The present invention realizes the accurate detection to subterrane form.

Description

Transient State Rayleigh Wave detects bearing calibration and device
Technical field
The present invention relates to geotechnical engineering investigation technical field more particularly to a kind of detection bearing calibration of Transient State Rayleigh Wave and dresses It sets.
Background technique
Rayleigh waves are a kind of surface waves propagated along some intracorporal interface of Rock And Soil Free Surface or ground, are along corresponding boundary A kind of secondary wave that the nonuniform plane P wave and nonuniform plane SV wave both bulk waves that face is propagated are formed by stacking, Wave energy It is concentrated mainly near the interface of wave propagation, caused particle vibration track is interface normal and direction of wave travel near interface Planar inverse into ellipse, transverse direction is consistent with interface normal orientation.In layering uniform dielectric or uneven Jie The Rayleigh waves propagated in matter have Dispersion, that is, Rayleigh waves are made of a series of single-frequency fluctuation of different frequency ingredients, different The fluctuation of frequency content has different spread speeds (phase velocity);In addition, the Effective depth penetration of single-frequency Rayleigh waves (has foot The influence of fluctuations depth of enough vibrational energies) it is directly proportional to wavelength, that is, and high frequency (shortwave) the ingredient reflection superficial part ground of Rayleigh waves is situated between Matter property, low frequency (long wave) ingredient reflect deep rock soil medium property.
Transient State Rayleigh Wave detection is substantially the Dispersion propagated according to Rayleigh waves, is generated using man-made explosion excitation more The Rayleigh waves of kind frequency content, search out velocity of wave with the variation relation of frequency, to finally determine the Rayleigh waves of subterrane Speed with site coordinate variation relation, the problems such as earthquake engineering to solve shallow engineering geology and subterrane.
Existing Transient State Rayleigh Wave detection method usually requires to see the landform for exploring work area as level of approximation, Cai Nengfen The form of subterrane is precipitated.However, the accuracy of the fluctuating meeting severe jamming detection of actual landform.Therefore, one kind is now needed Corrected Transient State Rayleigh Wave detection method can be influenced on hypsography.
Summary of the invention
The present invention provides a kind of detection bearing calibration of Transient State Rayleigh Wave and device, to solve existing Transient State Rayleigh Wave detection Bearing calibration causes the problem of detection result inaccuracy due to the fluctuating of landform.
The present invention provides a kind of Transient State Rayleigh Wave detection bearing calibration, comprising:
Geophone group is moved according to pre-determined distance in the abscissa direction of datum level, is obtained in each moving process The base position coordinate of multiple wave detectors in geophone group, wherein the abscissa of the datum level is any horizontal linear, it is described The ordinate of datum level is any downward straight line perpendicular to the abscissa of the datum level;
According in each moving process in geophone group each wave detector base position coordinate, determine multiple measuring points respectively Base position coordinate respectively corresponds to the center position coordinates of geophone group in each moving process;
It carries out Transient State Rayleigh Wave detection processing respectively at the base position coordinate of each measuring point, determines each measuring point to bedrock surface Depth measurement reference coordinate;
According in each moving process in geophone group each wave detector actual position coordinate, obtain each measuring point respectively respectively Corresponding matched curve;
According to the corresponding matched curve of each measuring point, processing is corrected to the depth measurement reference coordinate of each measuring point respectively, Obtain the correction of soundings coordinate of each measuring point.
Optionally, the base position coordinate for obtaining multiple wave detectors in geophone group in each moving process, comprising:
According to the origin of the datum level, the position coordinates of the focal point, the focal point to being moved through every time With a distance from the nearest wave detector of origin in journey in geophone group from the datum level and in the geophone group each wave detector it Between distance, determine the base position coordinate of each wave detector in geophone group in each moving process.
Optionally, described to carry out Transient State Rayleigh Wave spy respectively at the base position coordinate of each measuring point for any measuring point Survey processing, determine each measuring point to bedrock surface depth measurement reference coordinate, comprising:
Determine the horizontal seat in the base position coordinate of the abscissa and the measuring point in the depth measurement reference coordinate of the measuring point It marks equal;
Acquire the vibration signal of each wave detector in the corresponding geophone group of the measuring point;
Frequency dispersion processing is carried out to the vibration signal of each wave detector, obtains the spectrum curve of the measuring point;
Inversion procedure is carried out to the spectrum curve of the measuring point, obtains the vertical seat in the depth measurement reference coordinate of the measuring point Mark.
Optionally, described according to the corresponding matched curve of each measuring point for any measuring point, respectively to the survey of each measuring point Deep reference coordinate is corrected processing, obtains the correction of soundings coordinate of each measuring point, comprising:
Determine that the abscissa in the correction of soundings coordinate of the measuring point is the horizontal seat in the depth measurement reference coordinate of the measuring point Mark;
According to the matched curve of abscissa and the measuring point in the depth measurement reference coordinate of the measuring point, the measuring point is obtained Correction of soundings coordinate in ordinate.
Optionally, in the geophone group according to each moving process each wave detector actual position coordinate, After obtaining the corresponding matched curve of each measuring point respectively, the method also includes:
When the corresponding matched curve of each measuring point be straight line when, according to the corresponding matched curve of each measuring point with it is described The angle of the abscissa direction of datum level, determine respectively each measuring point to the batholith face correction of soundings direction.
The present invention also provides a kind of Transient State Rayleigh Waves to detect means for correcting, comprising:
Module is obtained, for moving geophone group according to pre-determined distance in the abscissa direction of datum level, is obtained In each moving process in geophone group multiple wave detectors base position coordinate, wherein the abscissa of the datum level is any Horizontal linear, the ordinate of the datum level are any downward straight line perpendicular to the abscissa of the datum level;
Determining module, for according in each moving process in geophone group each wave detector base position coordinate, respectively Determine that the base position coordinate of multiple measuring points respectively corresponds to the center position coordinates of geophone group in each moving process;
Processing module is determined for carrying out Transient State Rayleigh Wave detection processing respectively at the base position coordinate of each measuring point Depth measurement reference coordinate of each measuring point to bedrock surface;
The acquisition module is also used to according to the physical location of each wave detector is sat in geophone group in each moving process Mark, obtains the corresponding matched curve of each measuring point respectively;
The processing module is also used to according to the corresponding matched curve of each measuring point, respectively to the depth measurement base of each measuring point Quasi coordinates is corrected processing, obtains the correction of soundings coordinate of each measuring point.
Optionally, the module that obtains is specifically used for position according to the origin of the datum level, the focal point With a distance from the nearest wave detector of the origin of coordinate, the focal point into each moving process in geophone group from the datum level The distance between wave detector each in the geophone group, determines the benchmark of each wave detector in geophone group in each moving process Position coordinates.
Optionally, for any measuring point, the processing module is specifically used in the depth measurement reference coordinate for determining the measuring point Abscissa it is equal with the abscissa in the base position coordinate of the measuring point;
Acquire the vibration signal of each wave detector in the corresponding geophone group of the measuring point;
Frequency dispersion processing is carried out to the vibration signal of each wave detector, obtains the spectrum curve of the measuring point;
Inversion procedure is carried out to the spectrum curve of the measuring point, obtains the vertical seat in the depth measurement reference coordinate of the measuring point Mark.
Optionally, for any measuring point, the processing module is specifically also used to determine the correction of soundings coordinate of the measuring point In abscissa be the measuring point depth measurement reference coordinate in abscissa;
According to the matched curve of abscissa and the measuring point in the depth measurement reference coordinate of the measuring point, the measuring point is obtained Correction of soundings coordinate in ordinate.
Optionally, the processing module is used for when the corresponding matched curve of each measuring point is straight line, according to each measuring point The angle of the abscissa direction of corresponding matched curve and the datum level, determines each measuring point to the batholith respectively The correction of soundings direction in face.
Transient State Rayleigh Wave detection bearing calibration provided by the invention and device, pass through the abscissa direction along datum level According to the mobile geophone group of pre-determined distance, the base position for obtaining multiple wave detectors in geophone group in each moving process is sat Mark, wherein the abscissa of datum level is any horizontal linear, and the ordinate of datum level is times perpendicular to the abscissa of datum level One downward straight line.Corresponding inspection is determined further according to the base position coordinate of multiple wave detectors in geophone group in each moving process The center of wave device group obtains the base position coordinate of each measuring point.Then, distinguish at the base position coordinate of each measuring point Transient State Rayleigh Wave detection processing is carried out, so that it is determined that depth measurement reference coordinate of each measuring point to bedrock surface, further according to being moved through every time In journey in geophone group each wave detector actual position coordinate, obtain the corresponding matched curve of each measuring point respectively.Finally, logical The corresponding matched curve of each measuring point is crossed, processing is corrected to the reference coordinate of each measuring point respectively, obtains the survey of each measuring point Deep calibration coordinate.The present invention can fully consider the influence of hypsography, solve existing Transient State Rayleigh Wave detection correction side Method is since the fluctuating of landform causes the problem of detection result inaccuracy, to conveniently realize to subterrane form Detection, effectively increases the accuracy of detection.
Detailed description of the invention
Fig. 1 is the flow chart one that Transient State Rayleigh Wave provided by the invention detects bearing calibration;
Fig. 2 is the schematic diagram of a scenario one that Transient State Rayleigh Wave provided by the invention detects bearing calibration;
Fig. 3 is the flowchart 2 that Transient State Rayleigh Wave provided by the invention detects bearing calibration;
Fig. 4 is the schematic diagram of the dispersion curve of any measuring point in the present invention;
Fig. 5 is the schematic diagram of the inversion result of any measuring point in the present invention;
Fig. 6 is the schematic diagram of a scenario two that Transient State Rayleigh Wave provided by the invention detects bearing calibration;
Fig. 7 is the schematic diagram of the calibration curve of any measuring point in the present invention;
Fig. 8 is the flow chart 3 that Transient State Rayleigh Wave provided by the invention detects bearing calibration;
Fig. 9 is the schematic diagram of a scenario three that Transient State Rayleigh Wave provided by the invention detects bearing calibration;
Figure 10 is the structural schematic diagram that Transient State Rayleigh Wave provided by the invention detects means for correcting.
Specific embodiment
Fig. 1 is the flow chart one that Transient State Rayleigh Wave provided by the invention detects bearing calibration, and Fig. 2 is wink provided by the invention The schematic diagram of a scenario one of state Rayleigh wave survey bearing calibration.As shown in Figure 1, the Transient State Rayleigh Wave of the present embodiment detects bearing calibration May include:
S101, according to pre-determined distance geophone group is moved in the abscissa direction of datum level, obtain every time mobile In the process in geophone group multiple wave detectors base position coordinate, wherein the abscissa of datum level be any horizontal linear, base The ordinate in quasi- face is any downward straight line perpendicular to the abscissa of datum level.
Specifically, due to being needed when detecting using Transient State Rayleigh Wave by ground to be approximately horizontal plane, the present embodiment It can be taken up an official post abscissa of the horizontal linear as datum level by selecting ground, and the abscissa perpendicular to datum level also may be selected Ordinate of any downward straight line as datum level, i.e. the vertical seat perpendicular to ground any straight line directed downwardly as datum level Mark to complete the acquisition process of datum level, and can satisfy the requirement of Transient State Rayleigh Wave detection.
Further, in this embodiment geophone group includes multiple wave detectors, number and type of the present embodiment to wave detector Number without limitation.In order to meet the industry regulation of Transient State Rayleigh Wave detection, the number of wave detector is usually taken in the present embodiment 12 or more.Wherein for ease of description, as shown in Fig. 2, 12 wave detectors can be selected in the present embodiment as a geophone group. And in the present embodiment can the multiple wave detectors of random alignment, the spacing between multiple wave detectors can be identical, can also be different, the present embodiment Also without limitation to this, need to only meet known to the distance between multiple wave detectors.
Further, in this embodiment the moving direction of geophone group can be set as to the abscissa place side of datum level To, and repeatedly the pre-determined distance of mobile geophone group can be equidistant, can also be Unequal distance, the present embodiment to move geophone group away from From size without limitation, the pre-determined distance that need to only meet mobile geophone group is known.
Further, be by the spacing between each wave detector in this present embodiment it is known, datum level be also it is determining, Therefore, the present embodiment can be used various ways and determine that the base position of multiple wave detectors in geophone group in each moving process is sat Mark, the present embodiment to the implementation of the base position coordinate for obtaining multiple wave detectors without limitation.
Optionally, it is examined according to the origin of datum level, the position coordinates of focal point, focal point into each moving process The distance between each wave detector with a distance from the nearest wave detector of origin in wave device group from datum level and in geophone group, determines every In secondary moving process in geophone group multiple wave detectors base position coordinate.
Specifically, the present embodiment is needed when carrying out Transient State Rayleigh Wave detection in focal point progress jarring, the present embodiment The position of random epicentre point can be chosen, the present embodiment without limitation to the position of focal point, sit by the position that need to only meet focal point With a distance from the wave detector that mark is known and origin into each moving process in geophone group from datum level is nearest it is also known that. Again due to datum level is determining and geophone group in the distance between each wave detector it is known that therefore, for any secondary moving process Speech, the present embodiment can be according to the origin of datum level, the position coordinates of focal point, focal point detection into the secondary moving process The distance between each wave detector with a distance from the nearest wave detector of origin in device group from datum level and in geophone group, can determine In each moving process in geophone group each wave detector base position coordinate.
S102, according in each moving process in geophone group each wave detector base position coordinate, determine respectively multiple The base position coordinate of measuring point respectively corresponds to the center position coordinates of geophone group in each moving process.
Specifically, as shown in Fig. 2, the present embodiment can choose the center of geophone group in a moving process in office Base position as measuring point.It is due to the spacing of each wave detector in the origin of datum level, geophone group and mobile every time Distance be it is known, therefore, the present embodiment can determine the base position coordinate of corresponding measuring point in each moving process.
S103, it carries out Transient State Rayleigh Wave detection processing respectively at the base position coordinate of each measuring point, determines that each measuring point arrives The depth measurement reference coordinate of bedrock surface.
Specifically, during each mobile geophone group, before the present embodiment is constant in the position for keeping focal point Put, Transient State Rayleigh Wave detection processing carried out to current corresponding measuring point, obtain each measuring point to bedrock surface depth measurement reference coordinate. The present embodiment can concrete mode to the depth measurement reference coordinate of each measuring point of determination to bedrock surface without limitation.
In conjunction with Fig. 3, for any measuring point, the specific implementation of the depth measurement reference coordinate to the determining measuring point to bedrock surface It is described in detail.Fig. 3 is the flowchart 2 that Transient State Rayleigh Wave provided by the invention detects bearing calibration, and Fig. 4 is in the present invention The schematic diagram of the dispersion curve of any measuring point, Fig. 5 are the schematic diagram of the inversion result of any measuring point in the present invention, and Fig. 6 is this hair The schematic diagram of a scenario two of the Transient State Rayleigh Wave detection bearing calibration of bright offer.As shown in figure 3, the Transient State Rayleigh Wave of the present embodiment is visited Survey bearing calibration further include:
S1031, abscissa in the base position coordinate of the abscissa and measuring point in the depth measurement reference coordinate of measuring point is determined It is equal.
Specifically, due to carrying out Transient State Rayleigh Wave detection processing to measuring point when the base position for keeping measuring point is constant, Therefore, the present embodiment can be using the abscissa in the base position coordinate of measuring point as the horizontal seat in the depth measurement reference coordinate of measuring point Mark.
The vibration signal of each wave detector in S1032, the corresponding geophone group of acquisition measuring point.
S1033, frequency dispersion processing is carried out to the vibration signal of each wave detector, obtains the spectrum curve of measuring point.
S1034, inversion procedure is carried out to the spectrum curve of measuring point, obtains the ordinate in the depth measurement reference coordinate of measuring point.
Specifically, the present embodiment can for focal point carry out jarring, the measuring point geophone group center position, by examining Each wave detector in wave device group acquires vibration signal, then carries out frequency dispersion processing to the collected vibration signal of each wave detector again, Obtain the spectrum curve of measuring point.
As shown in figure 4, specifically can be according to Rayleigh waves depth of exploration formula in the present embodimentWherein H is depth, vRFor speed, f is frequency, by determining the model of depth and speed, can obtain the spectrum curve of the measuring point.
Further, the present embodiment carries out inversion procedure to the spectrum curve of the measuring point, can obtain measuring point to bedrock surface Depth size, i.e. ordinate in the depth measurement reference coordinate of the measuring point indulges as shown in figure 5, wherein abscissa is shear wave velocity Coordinate is depth, and the corresponding depth of velocity jump point is exactly using datum level as the distance of interface measuring point to bedrock surface, the i.e. survey Ordinate in the depth measurement reference coordinate of point.
Further, in this embodiment aforesaid way can be used for any point obtains corresponding depth measurement benchmark Ordinate in coordinate, as shown in fig. 6, wherein multiple measuring points are maintained on datum level in the present embodiment, by each measuring point point Not carry out Transient State Rayleigh Wave detection processing can determine each measuring point to bedrock surface depth measurement reference coordinate, in this way, can be by each The size of ordinate value of depth measurement reference coordinate of measuring point to bedrock surface determines that measuring point, and then can to the depth size of bedrock surface Understand the concrete condition of basement rock form.
Fig. 7 is the schematic diagram of the calibration curve of any measuring point in the present invention.Due to Transient State Rayleigh Wave detection processing need by The interface of measuring point is similar to horizontal plane, therefore, this implementation can using datum level as interface obtain each measuring point to bedrock surface depth away from From, but actual landform but more or less will appear fluctuating.In order to accurately detect the actual form of basement rock, the present embodiment can be tied Fig. 7 is closed the specific implementation process of S104 is described in detail.
S104, according in each moving process in geophone group each wave detector actual position coordinate, obtain each survey respectively The corresponding matched curve of point.
Specifically, in the present embodiment during multiple mobile geophone group, the modes such as mapping measuring tool can be passed through Measure the actual position coordinate of each wave detector.The present embodiment does not do the implementation for the actual position coordinate for obtaining each wave detector It limits.For ease of description, it as shown in fig. 7, the present embodiment can choose any measuring point in multiple measuring points, first passes through and measures the survey The actual position coordinate of each wave detector in geophone group in the corresponding moving process of point, then use the method for fitting by each inspection The corresponding actual position coordinate of wave device is processed into a matched curve.Wherein, matched curve can become straight by linear fit Line, can also be fitted by other fit approach becomes bending line, and the present embodiment does not limit this.Multiple measuring points in the present embodiment Aforesaid way all can be used, obtain the corresponding matched curve of each measuring point, be not repeated herein.
S105, according to the corresponding matched curve of each measuring point, the depth measurement reference coordinate of each measuring point is corrected respectively Processing, obtains the correction of soundings coordinate of each measuring point.
Specifically, after determining the corresponding matched curve of each measuring point, for any measuring point, the present embodiment can pass through Corresponding matched curve is corrected processing to the depth measurement reference coordinate of the measuring point, can fully consider the shadow that actual landform rises and falls It rings, accurately detects the actual conditions of basement rock form.
In conjunction with Fig. 8, for any measuring point, processing is corrected to the reference coordinate of the measuring point, obtains the depth measurement of the measuring point The specific implementation of calibration coordinate is described in detail.Fig. 8 is that Transient State Rayleigh Wave provided by the invention detects bearing calibration Flow chart 3, Fig. 9 are the schematic diagram of a scenario three that Transient State Rayleigh Wave provided by the invention detects bearing calibration.As shown in figure 8, this reality Apply the Transient State Rayleigh Wave detection bearing calibration of example further include:
S1051, determine the abscissa in the correction of soundings coordinate of measuring point for the abscissa in the depth measurement reference coordinate of measuring point.
S1052, according to the matched curve of abscissa and measuring point in the depth measurement reference coordinate of measuring point, obtain the depth measurement of measuring point Ordinate in calibration coordinate.
Specifically, usually using the abscissa in the depth measurement reference coordinate of the measuring point as the depth measurement of the measuring point in the present embodiment Abscissa in calibration coordinate determines that the correction of soundings of the measuring point is sat according to the matched curve of the abscissa value and the measuring point jointly Ordinate in mark.For example, when the corresponding matched curve of the measuring point is Z=aX+b, it can be by the depth measurement reference coordinate of the measuring point In abscissa be brought into formula, the correction of soundings coordinate (X, aX+b) of the measuring point can be obtained.In this way, as shown in figure 9, passing through The fluctuating of landform has been depicted in the variation of the physical location of each geophone station, and the correction of soundings coordinate of measuring point is just it can be considered that practically The case where shape, the form of basement rock is accurately detected, correction course is quick, effective, it is easy to accomplish.
Transient State Rayleigh Wave provided in this embodiment detects bearing calibration, by the abscissa direction along datum level according to Pre-determined distance moves geophone group, obtains the base position coordinate of multiple wave detectors in geophone group in each moving process, The abscissa of middle datum level is any horizontal linear, and the ordinate of datum level is any downward of the abscissa perpendicular to datum level Straight line.Corresponding geophone group is determined further according to the base position coordinate of multiple wave detectors in geophone group in each moving process Center obtain the base position coordinate of each measuring point.Then, wink is carried out respectively at the base position coordinate of each measuring point State Rayleigh wave survey processing, so that it is determined that depth measurement reference coordinate of each measuring point to bedrock surface, further according to being examined in each moving process The actual position coordinate of each wave detector in wave device group, obtains the corresponding matched curve of each measuring point respectively.Finally, passing through each survey The corresponding matched curve of point, is corrected processing to the reference coordinate of each measuring point respectively, obtains the correction of soundings of each measuring point Coordinate.The present embodiment can fully consider the influence of hypsography, solve the detection bearing calibration of existing Transient State Rayleigh Wave by The problem of detection result inaccuracy is caused in the fluctuating of landform, to conveniently realize the spy to subterrane form It surveys, effectively increases the accuracy of detection.
After above-mentioned execution S104 process, continuing with Fig. 8, to measuring point each in the present embodiment from datum level to batholith face The determination process of detection direction be described in detail.Depth measurement of the matched curve in addition to can determine each measuring point in the present embodiment It, can also be bent according to the corresponding fitting of each measuring point when the corresponding matched curve of each measuring point is straight line outside calibration coordinate The angle of the abscissa direction of line and datum level, determine respectively each measuring point to batholith face correction of soundings direction.
Specifically, in existing Transient State Rayleigh Wave detection method its detection direction for the ordinate of depth measurement reference coordinate institute In direction, and it is directed to any measuring point, when the corresponding matched curve of each measuring point for using the mode of linear fit to obtain is straight It is that the corresponding fitting of the measuring point is bent that the measuring point, which deviates the angle of the direction of the ordinate of reference coordinate, when line, in the present embodiment The angle of the abscissa direction of line and datum level, and the depth measurement direction in the measuring point to batholith face can be identified as measuring point deviation The corresponding direction of the angle of the direction of the ordinate of depth measurement reference coordinate.When fitting a straight line non-horizontal surface, the survey is corrected The depth measurement direction in point to batholith face can misalignment angle α, α=arctana.If landform is raised, depth measurement of the measuring point to batholith face Direction can deviate to moving direction;If topographic depression, the depth measurement direction of the measuring point to batholith face can be to the reversed of moving direction Deviate.
Figure 10 is the structural schematic diagram that Transient State Rayleigh Wave provided by the invention detects means for correcting, as shown in Figure 10, this reality Applying a Transient State Rayleigh Wave detection means for correcting includes:
Module 10 is obtained to obtain for moving geophone group according to pre-determined distance in the abscissa direction of datum level The base position coordinate of multiple wave detectors in geophone group in each moving process is taken, wherein the abscissa of the datum level is to appoint One horizontal linear, the ordinate of the datum level are any downward straight line perpendicular to the abscissa of the datum level;
Determining module 20, for according in each moving process in geophone group each wave detector base position coordinate, point Do not determine that the base position coordinate of multiple measuring points respectively corresponds to the center position coordinates of geophone group in each moving process;
Processing module 30, for carrying out Transient State Rayleigh Wave detection processing respectively at the base position coordinate of each measuring point, really Depth measurement reference coordinate of the fixed each measuring point to bedrock surface;
The acquisition module 10 is also used to according to the physical location of each wave detector is sat in geophone group in each moving process Mark, obtains the corresponding matched curve of each measuring point respectively;
The processing module 30 is also used to according to the corresponding matched curve of each measuring point, respectively to the depth measurement of each measuring point Reference coordinate is corrected processing, obtains the correction of soundings coordinate of each measuring point.
Optionally, the module 10 that obtains is specifically used for position according to the origin of the datum level, the focal point Set the nearest wave detector of the origin of coordinate, the focal point into each moving process in geophone group from the datum level away from From the distance between with each wave detector in the geophone group, the base of each wave detector in geophone group in each moving process is determined Quasi- position coordinates.
Optionally, for any measuring point, the processing module 30 is specifically used for determining the depth measurement reference coordinate of the measuring point In abscissa it is equal with the abscissa in the base position coordinate of the measuring point;
Acquire the vibration signal of each wave detector in the corresponding geophone group of the measuring point;
Frequency dispersion processing is carried out to the vibration signal of each wave detector, obtains the spectrum curve of the measuring point;
Inversion procedure is carried out to the spectrum curve of the measuring point, obtains the vertical seat in the depth measurement reference coordinate of the measuring point Mark.
Optionally, for any measuring point, the processing module 30 specifically is also used to determine that the correction of soundings of the measuring point is sat Abscissa in mark is the abscissa in the depth measurement reference coordinate of the measuring point;
According to the matched curve of abscissa and the measuring point in the depth measurement reference coordinate of the measuring point, the measuring point is obtained Correction of soundings coordinate in ordinate.
Optionally, the processing module 30 is used for when the corresponding matched curve of each measuring point is straight line, according to each survey The angle of the abscissa direction of the corresponding matched curve of point and datum level, determines the survey of each measuring point to batholith face respectively Deep orientation.
Transient State Rayleigh Wave provided in an embodiment of the present invention detects means for correcting, and above method embodiment can be performed, specific Implementing principle and technical effect, reference can be made to above method embodiment, details are not described herein again for the present embodiment.
The present invention also provides a kind of readable storage medium storing program for executing, it is stored with computer program in readable storage medium storing program for executing, when transient state is auspicious When at least one processor of Leibo detection means for correcting executes the computer program, arouse Transient State Rayleigh Wave detection means for correcting It executes Transient State Rayleigh Wave and detects bearing calibration.
The present invention also provides a kind of program product, which includes computer program, which is stored in In readable storage medium storing program for executing.At least one processor that Transient State Rayleigh Wave detects means for correcting can be read from readable storage medium storing program for executing should Computer program, at least one processor execute the computer program and make Transient State Rayleigh Wave detection means for correcting implementation transient state auspicious Leibo detects bearing calibration.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above-mentioned each method embodiment can lead to The relevant hardware of program instruction is crossed to complete.Program above-mentioned can be stored in a computer readable storage medium.The journey When being executed, execution includes the steps that above-mentioned each method embodiment to sequence;And storage medium above-mentioned include: ROM, RAM, magnetic disk or The various media that can store program code such as person's CD.
Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations;To the greatest extent Pipe present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: its according to So be possible to modify the technical solutions described in the foregoing embodiments, or to some or all of the technical features into Row equivalent replacement;And these are modified or replaceed, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution The range of scheme.

Claims (10)

1. a kind of Transient State Rayleigh Wave detects bearing calibration characterized by comprising
Geophone group is moved according to pre-determined distance in the abscissa direction of datum level, obtains detection in each moving process The base position coordinate of multiple wave detectors in device group, wherein the abscissa of the datum level is any horizontal linear, the benchmark The ordinate in face is any downward straight line perpendicular to the abscissa of the datum level;
According in each moving process in geophone group each wave detector base position coordinate, determine the benchmark of multiple measuring points respectively Position coordinates respectively correspond to the center position coordinates of geophone group in each moving process;
It carries out Transient State Rayleigh Wave detection processing respectively at the base position coordinate of each measuring point, determines the survey of each measuring point to bedrock surface Deep reference coordinate;
According in each moving process in geophone group each wave detector actual position coordinate, obtain each measuring point respectively and respectively correspond to Matched curve;
According to the corresponding matched curve of each measuring point, processing is corrected to the depth measurement reference coordinate of each measuring point respectively, is obtained The correction of soundings coordinate of each measuring point.
2. the method according to claim 1, wherein multiple in geophone group in each moving process of acquisition The base position coordinate of wave detector, comprising:
According to the origin of the datum level, the position coordinates of focal point, focal point into each moving process geophone group In the distance between each wave detector with a distance from the nearest wave detector of origin from the datum level and in the geophone group, determine In each moving process in geophone group each wave detector base position coordinate.
3. the method according to claim 1, wherein being directed to any measuring point, the base position in each measuring point Transient State Rayleigh Wave detection processing is carried out at coordinate respectively, determine each measuring point to bedrock surface depth measurement reference coordinate, comprising:
Determine the abscissa phase in the base position coordinate of the abscissa and the measuring point in the depth measurement reference coordinate of the measuring point Deng;
Acquire the vibration signal of each wave detector in the corresponding geophone group of the measuring point;
Frequency dispersion processing is carried out to the vibration signal of each wave detector, obtains the spectrum curve of the measuring point;
Inversion procedure is carried out to the spectrum curve of the measuring point, obtains the ordinate in the depth measurement reference coordinate of the measuring point.
4. according to the method described in claim 3, it is characterized in that, be directed to any measuring point, it is described respectively to be corresponded to according to each measuring point Matched curve, processing is corrected to the depth measurement reference coordinate of each measuring point respectively, obtains the correction of soundings coordinate of each measuring point, wrap It includes:
Determine that the abscissa in the correction of soundings coordinate of the measuring point is the abscissa in the depth measurement reference coordinate of the measuring point;
According to the matched curve of abscissa and the measuring point in the depth measurement reference coordinate of the measuring point, the survey of the measuring point is obtained Ordinate in deep calibration coordinate.
5. the method according to claim 1, wherein in the geophone group according to each moving process In each wave detector actual position coordinate, after obtaining the corresponding matched curve of each measuring point respectively, the method also includes:
When the corresponding matched curve of each measuring point is straight line, according to the corresponding matched curve of each measuring point and the benchmark The angle of the abscissa direction in face, determine respectively each measuring point to batholith face correction of soundings direction.
6. a kind of Transient State Rayleigh Wave detects means for correcting characterized by comprising
Module is obtained, for moving geophone group according to pre-determined distance in the abscissa direction of datum level, is obtained each In moving process in geophone group multiple wave detectors base position coordinate, wherein the abscissa of the datum level be any level Straight line, the ordinate of the datum level are any downward straight line perpendicular to the abscissa of the datum level;
Determining module, for according in each moving process in geophone group each wave detector base position coordinate, determine respectively The base position coordinate of multiple measuring points respectively corresponds to the center position coordinates of geophone group in each moving process;
Processing module determines each survey for carrying out Transient State Rayleigh Wave detection processing respectively at the base position coordinate of each measuring point Point arrives the depth measurement reference coordinate of bedrock surface;
The acquisition module, be also used to according in each moving process in geophone group each wave detector actual position coordinate, point The corresponding matched curve of each measuring point is not obtained;
The processing module is also used to according to the corresponding matched curve of each measuring point, is sat respectively to the depth measurement benchmark of each measuring point Mark is corrected processing, obtains the correction of soundings coordinate of each measuring point.
7. device according to claim 6, which is characterized in that the acquisition module is specifically used for according to the datum level The origin of origin, the position coordinates of focal point, focal point into each moving process in geophone group from the datum level The distance between each wave detector, determines wave detector in each moving process in the distance of nearest wave detector and the geophone group The base position coordinate of each wave detector in group.
8. device according to claim 6, which is characterized in that be directed to any measuring point, the processing module is specifically used for true Abscissa in the depth measurement reference coordinate of the fixed measuring point is equal with the abscissa in the base position coordinate of the measuring point;
Acquire the vibration signal of each wave detector in the corresponding geophone group of the measuring point;
Frequency dispersion processing is carried out to the vibration signal of each wave detector, obtains the spectrum curve of the measuring point;
Inversion procedure is carried out to the spectrum curve of the measuring point, obtains the ordinate in the depth measurement reference coordinate of the measuring point.
9. device according to claim 8, which is characterized in that be directed to any measuring point, the processing module is specifically also used to Determine that the abscissa in the correction of soundings coordinate of the measuring point is the abscissa in the depth measurement reference coordinate of the measuring point;
According to the matched curve of abscissa and the measuring point in the depth measurement reference coordinate of the measuring point, the survey of the measuring point is obtained Ordinate in deep calibration coordinate.
10. device according to claim 6, which is characterized in that the processing module is used for when each measuring point is corresponding When matched curve is straight line, according to the folder of the corresponding matched curve of each measuring point and the abscissa direction of the datum level Angle, determine respectively each measuring point to batholith face correction of soundings direction.
CN201710897962.5A 2017-09-28 2017-09-28 Transient State Rayleigh Wave detects bearing calibration and device Active CN107728212B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710897962.5A CN107728212B (en) 2017-09-28 2017-09-28 Transient State Rayleigh Wave detects bearing calibration and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710897962.5A CN107728212B (en) 2017-09-28 2017-09-28 Transient State Rayleigh Wave detects bearing calibration and device

Publications (2)

Publication Number Publication Date
CN107728212A CN107728212A (en) 2018-02-23
CN107728212B true CN107728212B (en) 2019-07-02

Family

ID=61208390

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710897962.5A Active CN107728212B (en) 2017-09-28 2017-09-28 Transient State Rayleigh Wave detects bearing calibration and device

Country Status (1)

Country Link
CN (1) CN107728212B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110617404B (en) * 2019-09-09 2021-02-02 安徽理工大学 Gas pipeline leakage point positioning method based on sound wave frequency shift inversion

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009081210A1 (en) * 2007-12-20 2009-07-02 Statoilhydro Asa Method of and apparatus for exploring a region below a surface of the earth
CN102681009A (en) * 2011-12-07 2012-09-19 中冶成都勘察研究总院有限公司 Multi-channel transient Rayleigh wave test system and method
CN203480038U (en) * 2013-09-30 2014-03-12 李振存 Rayleigh-wave detector automatic vibration source device
CN104216003A (en) * 2014-09-20 2014-12-17 中国地质大学(北京) Multichannel transient Rayleigh wave detection method
CN106772554A (en) * 2016-11-28 2017-05-31 安徽理工大学 A kind of multi-channel transient surface wave exploration method under MODEL OVER COMPLEX TOPOGRAPHY

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009081210A1 (en) * 2007-12-20 2009-07-02 Statoilhydro Asa Method of and apparatus for exploring a region below a surface of the earth
CN102681009A (en) * 2011-12-07 2012-09-19 中冶成都勘察研究总院有限公司 Multi-channel transient Rayleigh wave test system and method
CN203480038U (en) * 2013-09-30 2014-03-12 李振存 Rayleigh-wave detector automatic vibration source device
CN104216003A (en) * 2014-09-20 2014-12-17 中国地质大学(北京) Multichannel transient Rayleigh wave detection method
CN106772554A (en) * 2016-11-28 2017-05-31 安徽理工大学 A kind of multi-channel transient surface wave exploration method under MODEL OVER COMPLEX TOPOGRAPHY

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
瞬态面波法在岩土风化分层中的应用;李建宁 等;《物探化探计算技术》;20170430;第32卷(第4期);第121-124页

Also Published As

Publication number Publication date
CN107728212A (en) 2018-02-23

Similar Documents

Publication Publication Date Title
Dinehart et al. Averaged indicators of secondary flow in repeated acoustic Doppler current profiler crossings of bends
Salem et al. Analytic signal approach and its applicability in environmental magnetic investigations
CN101036069B (en) Amplitude preserving prestack migration method
CN110133715B (en) Microseism seismic source positioning method based on first-arrival time difference and waveform superposition
CN105093299A (en) Observation system optimization method based on offset vector tile technology and apparatus thereof
CN109581495A (en) Mountainous region surface velocity model construction method and system
CN107728204A (en) Based on the anisotropic crack prediction method of prestack compressional wave and system
CN105158795A (en) Seam hole detection method by means of stratum pre-stack texture attribute value
US20120186342A1 (en) Method for Detecting Underground Tunnels
Růžek et al. Three-dimensional seismic velocity model of the West Bohemia/Vogtland seismoactive region
Hoover Geophysical choices for karst investigations
Hammer et al. Exploring for stratigraphic traps with gravity gradients
Guan et al. Using fast marching eikonal solver to compute 3‐D Pds traveltime for deep receiver‐function imaging
CN106443776B (en) A kind of submarine seismograph method for relocating based on isochronous surface method
CN107728212B (en) Transient State Rayleigh Wave detects bearing calibration and device
Xu et al. Systemic comparison of seismometer horizontal orientations based on teleseismic earthquakes and ambient‐noise data
Moret et al. VSP traveltime inversion: Near-surface issues
CN102236101A (en) Method and device for predicting multilayer cracks
CN105866839A (en) Static correction method and static correction device on basis of common-attitude trace gathers
CN104142522B (en) A kind of detection method of city buried faults
Karimi et al. Using the ratio of the magnetic field to the analytic signal of the magnetic gradient tensor in determining the position of simple shaped magnetic anomalies
CN112684498A (en) Reservoir fracture prediction method and system based on wide-azimuth seismic data
CN102939548A (en) Method of determining the position of a detector disposed at the bottom of the sea
CN107942373A (en) Coherent algorithm based on the detection of Fractured oil and gas reservoir fracture system
CN113126146A (en) Detection method, control device and storage medium for complex geology of canyon region

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200327

Address after: 350001 4, floor 3, zone F, Fuzhou Software Park, 89 software tower, Gulou District, Fuzhou, Fujian.

Patentee after: FUZHOU HUAHONG INTELLIGENT TECHNOLOGY Co.,Ltd.

Address before: 232001 Anhui province Huainan Taifeng Street No. 168 earth and Environmental Sciences, Anhui University Of Science And Technology

Patentee before: Anhui University of Science and Technology