875 results about "Wave impedance" patented technology

The invention provides a method for predicting a dolomite reservoir. The method includes acquiring seismic attributes of a target stratum, increasing the signal-to-noise ratio of the seismic attributes, adopting the seismic attributes with the signal-to-noise ratio increased to perform post-stack wave impedance inversion; adopting the seismic attributes with the signal-to-noise ratio increased to perform pre-stack elastic parameter inversion; adopting the results of the post-stack wave impedance inversion and / or the pre-stack elastic parameter inversion to predict the dolomite reservoir of the target stratum. By an effective predicting method technical system of seismic data quality evaluation, pre-stack depth processing, post-stack attribute analysis and post-stack wave impedance inversion lithologic identification, pre-stack attribute analysis and pre-stack elastic parameter inversion of a high-quality reservoir and hydrocarbon detection of the dolomite reservoir, the dolomite reservoir can be effectively predicted.

The invention discloses a time-domain failure distance measuring method with two asynchronous ends of a direct-current power line and self-adaptive parameters, which comprises the following steps of: setting up an asynchronous time-domain failure observation equation of the direct-current line, which includes quantities of fault distance, distributed resistances of the line, transmission wave velocity, wave impedance, asynchronous time difference of data at the two ends and the like to be observed; taking the time of initial failure traveling waves which arrive at the two ends of the line respectively as reference time of data at each end; introducing a traveling wave characteristic equation; and eliminating the quantities to be observed of asynchronous time difference in the time-domain failure observing equation, thus obtaining a novel time-domain failure distance measuring method of the direct-current line so as to realize self-adaption of the line parameter self-adapting without synchronizing with a clock. The time-domain failure distance measuring method has higher transition resistance durability and can effectively improve the accuracy and reliability of failure distance measurement at the two ends of the direct-current line.

The invention relates to the earthquake exploration position calibration technology based on the pre-stack wave field simulation, belonging to the seismic data processing and interpretation technology field in oil exploration. The invention considers the seismic information acquisition and processing impact, considers multiple waves and conversion waves action in the reflection waves, enhances the accuracy of integrated seismic records, and enhances the consistency of integrated seismic records and real seismic records, this approach including the following eight steps: filtering and editing of density and acoustic logging well curve; using density and acoustic logging well curve to calculate the wave impedance curve; using statistics sub wave to integrate the seismic records; comparing the integrated seismic records and well-side seismic road, properly stretching the logging well curve; using well-side seismic road and logging well reflection coefficient curve to extract certainty sub wave; using the reflection rate method to simulate the pre-stack common reflection points roads set; simulating the real seismic data processing, and obtaining the final integrated seismic records; using the final integrated seismic records and comparing with the well-side seismic road, to process seismic exploration position calibration.

The invention discloses a prediction and evaluation technology of 3D seismic optimum time window river channel sand body reservoirs, belongs to the technical field of the prediction and evaluation of the 3D seismic reservoirs, and aims at solving technical problem that the river channel predictive resolution of the traditional methods is not enough. The technical proposal is as follows: the 3D visualization scanning is performed on each reservoir at the interval of 1-2ms, the optimum time window is determined according to the scope shown by the target river channel, the corresponding subfield is cut out, clairvoyance and scanning are carried out on the time window properties, such as root-mean-square amplitude, wave impedance and the like, auto-tracing is performed, top surface and bottom surface are picked up, a time isopach map is calculated and converted into a sand body isopach map, a top surface structure diagram is formed by time-depth conversion, and the reservoir physical properties are evaluated by curve reconstruction, thus realizing the prediction and evaluation of the plane morphology, longitudinal thickness and the reservoir physical properties of the river channel sand body. The method adopts the optimum time window to effectively inhibit interference, is applicable to various data volumes, and can effectively predict and evaluate a thin river channel sand body with the thickness far less than 1 / 4 of a wavelength under the condition of frequent interbedding of sand and mudstone, the thin river channel sand body comprises the river channel sand body which is not corresponding to wave crest or wave trough, and the method has good application effect in petroleum exploration and development.

The invention discloses a method and a device for analyzing the dispersion and attenuation of unsaturated double-porosity medium earthquake waves. The method comprises the following steps of: 1, acquiring geological data including penetration rate, porosity, shale content and mineral components and generating a reservoir rock dry framework model according to the geological data; 2, acquiring measured data including drilling data, logging data, logging information and fluid experiment measured data and generating a reservoir fluid model according to the measured data and a Van derWaals equation; 3, solving plane waves according to the reservoir rock dry framework model, the reservoir fluid model and a Biot-Rayleigh equation for describing the spread of the unsaturated double-porosity medium earthquake waves, and generating phase speed and inverse quality factors of longitudinal waves and horizontal waves; and 4, generating relationships between frequency and each of speed, attenuation,wave impedance, AVO (Amplitude Versus Offset) response characteristic and the like according to the phase speed and inverse quality factors of the longitudinal waves and the horizontal waves.

The invention provides a method for modeling a velocity model of a geologic structure, which is characterized by comprising the steps of performing stacking velocity analysis on a common midpoint gather to obtain stacking velocity; obtaining a layer plane graph of every layer of a wave impedance interface according to the stacking velocity; generating a root-mean-square velocity body by means of the stacking velocity; obtaining root-mean-square velocity of every layer by means of the generated layer plane graphs and the generated root-mean-square velocity body; obtaining initial layer velocity of every layer according to the obtained root-mean-square velocity by controlling velocity inversion; converting a time domain of layer information of every layer into a depth domain by utilizing the initial layer velocity and the graph migration technique to obtain layer information of the deep domain of every layer; and obtaining a first velocity model by utilizing the initial layer velocity of every layer and the layer information of the deep domain of every layer.

This invention disclose a kind of high resolution nonlinear reservoir object inversion method, uses a kind of nonlinear reservoir lithologic object convolution module to build the function relation of wave impedance and porosity / argillaceous content; describes the physical relation of complex rock by analyzing multilevel structure and determining two-way edge wavelet; realizes the combination of the certainty inversion, statistic inversion and nonlinear theoretics; then realizes the high resolution reservoir object parameter inversion of composite geology, well logging and earthquake wave impedance information base on both inverse operator inversion and positive operator reconstructed algorithm, gets high resolution reservoir object cross section (such as porosity, argillaceous content).

The invention is a seismic wave impedance inversion method, determining amplitude and phase spectrums of pos-stack seismic data, and extracting seismic subwaves on a near-well seismic channel and determining amplitude and phase spectrums, spectrum widths and master frequencies of the seismic subwaves; picking up seismic stratigraphic horizons corresponding to geological stratigraphic horizon and checking and correcting as well as smoothing and interpolating the picked horizons; demarcating logging acoustic time difference and density curves, and demarcating depth-domain logging acoustic time difference and density curves into time-domain ones, consistent with horizons of stacked seismic data, generating in-well wave impedances; making wavelet odd analysis and wavelet characteristic extraction on the stacked seismic data; providing high resolution wave impedance parameters for reservoir description in the exploration, development and recovery of oil fields.

The invention relates to a simple multi-parameter seismic inversion method for improving the precision and the reliability of the conventional pre-stack seismic inversion with the help of the seismic wave elastic impedance, thereby simultaneously extracting a variety of formation parameters. The steps are as follows: full wave logging information is utilized; an acoustic impedance AI and an elastic wave impedance EI are determined; a CMP gather is converted into an angle gather, and part of the angles are stacked; the conventional pre-stack elastic impedance inversion process is utilized to respectively carry out the AI and the EI inversion of cross-well measuring line of small angle stacked data and part of the stacked data of the middle incident angles; the channel calculation of the EI and the AI data which are obtained by inversion is carried out, thereby obtaining the seismic attributes which are closely related to the elastic density ratio, the lithology changes and the pore fluid ingredients; the reservoir lithology and the properties of the pore fluid are determined. The simple multi-parameter seismic inversion method only needs the seismic data of near and middle angles, and other parameters are obtained through the channel calculation of the relationship between the elastic density ratio and the acoustic impedance, thereby improving the calculation efficiency and being used in the pre-stack inversion of the seismic data obtained by the vast majority of the collection conditions.

The invention provides a method and device for predicating sand body thicknesses through logging constraint wave impedance inversion. The method includes the steps of acquiring corresponding seismic data and geological logging data after carrying out seismic exploration drilling sampling in a target exploration area, carrying out logging constraint wave impedance inversion on the target exploration area according to the seismic data with the geological logging data as constraint conditions to obtain a logging constraint wave impedance inversion model of the target exploration area, selecting a plurality of sampling points respectively and evenly for each sand group in the target exploration area on an inversion section of the logging constraint wave impedance inversion model to obtain sand body total wave impedance values of corresponding sand groups and sand body thickness relation curves of the corresponding sand groups through fitting, and obtaining the sand body thicknesses of the corresponding sand groups according to the relation curves. The method and device for predicating the sand body thicknesses through logging constraint wave impedance inversion improve reliability of sand body distribution and thickness prediction in stratums with complicated geological conditions, different reservoir stratum thicknesses and large physical property changes.

The invention relates to a two-end distance measuring method of transmission line fault under a dynamic condition, comprising the following steps of: acquiring two-end voltage and current signals through two-end phasor measurement units of the transmission line; then acquiring dynamic positive sequence impedance, admittance, wave impedance and propagation coefficients according to line parameter estimation and a dynamic line parameter method; and finally solving a nonlinear equation related to a fault distance according to the fact that positive sequence voltages derived from two ends to a fault point are equal by applying a Newton iteration method so as to obtain a dynamic distance measuring result of the transmission line fault. The two-end distance measuring method can more effectivelysolve the problem of line parameter and fault distance estimation under the dynamic condition and has accurate and reliable fault distance measuring result.

The invention discloses a simulation method for calculating the lightning back-striking performance of UHV transmission lines. The simulation method comprises the following steps: dividing wires among long-span transmission towers into a plurality of line segments, establishing a distributed parameter circuit model for each line segment, establishing a multi-wave impedance model for each tower, then selecting the lightning parameters according to lightning activities, integrating all models and calculating the lightning withstand level, and calculating to obtain the back-striking trip rate of the UHV transmission lines according to the lightning withstand level and a series of formulae. The simulation method disclosed by the invention establishes the distributed parameter models of the line segments in consideration of high long-span transmission towers and long line spans, subdivides the wave impedance of different parts of the towers, and establishes the multi-wave impedance models in consideration of the refraction and reflection processes of lightning waves on the tower bodies so as to make the simulation results more precise and practical.

The invention relates to an inversion method for a stratum elasticity constant parameter of pre-stack seismic data in geophysical prospecting for petroleum. The method comprises the following steps of: acquiring seismic data, and processing the pre-stack seismic data; obtaining logging data and extracting angle gather seismic data to form all stacked data; performing horizon calibration on the seismic data and the logging data; inverting the seismic data which is stacked at different incident angles to obtain a wave impedance parameter; calculating a corresponding reflection coefficient; generating a logging and seismic data pair; calculating a weight function and a model parameter of a function mapping network model, an impedance parameter and an elasticity constant; and drawing an elasticity constant profile for reservoir lithologic identification, petroleum-gas prediction, oil-water interface determination and petroleum reservoir description. By the method, the elasticity constant can be calculated by the conventional inversion of the seismic data and the logging data; and the method has certain adaptability to the inversion of faults and pinch-out zones, a wide application range, high resolution, high calculation speed, high stability, high calculation accuracy and the characteristic of certain anti-noise property.

The invention provides a multi-resolution wave impedance inversion method based on model constraints, and belongs to the field of seismic exploration and development of oil gas and coal bed methane. The method comprises the steps that an initial model used for wave impedance inversion is established by utilizing seismic records, well-logging data and seismic interpretation horizon data, and scale factors are led in to speed up the computation of algorithms and to reduce the solution multiplicity of the inversion. According to the method, the ground seismic data, the well-logging data and the seismic interpretation horizon data are organically combined, and a certain adaptability on the inversion of faults, the inversion of wedge-shaped stratums and the inversion of thin stratums is achieved; furthermore, the method has the advantages of being high in computational speed, good in stability, high in computational accuracy and strong in anti-noise capacity.

The invention provides a method for calculating harmonic wave impedance of an alternating current system for direct current transmission. The method includes the following steps that a node admittance matrix is formed according to a harmonic wave model of a transmission line, a transformer, a generator, a load and a parallel capacitor; nodes to be calculated are placed at the bottom line of the admittance matrix, the current value at the bottom line is set to be 1, and the current values at other lines are set to be 0; node optimizing numbering is conducted on the first n-1 nodes; the node admittance matrix is stored through a sparse technology; the node admittance matrix is solved, and the back substitution does not need to be conducted; a plurality of input files are simultaneously calculated through a parallel calculation method; statistic analysis is conducted according to the batch processing calculation result. According to the method, specific to the characteristic that multiple times of scanning needs to be conducted for calculating the harmonic wave impedance, the improved node optimizing numbering technology and the improved sparse technology are adopted, and therefore the workload for calculating the harmonic wave impedance is reduced, the calculating speed is greatly increased, burdens of operators are relieved, and precious time can be saved.

Provided is a multi-mode multi-frequency power amplifier (100), comprising: a controller (120) for receiving and outputting control signals based on external signals; a broadband amplifier channel (110) for receiving single-frequency or multi-frequency radio-frequency (RF) signals via an input end, conducting power amplification on the RF signals, and outputting the signals via an output end; a fundamental wave impedance convertor (170) comprising a first segment (130) shared by RF signals in respective frequency bands, a second segment (150) dedicated for RF signals in each frequency band, and a switch circuit (140), the switch circuit (140) being controlled by the controller (120) to switchably separate the RF signals with amplified power to the second segment (150) for multiple outputs.

A system and method for transferring wave energy into or out of a periodic structure having a characteristic wave impedance profile at a prime frequency, the characteristic wave impedance profile comprising a real portion and an imaginary portion, comprising: locating one or more energy transfer elements each having a wave impedance at the prime frequency approximately equal to the real portion of the characteristic wave impedance at one or more points on the periodic structure with the imaginary portion approximately equaling zero; and employing the one or more energy transfer elements to transfer wave energy into or out of the periodic structure. The energy transfer may be repeaters. Quarter-wave transformers can be provided at one or more points on the periodic structure with the imaginary portion approximately equaling zero to transmit waves across one or more discontinuities. A terminator can be employed for cancellation of waves. The invention substantially eliminates reflections of the wave energy at the prime frequency by joints between sections of the periodic structure.

The invention discloses a single-ended transient signal-based failure direction judgment element and a single-ended transient signal-based failure direction judgment method for an ultrahigh speed direct current transmission line. The element comprises a failure component calculation unit, a forward wave and backward wave energy calculation unit, and a comparison and command transmission unit. A failure direction is judged by comparing a calculated forward wave energy value with a calculated backward wave energy value, a failure is determined as a positive failure if the forward wave energy value is lower than the backward wave energy value, and on the contrary, the failure is determined as a negative failure, and a corresponding command is transmitted. A test result shows that the failure direction judgment element for the direct current transmission line can be used for accurately distinguishing the positive and negative failures of the line under different direct current transmission projects and the operating conditions of the different direct current transmission projects only by utilizing single-ended voltage and current signals of the line without the influence of a failure position and transition resistance and that only a normalized wave impedance value is required to be calculated for the protection and setting of the different direct current transmission projects so as to reduce the workload of protection workers and field operators.

The invention discloses a method for improving offshore seismic data resolution by utilizing air gun seismic focuses of different depths. The method is characterized in that air guns are respectively arranged in positions of different depths to be excited, the near-field wavelet of each seismic focus is recorded, and a far-field wavelet is obtained through simulation; and because of the existence of an offshore strong wave impedance interface, trap points occur on the frequency spectrum of each far-field wavelet, the high-frequency information of air gun seismic focuses at superficial parts is rich, the low-frequency advantages of seismic focuses at deep parts are obvious, and periodic extension is formed. The spectral analysis is carried out on the far-field wavelets to determine the dominant frequency band of each wavelet, the high-frequency dominant section of the wavelet at a superficial part is kept unchanged, the far-field wavelet of the seismic focus at a deep part is expected to be output only at the low-frequency end, and a frequency division matched filter is obtained by calculation and then is applied to all seismic records acquired by the air gun seismic focuses at the superficial parts, thereby realizing the combination of advantages of wavelets of seismic focuses in two different depths, widening the frequency band of the offshore seismic data and improving the resolution.

The invention discloses a method for improving the bench blasting quality of a deep-hole of super-hard rocks. The method includes the following steps that blast holes are formed in a step of a blast area, hole mesh parameters are determined, and perpendicular drilling is conducted; explosive loading structures are determined, and the explosive loading structures for segmented explosive loading and mixed explosive loading are adopted in the blast holes at the same time; explosives are determined, wherein the explosives are selected according to according to matching coefficients of explosive wave impedance and rock wave impedance; blasting parameters are determined and include the length W of a resistance line of a base plate, the hole pitch a, the row distance b, the super depth h and the blast hole filling length l; and blasting is conducted, wherein the explosives are loaded, the blast holes are filled with stemming on the tops of the blast holes, reverse blasting is conducted from the bottoms of the blast holes, and blasting is conducted in an inter-hole elementary error sequence. According to the method for improving the bench blasting quality of the deep-hole of super-hard rocks, due to the fact that the hole mesh parameters, the explosive loading structures, the explosives, the blasting parameters, the blasting position and the blasting sequence are reasonably selected, multi-way comprehensive functions are achieved, the performance of the explosives is fully achieved, the sufficient reaction of the explosives is ensured, the utilization rate of the blast holes is increased, and the boulder frequency and bedrocks of rock blasting are reduced.

The invention discloses a high-efficiency class-F and inverse class-F power amplifier. A transistor, a parasitic compensation circuit, a harmonic control circuit and an output fundamental wave impedance matching circuit are included. The harmonic control circuit is located between the parasitic compensation circuit and the output fundamental wave impedance matching circuit, and the parasitic compensation circuit is located between the transistor and the harmonic control circuit. As for a fixed work frequency, the input end of the harmonic control circuit forms a second harmonic short dot and a third harmonic open-circuit dot, the harmonic control circuit is composed of three micro-strips, and the parasitic compensation circuit is composed of an L-type micro-strip structure. A parasitic component compensation function can be achieved by adjusting electric length parameters of a micro-strip in the parasitic compensation circuit, and meanwhile drain electrode impedance conditions of the class-F and inverse class-F power amplifier can be achieved. By means of the class-F and inverse class-F power amplifier, influences of the parasitic component on the harmonic control circuit can be effectively reduced, precise control of second harmonics and the third harmonics is achieved, and work efficiency of the power amplifier is improved.

The invention relates to a simulation calculation method of total lightning stroke trip rate, in particular to the simulation calculation method of the total lightning stroke trip rate of a long-span circuit. The simulation calculation method of the total lightning stroke trip rate of the long-span circuit is based on a formula that: N is equal to N counter plus N shield; wherein, the N is the total lightning stroke trip rate of the long-span circuit, the N counter is the countering trip rate of the long-span circuit, and the N shield is the shielding failure trip rate of the long-span circuit, that is, the invention includes the simulation calculation method of the countering trip rate of the long-span circuit and the simulation calculation method of the shielding failure trip rate of the long-span circuit, the calculation method of the countering trip rate of the long-span circuit is based on an EMTP procedure, and the invention provides a complete simulation method by taking into consideration of the impacts of various factors such as pole tower wave impedance and wave speed, induction voltage component, air clearance flashover criterion, working voltage phase and ground line arc altitude.

The invention discloses a prediction method for reservoir thickness which is based on dipole wavelet. Steps are as follows, parameters are extracted, time difference of sound wave between a reservoir and a non-essential reservoir as well as the density are resolved; the time difference of sound wave and the density are converted into a time domain by time-depth conversion, the reflectances of the top and bottom of the reservoir are computed, and a time domain model at the well point is set up, the velocity of longitudinal wave, the density, the wave impedance, and the reflectances are involved in the time domain model; the seismic wavelet is extracted and optimized; a seismic responsing pattern store of reservoir is built by the following steps, the seismic wavelet is taken as a basic wavelet function, the reflectances of the top and bottom of the reservoir are respectively coefficients c and d of the dipole wavelet function to be inputs of the dipole wavelet function, and a reservoir seismic responsing pattern store with a given reservoir section of the work area is set up, spectrum of reservoir thickness is set up; the distribution of the reservoir in three dimensions is predicated quantitatively.

The invention provides a stratum aperture pressure prediction method based on variety earthquake attributes, which indirectly constructs the non-linear relation between the stratum aperture pressure and the variety earthquake attributes and realizes the prediction of the stratum aperture pressure. The advantages of the invention reside in two aspects. On one hand, the invention comprehensively uses variety earthquake attributes, and deducts other influence factors, like porosity, shale content, etc, while utilizing the inversed longitudinal wave speed to carry out stratum aperture pressure prediction so as to improve the prediction accuracy; on the other hand, the invention creatively uses two calculation methods to obtain the longitudinal wave speed, obtain the target sand-mud stratum longitudinal wave speed through high precision wave impedance inversion and the obtains the background longitudinal wave speed by using the DIX formula inversion. The former can relatively accurately inverts the longitudinal wave speed inside the target stratum so as to relatively calculate the effective stress Pe inside the target layer; and the later can reflect the relation between the longitudinal wave speed and the density as a whole so as to accurately calculate the overlying strata pressure Pov.

The invention discloses a pole tower impact wave impedance detecting method and device based on a tower top reflected wave. The method comprises the following steps: enabling one end of a current leadto be connected with a positive electrode of an impact generator, and enabling the other end of the current lead to be connected with the top of a pole tower, wherein a negative electrode of the impact generator is connected to the ground; enabling one end of a voltage measuring probe to be connected with the current lead connected to the top of the pole tower, and enabling the output end of thevoltage measuring probe to be connected with an oscilloscope, wherein the other end of the voltage measuring probe is connected to the ground; acquiring an impact step wave according to the impact generator, thereby acquiring a voltage waveform of the tower top according to the impact step wave, wherein the voltage waveform comprises an incident voltage wave and a reflected voltage wave; and according to an amplitude value of the incident voltage wave, an amplitude value of the reflected voltage wave and wave impedance of the current lead, acquiring impact wave impedance of the pole tower. Themethod is capable of measuring the impact wave impedance of the pole tower, effectively improving accuracy and reliability of detection, small in error, and simple and easy to implement.

The invention provides a reservoir wave impedance prediction method based on dipole wave, which comprises the following steps that: logging data is analyzed to acquire the corresponding parameter data of the seismic response characteristics of a reservoir; well seismic calibration is carried out according to the logging and seismic data so as to extract wavelets and determine the time-depth relationship; the seismic response characteristic wave function of impedance changes is constructed according to the parameter data, the wavelets and the time-depth relationship; the way integral characteristics of the logging data are extracted according to the seismic response characteristic wave function to form a seismic response model database; the seismic data beside a well is analyzed according to the seismic response model database to work out an impedance interpretation spectrum of the seismic data beside the well; a seismic trace is analyzed according to the seismic response model database to work out a data body wave impedance spectrum; and the wave impedance is predicted according to the seismic response model database, the impedance interpretation spectrum and the data body wave impedance spectrum. According to the reservoir wave impedance prediction method based on the dipole wave, the wave impedance of the reservoir is predicted through the seismic characteristics and wave impedance change rules based on the analysis to the logging data, so that the precision in predicting the wave impedance of the reservoir and the resolution to thin reservoirs are improved.

The invention discloses a well drilling rod and a pulse plasma drilling machine system. The well drilling rod comprises at least one well drilling rod unit, wherein each well drilling rod unit comprises at least two open-close type support rings and a connecting rod; two adjacent support rings are fixedly connected through one connecting rod; the well drilling rod units are connected in sequence; and two adjacent support rings of two adjacent well drilling units are fixedly connected. The well drilling machine system comprises a well drilling mechanism comprising the well drilling rod, a high-pressure pulse power supply, a high-pressure low-wave impedance coaxial water cable, a water circulating device and a drill bit electrode for discharging towards rock stone, so that the drill bit electrode can be always contacted with the rock stone and keep drilling vertically downwards, and the problems of collapse of a water pipe and a cable and the like occurring during deep drilling in the conventional plasma drilling machine are avoided simultaneously; and moreover, the length of the well drilling rod is adjusted by increasing or decreasing the well drilling rod units, so that well drilling of different depths can be realized.

The invention relates to a well-to-seismic integration paleo-formation pressure prediction method. The method includes the steps that a velocity cube, a density cube and a porosity cube of a formation are obtained through the fine reservoir wave impedance inversion technology; time domain and depth domain conversion of a three-dimensional model is completed by means of the velocity cube; a geologic model of the geological time of a pool-forming period is obtained by subtracting the thickness of a formation needing to be eroded away in the pool-forming period from the thickness of the current formation; a paleo-porosity cube is obtained according to the erosion compaction restoration thicknesses of grids in the longitudinal direction of the geologic model of the pool-forming period and the porosity and depth relation; the paleo-density of the formation is restored by means of the density and depth linear relation; a paleo-formation interval velocity cube is obtained through the paleo-porosity cube; the rock matrix velocity and relations between pore fluid velocities with depths are obtained through logging information; ancient overlaying formation pressure is obtained through a paleo-density cube; a formation pressure cube of the ancient geological time is obtained through the ancient overlaying formation pressure, the rock matrix velocity, the relations between pore fluid velocities with depths and the paleo-formation interval velocity.

The present invention relates to a seismic facies control-based sedimentary microfacies recognition method. The method includes the following steps that: 1) logging information is normalized; 2) fine horizon calibration is performed; 3) a model framework is determined; 4) a low frequency model is established; 5) a seismic data impedance value and a formation reflection coefficient are obtained, if synthetic trace data are in optimal matching with an actual seismic data trace, and drilling impedance is in optimal matching with seismic data impedance, and the method shifts to step 6), otherwise, the method shifts to the step 3); 6) well seismic constraint inversion is performed, a wave impedance inversion data body is obtained; 7) the top and bottom of a reservoir are traced and explained, with the top reflection layer and bottom reflection layer of the reservoir adopted as control surfaces, interlayer attributes are extracted from the inversion data body; and 8) seismic attributes are utilized to constrain the boundaries of sedimentary microfacies, and the spread of the sedimentary microfacies is described. According to the method, the logging information is adopted to constrain seismic inversion, and therefore, the resolution of an inversion result can be improved, the sedimentary microfacies of a target interval can be well depicted, the planar distribution of a thin sand layer can be predicted, and a basis can be provided for well location deployment.