CA2964862C - Filtering microseismic events for updating and calibrating a fracture model - Google Patents

Filtering microseismic events for updating and calibrating a fracture model Download PDF

Info

Publication number
CA2964862C
CA2964862C CA2964862A CA2964862A CA2964862C CA 2964862 C CA2964862 C CA 2964862C CA 2964862 A CA2964862 A CA 2964862A CA 2964862 A CA2964862 A CA 2964862A CA 2964862 C CA2964862 C CA 2964862C
Authority
CA
Canada
Prior art keywords
fracture
observed
computing device
geometry
hydraulic fracture
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.)
Active
Application number
CA2964862A
Other languages
French (fr)
Other versions
CA2964862A1 (en
Inventor
Jason D. Dykstra
Zhijie Sun
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Halliburton Energy Services Inc
Original Assignee
Halliburton Energy Services Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Halliburton Energy Services Inc filed Critical Halliburton Energy Services Inc
Publication of CA2964862A1 publication Critical patent/CA2964862A1/en
Application granted granted Critical
Publication of CA2964862C publication Critical patent/CA2964862C/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/282Application of seismic models, synthetic seismograms
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/12Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/288Event detection in seismic signals, e.g. microseismics
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06GANALOGUE COMPUTERS
    • G06G7/00Devices in which the computing operation is performed by varying electric or magnetic quantities
    • G06G7/48Analogue computers for specific processes, systems or devices, e.g. simulators
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/646Fractures

Abstract

A fracture model for a hydraulic fracture in a wellbore can be updated and calibrated. Information about a microseismic event can be received from a sensor that is monitoring a subterranean formation. The information can be received subsequent to a fracking fluid being introduced into the formation. An observed geometry of a hydraulic fracture can be determined based on the information and a predicted geometry of the fracture can be determined based on properties of the fracking fluid and a fracture model. The fracture model can be updated using the information about the microseismic event where it is determined that an uncertainty value of the observed geometry does not exceed a pre-set maximum. The uncertainty value can be based on the predicted geometry of the hydraulic fracture.

Description

FILTERING MICROSEISMIC EVENTS FOR UPDATING AND
CALIBRATING A FRACTURE MODEL
Technical Field [0001] The present disclosure relates generally to systems and methods for analyzing a subterranean formation. More specifically, but not by way of limitation, this disclosure relates to improved hydraulic fracture modeling.
Background
[0002] A well system (e.g., oil or gas wells for extracting fluids from a subterranean formation) can include a wellbore drilled into a formation.
Hydraulic fracturing operations can be performed on the wellbore to increase production by pumping a hydraulic fluid down the cased wellbore into the formation at pressures and injection rates sufficient to cause the formation rock to initiate and propagate a hydraulic fracture (or induced fracture) into the subterranean formation.
Accurate estimation of the geometry of the hydraulic fracture can improve production economics by increasing reservoir productivity and reducing completion costs.
Summary
[0003] A fracture model for a hydraulic fracture in a wellbore can be updated and calibrated. Information about a microseismic event can be received from a sensor that is monitoring a subterranean formation. The information can be received subsequent to a fracking fluid being introduced into the formation. An observed geometry of a hydraulic fracture can be determined based on the information and a predicted geometry of the fracture can be determined based on properties of the fracking fluid and a fracture model. The fracture model can be updated using the information about the microseismic event where it is determined that an uncertainty value of the observed geometry does not exceed a pre-set maximum. The uncertainty value can be based on the predicted geometry of the hydraulic fracture.
Brief Description of the Drawings
[0004] FIG. 1 is a schematic illustration of an active wellbore and a microseismic monitoring system according to one aspect of the present disclosure.
[0005] FIG. 2 is a schematic illustration of an active wellbore and a microseismic monitoring system according to one aspect of the present disclosure.
[0006] FIG. 3 is a is a block diagram depicting an example of a computing device for updating and calibrating a fracture model using microseismic event data according to one aspect of the present disclosure.
[0007] FIG. 4 is an illustration of the Perkins-Kern-Nordgren model of a fracture according to one aspect of the present disclosure.
[0008] FIG. 5 is an example of a flow chart for a process of updating and calibrating a fracture model using microseismic event data according to one aspect of the present disclosure.
Detailed Description
[0009] Certain aspects and examples of the disclosure are directed to using data acquired from sensors, for example geophones or tiltmeters, to generate and filter a data set used to update, and in some aspects calibrate, a fracture model. A
well updated and calibrated fracture model can provide an accurate estimation of the geometry of a hydraulic fracture. Accurate estimation of the geometry of the hydraulic fracture can improve production economics by increasing reservoir productivity and reducing completion costs.
[0010]
Hydraulic fracturing operations can cause microseismic events to occur in the subterranean formation. Microseismic events in subterranean formations are shear events that are slippages along hydraulic fractures as well as natural fractures, bedding, faults, dewatering features, and other planes of weakness. During a hydraulic fracturing operation, the formation stress and the pore pressure in the formation surrounding a hydraulic fracture are increased. Both the pore pressure increase and formation can cause microseismic events to occur in the subterranean formation proximate to the tip of a hydraulic fracture. A computing device can receive data from the sensors that monitor these microseismic events and can calculate the location of the microseismic event using the data and additional information about the formation. The computing device can use the location of the microseismic event to determine an observed length of the hydraulic fracture.
For example, in a bi-wing fracture the length of the fracture is the distance from the main wellbore to the tip where the microseismic event is assumed to have occurred.
The computing device can use information about the hydraulic fracturing operation and a fracture model to predict the geometry of the hydraulic fracture, for example the length of the hydraulic fracture. The fracture model can be updated with observed measurements to stop error accumulation by updating the state vector. In some aspects, the fracture model can be calibrated when necessary to update additional vectors of the fracture model. The computing device can use the observed length of the hydraulic fracture to update the fracture model and in some aspects calibrate the fracture model.
[0011]
Hydraulic fracturing can also cause microseismic events in the subterranean formation at locations other than the tip of the hydraulic fracture. For example, microseismic events can occur at a natural fracture. When the location of the microseismic event monitored by the sensors is not at the tip of the hydraulic fracture the observed length, or other geometry of the hydraulic fracture determined by the computing device, can be inaccurate. When the computing device uses an observed length associated with a microseismic event that did not occur at the tip of the hydraulic fracture to update the fracture model it can decrease the accuracy of the fracture model.
[0012] The computing device can increase the accuracy of the fracture model by calibrating the fracture model using observed lengths of the hydraulic fracture that correspond to microseismic events that occur near the tip of the fracture. The computing device can determine if the microseismic event occurred near the tip of the fracture by comparing an uncertainty value of the observed length to a maximum uncertainty value. The uncertainty value of the observed length can be determined by comparing the observed length to a predicted length of the hydraulic fracture based on the fracture model.
[0013] FIG. 1 is a schematic illustration of a microseismic monitoring system 100 according to one aspect. An active wellbore 102 extends through various earth strata. The active wellbore has a substantially vertical section 104 and a substantially horizontal section 106. The substantially vertical section 104 may include a casing string 108 cemented at an upper portion of the substantially vertical section 104. The substantially horizontal section 106 extends through a hydrocarbon bearing subterranean formation 110. A tubing string 112 extends from the surface 114 into the active wellbore 102. The tubing string 112 can provide a conduit for pumping a fracking fluid into the active wellbore 102 to perform hydraulic fracking operations on the active wellbore 102. The active wellbore 102 includes a hydraulic fracture 116 (or induced fracture) that extends from the substantially horizontal section 106. Proppant materials can be entrained in the fracturing fluid and deposited in the hydraulic fracture 116 to keep the hydraulic fracture 116 open.
[0014] The system 100 can include an observation wellbore 118 extending through various earth strata. The observation wellbore 118 has a substantially vertical section 120. Multiple geophones 122 are positioned on a tool 121 deployed within the observation wellbore 118. In some aspects, transducers, tiltmeters or other suitable sensors may be used in place of the geophones 122. The geophones 122 can monitor microseismic events in the formation 110. For example, the geophones 122 can detect the arrival of the Primary (P) waves and the Secondary (S) waves of a microseismic event in the formation 110 during a hydraulic fracturing operation. In some aspects, tiltmeters may be positioned at the surface 114 proximate the horizontal location of the hydraulic fracture 116 or in shallow holes at the surface 114. The surface tiltmeters can be communicatively coupled to the computing device 124 and can transmit data sufficient to determine a location of a microseismic event in the formation 110.
[0015] The geophones 122 can be communicatively coupled to a computing device 124. The computing device 124 can be positioned at the surface 114 or at a separate location. The computing device 124 can include a processor interfaced with other hardware via a bus. A memory, which can include any suitable tangible (and non-transitory) computer-readable medium, such as random-access memory ("RAM"), ROM, electrically erasable and programmable read-only memory ("EEPROM")õ or the like, can embody program components that configure operation of the computing device 124. In some aspects, the computing device can include input/output interface components (e.g., a display, keyboard, touch-sensitive surface, and mouse) and additional storage. The computing device 124 can transmit data to and receive data from the geophones 122 via a communication link 123. The communication link 123 is wireless and can include wireless interfaces such as IEEE 802.11, Bluetooth, or radio interfaces for accessing cellular telephone networks (e.g., transceiver/antenna for accessing a CDMA, GSM, UMTS, or other mobile communications network). In other aspects, the communication link 123 can be wired and can include interfaces such as Ethernet, USB, IEEE 1394, or a fiber optic interface. An example of the computing device 124 is described in further detail with respect to FIG. 3.
[0016] The computing device 124 can use the data received from the array of geophones 122 to determine the location of a microseismic event that occurs in the formation 110 during a fracking operation. For example, the computing device can receive from geophones 122 raw signal data corresponding to the arrival times of the P-waves and S-waves associated with the microseismic event. The computing device 124 can calculate the difference between the travel times of the P-waves and S-waves using data associated with the formation 110 to determine a location of the microseismic event. The computing device 124 can determine an observed length of the hydraulic fracture 116 based on the location of the microseismic event and the assumption that the event occurred at the tip of the hydraulic fracture.
[0017] The computing device 124 can also determine an estimated or a predicted length of the hydraulic fracture 116 (hereinafter "predicted length") based on a fracture model. For example, a bi-wing fracture can be characterized by the Perkins-Kern-Nordgren (PKN) model. In some aspects, the computing device can determine an estimated or predicted width or other geometry of the hydraulic fracture 116. The computing device 124 can also determine an uncertainty value of the observed length of the hydraulic fracture 116 that was calculated using the data received from the geophones 122. The uncertainty value of the observed length of the hydraulic fracture 116 can be the difference between the observed length and the predicted length (the length predicted by the model) of the hydraulic fracture 116.
[0018] The computing device 124 can determine whether the uncertainty value of the observed length of the hydraulic fracture is less than or equal to a maximum uncertainty value. The maximum uncertainty value can be an uncertainty region corresponding to the location of the tip of the hydraulic fracture and can be based on an uncertainty of the predicted length of the hydraulic fracture 116.
The computing device 124 can determine the uncertainty of the predicted length of the hydraulic fracture 116 based on the model. For example, the uncertainty of the = CkP(k)CkT +Vk predicted length can be determined as ZL If 95%
certainty that the microseismic event measured occurred at or near the fracture tip is desired, then the maximum uncertainty value may be 2ZL . If a 99% certainty is desired, then the maximum uncertainty value may be3ZL . The computing device 124 can receive the maximum uncertainty value from an input by a user.
[0019] An observed length of the hydraulic fracture 116 that has an uncertainty value that exceeds maximum uncertainty value may be considered to be an inaccurate measurement of the length of the hydraulic fracture 116. An observed length of the hydraulic fracture 116 that has an uncertainty value that does not exceed the maximum uncertainty value can be considered to be an accurate measurement of the length of the hydraulic fracture 116 and can be used to update the model.
[0020] In some aspects, the computing device 124 can store all the observed lengths of the hydraulic fracture 116, regardless of their uncertainty value, in a data store. In some aspects, the computing device 124 can store only the observed lengths of the hydraulic fracture 116 having an uncertainty value that does not exceed the maximum uncertainty value in the data store.
[0021] The observed lengths of the hydraulic fracture 116 can also include a classification that indicates whether the observed length is associated with an uncertainty value that exceeds or does not exceed the maximum uncertainty value.
An observed length of the hydraulic fracture 116 having an uncertainty value that exceeds the maximum uncertainty value can be considered to be a length measurement calculated from microseismic event far from the tip of the hydraulic fracture. An observed fracture length of the hydraulic fracture 116 having an uncertainty value that does not exceed the maximum uncertainty value can be considered to be a length measurement calculated from a microseismic event near the tip of the hydraulic fracture 116.
[0022] If the data store includes more observed lengths of the hydraulic fracture 116 associated with an uncertainty value that exceeds the maximum uncertainty value then calibration or adaptation of the fracture model may be needed. The computing device 124 can calibrate the fracture model by running a clustering algorithm on the observed lengths in the data store and determining a cluster of observed lengths that has the greatest possibility of representing the fracture growth. The computing device 124 can re-estimate or calibrate the fracture model using the cluster of observed lengths determined to have the greatest possibility of representing the fracture growth.
[0023] FIG. 2 is a schematic illustration of a microseismic monitoring system 200 according to one aspect. The wellbore 102 extends through various earth strata and has a substantially vertical section 104 and a substantially horizontal section 106. The substantially vertical section 104 has the casing string 108 cemented at an upper portion of the substantially vertical section 104. The substantially horizontal section 206 extends through the hydrocarbon bearing subterranean formation 110.
Tubing string 112 extends from the surface 114 into the wellbore 102 and can provide a conduit for pumping a fracking fluid into the wellbore 202 to perform hydraulic fracking operations on the wellbore 102. The hydraulic fracture 116 extends from the substantially horizontal section 106. Proppant materials can be entrained in the fracturing fluid and deposited in the hydraulic fracture 116 to keep the hydraulic fracture 116 open.
[0024] The system 200 includes multiple tiltmeters 218 positioned at the surface 114. In some aspects, each of the tiltmeters 218 can be positioned in a shallow hole at the surface 114 proximate the lateral location of the hydraulic fracture 116. The tiltmeters 218 can monitor microseismic events in the formation 110.
For example, the tiltmeters 218 can measure hydraulic fracture-induced deformations or microdeformations at the surface 114 using orthogonal tilt sensors. In some aspects, the tiltmeters 218 may be measure microdeformations in the subterranean formation when positioned downhole in an observation wellbore.
[0025] The tiltmeters 218 are communicatively coupled to the computing device 124 via a wireless communication link 220. In some aspects, the communication link 220 can be wired. The computing device 124 can transmit data to and receive data from the tiltmeters 218 via the communication link 220. An example of the computing device 124 is described in further detail with respect to FIG. 3.
[0026] The computing device 124 can receive data from the tiltmeters 218 related to a microseismic event in the formation 210 during a hydraulic fracturing operation. For example, the computing device 124 can receive measurements of hydraulic fracture-induced deformations or microdeformations at the surface 214.
The computing device 124 can use the data received from the tiltmeters 218 to determine a length of a hydraulic fracture. For example, the length of a vertical hydraulic fracture can be determined by calculating the relative distance on the horizontal plane between the microseismic event locations and perforations at a wellbore.
[0027] As described with respect to the system 100 of FIG. 1, the computing device 124 can also determine a predicted length of the hydraulic fracture 116 based on a fracture growth model, for example, by using the PKN model where the hydraulic fracture 116 is a bi-wing fracture. The computing device 124 can determine an uncertainty value of the observed length of the hydraulic fracture 116 that was calculated using the data received from the tiltmeters 218. The uncertainty value of the observed length of the hydraulic fracture 116 can be the difference between the observed length and the predicted length (the length predicted by the model) of the hydraulic fracture 116.
[0028] The computing device 124 can determine an uncertainty value of the V
predicted length of the hydraulic fracture 116 is:
=CkP(k)Ckr + k . The computing device 124 can receive an input from a user that is a maximum uncertainty value based on the uncertainty of the predicted length of the hydraulic fracture 116. For Ii example, if 95% certainty that the microseismic event measured occurred at or near the fracture tip, then the maximum uncertainty value can be2IL . If a 99%
certainty is desired, the maximum uncertainty value can be3IL .
[0029] The computing device 124 can determine whether the uncertainty value of the observed length of the hydraulic fracture 116 is less than or equal to the maximum uncertainty value. An observed length of the hydraulic fracture 116 that has an uncertainty value that is less than or equal to the maximum uncertainty value can be used to update the fracture model. In some aspects, the computing device 124 can store all the observed lengths of the hydraulic fracture 116 in a data store.
In some aspects, the computing device 124 can store only the observed lengths of the hydraulic fracture 116 having an uncertainty value that does not exceed the maximum uncertainty value in the data store.
[0030] An observed length of the hydraulic fracture 116 that has an uncertainty value that is greater than the maximum uncertainty value may be considered an inaccurate calculation of the length of the hydraulic fracture 116 and may not be used to update the model. The computing device 124 can increase the accuracy of the fracture model by calibrating the fracture model with the observed lengths of the hydraulic fracture 116 that have an uncertainty value that does not exceed the maximum uncertainty value.
[0031] FIG. 3 is a block diagram depicting an example of a computing device 124 for determining an observed length of a hydraulic fracture and calibrating the fracture model according to one aspect of the present disclosure. The computing device 124 includes a processing device 302, a memory device 306, and a bus 304.
[0032] The processing device 302 can execute one or more operations for calibrating a fracture model. The processing device 302 can execute instructions 308 stored in the memory device 306 to perform the operations. The processing device 302 can include one processing device or multiple processing devices.
Non-limiting examples of the processing device 302 include a Field-Programmable Gate Array ("FPGA"), an application-specific integrated circuit ("ASIC"), a microprocessor, etc.
[0033] The processing device 302 can be communicatively coupled to the memory device 306 via the bus 304. The non-volatile memory device 306 may include any type of memory device that retains stored information when powered off.
Non-limiting examples of the memory device 306 include EEPROM, flash memory, or any other type of non-volatile memory. In some aspects, at least some of the memory device 306 can include a medium from which the processing device 302 can read the instructions 308. A computer-readable medium can include electronic, optical, magnetic, or other storage devices capable of providing the processing device 302 with computer-readable instructions or other program code. Non-limiting examples of a computer-readable medium include (but are not limited to) magnetic disk(s), memory chip(s), ROM, RAM, an ASIC, a configured processor, optical storage, or any other medium from which a computer processor can read instructions. The instructions may include processor-specific instructions generated by a compiler or an interpreter from code written in any suitable computer-programming language, including, for example, C, C++, C#, etc.
[0034] FIG. 4 is an illustration of the PKN model 400 of a bi-wing hydraulic fracture where the fracture height H is fixed along the y-axis. The computing device 124 can use the PKN model to estimate a hydraulic fracture length L along the x-axis. The width W of the PKN model of the hydraulic fracture is depicted along the z-axis. In some aspects, the computing device 124 can use a different fracture model to predict the geometry of a different type of hydraulic fracture. For example, the Kristonovich-Geertsma-Daneshy (KGD) model or a discrete fracture network model can be used in some aspects. The computing device 124 can receive additional data to determine the predicted geometry of the hydraulic fracture when using some fracture models. For example, the computing device 124 may receive a proppant concentration value from an input by a user for use in predicting the geometry of the hydraulic fracture using a fracture model.
[0035] FIG. 5 is an example of a flow chart of a process 500 calibrating a fracture model using microseismic event data according to one aspect of the present disclosure. The process 500 below is described with reference to components described above with regard to the computing device 124 shown in FIG. 3. At block 502 the computing device receives microseismic event data from an array of geophones during a hydraulic fracturing operation. In some aspects, the computing device may receive the microseismic event data from other suitable sensors used to monitor microseismic events in a formation during a hydraulic fracture operation, for example tiltmeters.
[0036] At block 504 the computing device determines an observed length of the hydraulic fracture using the microseismic event data received from the geophones. For example, the computing device can determine the location of the microseismic event by calculating the difference between the travel times of the P-waves and S-waves of the microseismic event using additional information about the formation. The additional information can include, for example, the number of rock layers in the formation, the velocity of the P-waves and S- waves through each of the rock layers, and the height of each of the rock layers in the formation.
[0037] At block 506 the computing device determines a predicted length of the hydraulic fracture using a fracture model. For example, the PKN model can be used to determine a predicted fracture length of a bi-wing hydraulic fracture. The computing device can use inputs including the injection rate of the fracking fluid and the viscosity of the fracking fluid to determine the predicted length of the hydraulic fracture based on the fracture model. For example, a bi-wing fracture can be characterized by the Perkins-Kern-Nordgren (PKN) model:
aw(x, t) aq(x, t) +u/(x,t)= 0 at ax
[0038] q(0,t)= q0(t)
[0039] Where w(x,i) is the width of the fracture, q(x,t) is the unit-height injection rate, and ui(xn is the unit-height, unit-length, leak-off rate of the fracturing fluid. The injection rate of the fracking fluid at a time t is go(t) and can be received by the computing device 124 from a user input. The computing device 124 can use Carter's Theory to determine the leak-off rate 14/(x't), which in the PKN
model can be expressed as:

u (x,t) = , ___________
[0040] - r(x)
[0041] Where T(x) is a function of fracture arrival time (i.e., the time location x is first exposed to fluid) and Ci is the leak-off coefficient.
[0042] The computing device 124 can determine the predicted width of the hydraulic fracture in the PKN model by the pressure profile:
ap 2Kq q a " w2"+1
[0043] v v= Mei)
[0044] Where P is the net pressure, n is the flow behavior index of power-law fluid, and w is the shape factor of the fracture, and is the fracture compliance.
The fracture compliance Mc can be determined by known rock mechanics. In some aspects, where n =1 the fluid is Newtonian and the term K becomes the viscosity of the fluid. The computing device can receive the viscosity of the fracking fluid from an input by a user.
[0045] The computing device can determine the predicted length of the hydraulic fracture based on the fracture model using the following global material balance equation:
L(t) r L(ti) W(X,t)dx + u L(x,t1)dx dt, = q0(tI)dti
[0046] o 0 0 0
[0047] The computing device can derive a dynamic fracture model by combining the above equations and using the injection rate go and K . The dynamic fracture model can describe the relationship between g. and K and the predicted length of the hydraulic fracture as the following linear discrete-time-state-space form:
x(k +1) = Akx(k) + B ku(k)
[0048] L(k) = C kx(k)
[0049] This model can be rewritten with a Kalman filter as:
i(k +1) = Ak;(k) + Bku(k) + Jk(y(k) ¨ L(k))
[0050] L(k) = C (k)
[0051] Where u(k) is an input vector that stacks the known input variables together, for example, the injection rate and the fracking fluid viscosity.
Specifically u(k)=[q0(k) K(k)]T The ¶k) is the estimate of the state vector because the state vector x(k) cannot be measured directly. The L(k) is the predicted length of the hydraulic fracture at a time k given the information known about the hydraulic fracture and hydraulic fracturing operations, for example, the injection rate, viscosity of the fracking fluid, and previous observed lengths of the hydraulic fracture up to a time instant k -1. The Y(k) is the observed length of the hydraulic fracture at the time instant k as determined by the computing device using the data received from the sensors monitoring the microseismic event.
[0052] The vector Jk is the Kalman gain and adjusts how much prediction error the computing device uses to correct the fracture model. Specifically, the prediction error is a function 0fY(k)-11(k). In an aspect where the mechanical properties of the fracture are known perfectly, then the only errors present in the fracture model would be from the microseismic monitoring and Jk should be zero to prevent the microseismic measurement errors from be used to update the fracture model. In an aspect where the rock mechanics of the formation are not well known or the fracture propagation is random, A should be a large value since the microseismic results are a source of less error than the fracturing process itself. The value of Jk is given by the following Riccati Equation:
P(k +1) = A [P(k)- P(k)CkT (CkP(k)CkT +Vk)-1 C kP(k)]Akr +Wk
[0053] Jk AkP(k)C kT (C kP(k)C kr +Vk)-1 Where Wk is the covariance matrix of process noise at a time k (i.e., an estimate of randomness and model errors of the fracturing process) and Vk is the covariance matrix of error of the microseismic event location calculated by the microseismic monitoring at time k .
[0054] At block 508 the computing device can determine an uncertainty value of the observed length of the hydraulic fracture. The computing device can determine the uncertainty value of the observed length of the hydraulic fracture by comparing the observed length to the predicted length: AO l'(k). In some aspects, the uncertainty value can be a lateral error of the observed length of the hydraulic fracture (or observed microseismic event location) compared to the predicted length of the hydraulic fracture.
[0055] At block 510 the computing device can determine if the uncertainty value is greater than a maximum uncertainty value. The maximum uncertainty value can be based on an uncertainty of the predicted length of the hydraulic fracture. For example, the computing device can determine the uncertainty of the predicted length E Vk of the hydraulic fracturing using the Riccati Equation: L =
kP(k)CkT + For example, if 95% certainty that the microseismic event measured occurred at or near the fracture tip, then the maximum uncertainty value may be 2IL . If a 99%
certainty is desired, then the maximum uncertainty value may be 31/ .The computing device can receive the maximum uncertainty value from an input by a user.
[0056] The computing device can classify the observed lengths of the hydraulic fracture having an uncertainty value that exceeds the maximum uncertainty value as being associated with a microseismic event that occurred far from the tip of the hydraulic fracture. The computing device can also classify observed lengths have an uncertainty value that does not exceed the maximum uncertainty value as being associated with a microseismic event that occurred near the tip of the hydraulic fracture.
[0057] If at block 510 the computing device determines that the uncertainty value of the observed length of the hydraulic fracture does not exceed the maximum uncertainty value then at block 512 the computing device can update the fracture model using the observed length of the hydraulic fracture. In some aspects, the computing device can update the fracture model using the microseismic event data, from which the observed length of the hydraulic fracture can be inferred.
[0058] At block 514 the computing device can store the observed length of the hydraulic fracture, or in some aspects other data associated with the microseismic event, in a data store. The computing device can store the observed length of the hydraulic fracture with a classification indicating the observed length corresponds to an event near the fracture tip.
[0059] In response to determining the uncertainty value of the observed fracture length exceeds the maximum uncertainty value at block 510 the computing device can store the observed length of the hydraulic fracture in the data store at block 514 and can include a classification indicating the observed length corresponds to a microseismic event that occurred far from the fracture tip.
[0060] At block 516 the computing device determines if the number of observed lengths of the hydraulic fracture stored in the data store corresponding to a microseismic event that occurred far from the fracture tip exceeds a threshold number. The computing device can receive the threshold number from an input by a user. The threshold number can be in relation to the total number of observed lengths stored in the data store.
[0061] If at block 516 the computing device determines the number of stored observed lengths of the hydraulic fracture in the data store that correspond to a microseismic event far from the fracture tip exceeds the threshold value, then at block 518 the computing device can calibrate the fracture model by running a clustering algorithm on the observed lengths in the data store and determining a cluster of observed lengths that has the greatest possibility of representing the fracture growth. For example, the computing device can select the cluster of observed lengths that depict an increase in the fracture length over a period of time.
The computing device can at block 520 calibrate the fracture model using the cluster of observed lengths determined to have the greatest possibility of representing the fracture growth. The calibrated Kalman filter of the model will have new A, B, and C, matrices (k +1) = Ak.i(k)+ Bku(k)+ Jk(y(k)- i(k))
[0062] 1(k) = Ck.Z(k)
[0063] The computing device returns to block 502 after re-estimating the fracture model at block 520.
[0064] In some aspects, the computing device can instead determine an observed height or other geometry of the hydraulic fracture using the microseismic event data and update the fracture model with the observed geometry. For example, when a hydraulic fracture that only grows in a height direction. The computing device can then also determine a predicted or estimated height or other geometry of the hydraulic fracture using the fracture model.
In some aspects, a fracture model for a hydraulic fracture in a wellbore can be updated and calibrated. Information about a microseismic event can be received from sensors that are monitoring a subterranean formation. The information can be received subsequent to a fracking fluid being introduced into the formation.
An observed geometry of a hydraulic fracture can be determined based on the information and a predicted geometry of the fracture can be determined based on properties of the fracking fluid and a fracture model. The fracture model can be updated using the information where it is determined that an uncertainty value of the observed geometry does not exceed a pre-set maximum. The uncertainty value can be based on the predicted geometry of the hydraulic fracture.
[0065] In some aspects, a system can include a computing device that has a non-transitory computer-readable medium that has code executable for causing the computing device to receive information about a microseismic event in a subterranean formation from a plurality of sensors that is monitoring the formation.
The computing device can receive the information subsequent to a fracking fluid being introduced into the formation. The code can also be executable for causing the computing device to determine an observed geometry of a hydraulic fracture based on the information and can predict the geometry of the fracture, using a fracture model, based on properties of the fracking fluid. The code can also be executable for causing the computing device to determine an uncertainty value of the observed geometry based on the predicted geometry and can update the fracture model using the information when the uncertainty value is less than or equal to a pre-set maximum.
[0066] In some aspects, a system can include a plurality of sensors and a computing device. The plurality of sensors can transmit information about a microseismic event in a subterranean formation, the information can be transmitted subsequent to a fracking fluid being introduced into the formation. The computing device determine an observed geometry of a hydraulic fracture based on the information and can predict the geometry of the fracture, using a fracture model, based on properties of the fracking fluid. The computing device can also determine an uncertainty value of the observed geometry based on the predicted geometry and can update the fracture model using the information when the uncertainty value is less than or equal to a pre-set maximum.
[0067] The foregoing description of certain examples, including illustrated examples, has been presented only for the purpose of illustration and description and is not intended to be exhaustive or to limit the disclosure to the precise forms disclosed. Numerous modifications, adaptations, and uses thereof will be apparent to those skilled in the art without departing from the scope of the disclosure.

Claims (20)

Claims What is claimed is:
1. A method, comprising:
introducing a fracturing fluid into a subterranean formation;
receiving, by the computing device, information about a microseismic event in the subterranean formation from at least one sensor that is monitoring the subterranean formation;
determining, by the computing device, an observed geometry of a hydraulic fracture based on the information about the microseismic event;
determining, by the computing device, a predicted geometry of the hydraulic fracture based on properties of the fracking fluid using a fracture model;
determining, by the computing device, an uncertainty value of the observed geometry based on the predicted geometry;
determining, by the computing device, the uncertainty value is less than or equal to a pre-set maximum; and updating, by the computing device, the fracture model using the information about the microseismic event.
2. The method of claim 1, further comprising:
storing, by the computing device, the observed geometry and a classification of the observed geometry in a data store, the classification indicating if the uncertainty value associated with the observed geometry exceeds a pre-set maximum;

determining, by the computing device, a cluster of observed geometries that represent a growth of the hydraulic fracture; and calibrating, by the computing device, the fracture model using the cluster of observed geometries.
3. The method of claim 1, further comprising:
generating, by the computing device, a visual representation of the hydraulic fracture based on the fracture model.
4. The method of claim 1, wherein the fracture model is the Perkins-Kern-Nordgren model.
5. The method of claim 1, wherein the observed geometry is an observed length of the hydraulic fracture and the predicted geometry is a predicted length of the hydraulic fracture.
6. The method of claim 1, wherein the at least one sensor is positionable within a wellbore.
7. The method of claim 1, wherein the properties of the fracking fluid include at least one of an injection rate of the fracking fluid and a viscosity of the fracking fluid.
8. A system comprising:
a computing device having a non-transitory computer-readable medium on which is code that is executable for causing the computing device to:

receive information about a microseismic event in a subterranean formation from at least one sensor that is monitoring the subterranean formation, the information about the microseismic event being received subsequent to a fracking fluid being introduced into the subterranean formation;
determine an observed geometry of a hydraulic fracture based on the information about the microseismic event;
determine a predicted geometry of the hydraulic fracture based on properties of the fracking fluid using a fracture model;
determine an uncertainty value of the observed geometry based on the predicted geometry;
determine the uncertainty value is less than or equal to a pre-set maximum; and update the fracture model using the information about the microseismic event.
9. The system of claim 8, wherein the code is executable for causing the computing device to:
store the observed geometry and a classification of the observed geometry in a data store, the classification indicating if the uncertainty value associated with the observed geometry exceeds the pre-set maximum;
determine a cluster of observed geometries that represent a growth of the hydraulic fracture; and calibrate the fracture model using the cluster of observed geometries.
10. The system of claim 8, wherein the code is executable for causing the computing device to:

generate a visual representation of the hydraulic fracture based on the fracture model.
11. The system of claim 8, wherein the properties of the fracking fluid include at least one of an injection rate of the fracking fluid and a viscosity of the fracking fluid.
12. The system of claim 8, wherein the observed geometry is an observed length of the hydraulic fracture and the predicted geometry is a predicted length of the hydraulic fracture.
13. The system of claim 8, wherein the at least one sensor is positionable within a wellbore.
14. The system of claim 8, wherein the at least one sensor is positionable at a surface of the subterranean formation.
15. A system comprising:
at least one sensor for transmitting information about a microseismic event in a subterranean formation, the information about the microseismic event being transmitted subsequent to a fracking fluid being introduced into the subterranean formation; and a computing device for:
determining an observed geometry of a hydraulic fracture based on the information about the microseismic event, a predicted geometry of the hydraulic fracture based on properties of the fracking fluid and a fracture model, and an uncertainty value of the observed geometry based on the predicted geometry;
determining that the uncertainty value of the observed geometry is less than or equal to a pre-set maximum; and updating the fracture model using the information about the microseismic event.
16. The system of claim 15, the computing device also for:
storing the observed geometry and a classification of the observed geometry in a data store, the classification indicating if the uncertainty value associated with the observed geometry exceeds the pre-set maximum;
determining a cluster of observed geometries that represent a growth of the hydraulic fracture; and calibrating the fracture model using the cluster of observed geometries.
17. The system of claim 15, the computing device also for:
generating visual representation of the hydraulic fracture based on the fracture model.
18. The system of claim 15, wherein the observed geometry is an observed length of the hydraulic fracture and the predicted geometry is a predicted length of the hydraulic fracture.
19. The system of claim 15, wherein the at least one sensor is a tiltmeter.
20. The system of claim 15, wherein the at least one sensor is a geophone.
CA2964862A 2014-11-19 2014-11-19 Filtering microseismic events for updating and calibrating a fracture model Active CA2964862C (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2014/066385 WO2016080980A1 (en) 2014-11-19 2014-11-19 Filtering microseismic events for updating and calibrating a fracture model

Publications (2)

Publication Number Publication Date
CA2964862A1 CA2964862A1 (en) 2016-05-26
CA2964862C true CA2964862C (en) 2019-11-19

Family

ID=56014336

Family Applications (1)

Application Number Title Priority Date Filing Date
CA2964862A Active CA2964862C (en) 2014-11-19 2014-11-19 Filtering microseismic events for updating and calibrating a fracture model

Country Status (3)

Country Link
US (1) US10408955B2 (en)
CA (1) CA2964862C (en)
WO (1) WO2016080980A1 (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11125912B2 (en) * 2013-11-25 2021-09-21 Schlumberger Technology Corporation Geologic feature splitting
US10429528B2 (en) 2014-11-19 2019-10-01 Halliburton Energy Services, Inc. Reducing microseismic monitoring uncertainty
WO2016105351A1 (en) 2014-12-23 2016-06-30 Halliburton Energy Services, Inc. Microseismic monitoring sensor uncertainty reduction
US10509141B2 (en) * 2015-08-17 2019-12-17 Schlumberger Technology Corporation Method and apparatus for determining a fracture aperture in a wellbore
US11572777B2 (en) * 2019-01-28 2023-02-07 Landmark Graphics Corporation Constructing digital twins for oil and gas recovery using Ensemble Kalman Filter
CN110671089B (en) * 2019-10-22 2021-09-10 中海油田服务股份有限公司 Fracturing filling process parameter optimization design method
CN110910613B (en) * 2019-12-10 2022-04-05 大连理工大学 Rock micro-seismic wireless monitoring, receiving and early warning system
CN110942610A (en) * 2019-12-10 2020-03-31 大连理工大学 Intelligent rock microseismic acquisition and data wireless transmission system
US11143019B2 (en) 2020-03-03 2021-10-12 Halliburton Energy Services, Inc. Real time estimation of fracture geometry from the poro-elastic response measurements
CN111720104B (en) * 2020-08-04 2022-03-11 西南石油大学 Method for predicting fracture morphology of fractured reservoir multistage fracturing
CN112127882B (en) * 2020-11-02 2021-05-25 西南石油大学 Method for calculating dynamic fracture width of drilling fluid leakage of fractured formation

Family Cites Families (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2747172A (en) 1954-11-12 1956-05-22 Exxon Research Engineering Co Geophone arrangement for seismic prospecting
US6814160B1 (en) 1999-07-03 2004-11-09 Westerngeco, Llc Automated seismic geophone deployment
US6795773B2 (en) 2001-09-07 2004-09-21 Halliburton Energy Services, Inc. Well completion method, including integrated approach for fracture optimization
US7272973B2 (en) * 2005-10-07 2007-09-25 Halliburton Energy Services, Inc. Methods and systems for determining reservoir properties of subterranean formations
RU2324810C2 (en) 2006-05-31 2008-05-20 Шлюмберже Текнолоджи Б.В. Method for determining dimensions of formation hydraulic fracture
GB2439571B (en) * 2006-06-28 2008-11-12 Schlumberger Holdings Method for updating a model of the earth using microseismic measurements
US7451812B2 (en) * 2006-12-20 2008-11-18 Schlumberger Technology Corporation Real-time automated heterogeneous proppant placement
US7516793B2 (en) * 2007-01-10 2009-04-14 Halliburton Energy Service, Inc. Methods and systems for fracturing subterranean wells
US8148990B2 (en) 2007-04-30 2012-04-03 Kjt Enterprises, Inc. Marine electromagnetic acquisition apparatus with foldable sensor arm assembly
US7724608B2 (en) 2007-07-20 2010-05-25 Wayne Simon Passive reflective imaging for visualizing subsurface structures in earth and water
US20090316524A1 (en) 2008-06-23 2009-12-24 Stig Rune Tenghamn Flexible seismic data acquisition system for use in a marine environment
US20100252268A1 (en) 2009-04-03 2010-10-07 Hongren Gu Use of calibration injections with microseismic monitoring
US8612155B2 (en) 2009-04-08 2013-12-17 Schlumberger Technology Corporation Methods and systems for microseismic mapping
US8498852B2 (en) 2009-06-05 2013-07-30 Schlumberger Tehcnology Corporation Method and apparatus for efficient real-time characterization of hydraulic fractures and fracturing optimization based thereon
US8437962B2 (en) 2009-11-25 2013-05-07 Halliburton Energy Services, Inc. Generating probabilistic information on subterranean fractures
US9176245B2 (en) 2009-11-25 2015-11-03 Halliburton Energy Services, Inc. Refining information on subterranean fractures
US9410421B2 (en) 2009-12-21 2016-08-09 Schlumberger Technology Corporation System and method for microseismic analysis
RU2455665C2 (en) * 2010-05-21 2012-07-10 Шлюмбергер Текнолоджи Б.В. Method of diagnostics of formation hydraulic fracturing processes on-line using combination of tube waves and microseismic monitoring
US9305121B2 (en) 2010-06-28 2016-04-05 Exxonmobil Upstream Research Company Method and system for modeling fractures in ductile rock
US9618652B2 (en) 2011-11-04 2017-04-11 Schlumberger Technology Corporation Method of calibrating fracture geometry to microseismic events
CA2915625C (en) 2011-03-11 2021-08-03 Schlumberger Canada Limited Method of calibrating fracture geometry to microseismic events
FR2975782B1 (en) 2011-05-26 2014-02-21 Cggveritas Services Sa SYSTEM AND METHOD FOR DOPPLER POSITIONING OF SEISMIC SENSORS
WO2013134427A2 (en) * 2012-03-06 2013-09-12 Ion Geophysical Corporation Model predicting fracturing of shale
EP2831647A4 (en) 2012-03-28 2016-02-24 Geco Technology Bv Providing an objective function based on variation in predicted data
US9417348B2 (en) * 2012-10-05 2016-08-16 Halliburton Energy Services, Inc. Updating microseismic histogram data
US10429528B2 (en) 2014-11-19 2019-10-01 Halliburton Energy Services, Inc. Reducing microseismic monitoring uncertainty
WO2016105351A1 (en) 2014-12-23 2016-06-30 Halliburton Energy Services, Inc. Microseismic monitoring sensor uncertainty reduction

Also Published As

Publication number Publication date
US10408955B2 (en) 2019-09-10
US20170343689A1 (en) 2017-11-30
CA2964862A1 (en) 2016-05-26
WO2016080980A1 (en) 2016-05-26

Similar Documents

Publication Publication Date Title
CA2964862C (en) Filtering microseismic events for updating and calibrating a fracture model
CA2964863C (en) Reducing microseismic monitoring uncertainty
CN112593907B (en) System and method for calculating fracture area, volume and fluid loss rate, and program product
US9612359B2 (en) Generation of fracture networks using seismic data
US10280731B2 (en) Energy industry operation characterization and/or optimization
NO20150147A1 (en) Method to predict overpressure uncertainty from normal campaction trendline uncertailty
CA3080938C (en) Fracture length and fracture complexity determination using fluid pressure waves
AU2013293261B2 (en) Stratigraphic modeling using production data density profiles
US11879317B2 (en) Flow rate optimization during simultaneous multi-well stimulation treatments
US20200291774A1 (en) Determining Fracture Surface Area in a Well
AU2015406994B2 (en) Predicting wellbore operation parameters
US10338247B2 (en) Microseismic monitoring sensor uncertainty reduction
US20150331976A1 (en) Forecasting oil production for multiply-fractured horizontal wells
CA3097884C (en) Real time parent child well interference control
US11726230B2 (en) Subsurface strain estimation using fiber optic measurement
AU2021200537B2 (en) Subsurface strain estimation using fiber optic measurement
CA3107445C (en) Subsurface strain estimation using fiber optic measurement
Khasanov et al. A geomechanical approach to minimising sanding risk (Russian)

Legal Events

Date Code Title Description
EEER Examination request

Effective date: 20170418