Theory,technology and practice of shale gas three-dimensional development:A case study of Fuling shale gas field in Sichuan Basin,SW China

In the Jiaoshiba block of the Fuling shale gas field,the employed reserves and recovery factor by primary well pattern are low,no obvious barrier is found in the development layer series,and layered development is difficult.Based on the understanding of the main factors controlling shale gas enrichment and high production,the theory and technology of shale gas three-dimensional development,such as fine description and modeling of shale gas reservoir,optimization of three-dimensional development strategy,highly efficient drilling with dense well pattern,precision fracturing and real-time control,are discussed.Three-dimensional development refers to the application of optimal and fast drilling and volume fracturing technologies,depending upon the sedimentary characteristics,reservoir characteristics and sweet spot distribution of shale gas,to form"artificial gas reservoir"in a multidimensional space,so as to maximize the employed reserves,recovery factor and yield rate of shale gas development.In the research on shale gas three-dimensional development,the geological+engineering sweet spot description is fundamental,the collaborative optimization of natural fractures and artificial fractures is critical,and the improvement of speed and efficiency in drilling and fracturing engineering is the guarantee.Through the implementation of three-dimensional development,the overall recovery factor in the Jiaoshiba block has increased from 12.6%to 23.3%,providing an important support for the continuous and stable production of the Fuling shale gas field.

Integration system research and development for three-dimensional laser scanning information visualization in goaf

An integration processing system of three-dimensional laser scanning information visualization in goaf was developed. It is provided with multiple functions, such as laser scanning information management for goaf, cloud data de-noising optimization, construction, display and operation of three-dimensional model, model editing, profile generation, calculation of goaf volume and roof area, Boolean calculation among models and interaction with the third party soft ware. Concerning this system with a concise interface, plentiful data input/output interfaces, it is featured with high integration, simple and convenient operations of applications. According to practice, in addition to being well-adapted, this system is favorably reliable and stable.

Three-dimensional stability calculation method for high and large composite slopes formed by mining stope and inner dump in adjacent open pits

The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of significant importance to develop an effective 3D stability calculation method for composite slopes to enhance the design and stability control of open-pit slope engineering.Using the composite slope formed by the mining stope and inner dump in Baiyinhua No.1 and No.2 open-pit coal mine as a case study,this research investigates the failure mode of composite slopes and establishes spatial shape equations for the sliding mass.By integrating the shear resistance and sliding force of each row of microstrip columns onto the bottom surface of the strip corresponding to the main sliding surface,a novel 2D equivalent physical and mechanical parameters analysis method for the strips on the main sliding surface of 3D sliding masses is proposed.Subsequently,a comprehensive 3D stability calculation method for composite slopes is developed,and the quantitative relationship between the coordinated development distance and its 3D stability coefficients is examined.The analysis reveals that the failure mode of the composite slope is characterized by cutting-bedding sliding,with the arc serving as the side interface and the weak layer as the bottom interface,while the destabilization mechanism primarily involves shear failure.The spatial form equation of the sliding mass comprises an ellipsoid and weak plane equation.The analysis revealed that when the coordinated development distance is 1500 m,the error rate between the 3D stability calculation result and the 2D stability calculation result of the composite slope is less than 8%,thereby verifying the proposed analytical method of equivalent physical and mechanical parameters and the 3D stability calculation method for composite slopes.Furthermore,the3D stability coefficient of the composite slope exhibits an exponential correlation with the coordinated development distance,with the coefficient gradually decreasing as the coordinated development distance increases.These findings provide a theoretical guideline for designing similar slope shape parameters and conducting stability analysis.

Fracture of two three-dimensional parallel internal cracks in brittle solid under ultrasonic fracturing

Similar to hydraulic fracturing(HF), the coalescence and fracture of cracks are induced within a rock under the action of an ultrasonic field, known as ultrasonic fracturing(UF). Investigating UF is important in both hard rock drilling and oil and gas recovery. A three-dimensional internal laser-engraved crack(3D-ILC) method was introduced to prefabricate two parallel internal cracks within the samples without any damage to the surface. The samples were subjected to UF. The mechanism of UF was elucidated by analyzing the characteristics of fracture surfaces. The crack propagation path under different ultrasonic parameters was obtained by numerical simulation based on the Paris fatigue model and compared to the experimental results of UF. The results show that the 3D-ILC method is a powerful tool for UF research.Under the action of an ultrasonic field, the fracture surface shows the characteristics of beach marks and contains powder locally, indicating that the UF mechanism includes high-cycle fatigue fracture, shear and friction, and temperature load. The two internal cracks become close under UF. The numerical result obtained by the Paris fatigue model also shows the attraction of the two cracks, consistent with the test results. The 3D-ILC method provides a new tool for the experimental study of UF. Compared to the conventional numerical methods based on the analysis of stress-strain and plastic zone, numerical simulation can be a good alternative method to obtain the crack path under UF.

Discrete Element Modeling of Asphalt Concrete Cracking Using a User-def ined Three-dimensional Micromechanical Approach

We established a user-defined micromechanical model using discrete element method （DEM） to investigate the cracking behavior of asphalt concrete （AC）. Using the ＂Fish＂ language provided in the particle flow code in 3-Demensions （PFC3D）, the air voids and mastics in asphalt concrete were realistically built as two distinct phases. With the irregular shape of individual aggregate particles modeled using a clump of spheres of different sizes, the three-dimensional （3D） discrete element model was able to account for aggregate gradation and fraction. Laboratory uniaxial complex modulus test and indirect tensile strength test were performed to obtain input material parameters for the numerical simulation. A set of the indirect tensile test were simulated to study the cracking behavior of AC at two levels of temperature, i e, -10 ℃ and 15 ℃. The predicted results of the numerical simulation were compared with laboratory experimental measurements. Results show that the 3D DEM model is able to predict accurately the fracture pattern of different asphalt mixtures. Based on the DEM model, the effects of air void content and aggregate volumetric fraction on the cracking behavior of asphalt concrete were evaluated.

Simulation of Morphological Development of Soil Cracks in Yuanmou Dry-hot Valley Region, Southwest China

Soil cracking is an important process influencing water and solutes transport in the Yuanmou Dry-hot Valley region of Southwest China. Studying the morphological development of soil cracks helps to further reveal the close relationship between the soil cracking process and water movement in such semi-arid regions. Here we report regular changes on surface morphology of soil cracks with decreasing water in four different soils (Typ-Ustic Ferrisols,Ver-Ustic Ferrisols,Tru-Ustic Vertisols and Typ-Ustic Vertisols) through simulation experiments. Our results indicate the following: 1) Different soils ultimately have different development degrees of soil cracks,according to their various values of crack area density. Soil cracks in Typ-Ustic Ferrisols can only develop to the feeble degree,while those in the other three soils are capable of developing into the intensive degree,and even into the extremely intensive degree. 2) Soil crack complexity,as expressed by the value of the area-weighted mean of crack fractal dimension (AWMFRAC),is found to continuously decrease as a whole through the whole cracking process in all the studied soils. 3) Soil crack connectivity shows a uniform trend in the studied soils,that is to say,connectivity gradually increases with soil crack development.

THREE-DIMENSIONAL ELLIPTIC CRACK UNDER IMPACT LOADING

The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.

Interactions of penny-shaped cracks in three-dimensional solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.

Simulation of three-dimensional tension-induced cracks based on cracking potential function-incorporated extended finite element method

In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.

Analysis of Three-dimensional Crack Propagation by Using Displacement Discontinuity Method

A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.

A modified maximum tangential tensile stress criterion for three-dimensional crack propagation

The three-dimensional （3D） crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterion, it is supposed that cracks propagate only at crack front in the principal normal plane. The tangential tensile stress at crack front in the principal normal plane in local coordinates is employed to determine crack propagation, which is calculated through coordinate transformation from global to local coordinates. New cracks will propagate when the maximum tangential tensile stress at crack front in the principal normal plane reaches the tensile strength of rock-like materials. Compared with the previous crack propagation criteria, the modified crack propagation criterion is helpful in calculating 3D crack stress intensity factor, and can overcome the limitations of propagation step determined by individual experiences in previous studies. Finally, the 3D crack propagation process is traced by element-free Galerkin method. The numerical results agree well with the experimental ones for a frozen resin sample with prefabricated 3D cracks.

INTERFACIAL CRACK ANALYSIS IN THREE-DIMENSIONAL TRANSVERSELY ISOTROPIC BI-MATERIALS BY BOUNDARY INTEGRAL EQUATION METHOD

The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.

A STUDY ON STRESS INTENSITY FACTORS AND SINGULAR STRESS FIELDS OF THREE-DIMENSIONAL INTERFACE CRACK

By using the finite-part integral concepts and limit technique,the hypersingular inte- grodifferential equations ofthree-dimensional(3D)planar interface crack were obtained; then thedominant-part analysis of 2D hypersingular integral was further usedto investigate the stress fields near the crack front theoretically,and the accurate formulae were obtained for the singular stressfields and the complex stress intensity factors.

CONVOLUTION OF THE IMPACT THREE-DIMENSIONAL ELASTO-DYNAMICS AND DYNAMIC STRESS INTENSITY FACTOR FOR AN ELLIPTIC CRACK

This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method.

Discussion on Green Development of Fenlong for Yield Increase, Quality Enhancing, Water Retaining and Multiple Use of Natural Resources

Fenlong green ecological agriculture technology （Fenlong technology）, a new smash ridging farming method developed by Guangxi Academy of Agdcultural Sciences, has been elected as the recommended cultivation technique by the Ministry of Agriculture of China. It replaces the traditional plowshare with spiral drill, and its tilth depth is twice deeper than that by tractor tilthing. It also extends soil nutrient, moisture, oxygen and microorganism, the so-called ＂Four pools＂. Soil nutrient, oxygen, microorganism, light and rainfall use ratio is increased by 10%-100%, creating a platform for natural increase of more than 10% of crop yield. Its application to over 20 kinds of crops in 21 provinces has proved that the yield increases 10-30% with quality enhancing 5% and double water retaining capacity but no more input. When the application area of Fenlong could reach 67 million hm2, the amount of fertilizer can be reduced by 40-50 billion kg, saving 120-150 billion Yuan. In this paper, we put forward the strategy of ＂4＋1＂ （arable, saline-alkali soil, grasslands, Sponge City ＋ rivers） green development in China, and deepened the Fenlong cultivated tilled layer from 16.5 cm to 35 cm for 67 million hm2 arable land, ridged 13.3 million hm2 of saline-alkali soil for 35 cm, and also 35 cm for 67 million hm2 degraded steppe, which could have the following 3 effects： first, the 147 million hm2 of land with Fenlong cultivation could increase loosing soil to 315.491 billion m3, in＊ creasing by 159.26% for 120 million hm2 of arable land with the average tilled layer of 16.5 cm, which has loosing soil of only 198.1 billion m3, that is, the space of the land increases 1.6 times. Second, every hectare of plowland could store up to 450 m3/hm2 of natural rainfall, and the unused 60 m3 of saline-alkali soil and grasslands could store water of 102 billion m3, showing an increase of over 88.89% for the current plowland storage of 54 billion m3 at now, that is, double the natural rainfall storage capacity. Third, the two multiple increase of natural resources application can bring trillions of resource activation, environmental cleaning, food security, citizens, health, economic, ecological and social benefits, and makes the Chinese nation move forward in green development. Its application in ＂big scientific research＂ and ＂One Belt And One Road＂ will contribute Chinese strength to the world.

Development of RSDS-III Technology for Ultra-Low-Sulfur Gasoline Production

The 3rd generation catalytic cracking naphtha selective hydrodesulfurization(RSDS-III) technology developed by RIPP included the catalysts selective adjusting(RSAT) technology, the development of new catalysts and optimized process conditions. The pilot plant test results showed that the RSDS-III technology could be adapted to different feedstocks. The sulfur content dropped from 600 μg/g and 631 μg/g to 7 μg/g and 9 μg/g, respectively, by RSDS-III technology when feed A and feed B were processed to meet China national V gasoline standard, with the RON loss of products equating to 0.9 units and 1.0 unit, respectively. While the feed C with a medium sulfur content was processed according to the full-range naphtha hydrotreating technology, the sulfur content dropped from 357 μg/g in the feed to 10 μg/g in gasoline, with the RON loss of product decreased by only 0.6 units. Thanks to the high HDS activity and good selectivity of RSDS-III technology, the ultra-low-sulfur gasoline meeting China V standard could be produced by the RSDS-III technology with little RON loss.

Three-dimensional Operation Design of Forest Engineering——A Case Study of Bagua Zhou in Nanjing City

[Objective] The aim was to construct a virtuous eco-system covering multi-populations, multi-layers and multi-products and to explore timber production of physical and non-physical products, in order to extend products operated in forestry engineering and seek a novel model for forestry engineering. [Method] A three-dimensional operation project was designed based on forestry lands in Bagua Zhou in Nanjing City. [Result] In the project, timber products, by-products in forest and relaxing products supplimented and supported each other, and a virtuous circle has been achieved. [Conclusion] The novel model makes simultaneous relaxing and working possible, which creat vaule together.

配筋率对钢箱-UHPC组合梁裂缝特征的影响

为研究钢箱-超高性能混凝土(UHPC)组合梁负弯矩区裂缝特征,设计制作了3根配筋率分别为1%、2%和3%的UHPC翼板部分充填砼窄幅钢箱组合梁,通过对不同配筋率的试验梁进行反向加载试验研究组合梁的裂缝发展特征,并对最大裂缝宽度、裂缝间距及开裂荷载进行对比分析。研究结果表明:试验梁在不同配筋率条件下,UHPC面层裂缝发展明显不同,配筋率越高最大裂缝宽度越小,裂缝越密集,开裂荷载越大;与配筋率为1%的试验梁相比,配筋率为2%和3%的试验梁的开裂弯矩分别提高了30.3%和65.2%;可根据《公路钢筋混凝土及预应力混凝土桥涵设计规范》(JTG 3362—2018)中的裂缝宽度计算公式乘以0.78的修正系数,来计算窄幅钢箱-UHPC组合梁负弯矩区的裂缝宽度。

滑坡诱发农村山区砖混结构变形及演化规律研究

作为西南山区最为频发的地质灾害之一,滑坡对山区中大量存在的砖混房屋造成严重影响,尤其农村山区中受灾最为严重,给人民群众生产生活造成巨大安全隐患.研究滑坡引起的上部砖混结构变形特征,对科学指导房屋设计防护及灾后监测点布置具有重要现实意义.本文以酉阳高园子滑坡为例,通过现场调查和裂缝参数拟合研究了滑坡区砖混结构变形破坏特征.同时基于PFC3D,构建了有限体积-离散元(Finite Volume Method-Discrete Element Method,FVM-DEM)滑坡-房屋单向耦合模型,模拟并分析了在不同工况下滑坡动态发育过程,进一步揭示了房屋裂缝产生及演化过程.结果表明:相比于普通无柱砖混结构,底框结构抵抗滑坡诱发变形破坏的能力更强,并且处于前后缘段的房屋相对中段房屋更易出现损伤破坏;房屋的长宽比不宜过小,布置房屋长边走向平行斜坡走向有利于减少滑坡对房屋造成的破坏.本文总结了在滑坡作用下,上部砖混结构的变形特征和裂缝开展规律,可为农村山区砖混结构设计和布置提供参考.

沥青路面反射裂缝影响因素分析及抗裂性能提升技术研究

半刚性基层开裂一直是世界性难题,在荷载、温度、湿度等耦合作用下逐步演变为沥青路面裂缝。路面开裂已经成为我国最为典型的路面病害类型,严重影响路面使用的耐久性。本研究分析了荷载场、温度场以及层间接触状态对裂缝产生与发展的影响,通过对沥青路面混合料类型、胶结料种类进行比选研究,以及层间接触及粘结状态对抗裂影响分析,提出消能层、材料抗裂功能提升、高粘层间结合状态抗裂技术措施,将抗裂技术效果与沥青路面结构层相匹配,进而大幅度提升路面耐久性。