Unconventional Well Test Analysis for Assessing Individual Fracture Stages through Post-Treatment Pressure Falloffs: Case Study

: Researchers and operators have recently become interested in the individual stage optimization of unconventional reservoir hydraulic fracture. These professionals aim to maximize well performance during an unconventional well’s early-stage and potential Enhanced Oil Recovery (EOR) lifespan. Although there have been advances in hydraulic fracturing technology that allow for the creation of large stimulated reservoir volumes (SRVs), it may not be optimal to use the same treatment design for all stages of a well or many wells in an area. We present a comprehensive review of the main approaches used to discuss applicability, pros and cons, and a detailed comparison between different methodologies. Our research outlines a combination of the Diagnostic Fracture Injection Test (DFIT) and falloff pressure analysis, which can help to design intelligent production and improve well performance. Our field study presents an unconventional well to explain the objective optimization workflow. The analysis indicates that most of the fracturing fluid was leaked off through natural fracture surface area and resulted in the estimation of larger values compared to the hydraulic fracture calculated area. These phenomena might represent a secondary fracture set with a high fracture closure stress activated in neighbor stages that was not well-developed in other sections. The falloff pressure analysis provides significant and vital information, assisting operators in fully understanding models for fracture network characterization.


Introduction
The advent of unconventional reservoir development is a turning point in the global oil and gas industry, since these resources contain massive hydrocarbon reserves larger than those found in conventional formations.Domestic oil production from liquid-rich shale (LRS) reservoirs in North America has seen significant development, according to the US Energy Information and Administration (EIA), with production dramatically increasing in the 'top producing' American oil shale plays: the Bakken, Eagle Ford, and Permian Basin.This hydrocarbon production improvement has been driven by the application of modern horizontal drilling and multi-stage hydraulic fracturing (MSHF) techniques, which makes it possible to access low-porosity (<10%) and low-permeability (<0.1 mD) formations [1,2].The creation of large stimulated reservoir volumes (SRVs) has been achieved through breakthroughs in hydraulic fracturing technology; however, fracture treatment is not necessarily effective [3].Operators have started utilizing tighter spaced clusters, longer stage lengths, and greater proppant volumes in their stimulation designs [4,5].The ultimate oil recovery reported by several studies is less than 8% due to a rapid decline in unconventional well performance, and by approximately 75% within the first two years of well production (Figure 1).The decline is due to a low to no hydrocarbon recharge from the ultra-tight matrix blocks since the natural and induced fractures close, and there is a high flow resistance at the matrix-fracture interface; therefore, the increase in net stress leads to a zero-pressure gradient, which obstructs the fluid flow from the rock matrix into the fracture [3,6,7].Proppant choice plays a significant role in fracture conductivity.An inappropriate choice may cause proppant deformation and/or crush under closure pressures.High-pressure drawdown may also cause formation rock compression, which leads to a reduction in matrix permeability with changes in reservoir pressure or stress [3,4,8].The productivity of unconventional wells heavily relies on the effective fracture contact area or the propped fracture area per cluster, which is critical in evaluating hydraulic fracture treatment performance [9].Performing a successful treatment design depends primarily on evaluating the stimulated formation before and during fracture treatment.Recommendations for the proper selection of slurry type and amount are provided to produce the optimal fracture geometry using a high-efficiency assessment [11].Current literature suggests that the use of diagnostic tools is critical when assessing unconventional plays since they can reveal opportunities for future exploration, evaluation, delineation, and development [12]; however, currently available traditional methods are not feasible for decision-makers, even with the appropriate data.The majority of the previously published research is associated with high uncertainty and lacks a thorough discussion of fracture network characterizations.Only long-term production data were utilized in the past, the application of which caused ambiguity in understanding the hydraulic fracture performance [13,14].This ambiguity was caused by the difficulty in utilizing conventional exploration and production techniques to establish commercial production rates [3,15].
Most of the previous post-treatment analysis studies of unconventional reservoirs do not provide quantitative discussions with detailed support case studies, and very few field studies have discussed the application of pressure falloff data [9].We present guidelines to better understand unconventional well test analysis through critical literature reviews and case studies of fracturing treatment analysis to address the lack of quantitative evaluation.
We have concluded that post-stimulation condition performance evaluation using indirect methods, such as pressure falloff data, is the most promising technology.This approach may provide a clear perspective about the created fracture's dimensions and properties; therefore, the effective fracture contact surface area for both natural and induced fractures can be determined from stage to stage during fracture treatment jobs.The application of pressure falloff data is a valuable tool that provides comprehensive information, such as the mechanics of the created open, closed, and propped hydraulic fractures, due to the tool's capability of reflecting the rock and fluid's physical behavior.This technology may overcome some limitations and weaknesses in most of the proposed techniques in the literature, such as production data analysis and micro-seismic methods.

Research Objectives
The main objective of diagnostic tool usage is to maximize well performance during both primary production and late-stage EOR processes by understanding the contribution of individual fracture stages [1,8,13].Important pressure falloff data gathered after fracture treatments were generally ignored in the past, even though that detailed information could reveal attributes necessary for successful fracturing evaluation.There were also no precise measurement technologies for recording production rate and pressure during each stage to evaluate each cluster or stage's contribution.The pressure falloff data tool has garnered significant interest in the oil and gas industry since it is a powerful method for defining prime fracture parameters to gain insight into the effectiveness of the treatment jobs' fracture [9].
Evaluating fracture designs can be implemented on a stage-by-stage basis to optimize the overall performance of a well, unlike common performance evaluation methods such as Rate Transient Analysis (RTA) and micro-seismic fracture imaging/mapping (MS).We will investigate the following key parameters to evaluate these designs: 1.The importance of closure stress, including closure behavior and geomechanics parameters, in fracture geometry and the impact of its variation from stage to stage.2. The variability of geology and mineralogy and their impacts on fracture propagation and geometry.3. The contribution of natural fractures in created well-reservoir contact areas and their impact on fracture geometry.4. The role of treatment designs on proppant distribution and characteristics, such as conductivity, crush, and embedment, which are related to closure stress and effective stress.5. Investigate if an optimal fracture contact area would exist for a specified treatment design.6. Identify optimal treatment design parameters for individual fracture stages, such as proppant volume and mesh size, fracturing fluid amount and type, and pump schedule.7. Prepare protocols for evaluating the success of individual fracture stages from the viewpoint of production performance, such as post-treatment analysis.8. Determine the contribution of individual fracture stages in the well's overall performance and determine the right spot for treatment execution.9. Evaluate the lessons learned from fracturing in the previous stages.Develop a workflow for the real-time treatment design optimization for next-stage application to address good or bad frac-hits, unsuccessful designs, and fracturing in undesired formations.10.Determine if the expected treatment design optimization by stage would justify the additional cost and treatment design adjustments through case studies.
We have compared common fracture diagnostic tools, discussed DFIT and pressure falloff data, and analyzed fracturing treatment case studies fracturing for unconventional wells.We have also presented the potential of combining DFIT and pressure falloff data in various pressure-time plots to identify fracture and reservoir behavior characteristics.Guidelines to better understand unconventional well test analysis that can lead to realtime optimization and adjustment of fracture job treatments have been provided as we proceed from one stage to another in an MSHF operation.

Principle of Fracture Diagnostic Tools
Monitoring the growth of fracture networks in the subsurface is the common process during stimulation treatments in unconventional wells.The created fractures are usually simulated by simple bi-wing and single planar fracture model definitions, such as Perkins-Kern-Nordgren (PKN) and Khristianovic-Geertsma-de Klerk (KGD) fracture geometry models.These models assume that the hydraulic fractures will primarily stay within the pay zone and extend significantly [16]; however, in reality, the fracture networks can grow in an asymmetrical shape due to variable confinement across the geologic interfaces and orientation changes.Fracture growth around natural fractures can add more complexities to the fracture system due to the interaction between the induced and fissure networks [17].An effective fracture network's morphology in shale resources is still largely unknown.As a result of this knowledge gap and the network's complexity, it is difficult to predict, obtain, and verify fracture geometry, such as fracture length, height, and containment.These challenges can lead to suboptimal outcomes when incorrect assumptions are used through fracture diagnostic applications [11].
There are several fracture diagnostic applications.These applications allow petroleum engineers to assess the success of the fracture stages and create optimal development strategies for effective reservoir drainage.These tools can also optimize the entire field development regarding well spacing and location, optimal design, and optimum interval/height coverage.Numerous fracture diagnostics are discussed by researchers that provide subsets of knowledge about treatment design optimization [13,18,19]: (a) direct farfield fracture diagnostic techniques, such as micro-seismic fracture mapping and tiltmeter, (b) direct near-wellbore fracture diagnostic techniques, such as production and temperature logging tools and radioactive tracers, and (c) indirect fracture-diagnostic techniques, such as transient pressure and rate transient analysis "PTA/RTA" and fracture modeling "net pressure analysis".
We have focused on indirect fracture-diagnostic techniques in this review paper.Table 1 provides a brief discussion and comparison reference for widely used diagnostic tools with their strengths and limitations explained.
The analysis of source rock and fluid behaviors is detected by the fracture diagnostic tools, highlighting the fracture and reservoir properties; therefore, the combination of diagnostic tools provides more confidence and allows fracture engineers to make decisions in real-time.DFIT, post-treatment pressure falloff, micro-seismic, flowback, and other data can be collected and interpreted in real-time to assess the created contact surface areas and evaluate fracture design properties, such as fracture half-length, number of clusters, proppant loading, and fracture complexity and direction.The pump schedule can be adjusted from stage to stage by assessing the proppant placement and injection volumes to ensure the maximum pay zone proppant coverage.An optimal well stimulation strategy should be established to avoid some far-field issues, such as well interventions and frac-hits [11].
Applying the DFIT methodology and using post-treatment pressure falloff data as a fracture-diagnostic tool may assist us in monitoring the pressure interference and offset well intervention based on the survey comparison mentioned in Table 1.This technology will provide valuable information needed to improve individual fracture stage treatment design and enhance the production of the propped fracture surface area in real-time.

Diagnostic Fracture Injection Test (DFIT)
DFITs have been widely used in the oil and gas industry over the last 20 years.This approach is based on the pressure transient data procured immediately after fracture treatments to obtain a reliable assessment of fracture and reservoir properties.A DFIT is a standard well testing technique for ultra-low-permeability formations, where a traditional pressure transient test, such as a buildup test, is impractical.A DFIT analysis does not require a long shut-in time to reach the radial flow regime.
The pressure falloff data are analyzed in a DFIT to estimate the in-situ stress, fluid efficiency, leak-off coefficient, reservoir properties, and net pressure, which are the critical factors used to design and implement a successful main fracture treatment.A DFIT provides the representative properties of an undamaged formation since the test creates a large area of investigation that can extend beyond the damaged near-wellbore zone.The results from a DFIT can be used in several ways: (a) characterizing in-situ stresses and fracture compliance [20], (b) modeling hydraulic fracture propagation [21], (c) designing fracture treatment jobs [22], (d) modeling reservoir simulation [23], and (e) post-fracture treatment analysis [24].This application can also be utilized in geologic carbon sequestration, nuclear waste repositories, and geothermal energy exploitation [17,25,26].
A DFIT implementation and interpretation should be studied in detail before any field execution or analysis.We have presented lessons learned and recommendations that may help operators design the optimum fracture jobs.

DFIT Design and Tactics
A typical DFIT operation pumps a small volume of the treatment fluid, such as water, without proppant at a constant rate for a short period of approximately 3-5 minutes.The injection pressure increases above the reservoir fracturing pressure, or breakdown pressure, creating a short artificial fracture in the target layer.The leak-off behavior is small during the injection period, and no filter cake forms on the fracture wall.The pressure falloff data is a recorded function during shut-in time immediately after the treatment.The injected fluid begins to leak-off into the formation until the fracture wall comes into contact, also called closure.The pressure falloff period recorded in the case studies we analyzed was extended for days or even weeks to observe the radial flow regime, depending on the reservoir characteristics.
Figure 2 displays a typical pressure profile of a DFIT in the absence of natural fractures and weak planes.Two distinct periods of before and after closure (BC and AC) were analyzed to characterize the properties of the created fractures and reservoirs.These two periods were separated by the fracture closure event, which is the primary outcome and supplies us with the fracture closure pressure, or minimum in-situ stress.The fracture and reservoir properties were obtained by analyzing the typical flow regimes observed during the BC, before closure for fracture-dominated, and AC, after closure for reservoir-dominated, respectively (Figure 3).
The conditions of modeling in unconventional reservoirs are more complex than in conventional formations.Significant conflict exists on how to model the fracture closure behavior since most of the models assume that the fracture surface is perfectly smooth; however, fractures exist everywhere in the subsurface in the form of small-scale cracks and fissures, and large-scale joints and faults.The mechanical resistance and fluid transport properties found in unconventional formations are complicated and controlled by several factors, such as in-situ stress, compliance or stiffness, rock mineralogy, fracturesurface roughness, treatment fluid pressure inside the fractures, and leak-off rate [11,20,27].A proper DFIT model must account for the effects of pressure-dependent leakoff and dynamic fracture compliance to precisely capture the fracture pressure response and obtain a realistic estimation of the fracture closure pressure and fluid leak-off behavior [20,28].These parameters are crucial factors necessary for calculating the fracture surface contact areas to obtain proper hydraulic fracture modeling and accurate post-treatment assessment.
Investigating compelling field evidence using a downhole measurement that indicates what exactly occurs in the subsurface is paramount.Figure 4 illustrates the tiltmeter measurements of well #2B from the Gas Research Institute/Department of Energy M-site, which is defined as an indication of fracture displacement or fracture width [29] (Table 1).The data are readily available at low costs, and their analysis is straightforward.

•
This application is critical for using historical production data for various purposes: a) Characterize reservoir and well stimulation properties, b) Predict production performance for development plans and reserve estimations.

•
Production data analysis is useful for checking the consistency of the data and identifying the flow regimes over time.

•
There are a couple of ways to apply production data techniques: a) Type curve analysis, b) Straight line (flow regime) analysis, c) Analytical and numerical simulation, d) Empirical methods to quantify the hydraulic fracture performance at different unconventional well life stages.

•
Production analysis is a tool to study linear flow regimes when assessing the productivity and effectiveness of completion designs.The impact of fracture-hit and offset stages can also be evaluated.

•
High-frequency measurements are required to obtain a reasonable analysis.

•
The results suffering from some degree of uncertainty; therefore, more information is required, such as geology, phase behavior, and completion practices.This information needed is more complex compared to other approaches.

•
Traditional production data analysis assumes that all clusters are similar, which appears to be an oversimplification of the problem.

•
All production data analysis methods are non-uniqueness associated with well and reservoir properties; therefore, the different methods must be validated and cross-checked.

•
The production data approaches for conventional rate transient methods need further modifications to analyze the MSHW data.

•
Production data analysis does not offer any procedures for quick and measured adjustments of completion designs and production operations.
Micro-Seismic Fracture Imaging/Mapping Pros Cons Results Authors • Micro-seismic imaging provides the best resolution and lowest uncertainties when characterizing fracture geometries in most cases.

•
Micro-seismic mapping can be coupled with real-time simulations to accurately predict fracture growth in the target zone.This combination can be utilized to assess the effectiveness of flushing an unexpected screen-out, imaging proppant placement, synthetic micro-seismic event prediction, and fracture geometry control.

•
This technology is not widely used due to its high cost.

•
This application is associated with a few uncertainties due to source mechanisms: a) Receiver-coupling resonances, such as improper sensor couplings with rock properties, including velocity-model limitations, formation anisotropy, noise, and mislocation.

b)
The interference of fluid leak-off and stress effects in some formations, such as shale plays.

•
Fracture geometries obtained from micro-seismic monitoring may not be accurate enough in certain situations; therefore, the geometries must be validated and cross-checked with other diagnostic tools.

•
The extent of fractures parallel to the lateral might be difficult to interpret from micro-seismicity alone, such as in the case of micro-seismicity vs. dimensions.

•
It is a simple approach that only requires shut-in pressure vs. time.

•
The data are recorded right after fracture treatments with no additional cost.

•
It is a unique method to assess the state of created fractures, whether they are open, closed, or propped.

•
It does not require long periods of shut-in data, as required in conventional buildup tests.

•
A half-hour period for pressure recordings after the main fracturing treatment would be enough to reliably obtain falloff pressure analysis and diagnostic plots to determine: a) Effective contact of natural fracture surface area, b) Effective contact of induced fractures surface area.

•
The application of fracture falloff pressure analysis is still in its infancy, and its theory is based on the assumptions of DFIT: a) The whole surface of created fractures/cracks would contribute to the fluid flow.Further research is necessary to improve the current effective well/fractures contact area workflow estimation fracture.

•
We recommend the inclusion of proppant-impact-factors (PIFs) in the estimation of the effective contact areas to improve the DFIT assumptions.

•
It may not be used independently as an accurate diagnostic tool to assess fracture network growth.It is recommended that this tool be integrated with other methods, such as: a) Micro-seismic, where the fracture height is required for area estimation, given by MS, b) Production logging data, c) Fiber-optic information, • This technique can be used at any time during the life of the well and can also be integrated with intelligent production studies.

•
This technique may facilitate the real-time evaluation of the created contact areas between the well and reservoir at each fracture stage, representing the success of a fracturing treatment job.

•
The method provides the flow regimes and fracture/reservoir behaviors before and after closure.

•
Evaluating individual fracture stages through this technique can assist in stage-by-stage treatment job optimization and further design improvements.

•
This technique has not yet been used for real-time fracturing treatment optimizations.

•
A low-pressure gauge resolution can cause inaccurate diagnostic plots and unreliable results.We recommend a measurement of the shut-in pressure period at a time interval of one second using high-resolution pressure gauges, which can assist in solving data quality issues and minimize the effects of water hammers.Wang and Sharma [20] used recorded data to explain the relationship between normalized tilt and formation pressure (Figure 4).The data for this downhole field measurement were recorded and gathered immediately after the end of the several week-long test.The Y-axis presents the normalized tilt calculated by dividing each tiltmeter measurement by the maximum fracture displacement.The data were plotted vs. the wellbore pressure to generate the diagnostic plot (Figure 4), demonstrating a direct measurement of the rock deformation during the fracture closure behavior.The results indicate that the pressure falloff data response on a normalized tilt vs. formation pressure plot (Figure 4) is proportional to the fracture compliance Equation (1), or inversely proportional to fracture stiffness, and the fracture closure behavior is a function of average displacement and fracture volume.
The fracture volume is proportional to the average fracture width as the pressure continues to decrease, and two distinct periods are indicated on the diagnostic plot: (1) The trend of the pressure falloff data is a linear decline until the point of measured pressure inside the fracture, or closure pressure, is greater than 21 MPa (3046 psi).At this point, the fractures are still open, the stiffness factor is constant, and the surface area remains constant until the fracture wall comes into contact, or closure.The fracture geometry can then be estimated directly using Table 2. (2) The pressure falloff data begin to deviate from a straight line at the inflection point on the plot, where the closure pressure is marked with a dashed green line; therefore, the fracture stiffness increases gradually, as a result of the fracture closure on the asperities of the fracture edges or tips.This information coincides with downhole measurements but is not a parameter in the available DFIT models; therefore, this application assists us with measuring the appropriate closure pressures and provides more information on the mechanics of the created hydraulic fractures, such as if they are open or closed.This information may assist us in developing DFIT assumptions by adding proppant-impact-factors (PIFs), which can be used in the post-treatment pressure falloff data analysis to estimate an effective contact fracture surface area, such as the propped fracture area per cluster

Fundamentals of DFIT
The leak-off behavior term was introduced in Nolte's work (1979,1986) [35,36], where he pioneered DFIT as a reliable test method before executing main fracture jobs.A poroelastic closure model was used to describe the pressure falloff behavior when fracturing fluid leak-off entered the fractures and formations.Equations ( 2)-( 10) are used for analyzing a DFIT to capture normal leak-off behavior.This analysis is based on the following assumptions, which are assumed in several models: (1) power-law fracture growth, (2) negligible spurt loss, (3) constant fracture surface area immediately after the end of the test if there is constant leak-off area and constant fracture compliance or stiffness, and (4) Carter's leak-off model, which defines one-dimensional fluid leak-off across a constant pressure boundary.The leak-off behavior is not pressure-dependent, and the solution to the diffusivity equation predicts that the leak-off rate will scale with the inverse of the square root of time.
These assumptions may be realistic due to the characterization of the fracturing fluids and unconventional formations; therefore, Nolte's technique may not work for unconventional formations and may yield overestimated results in parameters such as fluid efficiency, leak-off coefficient, and storage coefficient.This technique is a reliable application under some circumstances and is derived based on the G-function approach (Equation (2)).The pressure and G-function time are analyzed on log-log graphs to obtain the fracture closure point and other parameters, such as fluid efficiency and leak-off coefficient.
where,  is the pressure at the end of pumping,  is the pressure recorded at the surface during the falloff period, and  is the dimensionless time defined by: is the superposition time defined by: With pumping time  , falloff period , productive fracture ratio  = ℎ/ℎ , fracture height ℎ , propped height ℎ, leak-off coefficient  , and fracture compliance  ,  * :  / at the closure point.
Castillo [37] later introduced a new G-function plot to address the assumption of pressure-dependent behavior by linearizing the relationship between the pressure falloff and time during the closure behavior.The relation is the time, such as the G-function time and the square root of time, on the X-axis vs. falloff pressure and the pressure derivative on the Y-axis.This approach reduces the uncertainty of estimating fracture fluid efficiency and the leak-off coefficient, while overestimated outcomes have resulted from Nolte's method.The diagnostic plot estimates accurate pressure parameters, such as instantaneous shut-in pressure "ISIP," closure pressure "Pc," and Nolte match pressure " * ".
Barree and Mukherjee [38] presented several types of abnormal leak-off behaviors: (a) natural fracture opening or pressure-dependent leak-off, (b) fracture tip extension or recession, (c) height recession, (d) pressure-dependent fracture compliance, and (e) transient flow in the fracture.The authors developed Nolte's work (1979,1986) [35,36] for various closure behavior types and removed the ambiguity associated with understanding the complex fracture networks.The diagnostic plots allow us to predict an accurate estimation of the in-situ stress, fluid efficiency, leak-off coefficient, and pressure parameters.The G-function plot in the proposed models is the relationship between the following terms: falloff pressure,  , first pressure derivative,  /, and second pressure derivative,  /, on the Y-axis vs. the G-function time on the X-axis.The closure event is the point where the curve deviates from the straight line.Equations (11) through (13) are used to determine the primary outcomes, as listed in Table 3.This method enables us to identify the proper leak-off behavior for accurately estimating the hydraulic fracture geometries.
where,  , closure time, ′, Young's modulus,  , total pumping volume, and  , the ratio of the average net pressure inside fractures to the maximum net pressure at the wellbore during the shut-in time (Table 2).
A DFIT uses the basis of conventional mini-fracture treatments that focus on acquiring treatment design parameters, such as fluid efficiency and leak-off behavior; however, this application is subtly different for unconventional formation analysis.This approach is used to acquire significantly more information on the created fractures and formation properties, such as pore pressure, closure and fracture gradients [11,20,28], process zone stresses [6], transmissibility values [39], leak-off mechanisms [6], natural fracture properties [6], fracture stiffness and un-propped fracture conductivity as a function of closure stress [20], and stimulation complexity and net pressure [40].
We can evaluate the properties of the main hydraulic fractures and natural fractures in a fracture treatment job by adopting the DFIT analysis method with no proppant.This method may not be ideal due to its tendency to ignore the impact of the proppant; however, it can be used with caution to evaluate post-treatment production and unconventional well performance.

DFIT Models: Before-Closure Analysis
Several analytical or semi-analytical models have been proposed for before-closure analysis (BC) [35,36,[41][42][43], where all BC models were based on two fundamental concepts underlying the proposed methodology: (1) a material balance equation before fracture closure, and (2) the diffusive flow in the formation after closure.These models were founded based on two main conditions: (1) the total injection volume is equal to the sum of the fracture volume and cumulative leak-off volume, and (2) the fracture volume is estimated from the linear elasticity theory and a 2D fracture geometry during the pressure falloff period.
Cramer and Nguyen [44] reported that it would be rare to observe a normal leak-off behavior in the field, and closure behavior is commonly related to the abnormal leak-off concepts (Figure 5); therefore, the BC analysis must correctly address abnormal leak-off behaviors and near-wellbore friction losses.Liu and Ehlig-Economides [28] presented a model that was not limited to normal leak-off behavior compared to previous BC models [20,[41][42][43].These BC models relied on the assumptions of ideal leak-off behavior: (1) constant injection rate, (2) constant fracture surface area after shut-in, (3) creation of one main hydraulic fracture cluster without the effects of natural fractures, (4) constant fracture compliance during the operation, ( 5) assumption of similar fracture closure stresses for all stages, and ( 6) assumption that all injected fluid at the surface flows into the created fracture, meaning that the impact of the wellbore storage (WBS) is insignificant.Most current studies do not provide a quantitative discussion and do not address factors such as formation geology and mineralogy, resistance-dependent fluid distribution, and geomechanics parameters.
Table 4 summarizes DFIT analysis methods and the details of the fracture/reservoir property estimations from each period in a DFIT BC and AC.Table 5 lists the different leak-off models commonly used to analyze a DFIT, while Table 6 presents the outcomes from each model for several case studies.We have added our conclusions by summarizing the pros and cons of recently published models on abnormal leak-off behaviors (Table 7).
We suggest developing the concept of DFIT to address this lack of a quantitative evaluation, allowing us to evaluate the post-treatment pressure falloff analysis in realtime fracture treatment job optimization.The accurate estimations of closure pressure, ISIP, and perforation and tortuosity friction losses can be obtained to prevent some farfield issues, such as well interventions, frac-hits, high apparent net pressure, and stress shadow.An effective evaluation of a stimulated formation before and during fracture treatment can identify optimal treatment design parameters for individual fracture stages.

Type of Leak-Off Behavior Definition
Normal leak-off behavior • Two scenarios can be considered for normal leak-off behavior: a) All injected fluids leak-off through the fracture network.The fracturing fluid is accumulated and stored in the induced fractures, b) All injected fluids leak-off across the contact surface area of the hydraulic fracture and are then recovered during flowback.

•
The survey in this paper indicates that the assumptions mentioned above may not be appropriate for DFIT analysis in most field cases.
Abnormal leak-off behavior

•
In abnormal leak-off models, part of the injected fluids may leak-off into active secondary fractures.

•
This behavior can be detected on the diagnostic plots from the following signature: a) Higher stress on the diagnostic plot, compared to that of the main HF, may indicate the closure behavior of the secondary fracture network.The main hydraulic fracture propagates to the minimum principal stress.

Type of Leak-Off Behavior Duration Results
Normal leak-off behavior Fracture dimensions, such as fracture width and fracture half-length.

•
Permeability and skin.
Abnormal leak-off behavior BC

•
Friction losses in the wellbore.

•
Net pressure during and after shut-in.

•
Type of leak-off behavior and leak-off coefficient.

•
Extension of secondary fractures.

•
Tip extension distance after shut-in.

•
Fracture surface area for natural fractures.

•
Fracture geometry, such as fracture width and fracture half-length.

AC
• Permeability and skin.

•
The model does not account for fluid loss into natural fractures, such as diagnostic plots modeled based on normal leak-off behavior assumptions.The outcomes, such as the leak-off coefficient, fluid efficiency, and fracture surface area, are associated with high uncertainty. [28]

•
This study investigated the effect of changing fracture compliance on pressure transience for unconventional formations, previously neglected during the closure behavior by previous research; therefore, the model provides accurate estimates of closure pressure analysis.

•
The standard models, such as PKN and KGD, simulate the created fractures; however, the actual fracture geometry is more complex due to the interaction between induced and fissure networks.Newtonian injection fluids and a fracture with uniform leak-off along the fracture face were assumed to form fracture geometry equations. [9]

•
This model involves the impact of secondary fractures on fluids leak-off, where multiple closure events are observed on the G-function and Bourdet pressure derivative plots [45].

•
This approach reduces the need for a rate stepdown test since this model allows us to perform a DFIT analysis to determine friction losses from nearwellbore tortuosity.

•
This model can be used for the main fracture treatment design since it accounts for fluid leak-off into the natural and induced fractures.

•
The fluid efficiency estimated from the normal leak-off behavior is smaller than that of Nolte's model.The model indicates inconsistency with other approaches.

•
This work assumes that leak-off behaviors, such as closure behavior in secondary and created fractures, are governed by one constant leak-off coefficient for all apparent closure events.

•
The surface fracture area calculations are associated with uncertainty due to ignoring the impact of the proppant.We recommend the inclusion of PIFs to estimate the effective contact areas to improve the DFIT assumptions. [20]

•
This model accounts for fracture stiffness/compliance changes as the fracture closes; leak-off rate is a function of fracture pressure.

•
This new concept promotes an understanding of pressure falloff and coupled behaviors during a DFIT with detailed support case studies.

•
This model assumes the closed fracture still retains fracture conductivity.

•
Fluid leak-off remains constant through the fracture surface area before and after closure behavior.

Recommendations When Conducting and Interpreting DFIT
Table 8 provides suggestions and recommendations that interpreters and operators should follow before any field execution or analysis to correctly perform a DFIT and achieve the research goals.

Suggestions and Recommendations Field Execution
Data Analysis [46] • The authors suggest injecting fracturing fluids at a low pumping rate for a short period to create a short fracture.

•
Reduce the time to reach fracture closure and AC pseudo linear and pseudo-radial flow. [11]

•
The pumping rate for a DFIT should be set at a rate close to the planned rate, at least 75%, of the fracture treatment job.

•
The same fracturing fluids should be used for both DFIT and fracture treatment operations.

•
The range of the injection period should be between 3 to 5 minutes for a DFIT.The test should end by performing a rapid step-down test.

•
A successful DFIT design provides reliable closure identification and representative fracture and reservoir properties.

•
A major assumption may not always be true. [27]

•
The final execution relies on several factors, such as the design of the wellhead manifold, pressure gauge types and sampling time, stated objectives, time schedules, wellbore conditioning, pump rate, and pump volume.

•
There are three surface located pressure gauges that should adhere to the specifications of 0.02% full-scale accuracy and 0.01 psi resolution.

•
The test should be performed at the toe stage to fill the well with fracturing fluid.This process will activate pressure and bleeds back by circulating the treatment fluid until trapped gas bubbles are removed.

•
Completion engineers may not overtune or over-smooth the pressure Bourdet derivative when analyzing pressure falloff data since it can cause false flow regimes identification.

•
Distinguishing between actual flow regime behavior and false look-alikes is critical.

•
Reservoir engineers should link the fundamentals of fracture mechanisms with a physical response.

•
Testing should be performed on the well casing, packers, tubing, and wellhead at high pressure.

•
A diesel fluid should be injected on top of the wellbore water to avoid surface line and wellhead freezing.

•
The analysis of DFIT depends on the pressure accuracy; therefore, pressure sources can mitigate erroneous or uninterpretable DFIT pressure responses.

Lessons Learned from DFIT Operations
Table 9 provides lessons learned using a comparison study from DFIT field cases that have been reported in peer-reviewed journal articles.The objective is to explain the impact of operation conditions on DFIT interoperation and analysis of the outcomes.

•
Two similar DFITs were performed in the same formation: a) The first DFIT was performed with a small volume (20.8 bbl) and a low rate (≤5.5 bpm).b) The second DFIT was designed with a large volume (155 bbl) and a high rate (≤14 bpm).

•
Diagnostic plots of both DFITs presented abnormal leak-off behaviors, but the difference in the closure signature was unexpected.a) The diagnostic plots confirmed the transverse storage/height recession behavior for the first DFIT and the pressure-dependent leak-off (PDL) behavior for the second DFIT, b) The final closure events were chosen using a holistic methodology, and the outcomes exhibited consistent values from both DFITs.[48] NIMR/Oman • Two successive informative field DFITs were reported, where they were performed in the same zone in a tight gas formation.a) A small volume (24.8 bbl) at a low rate (0.94-1.38 bpm) was injected in the first DFIT, b) The second DFIT was designed to pump a larger volume (158 bbl) at a higher pumping rate in the range of 12.6 to 18.7 bpm.

•
The results demonstrated a different fracture closure behavior.a) The diagnostic plot of the first DFIT presented wellbore storage (WBS) followed by limited tip extension that defines a simple BC trend close to the behavior of normal leak-off, b) The second DFIT illustrated a complicated BC behavior of variable fracture compliance under different categories that could be pressure-dependent leak-off (PDL), apparent height recession, or transverse storage, c) The BC analysis was not consistent with the closure pressures quantified in both DFITs.[49] Eastern Alberta Shallow/Canada

•
The authors conducted the interpretation of two successive DFITs in a shallow gas-shale formation with a thrust fault setting.a) The first injection had a relatively large volume, at 33.2 bbl, and a high rate at 6.3 bpm compared to the second injection, with an ultra-small volume of 0.82 bbl.

•
The results demonstrated a different fracture closure behavior.a) The diagnostic plots of the first DFIT presented the behavior of transverse storage/height recession closure behavior, while the second DFIT illustrated a normal leak-off behavior, b) The final closure chosen from those two DFITs was consistent but with different closure behavior types.The created fractures did not reopen due to a very small injection volume in the second DFIT.

How DFIT Analysis Is Applicable for Pressure Falloff Data of Main Fracture Treatments
Both DFIT and fracture treatment operations in unconventional horizontal wells have some similarities, especially in the created fractures and behavior of the pressure falloff data immediately after the treatment; however, several dissimilarities in the primary assumptions and operation conditions exist between the two techniques that must be addressed, examined, and basic formulae must be modified to obtain a representative analysis for pressure falloff behavior on a stage-by-stage basis in real-time.Very few studies [31] have analyzed pressure falloff data applications on a stage-by-stage basis in a horizontal wellbore.We have addressed the main differences between the two tests: 1.A DFIT is performed by pumping a small volume (24.8 bbl) at a low rate (0.94-1.38 bpm), while the injection volume of the main fracture treatment in each stage is higher; therefore, the formation pore pressure may not be accounted for.We account for the formation pore pressure during the DFIT analysis.2. A DFIT only operates with water that does not include proppant, but the role of proppant exists in the main fracture treatment jobs.We recommend the inclusion of PIFs to improve the DFIT assumptions.We can apply the DFIT methodology to the main fracture treatment pressure falloff data to estimate an effective contact fracture surface area.3. A DFIT operation creates only one fracture, while hydraulic fracture treatments generate multiple perforation clusters in each stage.The fracture geometry is different from stage to stage on a well due to several factors, such as stress shadow effects, formation heterogeneity, formation lithology, and resistance-dependent fluid distribution.4. The falloff period immediately after DFIT is relatively long, up to a week compared to a short period of half an hour in the hydraulic fracture treatment, which is enough to perform the analysis.

Effective Fracture Surface Area Calculations
This section presents the workflow to analyze post-treatment pressure falloff data in order to estimate effective fracture surface areas for both natural fractures and hydraulic fractures on a stage-by-stage basis.Figure 6 illustrates our methodology, and Table 10 lists the main Equations ( 14) through (19) that are used to determine the primary outcomes.The same concept was used to analyze the pressure falloff data by Liu et al. [9], who determined the total fracture surface area for both natural fractures and hydraulic fractures on a stage-by-stage basis.In this paper, we modified the workflow, which enables us to calculate fracture half-length and identify the main fracture flow regimes after the treatment of the Meramec Formation, STACK Play of the Anadarko Basin, Oklahoma.

Case Study Description
The DFIT was performed at the toe stage of an unconventional horizontal well (W2) in the Meramec Formation, STACK Play of the Anadarko Basin, Oklahoma.The STACK Play is a multi-layered tight oil reservoir "Meramec and Woodford Formations" with porosity lying in the range of 3% to 10%, and permeability in the range of 0.0001 to 0.01 mD.The plug and perf stimulation technique was used to complete the well.The wellbore was prepared before testing by filling it with fracturing fluids, typically water, to pressurize the top of the wellbore up to the point of completion so the formation did not break down if there were no pre-existing fractures.The fracturing fluid was pumped into the formation through the first casing interval with an inner diameter (ID) of 4.67 inches, a total vertical depth (TVD) of 9648 ft, and a measured depth (MD) of 19,635 ft.The pressure rises linearly with the injection volume during the injection period, while the injection rate remains constant (Figure 7). Figure 7 presents the injection and pressure profiles for our DFIT case study, where the surface injection rate was maintained at 12 bpm for approximately 40 minutes.The first shut-in period was measured for the same injection period, then the injection test was performed again with a higher injection rate of approximately 13 bpm.The test was then completed, and the pressure falloff data were recorded during the next 4-5 days.
The main fracture treatment job was performed in the same well (W2) based on the evaluation report from the DFIT analysis.Figure 8 presents the main hydraulic fracturing operation, where the job consists of several pump schedule stages, as shown in the figure by green and blue colors, and the slurry rates and proppant concentrations were changed during the operation.The total job period is around 160 min at an average 100 bpm slurry rate.At the end of the job, the falloff pressure was recorded for a period of time of approximately 15 min on average for all fracture stages.The fracture treatment strategy was applied with a constant pump schedule to create 36 frac stages, where each stage consisted of four to five clusters with 50 ft spacing.This section presents a case study of the unconventional horizontal well, where Table 11 shows the main fracture treatment operation parameters as well as preprocessing data quality for each fracture stage.The combination of DFIT and pressure falloff data immediately after the treatment was used to identify fracture and reservoir behavior characteristics to assess the evaluation of the fracture stages.[45] (Figure 9), and (2) the plot of the G-function diagnostic analysis (Figure 10).Both plots indicate consistent results and the same height recession leak-off behavior.The mechanisms of this closure fracture behavior are defined as most of the fracturing fluid's leak-offs into neighboring layers, a common behavior for unconventional formations.The plots also indicate that natural fractures have contributed to the multiple fracture closure events observed before the closure analysis.
Different trends, such as circle 1, appear to be wellbore storage coupled with friction dissipation.The following flat trend, circle 2, appears to be tip extension with limited growth distance, as shown in Figure 11.Linear flow with ½ signature, circle 3, is indicated on the plot before the closure behavior, marked by the green dashed line.Two closure events overlay on a 3/2 slope, circle 4, that depicts the closure behavior of the natural and induced fractures.The reservoir dominated flow observed after the closure analysis with two flow regimes started by linear flow, a ½ slope, circle 5, was followed by radial flow, and a zero slope, circle 6, where formation permeability can be estimated.The main results from DFIT analysis are listed in Table 12, where the formation is classified with low permeability in the range of 0.004 md and closure pressure in 6100 psi.The same concept was used to analyze the pressure falloff data by Liu et al. [9], who determined the total fracture surface area for both natural fractures and hydraulic fractures on a stage-by-stage basis.However, we modified the workflow to improve the data quality, minimizing the effects of water hammers for the first 2 min in falloff pressure data and including the analysis of the G-function plot in order to validate fracture flow regimes, as well as considering variable parameters for each stage.This methodology has not been applied before in the literature, and it is more precise with low uncertainty compared to previous research, where constant values were used for all stages to estimate total fracture surface areas.Our case study shows three different behaviors after the main fracture treatment as follows: a. Case #1, ¼ slope, tip extension of the main fracture: Figure 12 presents the behavior of tip extension of the main hydraulic fracture, where the results showed that the fracture surface area for the main fracture is higher than the total fracture surface area for the natural fracture.The reason is that low total pressure losses and more fracturing fluids leak-off through main fractures instead of nearby neighbor layers.b.Case #2, -1/2 and -1 slope, fracture height recession: Figure 13  For more detailed results, Table 13 provides the analysis outcomes, including friction pressure losses and total fracture surface area.Very few studies [9] have investigated the pressure falloff behavior combined with DFIT analysis to evaluate the effectiveness of fracture treatment conditions stage-by-stage; however, the assessment of fracture design quality stage-by-stage must still be studied, especially after treatment.The principle of DFIT needs to be adjusted by counting PIFs to estimate the effective fracture contact area, or the propped fracture area per unit length of lateral, which is one key factor for evaluating post-treatment performance.An integral approach using several diagnostics tools is necessary to develop the technology for determining the contribution of individual fracturing treatments when multiple factors are considered.This approach can be accomplished using a series of neural network models to predict fracture geometry, fracture directions, number of clusters needed, proppant loading, fracture complexity, and treatment costs during fracture treatment.A data-driven model must be created to apply an integral approach in real-time, a challenge that must be addressed with thorough discussion, especially given the lack of literature on this subject.These combination methods involve machine learning tools that can assist us in understanding fracture treatment effectiveness, assess new completion technologies, and evaluate which formations are the most productive.This technology is feasible for realtime analysis to apply an optimum pump schedule for the current and next treatment stages on a well.

Summary
In this study, we employed an integral approach to identify fracture and reservoir behavior characteristics to assess the performance of hydraulic fractures stage-by-stage.The research findings highlight the following:

•
Diagnostic tools have become attractive to the oil and gas industry since they are powerful methods for identifying fracture and reservoir characteristics.The analysis of these approaches is based on data type, which guides the assessment of complex fracture networks generated by hydraulic fracture treatment operations.

•
The use of falloff pressure data is a simple approach since it does not require more information than shut-in pressure vs. time.The data are recorded immediately after shut-in at no additional cost; however, micro-seismic monitoring is expensive, and the production data analysis has high uncertainty.

•
The models proposed by Liu and Ehlig-Economides [28] allowed us to perform DFIT analysis to determine flow friction losses in the wellbore, perforation, and near-wellbore tortuosity separately without the need for rate step-down tests.These models are the primary means by which we can calculate the total fracture surface area for the secondary fractures and hydraulic fractures and gain insight into the effectiveness of the hydraulic fracturing process.

•
Barree et al. [11] and Liu and Ehlig-Economides [28] reported that their data supported a holistic methodology that allows them to pick the hydraulic fracture closure, which contradicts the variable fracture compliance approach from McClure et al. [50].

•
Rizwan [48] created an alternative methodology to apply a DFIT, where a change in the test operation strategy may cause a change in fracture closure behavior.This study provided insight into expected fracture closure behavior during the main fracture treatment.

•
We concluded that evaluating the performance of post-stimulation conditions on a stage-by-stage basis using indirect methods, such as pressure falloff data analysis, is the most promising technique for providing a wide range of information covering the mechanics of the hydraulic fracture, such as open, closed, and propped.This method may overcome limitations and weaknesses found in many of the proposed techniques reported in peer-reviewed journal articles, such as production data analysis and micro-seismic methods.

•
This technology is a critical factor in the economic development of unconventional reservoirs since the well completion cost is a significant portion of the capital cost compared to other expenses, and heavily influences production rate or ultimate recovery.

•
We suggest combining static and dynamic diagnostic methods to better estimate fracture geometry through pressure data, diagnostic fracture injection tests (DFIT), micro-seismic fracture mapping, distributed temperature sensing (DTS), production logs, and production data.The full suite of information can provide valuable evidence concerning the details of the treatment and well performance from complicated shale plays.• Three different cases were observed through diagnostic plots, where mainly the analysis indicates that most of the fracturing fluid was leaked off through the natural fracture surface area and resulted in the estimation of larger values compared to the hydraulic fracture calculated area.These phenomena might represent a secondary fracture set with a high fracture closure stress activated in neighbor stages that was not well-developed in other sections.

•
This conclusion fits with our discussion that provided detailed information with support case studies to apply the technology of post-treatment pressure analysis in realtime.

Figure 1 .
Figure 1.Average oil production per well in the Bakken Formation, Williston Basin, North Dakota [10].

Figure 2 .
Figure 2. Typical pressure behavior includes a sequence of main events observed in a DFIT.ISIP is instantaneous shut-in pressure.

Figure 3 .
Figure 3.Typical flow regimes before and after closure behaviors.

Figure 6 .
Figure 6.Workflow of research methodology to determine fracture stage performance.

Table 1 .
Comparison of different diagnostic methods for fracture treatment performance analysis.

a. Production Data Analysis as a Diagnostic Method: Pros Cons Results Authors •
• Micro-seismic imaging cannot provide accurate information on individual fractures and cracks, or whether they are open, closed, or propped/unpropped.
• It is not clear why micro-seismics may not detect the tensile failure at the advancing tip of a hydraulic fracture.

Table 3 .
Primary outcomes from a DFIT based on stipulated fracture geometry.

Table 4 .
Main results and methodologies of DFIT analysis.

Table 5 .
Definition of two main leak-off behaviors.

Table 7 .
Recently published DFIT models under the conditions of abnormal leak-off behaviors.

Table 8 .
Important points that extended in planning, executing, and achieving the successful DFIT.

Table 9 .
Field case studies of DFIT operations and lessons learned.

Table 11 .
Main hydraulic fracturing parameters and data quality check for post-treatment pressure analyses (W2).

Table 13 .
Outcomes of pressure falloff analyses for the main hydraulic fracture treatment of W2.