CN107587871A - Determine the method and device of horizontal fracture width - Google Patents
Determine the method and device of horizontal fracture width Download PDFInfo
- Publication number
- CN107587871A CN107587871A CN201710666181.5A CN201710666181A CN107587871A CN 107587871 A CN107587871 A CN 107587871A CN 201710666181 A CN201710666181 A CN 201710666181A CN 107587871 A CN107587871 A CN 107587871A
- Authority
- CN
- China
- Prior art keywords
- horizontal fracture
- stoneley wave
- fracture
- amplitude
- horizontal
- 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.)
- Granted
Links
Abstract
The embodiment of the present application provides a kind of method and device for determining horizontal fracture width, and this method includes:The numerical model of pending well is established, and simulates acoustic logging waveform of the numerical model under each bar horizontal fracture;According to the acoustic logging waveform under each bar horizontal fracture, the first Stoneley wave changes in amplitude amount of this horizontal fracture is determined;According to the first Stoneley wave changes in amplitude amount of each bar horizontal fracture, the relation curve between fracture width and Stoneley wave changes in amplitude amount is determined;According to the acoustic logging waveform of the pending well, the second Stoneley wave changes in amplitude amount of every horizontal fracture in the pending well is obtained;According to the second Stoneley wave changes in amplitude amount of every horizontal fracture in the relation curve and the pending well, the width of every horizontal fracture in the pending well is determined.The embodiment of the present application can improve the standard the degree of accuracy of calculatingcrackswidth.
Description
Technical field
The application is related to fracture-type reservoir oil-gas exploration technical field, more particularly, to a kind of determination horizontal fracture width
Method and device.
Background technology
Fracture-type reservoir is an important directions of current oil-gas exploration.In fracture-type reservoir, fracture width, crack length
Degree and fracture spacing are important reservoir evaluation parameter, and the quantitative calculating especially for fracture width is particularly important, because
The calculating of other fracture parameters will be directly affected for the degree of accuracy of its result of calculation.
Being currently used in the method for calculating fracture width mainly includes acoustic logging method and resistivity imaging method.Cement bond logging
Well refer to arrange sound source in the deep liquid-filled well bore of several kms and produce wellhole in various modes ripple, then utilize various patterns
The acoustic information of ripple evaluates a kind of technology of formation properties by well.Wherein monopole source emission, the ripple that monopole receiver receives
Mainly including these four mode waves of compressional wave, shear wave, pseudo-Rayleigh wave and Stoneley wave in shape.After low-pass filtered, just it can obtain pure
Net Stoneley wave signal.There is wellhole acoustic logging waveform during horizontal fracture in nineteen ninety Tang Xiaoming theory deduction, he sends out
Existing Stoneley wave back wave coefficient is sensitive to horizontal fracture width, further established Stoneley wave reflectance factor and fracture width
Relation.
Electric imaging logging is a kind of using can cover multiple rows of button-type electrode of most of borehole wall area on pole plate to the borehole wall
Stratum emission current, surrounding well wall resistivity everywhere is measured, so as to form a kind of technology of borehole wall image.In the 1990s
Phase is micro-resisitivity image fast-developing period.The micro-resisitivity image research of major Logging Company
Start to be attended by the most momentous results successively in this period and come into operation in oil field.The presence in crack can make to go out in resistivity imaging figure
An existing black low-resistance band, and fracture width is bigger, and the resistivity of its corresponding band is lower, and band color is deeper.1990
Year Luthi have studied some parameters such as formation resistivity, mud resistivity to resistivity by way of finite element numerical simulation
The influence of the black stripe of cracking initiation in image, and then give and fracture width is calculated in a set of figure from resistivity imaging
Method.
For Tang Xiaoming et al. propose using Stoneley wave back wave calculate fracture width method, when fracture width compared with
Hour, reflection wave amplitude will be very small, and the fracture width degree of accuracy now calculated is relatively low;When many cracks in stratum being present
When, the Stoneley wave back wave of different cracking initiations can produce again to be interfered, so as to influence the degree of accuracy of its result of calculation.And
Method for being calculated fracture width using Electrical imaging figure, resistivity are not only influenceed by fracture width, also receive stratum
The combined influence of many kinds of parameters such as resistivity, mud resistivity, instrument pole plate and borehole wall distance, therefore, utilizes Electrical imaging figure meter
Fracture width is calculated to be difficult to accurately calculate horizontal fracture width.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device for determining horizontal fracture width, to obtain more
Accurate horizontal fracture width.
To reach above-mentioned purpose, on the one hand, the embodiment of the present application provides a kind of method for determining horizontal fracture width, bag
Include:
The numerical model of pending well is established, and simulates acoustic logging ripple of the numerical model under each bar horizontal fracture
Shape;
According to the acoustic logging waveform under each bar horizontal fracture, the first Stoneley wave width of this horizontal fracture is determined
Spend variable quantity;
According to the first Stoneley wave changes in amplitude amount of each bar horizontal fracture, fracture width and Stoneley wave width are determined
The relation curve spent between variable quantity;
According to the acoustic logging waveform of the pending well, obtain every horizontal fracture in the pending well second this
Tonneau wave amplitude variable quantity;
According to the second Stoneley wave changes in amplitude of every horizontal fracture in the relation curve and the pending well
Amount, determine the width of every horizontal fracture in the pending well.
The method of the determination horizontal fracture width of the embodiment of the present application, the numerical model for establishing pending well, and mould
Intend acoustic logging waveform of the numerical model under each bar horizontal fracture, including:
The numerical model of pending well is established according to default method for numerical simulation, and simulates the numerical model in each bar
Acoustic logging waveform under horizontal fracture.
The method of the determination horizontal fracture width of the embodiment of the present application, the sound wave according under each bar horizontal fracture
Well logging waveform, the first Stoneley wave changes in amplitude amount of this horizontal fracture is determined, including:
By the mean amplitude of tide of the Mintrop wave ripple bag in the acoustic logging waveform under each bar horizontal fracture, it is defined as this water
First Stoneley wave amplitude corresponding to plane fracture seam;
For every horizontal fracture, it is determined that maximum by the first Stoneley wave amplitude before and after this horizontal fracture with
The Magnitude Difference of minimum value, and determine first ratio of the Magnitude Difference and the maximum;
First ratio is defined as the first Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
The method of the determination horizontal fracture width of the embodiment of the present application, it is described according to the first of each bar horizontal fracture this
Tonneau wave amplitude variable quantity, the relation curve between fracture width and Stoneley wave changes in amplitude amount is determined, including:
First Stoneley wave changes in amplitude amount of each bar horizontal fracture is fitted, obtains fracture width and stone
Relation curve between sharp wave amplitude variable quantity.
The method of the determination horizontal fracture width of the embodiment of the present application, the relation curve include:
FW=ea·VA+b+c;
Wherein, FW is fracture width;VA is Stoneley wave changes in amplitude amount;A, b and c is coefficient.
The method of the determination horizontal fracture width of the embodiment of the present application, the acoustic logging ripple according to the pending well
Shape, the second Stoneley wave changes in amplitude amount of every horizontal fracture in the pending well is obtained, including:
LPF is carried out to the acoustic logging waveform of the pending well, obtains single Stoneley wave;
By way of ETAD expected time of arrival and departure window, every water in the pending well is obtained from the single Stoneley wave
The Stoneley wave amplitude curve of plane fracture seam;
For the Stoneley wave amplitude curve of every horizontal fracture in the pending well, determine its maximum with amplitude most
The Magnitude Difference of small value, and determine second ratio of the Magnitude Difference and the maximum;
Second ratio is defined as the second Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
The method of the determination horizontal fracture width of the embodiment of the present application, in addition to:
According to the Electric Log Data of the pending well, the position of every horizontal fracture in the pending well is obtained
Put.
On the other hand, the embodiment of the present application provides a kind of device for determining horizontal fracture width, including:
Numerical model establishes module, for establishing the numerical model of pending well, and simulates the numerical model in each bar
Acoustic logging waveform under horizontal fracture;
First changes in amplitude determining module, for according to the acoustic logging waveform under each bar horizontal fracture, it is determined that should
First Stoneley wave changes in amplitude amount of bar horizontal fracture;
Relation curve determining module, for the first Stoneley wave changes in amplitude amount according to each bar horizontal fracture, really
Determine the relation curve between fracture width and Stoneley wave changes in amplitude amount;
Second changes in amplitude determining module, for the acoustic logging waveform according to the pending well, wait to locate described in acquisition
Manage the second Stoneley wave changes in amplitude amount of every horizontal fracture in well;
Fracture width determining module, for according to every horizontal fracture in the relation curve and the pending well
Second Stoneley wave changes in amplitude amount, determine the width of every horizontal fracture in the pending well.
The device of the determination horizontal fracture width of the embodiment of the present application, the numerical model for establishing pending well, and mould
Intend acoustic logging waveform of the numerical model under each bar horizontal fracture, including:
The numerical model of pending well is established according to default method for numerical simulation, and simulates the numerical model in each bar
Acoustic logging waveform under horizontal fracture.
The device of the determination horizontal fracture width of the embodiment of the present application, the sound wave according under each bar horizontal fracture
Well logging waveform, the first Stoneley wave changes in amplitude amount of this horizontal fracture is determined, including:
By the mean amplitude of tide of the Mintrop wave ripple bag in the acoustic logging waveform under each bar horizontal fracture, it is defined as this water
First Stoneley wave amplitude corresponding to plane fracture seam;
For every horizontal fracture, it is determined that maximum by the first Stoneley wave amplitude before and after this horizontal fracture with
The Magnitude Difference of minimum value, and determine first ratio of the Magnitude Difference and the maximum;
First ratio is defined as the first Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
The device of the determination horizontal fracture width of the embodiment of the present application, it is described according to the first of each bar horizontal fracture this
Tonneau wave amplitude variable quantity, the relation curve between fracture width and Stoneley wave changes in amplitude amount is determined, including:
First Stoneley wave changes in amplitude amount of each bar horizontal fracture is fitted, obtains fracture width and stone
Relation curve between sharp wave amplitude variable quantity.
The device of the determination horizontal fracture width of the embodiment of the present application, the relation curve include:
FW=ea·VA+b+c;
Wherein, FW is fracture width;VA is Stoneley wave changes in amplitude amount;A, b and c is coefficient.
The device of the determination horizontal fracture width of the embodiment of the present application, the acoustic logging ripple according to the pending well
Shape, the second Stoneley wave changes in amplitude amount of every horizontal fracture in the pending well is obtained, including:
LPF is carried out to the acoustic logging waveform of the pending well, obtains single Stoneley wave;
By way of ETAD expected time of arrival and departure window, every water in the pending well is obtained from the single Stoneley wave
The Stoneley wave amplitude curve of plane fracture seam;
For the Stoneley wave amplitude curve of every horizontal fracture in the pending well, determine its maximum with amplitude most
The Magnitude Difference of small value, and determine second ratio of the Magnitude Difference and the maximum;
Second ratio is defined as the second Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
The device of the determination horizontal fracture width of the embodiment of the present application, in addition to:
Crack location determining module, for the Electric Log Data according to the pending well, obtain described pending
The position of every horizontal fracture in well.
Another further aspect, the embodiment of the present application additionally provide another device for determining horizontal fracture width, including memory,
Processor and the computer program being stored on the memory, the computer program are held when being run by the processor
Row following steps:
The numerical model of pending well is established, and simulates acoustic logging ripple of the numerical model under each bar horizontal fracture
Shape;
According to the acoustic logging waveform under each bar horizontal fracture, the first Stoneley wave width of this horizontal fracture is determined
Spend variable quantity;
According to the first Stoneley wave changes in amplitude amount of each bar horizontal fracture, fracture width and Stoneley wave width are determined
The relation curve spent between variable quantity;
According to the acoustic logging waveform of the pending well, obtain every horizontal fracture in the pending well second this
Tonneau wave amplitude variable quantity;
According to the second Stoneley wave changes in amplitude of every horizontal fracture in the relation curve and the pending well
Amount, determine the width of every horizontal fracture in the pending well.
The technical scheme provided from above the embodiment of the present application, the embodiment of the present application first establish the numerical value of pending well
Model, and acoustic logging waveform of the simulation value model under each bar horizontal fracture;Secondly according to the sound under each bar horizontal fracture
Ripple well logging waveform, determine the first Stoneley wave changes in amplitude amount of this horizontal fracture;Then according to the of each article of horizontal fracture
One Stoneley wave changes in amplitude amount, determines the relation curve between fracture width and Stoneley wave changes in amplitude amount;Then according to
The acoustic logging waveform of pending well, obtain the second Stoneley wave changes in amplitude amount of every horizontal fracture in pending well;Most
Afterwards according to the second Stoneley wave changes in amplitude amount of every horizontal fracture in relation curve and pending well, pending well is determined
In every horizontal fracture width, respond section because the embodiment of the present application can obtain crack by above-mentioned numerical simulation technology
Characteristic information, thus it effectively reduces multi-solution, so as to correct formation lithology and other factors to a certain extent
Interference, thus improve the degree of accuracy of horizontal fracture width calculation.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, do not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.In the accompanying drawings:
Fig. 1 is the method flow diagram for the method that horizontal fracture width is determined in the embodiment of the application one;
Fig. 2 is the stone profit in the wellhole of the horizontal fracture of the presence 1mm width of numerical simulation in the embodiment of the application one
Wave amplitude curve synoptic diagram;
Fig. 3 is the Representation Level fracture width and stone drawn in the embodiment of the application one by handling numerical simulation data
The exponential fitting curve of sharp wave amplitude variation relation;
Fig. 4 is certain processing well acoustic logging and the electric imaging logging result for existing in the embodiment of the application one horizontal fracture;
Fig. 5 is the structured flowchart for the device that horizontal fracture width is determined in the embodiment of the application one;
Fig. 6 is the structured flowchart for the device that horizontal fracture width is determined in another embodiment of the application.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.It is common based on the embodiment in the application, this area
The every other embodiment that technical staff is obtained under the premise of creative work is not made, it should all belong to the application protection
Scope.
With reference to shown in figure 1, the method for the determination horizontal fracture width of the embodiment of the present application can include:
S101, the numerical model for establishing pending well, and simulate sound wave of the numerical model under each bar horizontal fracture
Well logging waveform.
In the embodiment of the application one, the pending well is a well point, corresponding to the well point place included in reservoir
The fracture width of horizontal fracture is to be determined.
In the embodiment of the application one, the numerical model of pending well, and mould can be established according to default method for numerical simulation
Intend acoustic logging waveform of the numerical model under each bar horizontal fracture.Such as can be according to default in the embodiment of the application one
Method for numerical simulation establish the numerical model of pending well, and simulate sound wave of the numerical model under each bar horizontal fracture
Well logging waveform.In the embodiment of the application one, the numerical model of foundation can be consistent with the real well condition of the pending well,
For example use the parameter such as identical hole diameter, mud density and the velocity of sound and stratum p-and s-wave velocity.
S102, according to the acoustic logging waveform under each bar horizontal fracture, determine the first stone of this horizontal fracture
Sharp wave amplitude variable quantity.
In the embodiment of the application one, the acoustic logging waveform according under each bar horizontal fracture, this water is determined
First Stoneley wave changes in amplitude amount of plane fracture seam can include:
First by the mean amplitude of tide of the Mintrop wave ripple bag in the acoustic logging waveform under each bar horizontal fracture, it is defined as this
First Stoneley wave amplitude corresponding to horizontal fracture;Then for every horizontal fracture, it is determined that before and after by this horizontal fracture
The first Stoneley wave amplitude maxima and minima Magnitude Difference, and determine the first of the Magnitude Difference and the maximum
Ratio;First ratio is finally defined as the first Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
Such as the exemplary embodiment shown in Fig. 2 has been suffered and used and numerical simulation under the conditions of the pending equal hole condition of well
The Stoneley wave amplitude curve in the wellhole of the horizontal fracture of 1mm width be present.Amplitude is maximum in the Stoneley wave amplitude curve
It is worth for 0.01269mV, amplitude min value 0.01224mV, the ratio of both differences and Amplitude maxima is 3.55%, the ratio
Stoneley wave changes in amplitude amount as corresponding to the crack.
S103, the first Stoneley wave changes in amplitude amount according to each bar horizontal fracture, determine fracture width and stone
Relation curve between sharp wave amplitude variable quantity.
In the embodiment of the application one, in above-mentioned steps S101, in order to obtain accurate fitting formula, generally require
Simulate the horizontal fracture (situation for for example, at least simulating the horizontal fracture of more than 5 different in width) of a plurality of different in width.At this
Apply in an exemplary embodiment, exemplified by simulating the horizontal fracture of 8 different in width, the width of 8 horizontal fractures can be distinguished
Select as 0mm, 1mm, 5mm, 10mm, 20mm, 40mm, 60mm and 80mm.
In the embodiment of the application one, the first Stoneley wave changes in amplitude amount of each bar horizontal fracture can be fitted, obtained
Obtain the relation curve between fracture width and Stoneley wave changes in amplitude amount.In the embodiment of the application one, by taking exponential fitting as an example,
First Stoneley wave changes in amplitude amount of each bar horizontal fracture is carried out into exponential fitting can obtain relation curve as follows:
FW=ea·VA+b+c;
Wherein, FW is fracture width;VA is Stoneley wave changes in amplitude amount;A, b and c is coefficient.
In the exemplary embodiment of the application one, by taking the horizontal fracture of 8 above-mentioned different in width as an example, by above-mentioned 8
First Stoneley wave changes in amplitude amount of the horizontal fracture of different in width carries out exponential fitting, can obtain relation as shown in Figure 3
Curve FW=e0.0328VA+2.6-13.838。
S104, the acoustic logging waveform according to the pending well, obtain in the pending well every horizontal fracture
Second Stoneley wave changes in amplitude amount.
In the embodiment of the application one, the acoustic logging waveform according to the pending well, the pending well is obtained
In the second Stoneley wave changes in amplitude amount of every horizontal fracture can include:
LPF first can be carried out to the acoustic logging waveform of the pending well using digital filter, obtained single
Stoneley wave;
Secondly, by way of ETAD expected time of arrival and departure window, obtained from the single Stoneley wave in the pending well
The Stoneley wave amplitude curve of every horizontal fracture;Described ETAD expected time of arrival and departure window primarily directed to Stoneley wave first ripple bag and
Second ripple unwraps window, and time window length is typically chosen in 1ms.
Then for the Stoneley wave amplitude curve of every horizontal fracture in the pending well, its maximum and width are determined
The Magnitude Difference of minimum value is spent, and determines second ratio of the Magnitude Difference and the maximum;
Second ratio is finally defined as the second Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
It is different from the acoustic logging waveform of above-mentioned S101 simulations, in this step, the acoustic logging waveform of the pending well
To carry out acoustic logging to the pending well, and obtained actual sound wave measuring well curve.The exemplary embodiment of the application one
In, as shown in figure 4, by carrying out acoustic logging to certain well, obtaining monopole measured waveform, (first of ripple is played in left side in such as Fig. 4
Shape figure), after low-pass filtered, obtaining single Stoneley wave signal, (second oscillogram is played in left side in such as Fig. 4.Count successively
The Stoneley wave in (region between two vertical lines in Stoneley wave as shown in Figure 4) is calculated in the time window at each depth point
Average amplitude value, just obtain Stoneley wave amplitude curve (the 3rd oscillogram is played in left side in such as Fig. 4).
S105, the second Stoneley wave amplitude according to every horizontal fracture in the relation curve and the pending well
Variable quantity, determine the width of every horizontal fracture in the pending well.
In the embodiment of the present application, the second Stoneley wave changes in amplitude of every horizontal fracture in the pending well is obtained
After amount, the relation curve is substituted into, you can obtain the width of every horizontal fracture in surely described pending well.
In the embodiment of the application one, such as shown in Fig. 4, horizontal fracture is located near XXX5m, the 0.6m scopes above and below crack
Interior Amplitude maxima is 448.54mV, minimum value 218.58mV, Stoneley wave corresponding to the horizontal fracture thus calculated
Changes in amplitude amount is 51.27%, substitutes into the above-mentioned relation curve fitted, it is 58.5mm that can draw fracture width.It is worth noting
, under actual conditions, typically a plurality of horizontal fracture is arranged together so as to form a slit band, just as shown in this example,
The fracture width now calculated can be the width sum of a plurality of horizontal fracture.
, can also be according to the resistance of the pending well before fracture width is calculated in the application other embodiment
Rate log data, obtain the position of every horizontal fracture in the pending well, you can resistivity survey is carried out to the pending well
Well, such as resistivity logging is carried out to the pending well by resistivity imaging tool device, resistivity well week image is obtained,
So as to determine the approximate location of horizontal fracture in figure.It so can further eliminate multi-solution so that the fracture width of calculating
Really the truth on stratum is reflected.
With reference to shown in figure 5, the device of the determination horizontal fracture width of the embodiment of the present application can include:
Numerical model establishes module 51, for establishing the numerical model of pending well, and simulates the numerical model each
Acoustic logging waveform under bar horizontal fracture;
First changes in amplitude determining module 52, for according to the acoustic logging waveform under each bar horizontal fracture, it is determined that
First Stoneley wave changes in amplitude amount of this horizontal fracture;
Relation curve determining module 53, for the first Stoneley wave changes in amplitude amount according to each bar horizontal fracture,
Determine the relation curve between fracture width and Stoneley wave changes in amplitude amount;
Second changes in amplitude determining module 54, for the acoustic logging waveform according to the pending well, treated described in acquisition
Second Stoneley wave changes in amplitude amount of every horizontal fracture in processing well;
Fracture width determining module 55, for according to every horizontal fracture in the relation curve and the pending well
The second Stoneley wave changes in amplitude amount, determine the width of every horizontal fracture in the pending well.
In the application other embodiment, the device for determining horizontal fracture width can also include:
Crack location determining module 56, for the Electric Log Data according to the pending well, wait to locate described in acquisition
Manage the position of every horizontal fracture in well.
The device of the embodiment of the present application is corresponding with the method for above-described embodiment, therefore, is related to the device details of the application,
The method for referring to above-described embodiment, will not be repeated here.
With reference to shown in figure 6, the device of the determination horizontal fracture width of the embodiment of the present application can include memory, processing
Device and the computer program being stored on the memory, performed such as when the computer program is run by the processor
Lower step:
The numerical model of pending well is established, and simulates acoustic logging ripple of the numerical model under each bar horizontal fracture
Shape;
According to the acoustic logging waveform under each bar horizontal fracture, the first Stoneley wave width of this horizontal fracture is determined
Spend variable quantity;
According to the first Stoneley wave changes in amplitude amount of each bar horizontal fracture, fracture width and Stoneley wave width are determined
The relation curve spent between variable quantity;
According to the acoustic logging waveform of the pending well, obtain every horizontal fracture in the pending well second this
Tonneau wave amplitude variable quantity;
According to the second Stoneley wave changes in amplitude of every horizontal fracture in the relation curve and the pending well
Amount, determine the width of every horizontal fracture in the pending well.
The device of the embodiment of the present application is corresponding with the method for above-described embodiment, therefore, is related to the device details of the application,
The method for referring to above-described embodiment, will not be repeated here.
For convenience of description, it is divided into various units during description apparatus above with function to describe respectively.Certainly, this is being implemented
The function of each unit can be realized in same or multiple softwares and/or hardware during application.
The present invention is the flow with reference to method according to embodiments of the present invention, equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that can be by every first-class in computer program instructions implementation process figure and/or block diagram
Journey and/or the flow in square frame and flow chart and/or block diagram and/or the combination of square frame.These computer programs can be provided
The processors of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices is instructed to produce
A raw machine so that produced by the instruction of computer or the computing device of other programmable data processing devices for real
The device for the function of being specified in present one flow of flow chart or one square frame of multiple flows and/or block diagram or multiple square frames.
These computer program instructions, which may be alternatively stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that the instruction being stored in the computer-readable memory, which produces, to be included referring to
Make the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one square frame of block diagram or
The function of being specified in multiple square frames.
These computer program instructions can be also loaded into computer or other programmable data processing devices so that counted
Series of operation steps is performed on calculation machine or other programmable devices to produce computer implemented processing, so as in computer or
The instruction performed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in individual square frame or multiple square frames.
In a typical configuration, computing device includes one or more processors (CPU), input/output interface, net
Network interface and internal memory.
Internal memory may include computer-readable medium in volatile memory, random access memory (RAM) and/or
The forms such as Nonvolatile memory, such as read-only storage (ROM) or flash memory (flash RAM).Internal memory is computer-readable medium
Example.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology come realize information store.Information can be computer-readable instruction, data structure, the module of program or other data.
The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM (SRAM), moved
State random access memory (DRAM), other kinds of random access memory (RAM), read-only storage (ROM), electric erasable
Programmable read only memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc read-only storage (CD-ROM),
Digital versatile disc (DVD) or other optical storages, magnetic cassette tape, the storage of tape magnetic rigid disk or other magnetic storage apparatus
Or any other non-transmission medium, the information that can be accessed by a computing device available for storage.Define, calculate according to herein
Machine computer-readable recording medium does not include temporary computer readable media (transitory media), such as data-signal and carrier wave of modulation.
It should also be noted that, term " comprising ", "comprising" or its any other variant are intended to nonexcludability
Comprising so that process, method, commodity or equipment including a series of elements not only include those key elements, but also wrapping
Include the other element being not expressly set out, or also include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that wanted including described
Other identical element also be present in the process of element, method, commodity or equipment.
It will be understood by those skilled in the art that embodiments herein can be provided as method, system or computer program product.
Therefore, the application can be using the embodiment in terms of complete hardware embodiment, complete software embodiment or combination software and hardware
Form.Deposited moreover, the application can use to can use in one or more computers for wherein including computer usable program code
The shape for the computer program product that storage media is implemented on (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.)
Formula.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Task is performed and connected remote processing devices by communication network.In a distributed computing environment, program module can be with
In the local and remote computer-readable storage medium including storage device.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment stressed is the difference with other embodiment.It is real especially for system
For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
Embodiments herein is the foregoing is only, is not limited to the application.For those skilled in the art
For, the application can have various modifications and variations.All any modifications made within spirit herein and principle, it is equal
Replace, improve etc., it should be included within the scope of claims hereof.
Claims (15)
- A kind of 1. method for determining horizontal fracture width, it is characterised in that including:The numerical model of pending well is established, and simulates acoustic logging waveform of the numerical model under each bar horizontal fracture;According to the acoustic logging waveform under each bar horizontal fracture, determine that the first Stoneley wave amplitude of this horizontal fracture becomes Change amount;According to the first Stoneley wave changes in amplitude amount of each bar horizontal fracture, determine that fracture width becomes with Stoneley wave amplitude Relation curve between change amount;According to the acoustic logging waveform of the pending well, the second stone for obtaining every horizontal fracture in the pending well is sharp Wave amplitude variable quantity;According to the second Stoneley wave changes in amplitude amount of every horizontal fracture in the relation curve and the pending well, really The width of every horizontal fracture in the fixed pending well.
- 2. the method for horizontal fracture width is determined as claimed in claim 1, it is characterised in that the number for establishing pending well It is worth model, and simulates acoustic logging waveform of the numerical model under each bar horizontal fracture, including:The numerical model of pending well is established according to default method for numerical simulation, and it is horizontal in each bar to simulate the numerical model Acoustic logging waveform under crack.
- 3. the method for horizontal fracture width is determined as claimed in claim 1, it is characterised in that described horizontal according to each bar Acoustic logging waveform under crack, the first Stoneley wave changes in amplitude amount of this horizontal fracture is determined, including:By the mean amplitude of tide of the Mintrop wave ripple bag in the acoustic logging waveform under each bar horizontal fracture, it is defined as this level and splits First Stoneley wave amplitude corresponding to seam;For every horizontal fracture, it is determined that maximum and minimum by the first Stoneley wave amplitude before and after this horizontal fracture The Magnitude Difference of value, and determine first ratio of the Magnitude Difference and the maximum;First ratio is defined as the first Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
- 4. the method for horizontal fracture width is determined as claimed in claim 1, it is characterised in that described horizontal according to each bar The first Stoneley wave changes in amplitude amount in crack, determines the relation curve between fracture width and Stoneley wave changes in amplitude amount, Including:First Stoneley wave changes in amplitude amount of each bar horizontal fracture is fitted, obtains fracture width and Stoneley wave Relation curve between changes in amplitude amount.
- 5. the method for horizontal fracture width is determined as claimed in claim 1, it is characterised in that the relation curve includes:FW=ea·VA+b+c;Wherein, FW is fracture width;VA is Stoneley wave changes in amplitude amount;A, b and c is coefficient.
- 6. the method for horizontal fracture width is determined as claimed in claim 1, it is characterised in that described according to the pending well Acoustic logging waveform, obtain the second Stoneley wave changes in amplitude amount of every horizontal fracture in the pending well, including:LPF is carried out to the acoustic logging waveform of the pending well, obtains single Stoneley wave;By way of ETAD expected time of arrival and departure window, the every level in the pending well that obtained from the single Stoneley wave is split The Stoneley wave amplitude curve of seam;For the Stoneley wave amplitude curve of every horizontal fracture in the pending well, its maximum and amplitude min value are determined Magnitude Difference, and determine second ratio of the Magnitude Difference and the maximum;Second ratio is defined as the second Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
- 7. the method for horizontal fracture width is determined as claimed in claim 1, it is characterised in that also include:According to the Electric Log Data of the pending well, the position of every horizontal fracture in the pending well is obtained.
- A kind of 8. device for determining horizontal fracture width, it is characterised in that including:Numerical model establishes module, and for establishing the numerical model of pending well, and it is horizontal in each bar to simulate the numerical model Acoustic logging waveform under crack;First changes in amplitude determining module, for according to the acoustic logging waveform under each bar horizontal fracture, determining this water First Stoneley wave changes in amplitude amount of plane fracture seam;Relation curve determining module, for the first Stoneley wave changes in amplitude amount according to each bar horizontal fracture, it is determined that splitting Relation curve between slit width degree and Stoneley wave changes in amplitude amount;Second changes in amplitude determining module, for the acoustic logging waveform according to the pending well, obtain the pending well In every horizontal fracture the second Stoneley wave changes in amplitude amount;Fracture width determining module, for second according to every horizontal fracture in the relation curve and the pending well Stoneley wave changes in amplitude amount, determine the width of every horizontal fracture in the pending well.
- 9. the device of horizontal fracture width is determined as claimed in claim 8, it is characterised in that the number for establishing pending well It is worth model, and simulates acoustic logging waveform of the numerical model under each bar horizontal fracture, including:The numerical model of pending well is established according to default method for numerical simulation, and it is horizontal in each bar to simulate the numerical model Acoustic logging waveform under crack.
- 10. the device of horizontal fracture width is determined as claimed in claim 8, it is characterised in that described according to each bar water Acoustic logging waveform under plane fracture seam, the first Stoneley wave changes in amplitude amount of this horizontal fracture is determined, including:By the mean amplitude of tide of the Mintrop wave ripple bag in the acoustic logging waveform under each bar horizontal fracture, it is defined as this level and splits First Stoneley wave amplitude corresponding to seam;For every horizontal fracture, it is determined that maximum and minimum by the first Stoneley wave amplitude before and after this horizontal fracture The Magnitude Difference of value, and determine first ratio of the Magnitude Difference and the maximum;First ratio is defined as the first Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
- 11. the device of horizontal fracture width is determined as claimed in claim 8, it is characterised in that described according to each bar water First Stoneley wave changes in amplitude amount of plane fracture seam, determine that the relation between fracture width and Stoneley wave changes in amplitude amount is bent Line, including:First Stoneley wave changes in amplitude amount of each bar horizontal fracture is fitted, obtains fracture width and Stoneley wave Relation curve between changes in amplitude amount.
- 12. the device of horizontal fracture width is determined as claimed in claim 8, it is characterised in that the relation curve includes:FW=ea·VA+b+c;Wherein, FW is fracture width;VA is Stoneley wave changes in amplitude amount;A, b and c is coefficient.
- 13. the device of horizontal fracture width is determined as claimed in claim 8, it is characterised in that described according to described pending The acoustic logging waveform of well, the second Stoneley wave changes in amplitude amount of every horizontal fracture in the pending well is obtained, including:LPF is carried out to the acoustic logging waveform of the pending well, obtains single Stoneley wave;By way of ETAD expected time of arrival and departure window, the every level in the pending well that obtained from the single Stoneley wave is split The Stoneley wave amplitude curve of seam;For the Stoneley wave amplitude curve of every horizontal fracture in the pending well, its maximum and amplitude min value are determined Magnitude Difference, and determine second ratio of the Magnitude Difference and the maximum;Second ratio is defined as the second Stoneley wave changes in amplitude amount corresponding to this horizontal fracture.
- 14. the device of horizontal fracture width is determined as claimed in claim 8, it is characterised in that also include:Crack location determining module, for the Electric Log Data according to the pending well, obtain in the pending well The position of every horizontal fracture.
- 15. the device of horizontal fracture width is determined a kind of, including memory, processor and is stored on the memory Computer program, it is characterised in that the computer program performs following steps when being run by the processor:The numerical model of pending well is established, and simulates acoustic logging waveform of the numerical model under each bar horizontal fracture;According to the acoustic logging waveform under each bar horizontal fracture, determine that the first Stoneley wave amplitude of this horizontal fracture becomes Change amount;According to the first Stoneley wave changes in amplitude amount of each bar horizontal fracture, determine that fracture width becomes with Stoneley wave amplitude Relation curve between change amount;According to the acoustic logging waveform of the pending well, the second stone for obtaining every horizontal fracture in the pending well is sharp Wave amplitude variable quantity;According to the second Stoneley wave changes in amplitude amount of every horizontal fracture in the relation curve and the pending well, really The width of every horizontal fracture in the fixed pending well.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710666181.5A CN107587871B (en) | 2017-08-07 | 2017-08-07 | Method and device for determining horizontal crack width |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710666181.5A CN107587871B (en) | 2017-08-07 | 2017-08-07 | Method and device for determining horizontal crack width |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107587871A true CN107587871A (en) | 2018-01-16 |
CN107587871B CN107587871B (en) | 2020-05-08 |
Family
ID=61042114
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710666181.5A Active CN107587871B (en) | 2017-08-07 | 2017-08-07 | Method and device for determining horizontal crack width |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107587871B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375785A (en) * | 2018-01-12 | 2018-08-07 | 中国石油天然气股份有限公司 | Slit band method for correcting position and device |
CN108490490A (en) * | 2018-02-08 | 2018-09-04 | 中国石油天然气股份有限公司 | Slit band characterizing method, device, electronic equipment and computer storage media |
CN111140263A (en) * | 2020-04-02 | 2020-05-12 | 中铁工程设计咨询集团有限公司 | Tunnel concrete lining crack width calculation method and crack treatment method |
CN112666554A (en) * | 2020-12-17 | 2021-04-16 | 江苏中路工程技术研究院有限公司 | Method for identifying radar amplitude characteristic crack width of asphalt pavement |
CN113123784A (en) * | 2020-01-10 | 2021-07-16 | 中国石油天然气股份有限公司 | Method and device for obtaining permeability of low-pore fractured sandstone reservoir fractures |
CN113533528A (en) * | 2020-04-22 | 2021-10-22 | 中国石油天然气股份有限公司 | Stoneley wave experiment method for shock tube of rock core model containing cracks |
CN114151072A (en) * | 2020-09-08 | 2022-03-08 | 中国石油天然气股份有限公司 | Crack recognition device |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064118A (en) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | Method of acoustic logging and quantifying cavern filling degree |
CN104345346A (en) * | 2014-10-30 | 2015-02-11 | 中国石油天然气股份有限公司 | Method for obtaining crack widths |
-
2017
- 2017-08-07 CN CN201710666181.5A patent/CN107587871B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064118A (en) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | Method of acoustic logging and quantifying cavern filling degree |
CN104345346A (en) * | 2014-10-30 | 2015-02-11 | 中国石油天然气股份有限公司 | Method for obtaining crack widths |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108375785A (en) * | 2018-01-12 | 2018-08-07 | 中国石油天然气股份有限公司 | Slit band method for correcting position and device |
CN108490490A (en) * | 2018-02-08 | 2018-09-04 | 中国石油天然气股份有限公司 | Slit band characterizing method, device, electronic equipment and computer storage media |
CN113123784A (en) * | 2020-01-10 | 2021-07-16 | 中国石油天然气股份有限公司 | Method and device for obtaining permeability of low-pore fractured sandstone reservoir fractures |
CN113123784B (en) * | 2020-01-10 | 2023-12-22 | 中国石油天然气股份有限公司 | Method and device for acquiring fracture permeability of low-hole fractured sandstone reservoir |
CN111140263A (en) * | 2020-04-02 | 2020-05-12 | 中铁工程设计咨询集团有限公司 | Tunnel concrete lining crack width calculation method and crack treatment method |
CN113533528A (en) * | 2020-04-22 | 2021-10-22 | 中国石油天然气股份有限公司 | Stoneley wave experiment method for shock tube of rock core model containing cracks |
CN113533528B (en) * | 2020-04-22 | 2023-12-26 | 中国石油天然气股份有限公司 | Shock tube Stoneley wave experimental method for core model containing cracks |
CN114151072A (en) * | 2020-09-08 | 2022-03-08 | 中国石油天然气股份有限公司 | Crack recognition device |
CN114151072B (en) * | 2020-09-08 | 2023-12-22 | 中国石油天然气股份有限公司 | Crack recognition device |
CN112666554A (en) * | 2020-12-17 | 2021-04-16 | 江苏中路工程技术研究院有限公司 | Method for identifying radar amplitude characteristic crack width of asphalt pavement |
Also Published As
Publication number | Publication date |
---|---|
CN107587871B (en) | 2020-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107587871A (en) | Determine the method and device of horizontal fracture width | |
CN104375188B (en) | Seismic wave transmission attenuation compensation method and device | |
CN104793247B (en) | Near-surface velocity modeling method and model building device | |
CN105158796B (en) | The method and apparatus for determining TOC contents | |
CN113919219A (en) | Stratum evaluation method and system based on logging big data | |
US9938818B2 (en) | Hydraulic fracture permeability characterization from electromagnetic measurements | |
CN106896406A (en) | The method and apparatus that High-quality Reservoir is predicted based on impedance Domain Properties | |
AU2012385250A1 (en) | Anisotropy parameter estimation | |
CN107956465A (en) | The more borehole logging tool log normalization method and devices of the whole district based on association well | |
CN107422380A (en) | Carbonate rock fractured cave type Reservoir Body partition of the scale and quantization method | |
CN108873065A (en) | Sandstone High-quality Reservoir prediction technique and device | |
US11209559B2 (en) | Method and system for analyzing a borehole using passive acoustic logging | |
CN104267434B (en) | Three-dimensional multi-component earthquake observation system geophone offset distribution obtaining method and device | |
CN106125133B (en) | It is a kind of based on gas cloud area constrain under fine velocity modeling method | |
US10379257B2 (en) | Distributing petrofacies using analytical modeling | |
CN111239809B (en) | Dominant reservoir determination method and system based on attribute fusion principle | |
CN109991663B (en) | Work area seismic velocity body correction method and device | |
CN111060416A (en) | Method and system for dividing surface outcrop and rock core mechanical layer | |
CN108490490A (en) | Slit band characterizing method, device, electronic equipment and computer storage media | |
EP2751601B1 (en) | Measurement of permeability in borehole in the presence of mudcake | |
CN111812716A (en) | Pre-stack quantitative prediction method, device and equipment for total organic carbon content of shale gas reservoir | |
CN112285782B (en) | Near-surface seismic wave absorption attenuation investigation method and device | |
AU2014394076A1 (en) | Methods and systems for identifying and plugging subterranean conduits | |
CN113267810B (en) | Seismic exploration full-depth velocity modeling method and device | |
WO2023123952A1 (en) | Method and apparatus for inverting formation wave impedance using das well seismic data |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |