57 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 an ultra-deep reflection imaging method based on passive source seismic data and a related device, and relates to the technical field of earth depth structure imaging, and the method comprises the following steps: firstly, for original continuous seismic waveform data, intercepting seismic waveform data with a first arrival of a P wave and a non-arrival of an S wave, and carrying out the preprocessing of the data, so as to obtain a far earthquake event data segment; the method comprises the following steps: firstly, analyzing and identifying a far earthquake event with a high signal-to-noise ratio by using a time-frequency-wavenumber domain, removing data which are greatly influenced by noise, and then, reconstructing a reflection wave field of a high-quality far earthquake data fragment by using a multi-dimensional deconvolution technology, and carrying out convolution on a seismic source wavelet to obtain ultra-deep reflection wave field data. And finally, performing reflection imaging processing in combination with velocity analysis and superposition imaging technologies, optimizing the analysis capability of a deep reflection interface, and finally obtaining an earth deep reflection imaging profile, so that a continuous structure interface extending from the shallow earth crust to the deep mantle can be clearly revealed, and coupling characteristics between shallow and deep earth structures can be effectively captured.

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 present invention discloses a method for determining the equivalent design wave of a ship's wave load under level 6 sea conditions, relating to the field of ship structural safety. The method comprises: treating the hull as a variable-section beam, and for a ship at a constant speed, calculating the wave load transfer function at the hull beam section under different wave conditions, wherein the wave conditions include different wave direction angles and different circular frequencies of the incident wave; based on the wave load transfer function at the hull beam section, calculating the short-term predicted extreme values ​​of the wave load under different wave direction angles and different zero-crossing periods within the significant wave height range of level 6 sea conditions; selecting a maximum value from all the short-term predicted extreme values ​​of the wave load, and determining the equivalent design wave parameters of the wave load under level 6 sea conditions based on the maximum value. The equivalent design wave is a sinusoidal regular wave, and its parameters include wave amplitude, circular frequency, and wave direction. By converting the study of the load on the hull under random waves into the study of the hull load under regular waves, this method is of great significance for subsequent research on hull load levels and structural safety.

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 a wave dissipation structure and a revetment, and relates to the field of river revetments. The wave dissipation structure comprises a bank protection retaining wall and a wave dissipation water tank, when the water level of a river channel changes, water flows in and out of the containing space through the water permeable holes, the wave dissipation water tank conducts buoyancy self-adaptive lifting motion, the wave dissipation channel continuously aligns to the wave face, wave energy is dispersed and consumed by floating damping of the wave dissipation channel and the wave dissipation water tank, and therefore the purpose of reducing bank protection impact is achieved. Water flow enters the containing space through the water permeable holes to consume energy and is secondarily broken by the wave dissipation channels of the floating water tank, and the scouring force of ship traveling waves is remarkably weakened. The self-adaptive water level fluctuation design solves the problems that a traditional retaining wall with the fixed height is insufficient in wave dissipation at a high water level and shields sight at a low water level; the water area and the land landscape plant area are communicated through the wave dissipation channel, the ecological connectivity is improved, water exchange and biological migration are promoted, and the ecological blocking problem of a hard structure is solved.

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 invention provides a high-speed wave head penetrating ship, and relates to the field of ship engineering, the high-speed wave head penetrating ship comprises a high-speed wave head penetrating ship and a stem post, the bottom edge is in a micro-oval shape, the top is in an inverted triangle shape, the high-speed wave head penetrating ship protrudes towards the head, the side view geometric feature is in a cutting edge shape, and the stem post is of a forward-inclining type structure, a backward-inclining type structure and a vertical type structure. According to the method, the high-speed wave piercing head is designed into a specific geometrical shape, the size of the high-speed wave piercing head meets the cross section area curve distribution requirement of the whole ship, the head protruding length is related to the wavelength and the speed, and the cross section fat and thin degree meets the design parameters of the molded line inflow section, so that the wave-making resistance and the shape resistance during ship navigation are effectively reduced; meanwhile, the seakeeping performance of the ship in the waves is improved, the pitching amplitude, wave attack and splashing are reduced, stall in the waves is reduced, the excellent and comprehensive sailing performance is achieved, the dual technical advantages of high sailing speed and high seakeeping performance are achieved, and the ship is particularly suitable for the high-speed ship in the transition state.

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 provides a microseismic event S-wave first arrival pickup method, which belongs to the technical field of seismic waveform data processing, and comprises the following steps: carrying out detrending and band-pass filtering preprocessing on three-component original data; calculating a covariance matrix, solving a characteristic value and a characteristic vector, and carrying out polarization correction; a Bezier curve is adopted to accurately position waveform cross points; and determining the S-wave first arrival by combining the waveform cross point and the amplitude characteristics. According to the method, the dielectric anisotropy influence is eliminated through polarization correction, accurate intersection point detection is achieved through the Bezier curve, the amplitude judgment criterion is combined, the accuracy and reliability of S-wave first arrival recognition are remarkably improved, and environmental noise interference is effectively restrained. The method realizes full-process automatic processing, and is suitable for the fields of micro-seismic monitoring, seismic exploration and the like.