CN108593771A - Damage strength computational methods and damage strength computing device - Google Patents
Damage strength computational methods and damage strength computing device Download PDFInfo
- Publication number
- CN108593771A CN108593771A CN201810392263.XA CN201810392263A CN108593771A CN 108593771 A CN108593771 A CN 108593771A CN 201810392263 A CN201810392263 A CN 201810392263A CN 108593771 A CN108593771 A CN 108593771A
- Authority
- CN
- China
- Prior art keywords
- shale
- damage
- coefficient
- wave
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/11—Analysing solids by measuring attenuation of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/12—Analysing solids by measuring frequency or resonance of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/015—Attenuation, scattering
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Damage strength computational methods provided by the invention and damage strength computing device are related to Shale Hydration Damage Evaluation technical field.Wherein, damage strength computational methods include：Obtain the first interval transit time, the first attenuation coefficient, the first peak swing and the first dominant frequency of shale；Obtain the second interval transit time, the second attenuation coefficient, the second peak swing and the second dominant frequency of the shale after hydration process；The damage strength after Shale Hydration is calculated according to the first interval transit time and the second interval transit time, the first attenuation coefficient and the second attenuation coefficient, the first peak swing and the second peak swing, the first dominant frequency and the second dominant frequency.By the above method, the relatively low problem of the damage strength reliability after the Shale Hydration being calculated in the prior art can be improved.
Description
Technical field
The present invention relates to Shale Hydration Damage Evaluation technical fields, in particular to a kind of damage strength computational methods
With damage strength computing device.
Background technology
Borehole well instability problem is always a weight difficulties of the industries such as oil, natural gas.Wherein, in Fissile Shale
During strata drilling, drilling fluid can cause aquation to damage in shale, borehole wall stability be reduced, to cause down hole problem
It takes place frequently, influences safety drilling.Numerous studies are carried out in drilling applications for the aquation Damage Evaluation of Fissile Shale, have examined
Consider drillng operation and be located at deep formation, be typically based on petrophysics property, formation rock information, evaluation are obtained using well logging information
Drilling fluid damages the aquation of shale formation.In numerous well logging informations, sound logging is most widely used.Sound wave is in rock mass
In communication process, extremely abundant acoustic intelligence is carried.Especially after hydration, rock interior structure, which changes, can lead to rock
Stone acoustic characteristic changes, therefore can be damaged to aquation according to rock acoustic property and carry out quantitative assessment.
Through inventor the study found that currently based on the aquation evaluation method of acoustic characteristic mainly using interval transit time as index,
Shale other acoustic characteristics are had ignored, to cannot all-sidedly and accurately embody Shale Hydration damage strength, cause to be calculated
The relatively low problem of damage strength reliability after Shale Hydration.
Invention content
In view of this, the purpose of the present invention is to provide a kind of damage strength computational methods and damage strength computing device,
To improve the relatively low problem of the damage strength reliability after the Shale Hydration being calculated in the prior art.
To achieve the above object, the embodiment of the present invention adopts the following technical scheme that：
A kind of damage strength computational methods, for evaluating the aquation damage of shale, the method includes：
The first interval transit time, the first attenuation coefficient, the first peak swing and the first dominant frequency of the shale are obtained,
In, first interval transit time, the first attenuation coefficient, the first peak swing and the first dominant frequency are based on surpassing the shale
Acoustic wave transmission processing generates；
The second interval transit time, the second attenuation coefficient, the second peak swing and the second dominant frequency of the shale are obtained,
In, second interval transit time, the second attenuation coefficient, the second peak swing and the second dominant frequency are based on to after hydration process
Shale carry out transmission ultrasonic wave processing generate；
According to first interval transit time and second interval transit time, first attenuation coefficient and second decaying
Coefficient, first peak swing and second peak swing, first dominant frequency and second dominant frequency are respectively according to pre-
If the first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th Damage coefficient is calculated in formula；
According to first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th Damage coefficient according to pre-
If the damage strength after the Shale Hydration is calculated in formula.
In the embodiment of the present invention preferably selects, in above-mentioned damage strength computational methods, first sound wave is obtained
The step of time difference includes：
The length of initial Mintrop wave take-off time, shale Mintrop wave take-off time and the shale is obtained respectively, wherein described
Initial Mintrop wave take-off time after docking the excitation probe and receiving transducer of transmission ultrasonic wave instrument based on carrying out at transmission ultrasonic wave
Reason generates, and the shale Mintrop wave take-off time is based on excitation probe and receiving transducer being respectively arranged at the shale length
Transmission ultrasonic wave processing is carried out behind both ends on direction to generate；
According to the length of the initial Mintrop wave take-off time, the shale Mintrop wave take-off time and the shale according to pre-
If first interval transit time is calculated in rule.
In the embodiment of the present invention preferably selects, in above-mentioned damage strength computational methods, first decaying is obtained
The step of coefficient includes：
The length of initial first wave amplitude, shale first wave amplitude and the shale is obtained respectively, wherein the initial Mintrop wave
Amplitude is generated based on transmission ultrasonic wave processing is carried out after docking the excitation probe and receiving transducer of transmission ultrasonic wave instrument, the page
Rock first wave amplitude is laggard based on the both ends that the excitation probe and receiving transducer are respectively arranged on the shale length direction
The processing of row transmission ultrasonic wave generates；
According to the length of the initial first wave amplitude, the shale first wave amplitude and the shale according to preset rules meter
Calculation obtains first attenuation coefficient.
In the embodiment of the present invention preferably selects, in above-mentioned damage strength computational methods, it is maximum to obtain described first
The step of amplitude and first dominant frequency includes：
It obtains based on carrying out the sound wave time-domain signal that transmission ultrasonic wave processing generates to the shale, and from the sound wave time domain
The first peak swing is obtained in signal；
The sound wave time-domain signal is subjected to Fourier transformation processing to obtain sound wave frequency-region signal, and from the sound wave frequency domain
The first dominant frequency is obtained in signal.
In the embodiment of the present invention preferably selects, in above-mentioned damage strength computational methods, first damage is calculated
Coefficient, the second Damage coefficient, third Damage coefficient and the 4th Damage coefficient preset formula include：
Wherein, St、Sa、Sr、SfRespectively the first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th damage
Hinder coefficient, Δ t (0) and Δ t (s) are respectively the first interval transit time and the second interval transit time, and At (0) and At (s) are respectively first
Attenuation coefficient and the second attenuation coefficient, r (0) and r (s) are respectively the first peak swing and the second peak swing, fre (0) and fre
(s) it is respectively the first dominant frequency and the second dominant frequency.
In the embodiment of the present invention preferably selects, in above-mentioned damage strength computational methods, the Shale Hydration is calculated
The preset formula of damage strength afterwards includes：
Wherein, S is the damage strength after the Shale Hydration, St、Sa、Sr、SfRespectively the first Damage coefficient, the second damage
Hinder coefficient, third Damage coefficient and the 4th Damage coefficient.
The embodiment of the present invention additionally provides a kind of damage strength computing device, for commenting the aquation damage of shale
Valence, described device include：
First parameter acquisition module, the first interval transit time, the first attenuation coefficient, the first maximum for obtaining the shale
Amplitude and the first dominant frequency, wherein first interval transit time, the first attenuation coefficient, the first peak swing and the first dominant frequency
It is generated based on transmission ultrasonic wave processing is carried out to the shale；
Second parameter acquisition module, the second interval transit time, the second attenuation coefficient, the second maximum for obtaining the shale
Amplitude and the second dominant frequency, wherein second interval transit time, the second attenuation coefficient, the second peak swing and the second dominant frequency
It is generated based on transmission ultrasonic wave processing is carried out to the shale after hydration process；
Damage coefficient computing module, for according to first interval transit time and second interval transit time, described first
Attenuation coefficient and second attenuation coefficient, first peak swing and second peak swing, first dominant frequency and
Second dominant frequency respectively according to preset formula be calculated the first Damage coefficient, the second Damage coefficient, third Damage coefficient with
And the 4th Damage coefficient；
Damage strength computing module, for according to first Damage coefficient, the second Damage coefficient, third Damage coefficient with
And the 4th Damage coefficient the damage strength after the Shale Hydration is calculated according to preset formula.
In the embodiment of the present invention preferably selects, in above-mentioned damage strength computing device, first parameter acquiring
Module includes：
Time parameter acquisition submodule, for obtain respectively initial Mintrop wave take-off time, shale Mintrop wave take-off time and
The length of the shale, wherein the initial Mintrop wave take-off time is visited based on popping one's head in and receiving the excitation of transmission ultrasonic wave instrument
It carries out transmission ultrasonic wave processing after head docking to generate, the shale Mintrop wave take-off time is based on that the excitation is popped one's head in and received to visit
Head carries out transmission ultrasonic wave processing generation after being respectively arranged at the both ends on the shale length direction；
Interval transit time computational submodule, for according to the initial Mintrop wave take-off time, the shale Mintrop wave take-off time
And first interval transit time is calculated according to preset rules in the length of the shale.
In the embodiment of the present invention preferably selects, in above-mentioned damage strength computing device, first parameter acquiring
Module further includes：
Range parameter acquisition submodule, for obtaining initial first wave amplitude, shale first wave amplitude and the shale respectively
Length, wherein the initial first wave amplitude after docking the excitation probe and receiving transducer of transmission ultrasonic wave instrument based on carrying out
Transmission ultrasonic wave processing generates, and the shale first wave amplitude is based on excitation probe and receiving transducer being respectively arranged at described
Transmission ultrasonic wave processing is carried out behind both ends on shale length direction to generate；
Attenuation coefficient computational submodule, for according to the initial first wave amplitude, the shale first wave amplitude and described
First attenuation coefficient is calculated according to preset rules in the length of shale.
In the embodiment of the present invention preferably selects, in above-mentioned damage strength computing device, first parameter acquiring
Module further includes：
Amplitude parameter acquisition submodule, for obtaining based on the sound wave for carrying out transmission ultrasonic wave processing generation to the shale
Time-domain signal, and obtain the first peak swing from the sound wave time-domain signal；
Dominant frequency parameter acquiring submodule, for the sound wave time-domain signal to be carried out Fourier transformation processing to obtain sound wave
Frequency-region signal, and obtain the first dominant frequency from the sound wave frequency-region signal.
Damage strength computational methods provided by the invention and damage strength computing device, by obtaining needs assessment respectively
Interval transit time, attenuation coefficient, peak swing and dominant frequency of the shale before and after hydration process, and based on each parameter to the damage of shale
Hinder intensity to be calculated, more comprehensively to evaluate the damage of the aquation of shale, and then improves and calculate in the prior art
The relatively low problem of damage strength reliability after the Shale Hydration arrived.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment cited below particularly, and coordinate
Appended attached drawing, is described in detail below.
Description of the drawings
Fig. 1 is the structure diagram of electronic equipment provided in an embodiment of the present invention.
Fig. 2 is the flow diagram of damage strength computational methods provided in an embodiment of the present invention.
Fig. 3 is the flow diagram of step S110 in Fig. 2.
Fig. 4 is another flow diagram of step S110 in Fig. 2.
Fig. 5 is another flow diagram of step S110 in Fig. 2.
Fig. 6 is the structure diagram of damage strength computing device provided in an embodiment of the present invention.
Fig. 7 is the structure diagram of the first parameter acquisition module provided in an embodiment of the present invention.
Fig. 8 is another structure diagram of the first parameter acquisition module provided in an embodiment of the present invention.
Fig. 9 is another structure diagram of the first parameter acquisition module provided in an embodiment of the present invention.
Icon：10- electronic equipments；12- memories；14- processors；100- damage strength computing devices；110- first joins
Number acquisition module；111- time parameter acquisition submodules；112- interval transit time computational submodules；113- range parameters obtain submodule
Block；114- attenuation coefficient computational submodules；115- amplitude parameter acquisition submodules；116- dominant frequency parameter acquiring submodules；130-
Second parameter acquisition module；150- Damage coefficient computing modules；170- damage strength computing modules.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment only
It is a part of the embodiment of the present invention, instead of all the embodiments.The present invention being usually described and illustrated herein in the accompanying drawings
The component of embodiment can be arranged and be designed with a variety of different configurations.
Therefore, below the detailed description of the embodiment of the present invention to providing in the accompanying drawings be not intended to limit it is claimed
The scope of the present invention, but be merely representative of the present invention selected embodiment.Based on the embodiments of the present invention, this field is common
The every other embodiment that technical staff is obtained without creative efforts belongs to the model that the present invention protects
It encloses.
It should be noted that：Similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined, then it further need not be defined and explained in subsequent attached drawing in a attached drawing.In description of the invention
In, term " first ", " second ", " third ", " the 4th " etc. are only used for distinguishing description, and should not be understood as only or imply opposite
Importance.
As shown in Figure 1, can be used for damaging the electronics evaluated to the aquation of shale an embodiment of the present invention provides a kind of
Equipment 10.Wherein, the electronic equipment 10 may include memory 12, processor 14 and damage strength computing device 100.
It is directly or indirectly electrically connected between the memory 12 and processor 14, to realize the transmission or friendship of data
Mutually.It is electrically connected for example, these elements can be realized between each other by one or more communication bus or signal wire.The damage
Strength co-mputation device 100 can be stored in the form of software or firmware (firmware) in the memory 12 including at least one
Software function module.The processor 14 is for executing the executable computer program stored in the memory 12, example
Such as, the software function module included by the damage strength computing device 100 and computer program etc., to realize in terms of damage strength
Calculation takes method.
Wherein, the memory 12 may be, but not limited to, random access memory (Random Access Memory,
RAM), read-only memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only
Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM),
Electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc..
Wherein, memory 12 is for storing program, and the processor 14 executes described program after receiving and executing instruction.
The processor 14 may be a kind of IC chip, the processing capacity with signal.Above-mentioned processor 14
Can be general processor, including central processing unit (Central Processing Unit, CPU), network processing unit
(Network Processor, NP) etc.；It can also be digital signal processor (DSP), application-specific integrated circuit (ASIC), scene
Programmable gate array (FPGA) either other programmable logic device, discrete gate or transistor logic, discrete hardware group
Part.It may be implemented or execute disclosed each method, step and the logic diagram in the embodiment of the present invention.General processor can be with
It is microprocessor or the processor can also be any conventional processor etc..
It is appreciated that structure shown in FIG. 1 is only to illustrate, the electronic equipment 10 may also include more than shown in Fig. 1
Either less component or with the configuration different from shown in Fig. 1.Hardware, software may be used in each component shown in Fig. 1
Or combinations thereof realize.
Optionally, the concrete type of the electronic equipment 10 is unrestricted, for example, it may be, but be not limited to, intelligent hand
Machine, PC (personal computer, PC), tablet computer, personal digital assistant (personal digital
Assistant, PDA), mobile internet surfing equipment (mobile Internet device, MID), web (website) server, data
The equipment that server, computer, mobile internet surfing equipment (mobile Internet device, MID) etc. have processing function.
In conjunction with Fig. 2, the embodiment of the present invention additionally provides a kind of damage strength calculating can be applied to above-mentioned electronic equipment 10
Method.Wherein, method and step defined in the related flow of the method can be realized by the processor 14.It below will be to figure
Detailed process shown in 2 is described in detail.
Step S110 obtains the first interval transit time of the shale, the first attenuation coefficient, the first peak swing and first
Dominant frequency.
In the present embodiment, first interval transit time, the first attenuation coefficient, the first peak swing and the first dominant frequency base
It is generated in carrying out transmission ultrasonic wave processing to the shale.
Wherein, to the shale carry out transmission ultrasonic wave processing before, for ensure obtain each parameter validity and
Reliability can be handled the shale of extraction, it is, for example, possible to use air bit drill through processing to obtain the shale
Rock core, and the both ends of the rock core are cut, so that two end faces are vertical with rock core axis.
Step S130 obtains the second interval transit time of the shale, the second attenuation coefficient, the second peak swing and second
Dominant frequency.
In the present embodiment, second interval transit time, the second attenuation coefficient, the second peak swing and the second dominant frequency base
It is generated in carrying out transmission ultrasonic wave processing to the shale after hydration process.
Wherein, the mode that hydration process is carried out to shale is unrestricted, for example, either passing through arbitrary liquid (example
Such as, water) it carries out, can also be to be carried out by drilling fluid.It in the present embodiment, can be more true anti-for guarantee hydration process
The aquation effect in practical application is reflected, drilling fluid may be used.
Step S150, according to first interval transit time and second interval transit time, first attenuation coefficient and institute
State the second attenuation coefficient, first peak swing and second peak swing, first dominant frequency and second dominant frequency
The first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th damage, which is calculated, according to preset formula respectively is
Number.
Wherein, first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th Damage coefficient are calculated
Preset formula is unrestricted, can be configured according to practical application request, for example, it is directed to interval transit time, it can be directly by the
The difference of one interval transit time and second interval transit time is as the first Damage coefficient.In the present embodiment, step is executed for guarantee
May include with higher consistency, the preset formula when being handled by each coefficient when rapid S170：
Wherein, St、Sa、Sr、SfRespectively the first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th damage
Hinder coefficient, Δ t (0) and Δ t (s) are respectively the first interval transit time and the second interval transit time, and At (0) and At (s) are respectively first
Attenuation coefficient and the second attenuation coefficient, r (0) and r (s) are respectively the first peak swing and the second peak swing, fre (0) and fre
(s) it is respectively the first dominant frequency and the second dominant frequency.
Step S170, according to first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th damage system
The damage strength after the Shale Hydration is calculated according to preset formula for number.
Wherein, the preset formula for calculating the damage strength is unrestricted, can be configured according to practical application request,
For example, it may be taking intermediate value to four coefficients, it can also be that mean value is taken to four coefficients, can also be other calculations.
In the present embodiment, the preset formula for calculating the damage strength may include：
Wherein, S is the damage strength after the Shale Hydration, St、Sa、Sr、SfRespectively the first Damage coefficient, the second damage
Hinder coefficient, third Damage coefficient and the 4th Damage coefficient.
Further, in the present embodiment, by step S110 to obtain the first interval transit time, the first attenuation coefficient,
The mode of one peak swing and the first dominant frequency is unrestricted, can be configured according to practical application request, in conjunction with Fig. 3, step
Rapid S110 may include step S111 and step S112 to obtain first interval transit time.
Step S111 obtains the length of initial Mintrop wave take-off time, shale Mintrop wave take-off time and the shale respectively.
In the present embodiment, the initial Mintrop wave take-off time is visited based on popping one's head in and receiving the excitation of transmission ultrasonic wave instrument
It carries out transmission ultrasonic wave processing after head docking to generate, the shale Mintrop wave take-off time is based on that the excitation is popped one's head in and received to visit
Head carries out transmission ultrasonic wave processing generation after being respectively arranged at the both ends on the shale length direction.
Step S112, according to the initial Mintrop wave take-off time, the shale Mintrop wave take-off time and the shale
First interval transit time is calculated according to preset rules in length.
In the present embodiment, the preset rules of calculating first interval transit time can be：
Wherein, Δ t is the first interval transit time, t1For initial Mintrop wave take-off time, t2For shale Mintrop wave take-off time, L is page
The length of rock.
In conjunction with Fig. 4, in the present embodiment, step S110 can also include step S113 and step S114, described to obtain
First attenuation coefficient.
Step S113 obtains the length of initial first wave amplitude, shale first wave amplitude and the shale respectively.
In the present embodiment, the initial first wave amplitude is based on by the excitation probe and receiving transducer pair of transmission ultrasonic wave instrument
It carries out transmission ultrasonic wave processing after connecing to generate, the shale first wave amplitude is based on respectively setting excitation probe and receiving transducer
Transmission ultrasonic wave processing generation is carried out after being placed in the both ends on the shale length direction.
Step S114, according to the length of the initial first wave amplitude, the shale first wave amplitude and the shale according to
First attenuation coefficient is calculated in preset rules.
In the present embodiment, the preset rules of calculating first attenuation coefficient can be：
Wherein, At is the first attenuation coefficient, A0For initial first wave amplitude, A is shale first wave amplitude, and L is the length of shale.
In conjunction with Fig. 5, in the present embodiment, step S110 can also include step S115 and step S116, described to obtain
First peak swing and first dominant frequency.
Step S115 is obtained based on the sound wave time-domain signal for carrying out transmission ultrasonic wave processing generation to the shale, and from
The first peak swing is obtained in the sound wave time-domain signal.
The sound wave time-domain signal is carried out Fourier transformation processing to obtain sound wave frequency-region signal by step S116, and from
The first dominant frequency is obtained in the sound wave frequency-region signal.
In the present embodiment, the time-domain signal of fluctuation can be unfolded by frequency order by Fourier transformation, to obtain
Sound wave frequency-region signal, and obtain the first dominant frequency from the sound wave frequency-region signal.Wherein, the principle of Fourier transformation can be：
Wherein, f is frequency, and x (t) is time-domain function, and X (f) is the complex function of sound wave frequency domain, and t is one in Acoustic Wave Propagation
A time point.
Wherein, step S130 is being executed to obtain the second interval transit time, the second attenuation coefficient, the second peak swing and the
The mode of two dominant frequency is unrestricted, can be configured according to practical application request, for example, either side with step S110
Formula is identical, can also be different from the mode of step S110.In the present embodiment, for ensure by first interval transit time with
Second interval transit time, first attenuation coefficient and second attenuation coefficient, first peak swing and described
Each Damage coefficient that two peak swings, first dominant frequency and second dominant frequency are calculated has higher reliability, step
Rapid S130 is identical as the acquisition modes for parameter in step S110, specifically, incorporated by reference to including to step S110 above
The explanation of each step.
Further, to ensure that the evaluation of the aquation damage to shale has higher reliability, in the present embodiment, also
It can be based on the different hydration times to same shale, to carry out the calculating of damage strength respectively, to obtain the damage of the shale
One correspondence of intensity and hydration time.
Also, the evaluation that aquation damage in the present embodiment, can also be carried out to different shale, with by will be different
Shale damage strength is compared, to obtain more structurally sound evaluation.
In conjunction with Fig. 6, the embodiment of the present invention also provides a kind of damage strength calculating dress can be applied to above-mentioned electronic equipment 10
Set 100.Wherein, the damage strength computing device 100 may include the first parameter acquisition module 110, the second parameter acquiring mould
Block 130, Damage coefficient computing module 150 and damage strength computing module 170.
First parameter acquisition module 110, for obtain the first interval transit time of the shale, the first attenuation coefficient,
First peak swing and the first dominant frequency, wherein first interval transit time, the first attenuation coefficient, the first peak swing and
First dominant frequency is based on carrying out transmission ultrasonic wave processing generation to the shale.In the present embodiment, the first parameter acquiring mould
Block 110 can be used for executing step S110 shown in Fig. 2, and the specific descriptions about first parameter acquisition module 110 can join
According to the description above to step S110.
Second parameter acquisition module 130, for obtain the second interval transit time of the shale, the second attenuation coefficient,
Second peak swing and the second dominant frequency, wherein second interval transit time, the second attenuation coefficient, the second peak swing and
Second dominant frequency is based on to the shale progress transmission ultrasonic wave processing generation after hydration process.In the present embodiment, described
Two parameter acquisition modules 130 can be used for executing step S130 shown in Fig. 2, the tool about second parameter acquisition module 130
Body description is referred to the description to step S130 above.
The Damage coefficient computing module 150, for according to first interval transit time and second interval transit time, institute
State the first attenuation coefficient and second attenuation coefficient, first peak swing and second peak swing, described first
Dominant frequency and second dominant frequency are calculated the first Damage coefficient, the second Damage coefficient, third according to preset formula respectively and damage
Coefficient and the 4th Damage coefficient.In the present embodiment, the Damage coefficient computing module 150 can be used for executing shown in Fig. 2
Step S150, the specific descriptions about the Damage coefficient computing module 150 are referred to the description to step S150 above.
The damage strength computing module 170, for being damaged according to first Damage coefficient, the second Damage coefficient, third
Hinder coefficient and the damage strength after the Shale Hydration is calculated according to preset formula for the 4th Damage coefficient.In the present embodiment
In, the damage strength computing module 170 can be used for executing step S170 shown in Fig. 2, and mould is calculated about the damage strength
The specific descriptions of block 170 are referred to the description to step S170 above.
In conjunction with Fig. 7, in the present embodiment, first parameter acquisition module 110 may include that time parameter obtains submodule
Block 111 and interval transit time computational submodule 112.
The time parameter acquisition submodule 111, when for obtaining initial Mintrop wave take-off time, shale Mintrop wave take-off respectively
Between and the shale length, wherein the initial Mintrop wave take-off time based on by the excitation of transmission ultrasonic wave instrument probe and
Receiving transducer docking after carry out transmission ultrasonic wave processing generate, the shale Mintrop wave take-off time be based on by the excitation probe with
Receiving transducer carries out transmission ultrasonic wave processing generation after being respectively arranged at the both ends on the shale length direction.In the present embodiment
In, the time parameter acquisition submodule 111 can be used for executing step S111 shown in Fig. 3, be obtained about the time parameter
The specific descriptions of submodule 111 are referred to the description to step S111 above.
The interval transit time computational submodule 112, for being risen according to the initial Mintrop wave take-off time, the shale Mintrop wave
First interval transit time is calculated according to preset rules in the length for jumping time and the shale.In the present embodiment, institute
It states interval transit time computational submodule 112 and can be used for executing step S112 shown in Fig. 3, about the interval transit time computational submodule
112 specific descriptions are referred to the description to step S112 above.
In conjunction with Fig. 8, in the present embodiment, first parameter acquisition module 110 can also include that range parameter obtains son
Module 113 and attenuation coefficient computational submodule 114.
The range parameter acquisition submodule 113, for obtaining initial first wave amplitude, shale first wave amplitude and institute respectively
State the length of shale, wherein the initial first wave amplitude is based on docking the excitation probe and receiving transducer of transmission ultrasonic wave instrument
It carries out transmission ultrasonic wave processing afterwards to generate, the shale first wave amplitude is based on excitation probe and receiving transducer being respectively set
Transmission ultrasonic wave processing is carried out behind the both ends on the shale length direction to generate.In the present embodiment, the range parameter
Acquisition submodule 113 can be used for executing step S113 shown in Fig. 4, about the specific of the range parameter acquisition submodule 113
Description is referred to the description to step S113 above.
The attenuation coefficient computational submodule 114, for according to the initial first wave amplitude, the shale first wave amplitude with
And first attenuation coefficient is calculated according to preset rules in the length of the shale.In the present embodiment, the decaying system
Number computational submodule 114 can be used for executing step S114 shown in Fig. 4, the tool about the attenuation coefficient computational submodule 114
Body description is referred to the description to step S114 above.
In conjunction with Fig. 9, in the present embodiment, first parameter acquisition module 110 can also include that amplitude parameter obtains son
Module 115 and dominant frequency parameter acquiring submodule 116.
The amplitude parameter acquisition submodule 115, for obtaining based on to shale progress transmission ultrasonic wave processing life
At sound wave time-domain signal, and obtain the first peak swing from the sound wave time-domain signal.In the present embodiment, the amplitude ginseng
Number acquisition submodule 115 can be used for executing step S115 shown in fig. 5, the tool about the amplitude parameter acquisition submodule 115
Body description is referred to the description to step S115 above.
The dominant frequency parameter acquiring submodule 116, for by the sound wave time-domain signal carry out Fourier transformation processing with
Sound wave frequency-region signal is obtained, and the first dominant frequency is obtained from the sound wave frequency-region signal.In the present embodiment, the dominant frequency parameter obtains
Submodule 116 is taken to can be used for executing step S116 shown in fig. 5, specific about the dominant frequency parameter acquiring submodule 116 is retouched
It states and is referred to the description to step S116 above.
In conclusion damage strength computational methods provided by the invention and damage strength computing device 100, by obtaining respectively
Interval transit time, attenuation coefficient, peak swing and dominant frequency of the shale of needs assessment before and after hydration process are taken, and is based on each ginseng
Several damage strengths to shale calculate, and more comprehensively to evaluate the damage of the aquation of shale, and then improve existing
The relatively low problem of damage strength reliability after the Shale Hydration being calculated in technology.
In several embodiments that the embodiment of the present invention is provided, it should be understood that disclosed device and method also may be used
To realize by another way.Device and method embodiment described above is only schematical, for example, in attached drawing
Flow chart and block diagram show the device of multiple embodiments according to the present invention, the possibility of method and computer program product is realized
Architecture, function and operation.In this regard, each box in flowchart or block diagram can represent module, a program
A part for a part for section or code, the module, section or code includes that one or more is patrolled for realizing defined
Collect the executable instruction of function.It should also be noted that at some as the function of in the realization method replaced, being marked in box
It can occur in a different order than that indicated in the drawings.For example, two continuous boxes can essentially be held substantially in parallel
Row, they can also be executed in the opposite order sometimes, this is depended on the functions involved.It is also noted that block diagram and/or
The combination of each box in flow chart and the box in block diagram and or flow chart can use function or dynamic as defined in executing
The dedicated hardware based system made is realized, or can be realized using a combination of dedicated hardware and computer instructions.
In addition, each function module in each embodiment of the present invention can integrate to form an independent portion
Point, can also be modules individualism, can also two or more modules be integrated to form an independent part.
It, can be with if the function is realized and when sold or used as an independent product in the form of software function module
It is stored in a computer read/write memory medium.Based on this understanding, technical scheme of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be expressed in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be
People's computer, electronic equipment or network equipment etc.) execute all or part of step of each embodiment the method for the present invention
Suddenly.And storage medium above-mentioned includes：USB flash disk, read-only memory (ROM, Read-Only Memory), is deposited mobile hard disk at random
The various media that can store program code such as access to memory (RAM, Random Access Memory), magnetic disc or CD.
It should be noted that herein, the terms "include", "comprise" or its any other variant are intended to the packet of nonexcludability
Contain, so that the process, method, article or equipment including a series of elements includes not only those elements, but also includes
Other elements that are not explicitly listed, or further include for elements inherent to such a process, method, article, or device.
In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including the element
Process, method, article or equipment in there is also other identical elements.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, any made by repair
Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.
Claims (10)
1. a kind of damage strength computational methods, for evaluating the aquation damage of shale, which is characterized in that the method packet
It includes：
Obtain the first interval transit time, the first attenuation coefficient, the first peak swing and the first dominant frequency of the shale, wherein institute
The first interval transit time, the first attenuation coefficient, the first peak swing and the first dominant frequency is stated to be based on carrying out ultrasonic wave to the shale
Transmission processing generates；
Obtain the second interval transit time, the second attenuation coefficient, the second peak swing and the second dominant frequency of the shale, wherein institute
The second interval transit time, the second attenuation coefficient, the second peak swing and the second dominant frequency are stated based on to the page after hydration process
Rock carries out transmission ultrasonic wave processing and generates；
According to first interval transit time and second interval transit time, first attenuation coefficient and second decaying system
Several, described first peak swing and second peak swing, first dominant frequency and second dominant frequency are respectively according to default
The first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th Damage coefficient is calculated in formula；
According to first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th Damage coefficient according to default public affairs
The damage strength after the Shale Hydration is calculated in formula.
2. damage strength computational methods according to claim 1, which is characterized in that obtain the step of first interval transit time
Suddenly include：
The length of initial Mintrop wave take-off time, shale Mintrop wave take-off time and the shale is obtained respectively, wherein described initial
Mintrop wave take-off time is given birth to based on transmission ultrasonic wave processing is carried out after docking the excitation probe and receiving transducer of transmission ultrasonic wave instrument
At the shale Mintrop wave take-off time is based on excitation probe and receiving transducer being respectively arranged at the shale length direction
On both ends after carry out transmission ultrasonic wave processing generate；
According to the length of the initial Mintrop wave take-off time, the shale Mintrop wave take-off time and the shale according to default rule
First interval transit time is then calculated.
3. damage strength computational methods according to claim 1, which is characterized in that obtain the step of first attenuation coefficient
Suddenly include：
The length of initial first wave amplitude, shale first wave amplitude and the shale is obtained respectively, wherein the initial first wave amplitude
It is generated based on transmission ultrasonic wave processing is carried out after docking the excitation probe and receiving transducer of transmission ultrasonic wave instrument, the shale is first
Wave amplitude behind the excitation probe and the both ends that are respectively arranged on the shale length direction of receiving transducer based on will surpass
Acoustic wave transmission processing generates；
It is calculated according to preset rules according to the length of the initial first wave amplitude, the shale first wave amplitude and the shale
To first attenuation coefficient.
4. damage strength computational methods according to claim 1, which is characterized in that obtain first peak swing and institute
The step of stating the first dominant frequency include：
It obtains based on carrying out the sound wave time-domain signal that transmission ultrasonic wave processing generates to the shale, and from the sound wave time-domain signal
In obtain the first peak swing；
The sound wave time-domain signal is subjected to Fourier transformation processing to obtain sound wave frequency-region signal, and from the sound wave frequency-region signal
In obtain the first dominant frequency.
5. damage strength computational methods according to any one of claims 1-4, which is characterized in that calculate first damage
The preset formula for hindering coefficient, the second Damage coefficient, third Damage coefficient and the 4th Damage coefficient includes：
Wherein, St、Sa、Sr、SfRespectively the first Damage coefficient, the second Damage coefficient, third Damage coefficient and the 4th damage system
Number, Δ t (0) and Δ t (s) are respectively the first interval transit time and the second interval transit time, and At (0) and At (s) are respectively the first decaying
Coefficient and the second attenuation coefficient, r (0) and r (s) are respectively the first peak swing and the second peak swing, fre (0) and fre (s)
Respectively the first dominant frequency and the second dominant frequency.
6. damage strength computational methods according to any one of claims 1-4, which is characterized in that calculate the shale water
The preset formula of damage strength after change includes：
Wherein, S is the damage strength after the Shale Hydration, St、Sa、Sr、SfRespectively the first Damage coefficient, the second damage system
Number, third Damage coefficient and the 4th Damage coefficient.
7. a kind of damage strength computing device, for evaluating the aquation damage of shale, which is characterized in that described device packet
It includes：
First parameter acquisition module, for obtaining the first interval transit time of the shale, the first attenuation coefficient, the first peak swing
And first dominant frequency, wherein first interval transit time, the first attenuation coefficient, the first peak swing and the first dominant frequency are based on
Transmission ultrasonic wave processing is carried out to the shale to generate；
Second parameter acquisition module, for obtaining the second interval transit time of the shale, the second attenuation coefficient, the second peak swing
And second dominant frequency, wherein second interval transit time, the second attenuation coefficient, the second peak swing and the second dominant frequency are based on
Transmission ultrasonic wave processing is carried out to the shale after hydration process to generate；
Damage coefficient computing module, for according to first interval transit time and second interval transit time, first decaying
Coefficient and second attenuation coefficient, first peak swing and second peak swing, first dominant frequency and described
The first Damage coefficient, the second Damage coefficient, third Damage coefficient and is calculated according to preset formula respectively in second dominant frequency
Four Damage coefficients；
Damage strength computing module, for according to first Damage coefficient, the second Damage coefficient, third Damage coefficient and the
The damage strength after the Shale Hydration is calculated according to preset formula for four Damage coefficients.
8. damage strength computing device according to claim 7, which is characterized in that the first parameter acquisition module packet
It includes：
Time parameter acquisition submodule, for obtaining initial Mintrop wave take-off time, shale Mintrop wave take-off time and described respectively
The length of shale, wherein the initial Mintrop wave take-off time is based on by the excitation probe and receiving transducer pair of transmission ultrasonic wave instrument
It carries out transmission ultrasonic wave processing after connecing to generate, the shale Mintrop wave take-off time is based on excitation probe and receiving transducer point
Transmission ultrasonic wave processing generation is carried out after not being set to the both ends on the shale length direction；
Interval transit time computational submodule, for according to the initial Mintrop wave take-off time, the shale Mintrop wave take-off time and
First interval transit time is calculated according to preset rules in the length of the shale.
9. damage strength computing device according to claim 7, which is characterized in that first parameter acquisition module is also wrapped
It includes：
Range parameter acquisition submodule, the length for obtaining initial first wave amplitude, shale first wave amplitude and the shale respectively
Degree, wherein the initial first wave amplitude is ultrasonic based on being carried out after docking the excitation probe and receiving transducer of transmission ultrasonic wave instrument
Wave transmission processing generates, and the shale first wave amplitude is based on excitation probe and receiving transducer being respectively arranged at the shale
Transmission ultrasonic wave processing is carried out behind both ends on length direction to generate；
Attenuation coefficient computational submodule, for according to the initial first wave amplitude, the shale first wave amplitude and the shale
Length first attenuation coefficient is calculated according to preset rules.
10. damage strength computing device according to claim 7, which is characterized in that first parameter acquisition module is also
Including：
Amplitude parameter acquisition submodule, for obtaining based on the sound wave time domain for carrying out transmission ultrasonic wave processing generation to the shale
Signal, and obtain the first peak swing from the sound wave time-domain signal；
Dominant frequency parameter acquiring submodule, for the sound wave time-domain signal to be carried out Fourier transformation processing to obtain sound wave frequency domain
Signal, and obtain the first dominant frequency from the sound wave frequency-region signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810392263.XA CN108593771A (en) | 2018-04-27 | 2018-04-27 | Damage strength computational methods and damage strength computing device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810392263.XA CN108593771A (en) | 2018-04-27 | 2018-04-27 | Damage strength computational methods and damage strength computing device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108593771A true CN108593771A (en) | 2018-09-28 |
Family
ID=63610762
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810392263.XA Pending CN108593771A (en) | 2018-04-27 | 2018-04-27 | Damage strength computational methods and damage strength computing device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108593771A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109541688A (en) * | 2018-11-30 | 2019-03-29 | 中铁第四勘察设计院集团有限公司 | A kind of structural plane flaw evaluation method based on energy feature |
| CN109596713A (en) * | 2018-12-29 | 2019-04-09 | 西南石油大学 | A kind of coal petrography ultrasonic wave attenuation coefficient evaluation method based on Fractal |
| CN109828031A (en) * | 2019-02-15 | 2019-05-31 | 西南石油大学 | Rock brittleness evaluation method and device |
| CN109883920A (en) * | 2019-03-08 | 2019-06-14 | 西南石油大学 | Characterize the method and device of rock thermal damage |
| CN111855484A (en) * | 2020-07-30 | 2020-10-30 | 西南石油大学 | Method for evaluating well wall capability of drilling fluid for stabilizing shale formation based on acoustoelectric response |
| CN113624847A (en) * | 2021-08-12 | 2021-11-09 | 西南石油大学 | Method for establishing prediction model of shale hydration damage coefficient and prediction method |
| CN113958315A (en) * | 2021-06-16 | 2022-01-21 | 西南石油大学 | Shale stratum collapse pressure prediction method based on self-absorption-constitutive model |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102183585A (en) * | 2011-03-09 | 2011-09-14 | 西南石油大学 | Rock core sampling method |
| CN102519784A (en) * | 2011-12-16 | 2012-06-27 | 武汉大学 | Method for determining rock conjugate damage strength through adopting supersonic waves |
| CN202305478U (en) * | 2011-08-10 | 2012-07-04 | 建维科技（深圳）有限公司 | Sensor and concrete structure hydration and damage monitoring system |
| CN102749386A (en) * | 2011-04-19 | 2012-10-24 | 香港科技大学 | System and method for in-situ hydration monitoring and damage detection of concrete structure and sensors used by system and method |
| CN104502452A (en) * | 2014-12-17 | 2015-04-08 | 江苏大学 | Structural damage assessment method based on ultrasonic guided-wave drive/sensing array |
| CN104675395A (en) * | 2015-02-12 | 2015-06-03 | 中国石油大学(北京) | Evaluation method for hydration characteristics of layered hard brittle mud shale |
| CN105863626A (en) * | 2016-05-10 | 2016-08-17 | 西南石油大学 | Evaluation method for physical and chemical action of drilling fluid and shale formation |
-
2018
- 2018-04-27 CN CN201810392263.XA patent/CN108593771A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102183585A (en) * | 2011-03-09 | 2011-09-14 | 西南石油大学 | Rock core sampling method |
| CN102749386A (en) * | 2011-04-19 | 2012-10-24 | 香港科技大学 | System and method for in-situ hydration monitoring and damage detection of concrete structure and sensors used by system and method |
| CN202305478U (en) * | 2011-08-10 | 2012-07-04 | 建维科技（深圳）有限公司 | Sensor and concrete structure hydration and damage monitoring system |
| CN102519784A (en) * | 2011-12-16 | 2012-06-27 | 武汉大学 | Method for determining rock conjugate damage strength through adopting supersonic waves |
| CN104502452A (en) * | 2014-12-17 | 2015-04-08 | 江苏大学 | Structural damage assessment method based on ultrasonic guided-wave drive/sensing array |
| CN104675395A (en) * | 2015-02-12 | 2015-06-03 | 中国石油大学(北京) | Evaluation method for hydration characteristics of layered hard brittle mud shale |
| CN105863626A (en) * | 2016-05-10 | 2016-08-17 | 西南石油大学 | Evaluation method for physical and chemical action of drilling fluid and shale formation |
Non-Patent Citations (1)
| Title |
|---|
| 王光兵 等: "硬脆性页岩水化的超声波透射实验研究", 《科学技术与工程》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109541688A (en) * | 2018-11-30 | 2019-03-29 | 中铁第四勘察设计院集团有限公司 | A kind of structural plane flaw evaluation method based on energy feature |
| CN109596713A (en) * | 2018-12-29 | 2019-04-09 | 西南石油大学 | A kind of coal petrography ultrasonic wave attenuation coefficient evaluation method based on Fractal |
| CN109828031A (en) * | 2019-02-15 | 2019-05-31 | 西南石油大学 | Rock brittleness evaluation method and device |
| CN109883920A (en) * | 2019-03-08 | 2019-06-14 | 西南石油大学 | Characterize the method and device of rock thermal damage |
| CN111855484A (en) * | 2020-07-30 | 2020-10-30 | 西南石油大学 | Method for evaluating well wall capability of drilling fluid for stabilizing shale formation based on acoustoelectric response |
| CN111855484B (en) * | 2020-07-30 | 2022-05-20 | 西南石油大学 | Method for evaluating drilling fluid stable shale stratum well wall capability based on acoustoelectric response |
| CN113958315A (en) * | 2021-06-16 | 2022-01-21 | 西南石油大学 | Shale stratum collapse pressure prediction method based on self-absorption-constitutive model |
| CN113958315B (en) * | 2021-06-16 | 2022-05-17 | 西南石油大学 | Shale stratum collapse pressure prediction method based on self-absorption-constitutive model |
| CN113624847A (en) * | 2021-08-12 | 2021-11-09 | 西南石油大学 | Method for establishing prediction model of shale hydration damage coefficient and prediction method |
| CN113624847B (en) * | 2021-08-12 | 2022-07-08 | 西南石油大学 | Method for establishing prediction model of shale hydration damage coefficient and prediction method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108593771A (en) | Damage strength computational methods and damage strength computing device | |
| US11098565B2 (en) | Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves | |
| US10995606B2 (en) | Well integrity analysis using sonic measurements over depth interval | |
| US11795817B2 (en) | System and method for determining formation characteristics using electrical arc modeling | |
| US9322807B2 (en) | Ultrasonic signal time-frequency decomposition for borehole evaluation or pipeline inspection | |
| CN105891089B (en) | Reservoir permeability determines method and device | |
| US20130204534A1 (en) | Method Of Estimating A Subterranean Formation Property | |
| CN110736400B (en) | Underwater drilling blasting vibration velocity calculation method considering internal rock structure | |
| CN106934183B (en) | Dispersion curve determines that method and apparatus and p-and s-wave velocity determine method and apparatus | |
| US20150260872A1 (en) | Method and system for realtime determination of formation anisotropy, dip, and strike with mci data | |
| CN110426739B (en) | Geological exploration detection method and device and storage medium | |
| CN107450108B (en) | The determination method and apparatus in dessert area | |
| CN108491646A (en) | Method for identifying vibration damage caused by blasting when tunnel passes through important building in short distance | |
| CN103852785B (en) | The evaluation method of formation anisotropy | |
| CN106501872B (en) | Method and device for calculating ground stress characteristics of fractured reservoir | |
| CN115992691A (en) | Well cementation quality detection method and device based on ultrasonic lamb waves | |
| CA2895025C (en) | Sourceless density determination apparatus, methods, and systems | |
| CN117270053B (en) | Stratum transverse wave slowness and attenuation calculation method based on dipole acoustic logging | |
| CN109828031A (en) | Rock brittleness evaluation method and device | |
| Godbeer et al. | Environment-induced dephasing versus von Neumann measurements in proton tunneling | |
| CN116556925A (en) | Method, device, computing equipment and storage medium for computing borehole ovality | |
| Movassagh et al. | A comparison of fractal methods for evaluation of hydraulic fracturing surface roughness | |
| DE112016000974T5 (en) | Method of assessing cement bonding | |
| US10295696B2 (en) | Multi-component induction logging data processing in non-circular boreholes | |
| CN111060964B (en) | Method and device for determining stratum elasticity parameters |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180928 |
|
| RJ01 | Rejection of invention patent application after publication |