Using fracture-based continuum modeling of coupled geomechanical-hydrological processes for numerical simulation of hydraulic fracturing

This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.

Fourth-order phase-field modeling for brittle fracture in piezoelectric materials

Failure analyses of piezoelectric structures and devices are of engineering and scientific significance.In this paper,a fourth-order phase-field fracture model for piezoelectric solids is developed based on the Hamilton principle.Three typical electric boundary conditions are involved in the present model to characterize the fracture behaviors in various physical situations.A staggered algorithm is used to simulate the crack propagation.The polynomial splines over hierarchical T-meshes(PHT-splines)are adopted as the basis function,which owns the C1continuity.Systematic numerical simulations are performed to study the influence of the electric boundary conditions and the applied electric field on the fracture behaviors of piezoelectric materials.The electric boundary conditions may influence crack paths and fracture loads significantly.The present research may be helpful for the reliability evaluation of the piezoelectric structure in the future applications.

Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.

Development of an experimental method for well-controlled blast induced traumatic limb fracture in rats

Heterotopic ossification(HO)is a consequence of traumatic bone and tissue damage,which occurs in 65%of military casualties with blast-associated amputations.However,the mechanisms behind blast-induced HO remain unclear.Animal models are used to study blast-induced HO,but developing such models is challenging,particularly in how to use a pure blast wave(primary blast)to induce limb fracture that then requires an amputation.Several studies,including our recent study,have developed platforms to induce limb fractures in rats with blast loading or a mixture of blast and impact loading.However,these models are limited by the survivability of the animal and repeatability of the model.In this study,we developed an improved platform,aiming to improve the animal's survivability and injury repeatability as well as focusing on primary blast only.The platform exposed only one limb of the rat to a blast wave while providing proper protection to the rest of the rat's body.We obtained very consistent fracture outcome in the tibia(location and pattern)in cadaveric rats with a large range of size and weight.Importantly,the rats did not obviously move during the test,where movement is a potential cause of uncontrolled injury.We further conducted parametric studies by varying the features of the design of the platform.These factors,such as how the limb is fixed and how the cavity through which the limb is placed is sealed,significantly affect the resulting injury.This platform and test setups enable well-controlled limb fracture induced directly by pure blast wave,which is the fundamental step towards a complete in vivo animal model for blast-induced HO induced by primary blast alone,excluding secondary and tertiary blast injury.In addition,the platform design and the findings presented here,particularly regarding the proper protection of the animal,have implications for future studies investigating localized blast injuries,such as blast induced brain and lung injuries.

Three predictive scores compared in a retrospective multicenter study of nonunion tibial shaft fracture

BACKGROUND Delayed union,malunion,and nonunion are serious complications in the healing of fractures.Predicting the risk of nonunion before or after surgery is challenging.AIM To compare the most prevalent predictive scores of nonunion used in clinical practice to determine the most accurate score for predicting nonunion.METHODS We collected data from patients with tibial shaft fractures undergoing surgery from January 2016 to December 2020 in three different trauma hospitals.In this retrospective multicenter study,we considered only fractures treated with intramedullary nailing.We calculated the tibia FRACTure prediction healING days(FRACTING)score,Nonunion Risk Determination score,and Leeds-Genoa Nonunion Index(LEG-NUI)score at the time of definitive fixation.RESULTS Of the 130 patients enrolled,89(68.4%)healed within 9 months and were classified as union.The remaining patients(n=41,31.5%)healed after more than 9 months or underwent other surgical procedures and were classified as nonunion.After calculation of the three scores,LEG-NUI and FRACTING were the most accurate at predicting healing.CONCLUSION LEG-NUI and FRACTING showed the best performances by accurately predicting union and nonunion.

Estimation of fracture size and azimuth in the universal elliptical disc model based on trace information

The geometric characteristics of fractures within a rock mass can be inferred by the data sampling from boreholes or exposed surfaces.Recently,the universal elliptical disc(UED)model was developed to represent natural fractures,where the fracture is assumed to be an elliptical disc and the fracture orientation,rotation angle,length of the long axis and ratio of short-long axis lengths are considered as variables.This paper aims to estimate the fracture size-and azimuth-related parameters in the UED model based on the trace information from sampling windows.The stereological relationship between the trace length,size-and azimuth-related parameters of the UED model was established,and the formulae of the mean value and standard deviation of trace length were proposed.The proposed formulae were validated via the Monte Carlo simulations with less than 5%of error rate between the calculated and true values.With respect to the estimation of the size-and azimuth-related parameters using the trace length,an optimization method was developed based on the pre-assumed size and azimuth distribution forms.A hypothetical case study was designed to illustrate and verify the parameter estimation method,where three combinations of the sampling windows were used to estimate the parameters,and the results showed that the estimated values could agree well with the true values.Furthermore,a hypothetical three-dimensional(3D)elliptical fracture network was constructed,and the circular disc,non-UED and UED models were used to represent it.The simulated trace information from different models was compared,and the results clearly illustrated the superiority of the proposed UED model over the existing circular disc and non-UED models。

Geological characteristics and models of fault-foldfracture body in deep tight sandstone of the second member of Upper Triassic Xujiahe Formation in Xinchang structural belt of Sichuan Basin,SW China

In the second member of the Upper Triassic Xujiahe Formation(T_(3)x_(2))in the Xinchang area,western Sichuan Basin,only a low percent of reserves has been recovered,and the geological model of gas reservoir sweet spot remains unclear.Based on a large number of core,field outcrop,test and logging-seismic data,the T_(3)x_(2) gas reservoir in the Xinchang area is examined.The concept of fault-fold-fracture body(FFFB)is proposed,and its types are recognized.The main factors controlling fracture development are identified,and the geological models of FFFB are established.FFFB refers to faults,folds and associated fractures reservoirs.According to the characteristics and genesis,FFFBs can be divided into three types:fault-fracture body,fold-fracture body,and fault-fold body.In the hanging wall of the fault,the closer to the fault,the more developed the effective fractures;the greater the fold amplitude and the closer to the fold hinge plane,the more developed the effective fractures.Two types of geological models of FFFB are established:fault-fold fracture,and matrix storage and permeability.The former can be divided into two subtypes:network fracture,and single structural fracture,and the later can be divided into three subtypes:bedding fracture,low permeability pore,and extremely low permeability pore.The process for evaluating favorable FFFB zones was formed to define favorable development targets and support the well deployment for purpose of high production.The study results provide a reference for the exploration and development of deep tight sandstone oil and gas reservoirs in China.

A fracture model for assessing tensile mode crack growth resistance of rocks

Evaluating the fracture resistance of rocks is essential for predicting and preventing catastrophic failure of cracked structures in rock engineering.This investigation developed a brittle fracture model to predict tensile mode(mode I)failure loads of cracked rocks.The basic principle of the model is to estimate the reference crack corresponding to the fracture process zone(FPZ)based on the maximum normal strain(MNSN)ahead of the crack tip,and then use the effective crack to calculate the fracture toughness.We emphasize that the non-singular stress/strain terms should be considered in the description of the MNSN.In this way,the FPZ,non-singular terms and the biaxial stress state at the crack tip are simul-taneously considered.The principle of the model is explicit and easy to apply.To verify the proposed model,laboratory experiments were performed on a rock material using six groups of specimens.The model predicted the specimen geometry dependence of the measured fracture toughness well.More-over,the potential of the model in analyzing the size effect of apparent fracture toughness was discussed and validated through experimental data reported in the literature.The model was demonstrated su-perior to some commonly used fracture models and is an excellent tool for the safety assessment of cracked rock structures.

Analytical model of 2D leakoff in waterflood-induced fractures

Waterflood-induced fractures,also known as self-induced fractures,spontaneously form at injection wells during waterflooding.These fractures propagate long distances through rock,allowing injected fluids to travel far away from a well,both within and outside the flooding layer.Essentially,the me-chanics of waterflood-induced fracture propagation is similar to that of hydraulic fractures,which are intentionally created for reservoir stimulation.Fracturing models developed for hydraulic fractures can also be applied to waterflood-induced fractures.However,waterflood-induced fractures are typically pumped with much larger volumes of water or brine and grow much longer in time.As a result,fluid leakoff from waterflood fractures into the formation is more extensive and two-dimensional(2D),a characteristic that is often ignored in a majority of modern fracturing simulators,making their appli-cation to waterflood fractures unreliable.In this work,we revisit the problem of leakoff for long-growing waterflood-induced fractures and develop a new analytical model for fluid leakoff that provides improved predictions of fracture geometry and can be easily implemented in fracturing simulators.We incorporate the developed solution into the classical Perkins-Kern-Nordgren(PKN)model of fracture growth,which shows that the choice of the Carter or a 2D leakoff model greatly impacts fracture ge-ometry at large time.The conducted parametric study shows while a toughness-dominated regime af-fects fracture evolution,most of fracture lifetime occurs in a viscosity-and-leakoff-dominated regime.We also develop an asymptotic solution for a leakoff profile in the limiting case of 2D leakoff domination(~~M and~~K).Finally,we study 3D fracture growth and out-of-zone injection with three layers and a complex structure of zones.The study shows that ignoring the 2D leakoff during simulation results in a significant overestimation of fracture geometry predictions.The present work,thus,plays an important role in improving waterflood fracture modelling,as it highlights the significance of 2D leakoff in waterflood-induced fractures and provides a reliable analytical model for fluid leakoff that can be incorporated into modern fracture simulators.

Fracture propagation laws of staged hydraulic fracture in fractured geothermal reservoir based on phase field model

Hydraulic fracturing is widely used in geothermal resource exploitation, and many natural fractures exist in hot dry rock reservoirs due to in-situ stress and faults. However, the infuence of natural fractures on hydraulic fracture propagation is not considered in the current study. In this paper, based on the phase feld model, a thermo-hydro-mechanical coupled hydraulic fracture propagation model was established to reveal the infuence of injection time, fracturing method, injection fow rate, and natural fracture distribution on the fracture propagation mechanism. The results show that fracture complexity increases with an increase in injection time. The stress disturbance causes the fracture initiation pressure of the second cluster signifcantly higher than that of the frst and third clusters. The zipper-type fracturing method can reduce the degree of stress disturbance and increase fracture complexity by 7.2% compared to simultaneous hydraulic fracturing. Both low and high injection fow rate lead to a decrease in fracture propagation time, which is not conducive to an increase in fracture complexity. An increase in the natural fracture angle leads to hydraulic fracture crossing natural fracture, but has a lesser efect on fracture complexity. In this paper, we analyzed the infuence of diferent factors on initiation pressure and fracture complexity, providing valuable guidance for the exploitation of geothermal resources.

Development of an improved three-dimensional rough discrete fracture network model:Method and application

Structure plane is one of the important factors affecting the stability and failure mode of rock mass engineering.Rock mass structure characterization is the basic work of rock mechanics research and the important content of numerical simulation.A new 3-dimensional rough discrete fracture network(RDFN3D)model and its modeling method based on the Weierstrass-Mandelbrot(W-M)function were presented in this paper.The RDFN3D model,which improves and unifies the modelling methods for the complex structural planes,has been realized.The influence of fractal dimension,amplitude,and surface precision on the modeling parameters of RDFN3D was discussed.The reasonable W-M parameters suitable for the roughness coefficient of JRC were proposed,and the relationship between the mathematical model and the joint characterization was established.The RDFN3D together with the smooth 3-dimensional discrete fracture network(DFN3D)models were successfully exported to the drawing exchange format,which will provide a wide application in numerous numerical simulation codes including both the continuous and discontinuous methods.The numerical models were discussed using the COMSOL Multiphysics code and the 3-dimensional particle flow code,respectively.The reliability of the RDFN3D model was preliminarily discussed and analyzed.The roughness and spatial connectivity of the fracture networks have a dominant effect on the fluid flow patterns.The research results can provide a new geological model and analysis model for numerical simulation and engineering analysis of jointed rock mass.

Revolutionizing Tight Reservoir Production: A Novel Dual-Medium Unsteady Seepage Model for Optimizing Volumetrically Fractured Horizontal Wells

This study presents an avant-garde approach for predicting and optimizing production in tight reservoirs,employing a dual-medium unsteady seepage model specifically fashioned for volumetrically fractured horizontal wells.Traditional models often fail to fully capture the complex dynamics associated with these unconventional reservoirs.In a significant departure from these models,our approach incorporates an initiation pressure gradient and a discrete fracture seepage network,providing a more realistic representation of the seepage process.The model also integrates an enhanced fluid-solid interaction,which allows for a more comprehensive understanding of the fluid-structure interactions in the reservoir.This is achieved through the incorporation of improved permeability and stress coupling,leading to more precise predictions of reservoir behavior.The numerical solutions derived from the model are obtained through the sophisticated finite element method,ensuring high accuracy and computational efficiency.To ensure the model’s reliability and accuracy,the outcomes were tested against a real-world case,with results demonstrating strong alignment.A key revelation from the study is the significant difference between uncoupled and fully coupled volumetrically fractured horizontal wells,challenging conventional wisdom in the field.Additionally,the study delves into the effects of stress,fracture length,and fracture number on reservoir production,contributing valuable insights for the design and optimization of tight reservoirs.The findings from this study have the potential to revolutionize the field of tight reservoir prediction and management,offering significant advancements in petroleum engineering.The proposed approach brings forth a more nuanced understanding of tight reservoir systems and opens up new avenues for optimizing reservoir management and production.

Calculation Model of Equivalent Strength for Induced Crack Based on Double-K Fracture Theory and Its Optimizing Setting in RCC Arch Dam

By means of fracture testing on roller-compacted concrete （RCC） three-point bending beams with two different specimen sizes, the P-CMOD complete curve for RCC was gained. Furthermore, by applying double-K fracture theory, KiniⅠC,KunⅠC, as well as the critical effective crack length and the critical crack tip opening displacement, were evaluated. Based on the double-K fracture parameters above, the calculation model of equivalent strength for induced crack was established, thus the calculation method on its initiation, stable propagation and unstable fracture was ascertained. Moreover, the finite element simulation analysis of stress field in ShaPai arch dam and the on-site observational splaying points of induced crack at different altitudes validated the reliability of the model. Finally, crack inducer′s optimal setting in RCC arch dam was studied. It improves the design level of induced crack in RCC arch dam and satisfies the necessity of engineering practice.

Probing the influence of secondary fracture connectivity on fracturing fluid flowback efficiency

A deep understanding of the geometric impacts of fracture on fracturing fluid flowback efficiency is essential for unconventional oil development. Using nuclear magnetic resonance and 2.5-dimensional matrix-fracture visualization microfluidic models, qualitative and quantitative descriptions of the influences of connectivity between primary fracture and secondary fracture on flowback were given from core scale to pore network scale. The flow patterns of oil-gel breaking fluid two-phase flow during flowback under different fracture connectivity were analyzed. We found some counterintuitive results that non-connected secondary fracture (NCSF, not connect with artificial primary fracture and embedded in the matrix) is detrimental to flowbackefficiency. The NCSF accelerates the formation of oil channeling during flowback, resulting in a large amount of fracturing fluid trapped in the matrix, which is not beneficial for flowback. Whereas the connected secondary fracture (CSF, connected with the artificial primary fracture) is conducive to flowback. The walls of CSF become part of primary fracture, which expands the drainage area with low resistance, and delays the formation of the oil flow channel. Thus, CSF increases the high-speed flowback stage duration, thereby enhancing the flowback efficiency. The fracturing fluid flowback efficiency investigated here follows the sequence of the connected secondary fracture model (72%) > the matrix model (66%) > the non-connected secondary fracture model (38%). Our results contribute to hydraulic fracturing design and the prediction of flowback efficiency.

Improved Staggered Algorithm for Phase-Field Brittle Fracture with the Local Arc-Length Method

The local arc-length method is employed to control the incremental loading procedure for phase-field brittle fracture modeling.An improved staggered algorithm with energy and damage iterative tolerance convergence criteria is developed based on the residuals of displacement and phase-field.The improved staggered solution scheme is implemented in the commercial software ABAQUS with user-defined element subroutines.The layered system of finite elements is utilized to solve the coupled elastic displacement and phase-field fracture problem.A one-element benchmark test compared with the analytical solution was conducted to validate the feasibility and accuracy of the developed method.Our study shows that the result calculated with the developed method does not depend on the selected size of loading increments.The results of several numerical experiments show that the improved staggered algorithm is efficient for solving the more complex brittle fracture problems.

Simulation study of supercritical carbon dioxide jet fracturing for carbonate geothermal reservoir based on fluid-thermo-mechanical coupling model

Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon dioxide(SC-CO_(2))jet fracturing is expected to efficiently stimulate the carbonate geothermal reservoirs and achieve the storage of CO_(2) simultaneously.In this paper,we established a transient seepage and fluid-thermo-mechanical coupled model to analyze the impact performance of sc-CO_(2) jet fracturing.The mesh-based parallel code coupling interface was employed to couple the fluid and solid domains by exchanging the data through the mesh interface.The physical properties change of sC-CO_(2) with temperature were considered in the numerical model.Results showed that SC-CO_(2) jet frac-turing is superior to water-jet fracturing with respect to jetting velocity,particle trajectory and pene-trability.Besides,stress distribution on the carbonate rock showed that the tensile and shear failure would more easily occur by SC-CO_(2) jet than that by water jet.Moreover,pressure and temperature control the jet field and seepage field of sC-CO_(2) simultaneously.Increasing the jet temperature can effectively enhance the impingement effect and seepage process by decreasing the viscosity and density of SC-CO_(2).The key findings are expected to provide a theoretical basis and design reference for applying SC-CO_(2) jet fracturing in carbonate geothermal reservoirs.

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.

A 3D analysis of the occurrence of fractures in hot dry rock reservoirs based on the spatial distribution of natural fractures

Hot dry rock(HDR)geothermal energy is a kind of widely distributed clean energy with huge reserves.However,its commercial development has been constrained by reservoir stimulation.In the early stage of HDR geothermal energy development,properly determining spatial distribution patterns of natural fractures in HDR reservoirs can effectively guide reservoir stimulation.This study analyzes the spatial distribution of natural fractures by using FracMan software based on the actual geological data and log data of well M-2 in the Matouying Uplift area,Hebei Province.The fracture parameters are counted and Monte Carlo simulation technique is introduced to optimize the parameters,which makes the natural fracture model more accurate and reliable.Furthermore,this study simulates hydraulic fracturing using the model combined with the actual in-situ stress parameters and the construction scheme.As verified by fitting the changes in simulated wellhead pressure during hydraulic fracturing with the actual wellhead pressure data detected during construction,the methods for natural fracture modeling used in this study are scientific and reasonable.The preliminary prediction results show that the displacement design scheme with a pump displacement of 2.0-3.0 m^(3)/min,4.0-5.5 m^(3)/min and 6-7 m^(3)/min in the early,middle and late stages,respectively,has good fracturing effect.The results of this study can be utilized as a reference for preparing development schemes for HDR reservoirs.

Numerical investigation of dual-porosity model with transient transfer function based on discrete-fracture model

Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.

The tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy: A case study in the Wangyao Oilfield of Ordos Basin, China

Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantitative information on a large number of fractures, based on which the accuracy of subsurface fracture characterization can be improved. Here we take the tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy, a case study of the Chang6t~ Formation of the Upper Triassic Yanchang Group of the Wangyao Oilfield in the Ordos Basin of China. An outcrop at the edge of the basin is a suitable analog for the reservoir, but the prerequisite is that they must have equivalent previous stress fields, similar final structural characteristics, relative timing and an identical depositional environment and diagenesis. The relationship among fracture density, rock type and bed thickness based on the outcrop is one of the most important fracture distribution models, and can be used to interpret fracture density in individual wells quantitatively. Fracture orientation, dip, geometry and scale, also should be described and measured in the outcrop, and can be used together with structure restoration and single well fracture density interpretation to guide fracture intensity prediction on bed surfaces and to constrain the construction of the 3D fracture geometry model of the subsurface reservoir. The application of the above principles shows the outcrop-based tectonic fracture models of the target ultra-low permeability sandstone reservoir are consistent with fractures inferred from microseismic interpretation and tracer tests. This illustrated that the fracture modeling based on the outcrop analogy is reliable and can reduce the uncertainty in stochastic fracture modeling.