50 results about "S-wave" patented technology

The invention discloses a multi-wave chromatography and reflection combined detection method and device for a thin-covering-layer urban active fault and a medium. The method comprises the steps that longitudinal wave seismic data and transverse wave seismic data are collected; recording vertical component data based on the longitudinal wave seismic data for first arrival travel time pickup and tomography inversion to obtain a longitudinal wave velocity structure model Vp; extracting Rayleigh surface wave frequency dispersion through the vertical component data and carrying out multichannel surface wave analysis inversion to obtain a transverse wave velocity structure model Vs; calculating a longitudinal and transverse wave velocity ratio parameter structure model Vp / Vs, and calculating to obtain a Poisson's ratio parameter structure model v; obtaining a longitudinal wave (P wave) reflection profile and a transverse wave (S wave) reflection profile; identifying and summarizing into an interpretation fracture set Fc, and recording key elements; and by taking the interpretation fracture set Fc as a structural skeleton, completing unified registration and superposition display on Vp, Vs, Vp / Vs and v, a longitudinal wave reflection profile and a transverse wave reflection profile, and outputting a shallow underground structure comprehensive interpretation map of the thin covering layer region.

The invention discloses a method and a system for calculating three-dimensional S wave velocity structure interpolation based on an OTG model, which are based on a hybrid interpolation model fusing OK, TPS and GCN technologies, dynamically adjust the weight of an S wave velocity interpolation predicted value calculated based on three methods of OK, TPS and GCN through an adaptive fusion mechanism, and calculate the S wave velocity interpolation predicted value based on the OTG model. And combining with a predicted point coordinate to obtain a final three-dimensional S-wave velocity structure interpolation. In the process of calculating the final S-wave velocity structure interpolation, the multi-dimensional composite loss function is fused, the deviation between the interpolation result and the real physical attribute of the underground medium is reduced, the optimized interpolation result is deeply coupled with the three-dimensional geologic modeling process, accurate mapping of the S-wave velocity in the three-dimensional space is achieved, and the accuracy of the S-wave velocity in the three-dimensional space is improved. And high-fidelity input is provided for subsequent tasks. The method provides an effective means for the three-dimensional fine modeling of the S-wave velocity, and has important reference value for improving the exploration efficiency of oil and gas resources and optimizing the early warning of engineering geological disasters.

The invention discloses an elastic reverse time migration imaging method based on deviatoric stress characterization, and relates to the technical field of exploration geophysics. Deviatoric stress is introduced into a conventional elastic wave equation, and a speed-deviatoric stress elastic wave equation describing stress propagation characteristics is obtained through derivation; based on Hooke's law, P wave contribution in total partial strain is removed through partial strain decomposition, decoupling of P-S wave stress components is achieved, and correctness of amplitude, phase and physical unit of a decoupling wave field is guaranteed; and then a joint imaging condition with the average stress as a seismic source wave field and the multi-stress component as a receiving wave field is constructed, PP waves are subjected to normalization superposition through the average stress and the deviatoric stress, and PS waves are subjected to shear stress imaging. According to the method, the sensitivity to heterogeneity and local stress variation of underground media can be remarkably improved, P-S wave crosstalk and imaging noise are effectively eliminated, the imaging precision and resolution of a complex geologic structure are greatly improved, and a reliable technical means is provided for high-precision seismic exploration of complex oil and gas reservoirs.

The application discloses a coal seam roof sandstone aquifer water enrichment evaluation method based on three-dimensional seismic, which comprises the following steps: S10, collecting drilling data, well logging data and high-precision three-dimensional seismic data of a research area; S20, inverting pre-stack seismic data to obtain a P-S wave velocity ratio of the roof sandstone aquifer; S30, performing time-frequency analysis on post-stack seismic data and performing spectral equalization on frequency domain seismic data; S40, inverting a longitudinal wave frequency dispersion attribute of a target layer according to a reflection coefficient approximate formula and the time-frequency spectrum of the post-stack seismic data after the spectral equalization; S50, using an entropy weight method to calculate weight values of the P-S wave velocity ratio and the longitudinal wave frequency dispersion attribute; and S60, constructing a water enrichment indication factor calculation formula and evaluating water enrichment of the sandstone aquifer according to the water enrichment indication factor. The method is based on the theoretical basis of seismic rock physics, has low cost, high operability, more statistical significance, can quickly and accurately evaluate water enrichment of the coal seam roof sandstone aquifer, and provides technical support for development and utilization of coal resources.

The invention discloses a novel multi-component conversion multiple imaging method for a deep-sea rugged seabed structure, and relates to the technical field of geophysical exploration for petroleum, and the method comprises the steps: inputting a longitudinal and transverse wave velocity field, a seismic source wavelet, observation system parameters, and a rugged seabed elevation; constructing a non-uniform curved grid of the acoustic-viscoelasticity model; calculating an acoustic-viscoelastic primary wave field continuation operator of vector wave separation under the curved coordinate system; calculating an acoustic-viscoelastic primary wave adjoint wave field of vector wave separation under the curved coordinate system based on an adjoint state theory; calculating an acoustic-viscoelastic n-order multiple wave field continuation operator of vector wave separation under the curved coordinate system; and calculating an acoustic-viscoelastic n-order multiple accompanying wave field of vector wave separation under the curved coordinate system, and generating a multiple imaging result by applying an elastic multiple imaging condition of vector wave separation. According to the invention, full-path compensation and longitudinal and transverse wave vector imaging of multiple waves can be realized, the imaging range is expanded, and the information amount of imaging is increased.

The invention provides a Sigmoid modulation phase correlation surface wave full waveform inversion method and system. The method comprises the following steps: acquiring a surface wave seismic record and constructing an initial S wave velocity model; performing forward modeling to generate a synthetic record; calculating instantaneous phase correlation of observation and synthetic records to construct an amplitude-independent target function; modulating the target function by using a Sigmoid function, and mapping the value of the target function to a predetermined interval; and calculating a gradient based on the modulated target function and iteratively updating the speed model until a termination condition is met, and outputting a final model. According to the method, the smooth saturation and bounded output characteristics of the Sigmoid function are utilized, the target function response is converted into a self-adaptive gradient scaling mechanism, a self-stabilization optimizer is introduced for the inversion process, and the convergence domain is widened by applying smooth constraint, so that the dependence on the initial model precision is reduced, and the first success rate and the result reliability of inversion under complex conditions are improved.

The invention relates to the technical field of seismic exploration, and particularly discloses a VVAZ and AVAZ combined inversion anisotropy parameter calculation method and device, and the method comprises the steps: carrying out the VVAZ inversion of a five-dimensional pre-stack gather, and obtaining an anisotropy parameter delta; carrying out residual time difference correction on the five-dimensional pre-stack gather to obtain a corrected five-dimensional pre-stack gather; aVAZ inversion is carried out on the corrected five-dimensional pre-stack gather, and anisotropic fracture parameters are obtained; based on logging information of the work area, elastic parameter inversion is carried out on the corrected five-dimensional pre-stack gather, and the longitudinal and transverse wave velocity ratio of the work area is obtained; and solving anisotropy parameters delta N and delta T based on the anisotropy parameter delta, the anisotropy crack parameter and the longitudinal and transverse wave velocity ratio. According to the method, rapid solving of fracture anisotropy parameters deltaN and deltaT is achieved by utilizing the solving relation of anisotropy parameters of VVAZ and AVAZ and the relation among the Thomsen anisotropy parameter delta, the Schobergn anisotropy parameters deltaN and deltaT and the anisotropy gradient Bani.

The application discloses a microseismic source positioning method, device, equipment and medium, and relates to the field of earthquake monitoring, and comprises the following steps: acquiring initial seismic data, converting the initial seismic data into an envelope signal; determining a reference channel, identifying P-wave peaks and S-wave peaks on the envelope signal of the reference channel to determine a reference P-S separation point; aligning the envelope signal by using a dynamic time warping algorithm to obtain an alignment path; mapping the reference P-S separation point to each seismic channel according to the alignment path to obtain an aligned P-S separation point, and determining a wave phase travel time offset according to the alignment path; based on the aligned P-S separation point, extracting P-wave signals and S-wave signals from the envelope signal respectively to determine a P / S wave double-difference correlation function; constructing a weighted imaging operator based on the P / S wave double-difference correlation function and the wave phase travel time offset, and scanning the weighted imaging operator to determine a microseismic source position in a three-dimensional grid map. The source position can be accurately positioned under the 3D grid.

The utility model provides a spiral conical wave-absorbing structure, which relates to the technical field of wave-absorbing, and comprises a base, a plurality of cone bodies which rotate clockwise to rise or anticlockwise to rise are arranged on the base, and the cross section of each cone body is conical. The spiral conical wave-absorbing structure comprises the base and the cone, a plurality of cones rising in a clockwise rotating mode or in an anticlockwise rotating mode are connected to the base in an inserted mode, electromagnetic waves are propagated in the cones in a spiral mode, and the spiral structures of the cones increase the transmission distance of the electromagnetic waves in materials. The retention time of electromagnetic waves in the material is prolonged, and the loss capability of the material to the electromagnetic waves is improved, so that the wave absorbing performance of the material, especially the absorbing performance of a low-frequency part, is improved.

The invention discloses an inversion method and device for a longitudinal-transverse wave velocity ratio of a coal measure strata, and the method comprises the steps: determining the longitudinal-transverse wave velocity ratio reflectivity of a target interval according to the prestack angle gather data of the target interval, and extracting the attribute data of the longitudinal-transverse wave velocity ratio reflectivity through employing the AVO technology; according to the longitudinal and transverse wave velocity ratio curve, the attribute data of the longitudinal and transverse wave velocity ratio reflectivity and preset wavelet extraction time window parameters, wavelet extraction is carried out on the target layer section through the wavelet extraction technology; the longitudinal and transverse wave velocity ratio curve is determined according to logging data; and according to the wavelet and the attribute data of the velocity ratio reflectivity of the longitudinal and transverse waves, performing post-stack compressed sensing reflection coefficient inversion to obtain an inversion result of the velocity ratio reflectivity of the longitudinal and transverse waves. According to the invention, a stable inversion result of the longitudinal and transverse wave velocity ratio can be effectively obtained.

The utility model relates to the technical field of geothermal wells, in particular to a damping ring for damping of a geothermal well, which comprises a circular annular shell, a plurality of damping mechanisms are arranged in the annular shell, and the damping mechanisms are uniformly distributed in the annular shell and are symmetrically arranged in pairs along the radial direction of the annular shell; the damping mechanism comprises a plurality of longitudinal damping assemblies which are longitudinally arranged and a plurality of transverse damping assemblies which are horizontally arranged, the longitudinal damping assemblies are used for conducting damping and energy dissipation on vibration generated by P waves in seismic waves, and the transverse damping assemblies are used for conducting damping and energy dissipation on vibration generated by S waves in the seismic waves. Vibration caused by P waves and S waves in an earthquake can be subjected to vibration reduction and energy dissipation treatment, and the anti-seismic performance of the geothermal well is effectively improved.

The present application provides a kind of ground shallow well microseismic monitoring event identification technology, device, equipment and medium.The identification method includes: based on the three-dimensional space scanning horizontal direction and depth direction of S wave, the event energy superposition value of each grid point at each time is obtained, and the multi-channel superposition value S (t) is calculated;The maximum value of the event energy superposition value of each grid point in space at each time is obtained, and the maximum value of the multi-channel superposition value S (t) of each point at each time is updated as M (t);The envelope B (t) of M (t) is calculated by Hilbert transform, and the maximum value of M (t) greater than B (t) in fixed time window is detected as the event in the time window.The present application establishes the microseismic judgment criterion of three-dimensional space P-S wave multi-channel superposition according to the ground shallow well microseismic monitoring mode and the three-component recording characteristics, realizes the automatic detection of microseismic event based on event similarity method, the operation speed is fast, the efficiency is high, and the microseismic event can be detected automatically and quickly.

The invention belongs to the field of rock-soil earthquake engineering, and particularly relates to a method for evaluating the anti-seismic performance of an end bearing pile in stratified soil under the action of earthquake S waves. According to the method, a pile-soil model is established on the basis of a continuous medium theory under a linear elastic frame, and displacement of a free field and a scattered field of a soil body is decoupled through a Helmholtz decomposition method; the soil layer boundary is processed by adopting the orthogonality of a modal function characteristic function, and the motion response of the pile foundation is described by adopting a strict closed form solution. Frequency domain amplitude-frequency characteristic analysis is adopted, a motion corresponding coefficient and a motion amplification coefficient are introduced to quantitatively evaluate the influence of the soil layer thickness and the soil layer rigidity ratio on the anti-seismic performance of the end bearing pile, and theoretical support is provided for pile foundation optimization design.

This invention relates to an intelligent P / S wave arrival time picking method based on structure and target synergistic enhancement, comprising the following steps: Step 1, constructing a bimodal joint probability modeling system and generating P / S wave joint probability distribution labels; Step 2, designing a HybridAttnUNet1D network architecture based on the task paradigm of Step 1; Step 3, executing a general method flow for model training based on the labels of Step 1 and the model of Step 2; Step 4, realizing P / S wave arrival time inference picking and performance verification based on the model trained in Step 3. This invention jointly predicts P and S waves, effectively preserving the relative temporal relationship between P and S waves, thereby improving the stability and noise resistance of arrival time picking. This invention also designs a HybridAttnUNet1D network model, which aims to take into account both the multi-scale local features and global temporal dependencies of microseismic signals. The network input is the original one-dimensional microseismic waveform time series, and the final output is a P / S wave bimodal joint probability distribution sequence of the same length as the input.

The utility model discloses a counter flow type cooling tower, which comprises a frame structure and a cooling tower body arranged on the frame structure, an air duct, a water collector, a water distribution system and a filler are sequentially arranged along the cooling tower body from top to bottom, a fan is arranged in the air duct, the filler comprises a filler layer with an S-wave-shaped surface, and the surface of the filler layer is provided with a water inlet and a water outlet. The wave crest and the wave trough of the S wave form are arc-shaped, and the extending direction of the S wave is consistent with the water flow direction and the air flow direction; according to the utility model, not only is the heat exchange efficiency higher, but also dirt, impurities or sediments are not easy to accumulate between wave crests and wave troughs of the filler, water flow blockage caused by blockage is avoided, and the cooling efficiency is ensured.

The application relates to a method for constructing a low-frequency model of a fluvial facies dense heterogeneous reservoir, and relates to the technical field of geophysical reservoir prediction and seismic inversion, and comprises the following steps: dividing a fluvial facies dense heterogeneous reservoir into lithofacies types and counting the proportion of each lithofacies; obtaining a P-wave impedance, a S-wave impedance and a density parameter inversion three-dimensional data body after first-time inversion; obtaining a P-wave impedance, a S-wave impedance and a density parameter low-frequency model after last-time iteration, and a P-wave impedance, a S-wave impedance and a density parameter three-dimensional data body after prestack inversion of last-time iteration; calculating a P-S wave velocity ratio data body to obtain a fluvial facies dense sandstone reservoir distribution rule. The method for constructing a low-frequency model of a fluvial facies dense heterogeneous reservoir disclosed by the application provides an accurate low-frequency model for reservoir prediction, and the prestack inversion result can better depict channel distribution characteristics.

This invention provides a method for evaluating the safety of climbing utility poles based on structural deformation and material properties. The method includes: collecting the basic parameters of the utility pole; calculating the pole's tilt angle based on the basic parameters and issuing a hazard warning based on the tilt angle; calculating geometric evaluation indicators based on the basic parameters and the tilt angle; calculating the pole's wave velocity based on the basic parameters, judging the material quality, and calculating material performance indicators; calculating the pole's burial depth based on the wave velocity, judging the burial depth quality, and calculating a burial depth safety index; calculating a comprehensive climbing safety index based on the geometric evaluation index, material performance index, and burial depth safety index, and classifying the utility pole's safety based on the comprehensive climbing safety index. This method quantifies latent deterioration, forming a standardized evaluation system, avoiding misjudgments, thereby improving the accuracy and objectivity of pole climbing safety evaluation and effectively reducing operational safety risks.

The invention relates to a seismic source parameter prediction method based on physical constraint and array geometric adaptation, which comprises the following steps of: determining input, output and initial prediction models, sequentially obtaining station characteristics and predicted seismic source parameters, comparing the predicted seismic source parameters with a seismic source true value to obtain regression loss (if the seismic source true value does not exist, the regression loss is not recorded); obtaining phase loss according to the predicted seismic source parameter and the seismic phase label, and then obtaining total loss; carrying out a training process on the model according to the previous steps, and storing an optimal model; inputting related data into the optimal model, and obtaining predicted seismic source parameters; the phase loss only needs the arrival time data of the P wave and the S wave, and is the signal feature which is the most easily obtained signal feature, so that the application range is relatively wide; phase loss takes the physical law of seismic wave propagation as a hard constraint to be embedded into a model optimization process, so that the interpretability of a prediction result is improved; therefore, the method is wide in application range and high in interpretability when used for predicting the seismic source parameters.

This invention discloses an in-situ protection and reinforcement device for the ancient embankment heritage of the Grand Canal, belonging to the field of ancient embankment protection technology. It includes a plate installation assembly, an impact dispersion assembly, and an impact conversion structure. The main body of the baffle is made of high-strength, corrosion-resistant composite material. The bottom conical insertion nails guide and embed into the underwater soil, ensuring precise initial positioning and strong anchoring force, enabling rapid vertical fixation and effectively preventing soil collapse and water erosion. The device is installed in situ without large-area excavation, minimizing disturbance to the ancient embankment heritage and its surrounding environment, adhering to the principle of minimal intervention in cultural heritage protection. The arc-shaped mesh plate on the wave-facing side uses its curved surface to redirect the flow and permeable holes to divert water, initially dispersing continuous surging waves and reducing impact intensity. Combined with connecting rods, sliding blocks, dampers, and buffer springs, a mechanical buffer circuit is formed, converting impact kinetic energy into elastic potential energy and damping dissipation energy. This provides efficient buffering for small waves and stable unloading for large waves, significantly improving the device's wave resistance and anti-overturning capabilities.

The invention relates to the technical field of ground pressure early warning, and discloses a mine broadband earthquake impact ground pressure early warning method and device, and the method comprises the steps: arranging shot points and a three-component detector on the same side wall of a roadway, exciting slot waves, collecting signals, extracting reflected waves through time-frequency, band-pass and median filtering, and positioning a broken structural surface; for a target area, collecting a microseismic original waveform, automatically scanning a potential event, picking up P / S wave arrival time, and inverting a wave velocity model to calibrate a seismic source position; dispersing the broken structural surface into grids, counting the number of micro-seismic events in a set time-space, and calculating a dynamic fracture probability; and constructing a gray model in combination with the stress sequence of the anchor rod strain gauge, generating a cumulative reduction stress prediction value, and constructing risk early warning classification through the stress prediction value, the crushing dynamic probability and the crushing frequency. According to the method, the stress spatio-temporal evolution law is captured through the gray scale prediction model, the migration path of the stress concentration area is speculated, and real-time early warning can be carried out on the roadway stress change.

The embodiment of the invention provides a high-precision seismic elastic wave field forward modeling method, device and equipment and a storage medium, and the method comprises the steps: carrying out the discretization of a wave equation in an elastic medium through the combined application of an improved time high-order implicit staggered grid difference method, and obtaining two types of time-space domain frequency dispersion relationships related to a P wave and an S wave based on the plane wave arrangement and simplification; solving the two types of time-space domain dispersion relations of the P wave and the S wave by combining a least square method and a Taylor series expansion method to obtain high-order difference coefficients corresponding to the P wave and the S wave; based on a wave field separation technology, a P wave field component is calculated by using a high-order difference coefficient corresponding to a P wave, an S wave field component is calculated by using a high-order difference coefficient corresponding to an S wave, and the components are superposed to obtain a high-precision full-elastic wave field. According to the technical scheme, under the same parameter condition, better simulation precision can be generated, so that an effective wave field continuation scheme is provided for subsequent elastic wave migration imaging and waveform inversion.

The invention discloses an undrilled formation velocity prediction method, system and product based on VSP longitudinal and transverse waves, the method is excited through VSP longitudinal and transverse wave double seismic sources and directly receives downlink waves and uplink waves through a detector arranged in a well, a high-resolution time-depth relation and interval velocity information can be obtained, the high-resolution time-depth relation and interval velocity information comprise underground deep-layer fidelity high-resolution reflection characteristics, and the high-resolution time-depth relation and the interval velocity information can be used for predicting the velocity of an undrilled formation. A reliable starting point is provided for speed prediction below the shaft bottom; secondly, according to the corresponding rule of the speed change and the high-resolution longitudinal and transverse wave data reflection characteristics, the relation between the reflection characteristics and the speed change is established, the speed of the underground undrilled formation is finely predicted through a joint inversion means, a more real and reliable data basis is provided for research of underground exploration and development and the like, and the method is suitable for popularization and application. The exploration and development risk and the drilling cost are greatly reduced.

The present specification relates to the technical field of exploration geophysics, and particularly relates to an anisotropic elastic wave decoupling method and device and computer equipment. The anisotropic elastic wave decoupling method comprises: determining a set of Thomson parameters included in an anisotropic model according to a received wave field decomposition request to be decomposed, the wave field decomposition request to be decomposed comprising a wave field to be decomposed; converting the set of Thomson parameters to obtain an initial set of elastic parameters; performing P-S wave velocity separation processing on the initial set of elastic parameters to obtain a target P-wave elastic parameter set and a target S-wave elastic parameter set; and substituting the target P-wave elastic parameter set and the target S-wave elastic parameter set into the anisotropic model to process the wave field to be decomposed to obtain a target P-wave matrix and a target S-wave matrix. The embodiments of the present specification simplify the P-S wave field decomposition process and reduce the calculation cost.

The present application relates to a seismic phase picking method based on time-frequency analysis and deep learning, belonging to the field of seismology, a deep neural network comprising an empirical mode decomposition module, a residual shrinkage module, an LSTM layer, a transformer layer, an up-sampling layer and a convolution layer is established, a linear rectifier function is used as an activation function, an Adam adaptive optimization function is used for back propagation, a binary cross-entropy loss function is used as a loss function, and the average absolute error function, error standard deviation, precision, recall and F1 score are used as evaluation indexes to obtain a deep neural network model. When there are differences in the frequencies of noise, P waves and S waves, empirical mode decomposition can effectively separate the three, thereby helping the network to better identify P waves and S waves. Even for wideband noise, the noise will be decomposed into different intrinsic mode functions, thereby reducing the noise content in the intrinsic mode function where the P wave and S wave are located. This time-frequency analysis method can effectively improve the noise immunity of the network.

The invention discloses a ship sound vibration characteristic identification and positioning method based on a DAS system. The method comprises the following steps: S1, obtaining monitoring data based on the DAS system; s2, classifying data based on two modes of frequency domain features and differential distance, and extracting underwater acoustic signals and vibration signals; s3, performing feature extraction and beam forming on the underwater acoustic signal; s4, carrying out beam forming on the vibration signal, and outputting a final beam forming result through combined judgment with beam forming of the underwater acoustic signal; obtaining the time difference between the P wave and the S wave based on the speed difference between the P wave and the S wave; and S5, finally realizing target identification based on feature extraction of the underwater acoustic signal, realizing target orientation estimation based on a combined comparison result in the step S4, and determining the distance between the target and the optical cable array based on the time difference between the P wave and the S wave. The method effectively solves the technical problems that in the prior art, accurate positioning is difficult due to single dependence on underwater acoustic signals, and target recognition cannot be achieved only according to vibration signals.

The invention belongs to the technical field of elastic wave reverse time migration, and particularly relates to a decoupling Poynting vector based on divergence and rotation operators and an application method of the decoupling Poynting vector in elastic wave reverse time migration, and the method comprises the following steps: P wave and S wave decoupling and Poynting vector calculation based on the divergence and rotation operators; the invention relates to an elastic wave RTM imaging condition based on decoupling Poynting vector constraint. According to the method, the Poynting vector calculation method compatible with divergence / rotation wave field decoupling is developed theoretically, the comprehensive advantages of the method in the aspects of improving the elastic wave RTM imaging quality (especially the reliability of PS wave imaging), effectively suppressing noise and improving the calculation efficiency are shown in the practical application level, and the method has important scientific value and wide application prospects and is worthy of popularization and application. And an effective new way is provided for further improving the practical level and the imaging precision of the elastic wave RTM under the existing computing resource condition.

The application provides an ice cover granular snow layer density estimation method and system based on a vehicle-mounted noise source, relates to the technical field of physical exploration and ice cover detection, and the method comprises the following steps: according to the exact start and end time period of each vehicle-mounted noise event, a corresponding multi-channel vehicle-mounted noise event data is obtained by intercepting from an original continuous noise record; the mutual coherence processing and linear superposition of the intercepted multi-channel vehicle-mounted noise event data are carried out to obtain a virtual shot record with improved signal-to-noise ratio; P wave and surface wave in the record are jointly inverted to obtain a P wave and S wave velocity structure, and on the basis, a double-parameter granular snow layer density formula fitting is carried out. The application is suitable for realizing the quantitative estimation of the regional scale ice cover granular snow layer density structure by using passive source seismic signals, and provides parameter support for ice cover engineering activity site evaluation and environmental evolution research.

The invention discloses a microseismic event positioning method based on microseismic signals, based on a microseismic acquisition device, the microseismic acquisition device comprises a plurality of three-component microseismic detectors arranged in a long drill hole at equal intervals, and the three-component microseismic detectors are used for acquiring the microseismic signals generated by excitation of a microseismic event in the deep position of a coal seam in the long drill hole. Comprising the following steps: S1, acquiring a microseismic signal by using a microseismic acquisition device; p waves and S waves of microseismic signals are simultaneously collected through a microseismic collection device, and an absolute distance to a seismic source is directly provided for each receiver by using the time difference between the P waves and the S waves, so that the positioning process is converted into a spherical surface intersection problem with each receiver as a sphere center and the calculated distance as a radius. According to the method, systematic errors caused by uncertainty of the earthquake generating moment are eliminated, and a large amount of redundant data is provided for cross validation, so that the precision and reliability of a positioning result are greatly improved. And the accuracy and comprehensiveness of signal acquisition are improved.

The invention provides a mine Ordovician limestone water pressure and top interface crack detection method based on a micro-seismic technology. The method comprises the following steps: collecting actual micro-seismic data; de-noising processing is carried out on the collected actual micro-seismic data by using a wavelet threshold de-noising algorithm; a P-S wave combined positioning algorithm based on P wave and S wave first arrival is adopted to position the denoised micro-seismic event; seismic source parameter calculation and seismic source mechanism inversion are carried out based on the positioned micro-seismic event waveform and the initial polarity; ground stress inversion work is carried out, and the change condition of Ordovician limestone water pressure is deduced; calculating advanced attributes representing stratum stress or fracture characteristics on the basis of micro-seismic event source parameters, determining an advanced attribute threshold range causing Ordovician limestone water damage, and judging possible crack conditions; the Ordovician limestone top crack condition and Ordovician limestone water pressure are comprehensively analyzed, and whether the risk that Ordovician limestone water protrudes along the crack exists or not is judged. The method can monitor the water pressure change of the Ordovician limestone in real time, accurately detect the water flowing fracture of the Ordovician limestone top interface, and provide a scientific basis for mine water disaster prevention and control.

The invention relates to a seismic data processing method and device, equipment, a storage medium and a program product, and relates to the technical field of geophysical exploration. The method comprises the following steps: acquiring seismic profile data of a target area; obtaining logging data corresponding to one or more drilling wells in the target area; determining well seismic calibration results corresponding to the one or more drilling wells according to data of the positions of the corresponding drilling wells in the seismic profile data and the logging data; the well-seismic calibration result is used for indicating amplitude data corresponding to the longitudinal wave and the converted shear wave in the synthesis channel and the seismic channel respectively; reconstructing the seismic profile data according to the well seismic calibration results corresponding to the one or more drilling wells; and obtaining a multi-wave fusion data volume of the target area according to the reconstructed seismic section data. When the fusion data volume of the scheme is used for carrying out subsequent depiction, the data of the longitudinal waves and the converted transverse waves do not need to be independently depicted and artificially fused.