43 results about "Exploration geophysics" patented technology

The invention discloses an elastic reverse time migration imaging method by combining seismic multi-component and belongs to the field of exploration geophysics. The method comprises the following steps: directly inputting seismic multi-component data; without performing the wave field separation on the input data, extrapolating forwards at inverse time based on an elastic wave equation; building an underground elastic vector seismic wave field by combining with multi-component; acquiring four elastic wave imaging results having specific physical significance, such as longitudinal wave, transverse wave, transition longitudinal wave and transition transverse wave, by applying an elastic wave cross-correlation imaging condition of lighting compensation; and suppressing the low wave number noise in the elastic reverse time migration by performing low-pass angle filtering, thereby acquiring an imaging result of final elastic wave migration. By directly inputting the seismic multi-component data and building the underground elastic vector seismic wave field by combining with multi-component, the method can be used for accurately imaging for complex earth medium.

The invention belongs to the field of exploration geophysics and particularly relates to an elastic seismic wave field forward modeling method. The seismic forward modeling method includes the steps of firstly, according to the principle of wave field vector decomposition, decomposing a seismic wave field into a longitudinal wave field and a transverse wave field and respectively calculating corresponding propagation matrixes of the longitudinal wave field and the transverse wave field; secondly, according to the principle of low-rank decomposition, decomposing the propagation matrixes to obtain propagators of vector longitudinal wave and vector transverse wave so as to reduce calculation amount of numerical simulation to improve calculation efficiency; and thirdly, coupling updated records of the vector longitudinal wave and the vector transverse wave acquired after the numerical simulation so as to obtain a forward modeling result of the seismic wave field. The calculation precision is high due to the fact that the propagation matrixes include compensations to medium parameter and numerical simulation parameter. Numerical dispersion of seismic wave field records acquired by the forward modeling method is nearly eliminated and stability of the numerical simulation is high; accordingly, the forward modeling method is a high-precision elastic wave field numerical simulation method.

The invention belongs to the field of exploration geophysics, and particularly relates to a carbonatite oil and gas reservoir crack earthquake detection method. The method comprises the steps of: carrying out amplitude versus offset and azimuth angle analysis, speed versus offset and azimuth angle analysis, frequency versus offset and azimuth angle analysis respectively based on pre-stack seismic azimuth angle dataset, further carrying out ellipse fitting respectively to obtain a crack development azimuth angle, and predicting a crack development azimuth through a seismic attribute comprehensive analysis way; carrying out quality factor parameter extraction, anisotropic parameter inversion and non-hyperbolic moveout velocity analysis respectively based on pre-stack seismic trace dataset and shear wave velocity data, and predicting a reservoir crack development region by comprehensively utilizing the parameters above through a seismic attribute comprehensive analysis way. The carbonatite oil and gas reservoir crack earthquake detection method can give more comprehensive consideration to the geology of the reservoir and the physic characteristics of the globe, and the estimated crack azimuth and the predicted crack development region are more exact.

The invention discloses a method for establishing an earthquake response mode based on a fracture equivalent medium model, which belongs to the field of exploration geophysics. The method is characterized by comprising the following steps of: determining a relationship between fracture characteristic parameters and equivalent medium elastic parameters by quantificationally applying a fracture equivalent medium theory, establishing a constitution relation of a fracture network anisotropic equivalent medium model based on a Bond conversion and superposition principle, and placing the constructed fracture network equivalent medium model in an actual space position of strata, and establishing an actual stratigraphic fracture equivalent medium model; and establishing an earthquake response model of different tectonic positions and different lithologic fractured reservoirs by use of an earthquake wave field numerical simulation method. The method for establishing the fractured reservoir earthquake response mode based on the fracture equivalent medium model has strong practicability, and provides the basis for identifying the fractured reservoir from the prestack earthquake data.

The invention provides an adaptive hybrid norm dictionary learning seismic wave impedance inversion method, and belongs to the field of seismic inversion. Through the combination of the post-stack wave impedance inversion in exploration geophysics with dictionary learning and denoising, the adaptive denoising and dictionary learning based post-stack wave impedance inversion method can be provided.The method adopts an adaptive hybrid norm objective function to replace an original 2-norm objective function, so that an objective function can better overcome non-gaussian noise; well logging information is introduced through a dictionary learning method, so that influences of different seismic noise can be effectively overcome, and the vertical resolution of inversion can be enhanced by fullyutilizing the well logging information; and compared with traditional methods of least squares, robustness can be stronger, the noise immunity and inversion result precision of the seismic wave impedance inversion can be enhanced, multiple solutions can be reduced, and convenient data processing in geological exploration can be achieved.

A plane seismic exploration signal decomposition method based on a Kriging method belongs to the field of exploration geophysics. The method comprises the following steps: acquiring a seismic exploration signal, setting the basic lag distance and the maximum lag distance, and calculating a two-dimensional experimental variation function of a seismic attribute slice to obtain a two-dimensional variation function graph; determining ellipse information corresponding to target information from the variation function graph, extracting experimental variation function curves of long-axis and short-axis directions (respectively corresponding to primary and secondary directions of the target information), and using a theoretical model to perform one-dimensional second-order nested fitting; determining the initial parameters, including sill, long-axis range, ellipse direction and the ratio of long axis to short axis, of the two-dimensional variation function according to the parameters of the separated variation function curves of the primary and secondary directions, tuning parameters, selecting different models (spherical model, exponential model, Gauss model), and making the fitted two-dimensional variation function more consistent with a target signal experimental variation function; and determining a two-dimensional theoretical variation function, substituting the function into a Kriging interpolation equation, and decomposing to obtain the target information.

The invention relates to the technical field of exploration geophysics, and in particular relates to a near-surface scattering surface wave removing method. According to the method provided by the invention, the problem that surface wave noise, direct waves and upstream reflected waves, which are caused by near-surface heterogeneity induction in seismic exploration, are affected by heterogeneous body scattering waves; scattering occurs in the transition from body waves to surface waves and from surface waves to body waves; scattering waves are usually close to or greater than reflection records; surface wave denoising is carried out on actual records by means of a time-frequency analysis technology; scattering surface waves are removed by a direction adjustable filter; and time-frequency analysis and the direction adjustable filter are combined to remove the surface wave information, which can effectively improve the signal to noise ratio of a signal.

The invention discloses a method extracting fault anomaly from geophysical prospecting gravity anomaly and belongs to exploration geophysics in the earth sciences. An existing method for separating regional anomaly and local anomaly (a trend analysis method, a frequency domain filtering method, a moving average method, upward extension and the like) has the defects that the boundary effect is obvious and spatial structural characteristics of the regional anomaly are not considered, and accordingly the separated and extracted regional anomaly is not high in accuracy. For obtaining a high-accuracy method for separating the regional anomaly and local anomaly of the fault, a structure function method is given. Under the condition that a spatial distribution geometrical characteristic, a covariance function and the regional anomaly form of gravity observation data are known, the fault region anomaly is obtained by means of the optimal, linear and unbiased interpolation estimation method. The trend analysis method is only a particular case of the structure function method. The structure function method is suitable for separation and extraction of the fault gravity anomaly with a linear structure. An experiment proves that the structure function method is effective and feasible on the aspect of implementation.

The invention discloses a transient electromagnetic real-time dynamic advanced detection method and system based on a magnetic gradient tensor, and belongs to the field of exploration geophysics. Themethod comprises the steps: enabling a transient electromagnetic field to be equivalent to a series of continuous static magnetic fields, enabling a water conduction structure to be equivalent to a magnetic dipole body, directly observing magnetic field intensity components of a secondary field in x, y and z directions through five observation points which are arranged in the center of a roadway space and form a cross-shaped measuring line, and performing pairwise differential combination to obtain a magnetic gradient tensor; performing inversion calculation through the obtained magnetic gradient tensor to realize positioning of the magnetic dipole. According to the invention, data processing is directly carried out based on the magnetic gradient tensor, the precision and the accuracy arehigher than those of a scalar measurement and inversion method, and the device is simple, has instantaneity and can be effectively applied to advanced detection and monitoring of hidden water inrush disaster-causing factors in the roadway tunneling and working face mining process.

Embodiments of the present application provide a process communication method and apparatus in a cluster system. The method includes: determining a startup relationship between process numbers, wherein the startup relationship includes a cell division type startup relationship; starting a process in a cell division manner according to the startup relationship between the process numbers; causing each process to send address information thereof to the 0 process during the process startup; performing inter-process communication; when first communication is made between the processes, the two communicating parties inquire about the address information of the other party from the 0 process. The method and the device can improve the parallel processing efficiency of a cluster system and have the advantages of high performance, high scalability and ease of use; the method and teh device can meet the needs of the development of a massive data parallel computing framework in computation-intensive and data-intensive exploration geophysics, and significantly improve the development efficiency and quality of large-scale parallel and distributed seismic applications.

The invention discloses a mixed domain seismic migration Hessian matrix estimation method, and relates to the technical field of exploration geophysics, and the method comprises the steps: firstly obtaining an initial imaging result through employing a conventional migration method, then taking the migration result as a reflectivity model, carrying out the linear forward modeling and migration, and obtaining a secondary imaging result, using the two offset results for estimating an unstable phase filter in a space-wavenumber domain for approximating a Hansen matrix, calculating the inverse of the Hansen matrix by using mixed domain deconvolution, and using the estimated inverse of the Hansen matrix as a functional gradient in least square offset to improve the convergence rate of inversion, so that the inversion accuracy is improved. And the least square inversion is converged to an accurate solution within 4-5 iterations. The method has the advantages that the method is high in calculation efficiency, the least square migration convergence speed can be increased, the calculation cost of least square reflectivity inversion is reduced, and the method has good application prospects.

The invention discloses a seismic wavelet signal extraction method based on the maximum joint entropy; a seismic wavelet estimate method in the exploration geophysics and the maximum joint entropy criterion in information theory learning are combined, and the maximum joint entropy criterion can replace the original least squares criterion, thus allowing an object function can better solve non-gauss noises. Compared with the conventional least squares criterion, the seismic wavelet signal extraction method based on the maximum joint entropy is better in robustness, higher in precision, and can better process the geology exploration data.

The invention belongs to the field of exploration geophysics, and relates to an elastic reverse time migration method based on an imaging energy flow vector. The method comprises the following steps:on the basis of an elastic wave forward continuation operator, obtaining an elastic wave mass point vibration velocity vector field at a seismic source end; on the basis of the elastic wave reverse continuation operator, obtaining an elastic wave stress tensor field at a detection end; carrying out dot product operation on the seismic source end mass point vibration velocity vector field and the detection end stress tensor field, and performing constructing to obtain an imaging energy flow vector of each grid node at a check point moment; and imaging the constructed imaging energy flow vectoraccording to the energy flow vector divergence imaging condition to obtain an elastic wave reflection energy imaging result. According to the method, the extraction of the standardized elastic wave reflection energy imaging result is realized from the vectorized elastic wave field, and the obtained imaging result is clear in physical significance; moreover, imaging is carried out only by using theelastic wave mass point vibration velocity vector field at the seismic source end and the elastic wave stress tensor field at the detection end, so that the storage consumption is low and the storageefficiency is high.

The invention discloses an angle domain least square reflectivity inversion method based on coherent superposition, which relates to the technical field of exploration geophysics, and comprises the following steps of: parameterizing and linearizing a wave equation by using an angle domain reflection coefficient as a reflectivity model, and constructing a linear positive operator based on Bern approximation; then, calculating the gradient of linear inversion by using an adjoint state method; and finally, updating the reflectivity model in an angle domain by using a precondition conjugate gradient method; and carrying out automatic event tracking on the reflection axis on the inverted angle domain reflectivity gather, and then carrying out coherent superposition along the tracked event to obtain a final superposition imaging result. The method has the beneficial effects that the problem of offset illusion caused by speed errors is avoided, and a deep high-quality imaging result is obtained; the sensitivity of the least square migration to the speed error is reduced, the deep-ultra deep oil gas target can be accurately imaged, and finally the practical process of the least square migration in deep oil gas exploration is promoted.

The invention discloses a method for estimating the velocity of transverse waves based on the inversion of an equivalent elastic modulus for self-adapting matrix minerals, which belongs to the field of exploration geophysics, and used for solving the problem that in a conventional method for estimating the velocity of transverse waves based on a Xu-White model, an elastic modulus for matrix minerals is difficult to set accurately. The method comprises the following steps of: (1) setting the Poisson ratio of a dry core sample at a single logging depth point and the value selection range and initial value of an equivalent bulk modulus for matrix minerals; (2) carrying out inversion solution on a same fluid item from two different angles by using two different fluid item calculation methods so as to invert an optimal Poisson ratio of the dry core sample and an optimal equivalent bulk modulus and an optimal equivalent shear modulus for matrix minerals; and (3) estimating the velocity of transverse waves by using a simplified Xu-White model. In the method disclosed by the invention, the equivalent elastic modulus for matrix minerals is taken as an object of study, through reducing the number of target parameters, the precision and reliability of the method for estimating the velocity of transverse waves are effectively improved.

The invention discloses an exploration geophysics teaching experimental data collection device based on a wireless network. A 24-bit ADS127L01AD conversion circuit with the collection rate of 512kSPS is used as the core of an analog circuit part for data collection. The analog part further comprises a front end amplifier circuit, a linear amplifier circuit, a filter forming circuit and a differential input circuit. According to the invention, the operation function of exploration geophysics experiment is realized through the wireless network; and geophysics experimental data collection and experimental teaching are well combined, which is convenient for students to understand and convenient for teachers to carry out experimental teaching.

The invention belongs to the field of exploration geophysics, and relates to a viscoelasticity parameter extraction method of a carbonatite fractured reservoir. The viscoelasticity of a rock matrix of the carbonatite fractured reservoir and the anisotropy induced by fractures are introduced into the method, and an equivalent medium model of the viscoelasticity of the fractured reservoir is established, wherein a quality factor, the fracture density, the fracture opening degree and the fracture occurrence are used as the parameters. A rotating staggered mesh high-order finite difference method is used for conducting seismic wave field simulation, and forward modeling records of various fracture combinations in different buried depths and at different construction positions are obtained. On the basis of seismic wave field transmission mechanism analysis, a fracture sensitivity parameter is obtained. Optimizing processing is conducted on high precision prestack earthquake location angle set extracted through a sticky sound wave reverse time migration method according to the fracture sensitivity parameter. The prestack earthquake location angle set after optimizing processing is used for extracting the viscoelasticity parameter of the carbonatite fractured reservoir. The viscoelasticity parameter extraction method of the carbonatite fractured reservoir is suitable for fracture detection of the carbonatite fractured reservoir.

The invention discloses a transient electromagnetic real-time dynamic advanced detection method and system based on magnetic gradient tensor, which belongs to the field of exploration geophysics. The invention equates the transient electromagnetic field into a series of continuous static magnetic fields, guides water-bearing structures, etc. The valence is a magnetic dipole body. Through the five observation points arranged in the center of the roadway space to form a "ten"-shaped survey line, the magnetic field intensity components of the secondary field in the x, y, and z directions are directly observed, and two-by-two The differential combination yields a magnetic gradient tensor. Then, the inverse calculation is performed through the obtained magnetic gradient tensor to realize the positioning of the magnetic dipole. The present invention processes data directly based on the magnetic gradient tensor, and its precision and accuracy are higher than those based on scalar measurement and inversion, and the device is simple and instant, and can be effectively applied to concealed detection during roadway excavation and working face mining. Advance detection and monitoring of water inrush disaster-causing factors.

The invention discloses a double-loop zero-flux transient electromagnetic detection device, and relates to the field of geology and exploration geophysics. The double-loop zero-flux transient electromagnetic detection device comprises a transmitting coil, a transmitter and a receiving coil, wherein the transmitting coil is used for generating a primary pulse magnetic field and adopts a double-loop structure; the transmitter is used for supplying power to the transmitting coil; and the receiving coil is used for receiving a secondary induction eddy current field caused in the underground medium during the intermittent period of the primary pulse magnetic field. Compared with the prior art, the device has the beneficial effects that the transmitting coil in the scheme adopts the double-loop transmitting coil, the receiving coil and the transmitting coil are vertically arranged and coincide in center, the magnetic flux of a primary field in the receiving coil is zero, the mutual inductance influence between the transmitting coil and the receiving coil is eliminated, the receiving coil observes a pure secondary field, the detection blind area is reduced, and the detection precision is improved; and the secondary field directions of geologic bodies on the two sides of the coil are opposite, so that geologic anomalous bodies on the two sides of the transmitting coil can be distinguished.

The invention relates to the field of exploration geophysics, and discloses a device for quickly positioning a ditch leakage point and a control method thereof. According to the device and the controlmethod thereof, a transmitting electrode and a receiving electrode are driven to move in water by utilizing an overwater mobile carrier, an abnormal point is determined according to an electric fielddifference in the transverse direction of a ditch, positioning information is acquired and recorded by means of a Beidou or GPS positioning device, and a leakage point is further determined in an abnormal region, thereby solving the problems that an existing geophysical method is not suitable for the ditch and the ditch leakage point cannot be rapidly, precisely and automatically positioned.

The invention belongs to the technical field of seismic physical models of oil-gas exploration geophysics, and discloses a stress-sensitive compact sandstone reservoir physical model material, a stress-sensitive compact sandstone reservoir physical model and a preparation method. The physical model material comprises quartz sand, kaolin, polyurethane, a curing agent and a diluent. The method comprises the following steps: carrying out heating treatment on polyurethane and a curing agent; performing curing pretreatment on the inner surface of the physical model curing mold; uniformly mixing the quartz sand and the kaolin, adding polyurethane, a curing agent and a diluent, and uniformly mixing to obtain a mixture; and putting the mixture into a physical model curing mold, and carrying out pressing treatment and curing treatment on the mixture to obtain the stress-sensitive compact sandstone reservoir physical model. According to the method, the soft organic material is used for cementing the quartz sand and the kaolin, so that physical simulation research on the ground stress-seismic response change rule of the compact sandstone reservoir is realized.

The invention discloses an elastic medium undulating surface Fresnel beam migration imaging method, and relates to the field of exploration geophysics. According to the method, the seismic wave beam has the effective half-width of 1 / 4 wavelength near the near surface, and the calculation precision of the Green function of the near-surface region can be remarkably improved. The novel seismic wave beam is applied to continuation of a multi-wave multi-component seismic wave field, an elastic medium Fresnel beam migration imaging method is developed, and the near-surface imaging quality of land data can be remarkably improved while accurate imaging of a middle-deep part is guaranteed.

The invention discloses a data driving method and system for seismic wave first arrival pickup, and relates to the technical field of exploration geophysics, and the method comprises the steps: carrying out the linear superposition summation of all first virtual seismic traces formed at the position of a stationary-phase seismic source, and obtaining a first linear superposition result; performing phase weighted stacking on the first linear stacking result by using a first weighting coefficient to obtain a second virtual seismic trace; obtaining a super virtual seismic trace according to the second virtual seismic trace; performing linear superposition summation on all the super virtual seismic traces generated at the stable-phase detection point position to obtain a second linear superposition result; and performing phase weighted stacking on the second linear stacking result by using a second weighting coefficient to obtain a final super virtual seismic trace of the target trace gather in the target shot gather. According to the invention, under the condition that the signal-to-noise ratio of the original data is very poor and the original data volume is limited, the signal-to-noise ratio of the seismic data first-arrival wave is obviously improved, so that high-quality first-arrival wave input data is provided for seismic wave first-arrival pickup.

The invention belongs to the field of exploration geophysics, and in particular relates to a forward modeling method for elastic seismic wave field. First, the method decomposes the elastic seismic wave field into a longitudinal wave field and a shear wave field based on the principle of wave field vector decomposition, and calculates the corresponding propagation matrices of the longitudinal and shear waves respectively. Parameter compensation, so it has high calculation accuracy; then, based on the principle of low-rank decomposition, the propagation matrix is ​​decomposed to obtain the propagation operators of the vector longitudinal wave and the vector shear wave, thereby reducing the calculation amount of numerical simulation and improving the calculation efficiency; finally, the The updated vector P-wave records and vector S-wave records obtained by numerical simulation are coupled to obtain the seismic wavefield forward modeling results. The seismic wave field record obtained by the forward modeling method constructed by the invention has almost no numerical dispersion and high numerical simulation stability, and is a high-precision elastic wave field numerical simulation method.

The present invention is applicable to the technical field of exploration geophysics, and provides a full-band magnetic sensor. The full-band magnetic sensor includes an induction coil, a voltage measurement module and an integrator. The full-frequency magnetic sensor detects the induced electromotive force generated by the induction coil through the voltage measurement module. Finally, the integrator used has the characteristics of large low-frequency magnification and small high-frequency magnification, so as to realize the expansion of the measurement bandwidth, and extend the measurement bandwidth of the magnetic sensor to the ultra-low frequency of 0.001Hz. Under the condition of no magnetic flux negative feedback Realize the measurement of the whole frequency band.

The invention is suitable for the technical field of exploration geophysics, and provides a full-band magnetic sensor. The magnetic sensor comprises an induction coil, a voltage measurement module andan integrator. The full-band magnetic sensor detects an induced electromotive force generated by the induction coil through the voltage measurement module, and then the integrated has the characteristics of a large low-frequency amplification factor and a small high-frequency amplification factor, thereby achieving the extension of measurement bandwidth, enabling the measurement bandwidth of themagnetic sensor to extend to the ultralow frequency of 0.001Hz, and achieving the full-band measurement under the condition that there is no flux negative feedback.

The invention discloses a micro-seismic first arrival intelligent pickup method, system and device based on an FFPN model, and a storage medium, and belongs to the field of exploration geophysics, and the model provided by the invention retains the multi-scale feature fusion capability of an original FPN network model, is suitable for a first arrival intelligent pickup problem, and can achieve a more refined first arrival output result. According to the network model, output of two scale features is selected as an optimization condition of the model, finer first-arrival pickup output is calculated through low-level features with rich detail features, dichotomous output is calculated through high-level features with a larger receptive field, and the output loss of the two parts is calculated. And back propagation is carried out in the training process to optimize the network model, so that a more accurate dichotomy result (a first arrival picking result) is obtained. Moreover, the point-aware loss function proposed by the model can calculate the loss of low-level feature output, and the low-level features are optimized to achieve the effect of optimizing high-level feature output.

The invention belongs to the technical field of seismic physical models of oil-gas exploration geophysics, and discloses a hole type reservoir seismic physical model material, a hole type reservoir seismic physical model and a manufacturing method. The physical model material comprises the following components in parts by weight: 15-25 parts of kaolin, 20-100 parts of water, 1000 parts of quartz sand, 60-100 parts of epoxy resin and 5-20 parts of a curing agent. Kaolin is used for simulating a collapse object, the kaolin and water are mixed in proportion and then quick-frozen into an ice-like solid, the ice-like solid is processed into hole models of different forms and sizes under the low-temperature condition, the hole models are embedded into a high-speed surrounding rock material to simulate a hole reservoir, and a hole space can be formed after the hole-type reservoir seismic physical model is solidified. According to the method, the seismic physical model with higher similarity can be obtained, so that more accurate experimental data can be obtained to be used for carrying out hole reservoir oil-gas exploration method research.

The invention discloses a seismic attenuation estimation method, system, equipment and storage medium based on sparse S transformation and adaptive logarithmic spectral ratio method, belonging to the research field of exploration geophysics. The invention uses sparse S transformation to solve the time spectrum of seismic data , compared with the traditional S transform, the time spectrum energy is more concentrated, which is beneficial to improve the accuracy of the Q value estimated by the logarithmic spectral ratio method; the present invention introduces an adaptive frequency band selection method, which avoids the artificial influence that may be caused by manually selecting the frequency band range. , which improves the reliability of the Q-value estimation results. Secondly, the invention uses sparse S transform, which has a certain inhibitory effect on noise, which makes it possible to avoid serious interference of noise when extracting the spectrum at the upper and lower layers of the target area and use the logarithmic spectral ratio method to estimate the Q value, and improves the Noise immunity of Q-value estimates.