• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Effect of proppant pumping schedule on the proppant placement for supercritical CO2 fracturing

    2022-06-02 05:00:16YongZhengHiZhuWngYingJieLiGngHuTinBingYngChengMingZhoMingShengLiu
    Petroleum Science 2022年2期

    Yong Zheng ,Hi-Zhu Wng ,*,Ying-Jie Li ,Gng-Hu Tin ,Bing Yng ,Cheng-Ming Zho ,Ming-Sheng Liu

    a State Key Laboratory of Petroleum Resources and Prospecting,China University of Petroleum (Beijing),Beijing,102249,China

    b CNPC Offshore Engineering Co.Ltd.,Beijing,100028,China

    Keywords:Supercritical CO2 fracturing Proppant placement CFD-DEM Unconventional natural gas CO2 utilization

    ABSTRACT Supercritical CO2 fracturing is a potential waterless fracturing technique which shows great merits in eliminating reservoir damage,improving shale gas recovery and storing CO2 underground.Deep insight into the proppant-transport behavior of CO2 is required to better apply this technique in the engineering field.In the present paper,we adopted a coupled Computational Fluid Dynamics and Discrete Element Method (CFD-DEM) approach to simulate the proppant transport in a fracking fracture with multiple perforation tunnels.Previous experiments were first simulated to benchmark the CFD-EDM approach,and then various pumping schedules and injection parameters (injection location,multi-concentration injection order,multi-density injection order and injection temperature) were investigated to determine the placement characteristics of proppant.Results indicate that the swirling vortex below the injection tunnels dominates the proppant diffusion in the fracture.The velocity of fluid flow across the proppant bank surface in multi-concentration injection shows a positive correlation with the proppant concentration.Injecting high-density proppant first can promote the transportation of low-density proppant injected later in the fracture to a certain extent.Decreasing the initial injection temperature of supercritical CO2 slurry helps enhance the particle-driving effect of fluid and improve the performance of supercritical CO2 in carrying proppant.

    1.Introduction

    Over the past few decades,as a result of rapid economic development,industrialization and urbanization,global CO2emissions have continued to increase,and the environmental challenges posed by rising levels of CO2in the atmosphere have become a matter of growing public concern.It is generally accepted that a reduction of the CO2emissions is an effective way to promote climate change.Hence,many CO2reduction technologies and methods have been proposed worldwide,among which carbon capture and storage (CCS) is one of the important strategies(Holloway,2005;Leung et al.,2014).However,due to concerns about the high cost of CCS,carbon capture utilization and storage(CCUS) technology has been developed,which enables resource utilization while storing CO2(Cuéllar-Franca and Azapagic,2015;Li et al.,2013;Norhasyima and Mahlia,2018),making it a strong prospect.

    Among various utilization ways of CO2,the application in oil and natural gas industry is a crucial component (Norhasyima and Mahlia,2018).Enhancing oil recovery by injecting CO2into depleted reservoirs (CO2-EOR) is a well-known practice in oil and natural gas industry and has developed into a mature technology(Alvarado and Manrique,2010).CO2-EOR technology is likely to be suitable for more than 90% of the world's estimated ultimately recoverable oil(Godec,2011),and it is also considered an option for permanent CO2storage(Dai et al.,2014).Thus,the use of CO2-EOR has the dual advantages of reducing CO2emissions and increasing oil production.Furthermore,CO2also plays an important role in the extraction of unconventional natural gas such as shale gas(Iddphonce et al.,2020).The application of hydraulic fracturing has led to a boom in shale gas production over the past few years,yet the significant environmental controversy it faces has prompted the industry and researchers to explore waterless fracturing fluids to reduce water use (Clark et al.,2013;Middleton et al.,2017;Vengosh et al.,2014).The potential effectiveness and unique advantages of CO2make it one of the most interested waterless fracturing fluids.

    Similar to hydraulic fracturing,liquid or supercritical CO2is injected at high pressure into the shale formation at a flow that exceeds the absorption capacity of the formation to open the existing fractures or initiate new fractures,thereby increasing conductivity of the shale formation.Under the temperaturepressure conditions of shale reservoirs,CO2can easily reach the supercritical state (pressure over 7.38 MPa,temperature over 304.1 K).Compared with water-based fracturing fluids,supercritical CO2as a working fluid has several significant advantages.Firstly,supercritical CO2fracturing is inclined to form a complex network of fractures in the reservoir due to its low viscosity,and enhance the reservoir stimulation effect(Chen et al.,2021b;Li et al.,2018;Song et al.,2019;Yang et al.,2021).At the same time,there is no water in the fluid,which will not cause water block and swelling of clay minerals that can impair the hydrocarbon migration pathway(Bahrami et al.,2012;Dehghanpour et al.,2012;Middleton et al.,2015;Naik et al.,2019).Moreover,CO2has an obvious tendency to adsorb to shale compared to methane (Hamza et al.,2021;Huo et al.,2017;Wang et al.,2021),which means that it can replace the methane adsorbed in shale and achieve permanent storage of CO2while increasing shale gas production and recovery.Based on these potential advantages,supercritical CO2fracturing is considered to be a waterless fracturing technology with strong growth prospects (Chen et al.,2021a;Middleton et al.,2015;Mojid et al.,2021;Wang et al.,2012,2019).

    In addition,fracturing the reservoir will require the pumping of supercritical CO2slurry,a mixture of supercritical CO2fluid and proppant particles,into the fracture,which will determine the future production performance of the well(Rightmire et al.,2005).A proppant is a solid particle used to support fractures created by fracturing to prevent them from closing again after the operation is completed.It works by forming a support fracture zone that is much higher than the reservoir permeability,providing a highly conductive channel for the flow of oil and gas from the reservoir to the wellbore.In order to obtain optimal fracture conductivity and to predict the productivity of the fractured well,the transport behavior and placement of proppant in the fracture has been studied extensively(Qu et al.,2021;Roostaei et al.,2020;Sahai and Moghanloo,2019).However,it is noteworthy that the majority of these existing studies are for conventional fracturing fluids,i.e.,gel and slickwater,and relatively little effort has been conducted in the past on the proppant carried by supercritical CO2.

    Carrying the proppant is a major challenge facing supercritical CO2fracturing technology currently(Mojid et al.,2021;Wang et al.,2020).The low viscosity characteristic of supercritical CO2,while favoring the induction of complex fracture networks,also limits its ability to carry proppant and thus largely hinders its widespread application.Although CO2fracturing has been tested and explored by the industry for a long time (Gupta and Bobier,1998),research into its ability to carry proppant for placement has only begun in recent years (Sun et al.,2018;Wang et al.,2018;Wang and Elsworth,2020;Zheng et al.,2020).What's more,it is surprising that little research has been reported on supercritical CO2slurry pumping schedule,and the effect of different pumping schedule on proppant placement is poorly understood.

    In the present study,we studied the proppant transport and placement in a planar fracture with perforation tunnel under different pumping schedule using developed CFD-DEM model.Further validation of the model is performed with experimental results.Based on the model,the effects of injection location,multiconcentration injection order,multi-density injection order and injection temperature on the proppant placement are investigated.The results of this study are to provide a better understanding of the proppant placement in supercritical CO2fracturing and can assist fracturing engineers solve problems related to proppant pumping design in fracturing.

    2.Numerical model

    2.1.Mathematics model

    A CFD-DEM based multiphase model was established to simulate the flow of supercritical CO2slurry.The fluid phase is described by Navier-Stokes equation based on an Eulerian approach in CFDDEM model.The DEM method is used to track the proppant particles in the system,and Newton's second law is applied to govern its movement.A detailed description of this model can be found in our previous work (Zheng et al.,2020).However,in the previous model,the heat transfer between fluid-matrix is neglected in simulation.In the present study,we further considered the heat transfer between the supercritical CO2fluid and the fracture wall.The energy equation for fluid flow and fracture wall are respectively written as:

    2.2.CO2 properties

    The physical properties of CO2fluid are very sensitive to the changes of temperature and pressure,and CO2will reach supercritical state under the reservoir temperature and pressure.Therefore,it is important to effectively calculate the physical properties of supercritical CO2to study the flow of supercritical CO2slurry within the fracture.The parameters of CO2properties used in the model include density,viscosity and thermal conductivity.The widely accepted Span-Wagner equation of state(Span and Wagner,1996) is used to calculate the fluid density,and the equation is adopted and recommended by the National Institute of Standards and Technology (NIST).The viscosity and thermal conductivity of the fluid are calculated using the explicit models developed by Heidaryan et al.(2011) and Jarrahian and Heidaryan (2012) for calculating CO2in the supercritical state,respectively.The average absolute relative error(AARE)of the two models are 1.82%and 2.4%,respectively,compared with NIST web book data,and the accuracy can meet the engineering calculation requirements.The use of explicit CO2physical property calculation models is more efficient because they do not require a large number of parameters and complex calculations.

    The viscosity can be formulated as (Heidaryan et al.,2011):

    where the viscosity is expressed in centipoise (mPa·s);and the temperature and pressure are expressed in K and bar,respectively.

    The thermal conductivity can be formulated as (Jarrahian and Heidaryan,2012):

    where the temperature and pressure are expressed in K and MPa,respectively;λ is the thermal conductivity,mW/(m·K);A and C are tuned coefficients,the values of which can be found in the corresponding references.

    2.3.Model geometry

    In reservoir fracturing,it is generally assumed that there are symmetrical plane fractures with two wings.Fracturing operations usually result in irregular and tortuous fractures.However,in previous numerical models of proppant transport,these fractures are usually represented as a vertical planar,which is a widely used and effective simplification (Sahai and Moghanloo,2019).Therefore,a rectangular channel with perforation tunnel is used to mimic the supercritical CO2slurry injection process in the vertical fracture of a perforation completion.Due to the large consumption of computational resources by CFD-DEM method,a small-scale fracture characterization unit is used in this paper to simulate the transport process of proppant within artificial fractures.The fracture geometry is 500 mm in length,150 mm in height,and 2 mm in width.Four injection points 3 mm high and 2 mm wide are on the left boundary,whose centers are separated by 12 mm,as shown in Fig.1.On the right side,the fluid can flow out along the height of the fracture,but proppant particles will be left in the fracture.Furthermore,the effects of fracture propagation and fluid leak-off are neglected in the model.

    2.4.Simulation conditions

    Fig.1.Sketch of the fracture geometry (not to scale).

    At the inlets,the supercritical CO2slurry enters the computational domain with a specified velocity and temperature.The velocity and temperature of the proppant particles at the initial injection are the same as those of the fluid.In a conventional hydraulic fracturing operation,the average velocity of slurry in the fracture is mostly below 0.25 m/s,and the flow velocity of slurry in such fractures has been discussed by others before (Alotaibi et al.,2019).We choose the higher value,0.16 m/s,as the average velocity of slurry in the fracture in this simulation.Therefore,for the fractures(150 mm high and 2 mm wide)and perforation clusters(4 perforations of 3 mm high and 2 mm wide)used in our simulation,the fluid velocity in each perforation is approximately 2.0 m/s.In this paper,uniform spheres are used to represent proppant particles.The fluid flows out under constant pressure boundary conditions.The external pressure at the outlet is set at 18 MPa.The proppant concentration used in the simulation is the ratio of proppant particle volume to fluid volume.In order to explore the effect of heat transfer between fluid and fracture wall on proppant placement,the initial temperature of fluid injection is lower than the wall temperature in the simulation,and this setting is also similar to the actual flow situation in operation(Meng et al.,2018).Therefore,the temperature of fracture wall is fixed at 358 K,which is higher than the initial temperature of supercritical CO2slurry.

    The fluid governing equation is solved by the finite volume method.The Phase Coupled SIMPLE scheme (Patankar,1980) is used to treat the coupling between pressure and momentum.The momentum,volume fraction and energy equations are discretized by the first order upwind difference scheme.The wall of the fracture model adopts the non-slip fluid boundary condition,that is,assuming that the velocity of the fluid on the fracture wall is zero.The time step for the fluid phase is 1 × 10-3s while the time step for the solid phase is 1 × 10-5s.The Young's modulus,restitution coefficient,Poisson's ratio and friction coefficient of the proppant used by Baldini et al.(2018)were used in the simulation.The elastic parameters of the reservoir rock are obtained from the work of Moghadam et al.(2019).The values of some main parameters used in simulations are listed in Table 1.

    3.Results and analysis

    3.1.Model validation

    This CFD-DEM model has been validated in our previous work,which has an efficient performance in simulating the proppant transport behavior within the fracture (Zheng et al.,2020).However,for the sake of completeness,we further compared the proppant bank profiles in the simulation with the experimental results of Tong and Mohanty (2016) on the basis of the previous verification.The length,height and width of the main fracture in the experiment are 381,76.2,and 2 mm,respectively.The anglebetween the secondary fracture and the main fracture is 90°,the length of the secondary fracture is 190.5 mm,and the height and width are the same as that of the main fracture.The slickwater slurries were injected from the right side of the fracture.The fracture size used in the simulation is the same as that in the laboratory experiment.The basic parameters used in model validation are listed in Table 2.

    Table 1Numerical parameters.

    Table 2Parameters used in model validation.

    Fig.2 shows the profiles of the proppant dune in the main fracture and bypass fracture from numerical results and experimental data under different time moments.It is clear that the proppant distribution obtained from the numerical simulation is in line with the experimental observation.The discrepancy is mainly due to the lack of understanding of the composition and properties of the proppant and fracture wall materials in the experiment,which is needed for CFD-DEM calculation.In addition,the proppant used in the simulation was a sphere of consistent diameter,whereas the experiments used 20/40 mesh sand,which again introduces some error.Considering the simulation bias caused by some parameters in the simulation that need to be estimated,the CFD-DEM model is reliable to study the fluid-solid flow characteristics of the proppant slurry in the fracture.

    3.2.Effect of injection location

    In fracturing,the location of the perforation is usually selected based on reservoir characteristics as well as the geomechanical properties of a formation in order to obtain better productivity.To study the effect of injection location on proppant placement,we simulated the supercritical CO2slurry injection at three different heights on the fracture inlet side.The distances from the center of the three injection clusters to the bottom are,respectively,112.5(top),74.5 (middle),and 37.5 mm (bottom).In each case,the injection temperature is 338 K,the proppant density is 1250 kg/m3,the injection proppant concentration is 0.04 and the injection time 20 s.Other parameters are shown in Table 1.

    Table 3Design of pumping schedules.

    Fig.3 shows a comparison of the proppant bank building process under the three injection locations.It is clear from the results shown in the figure that different perforation positions have a greater impact on the transport and placement of proppant.The top injection causes the proppant particles to settle in a relatively dispersed manner,resulting in a more uniform proppant bank distribution than middle and bottom injection.This is because the injected proppant slurry and the surrounding supercritical CO2fluid will generate a large number of vortices under the effect of viscous shear(Longmire and Eaton,1992),and then form a swirling vortex.The motion of swirling vortex promotes the diffusion of proppant particles in the fracture.This phenomenon was also observed by Baldini et al.(2018) in their study of slickwater transport proppant.However,since the viscosity of supercritical CO2is considerably lower than that of slickwater,the particles are not carried by the swirling vortex to be dispersed above the injected stream in our simulation.Therefore,the swirling vortex that form beneath the injected beam influence the development of the proppant bank.

    Fig.2.Comparison of proppant bank profiles in the main fracture (left) and bypass fracture (right) in simulation result and experimental data.

    Fig.3.The proppant placement process for different injection locations.

    Fig.4.The streamline of top injection (first column),middle injection (second column) and bottom injection (third column).

    To clarify the distribution of swirling vortices in the fracture at the three injection positions,we have given the corresponding streamlines on the center plane of the fracture,as shown in Fig.4.It can be observed that the swirling vortex is mainly below the injection point in the top injection(first column),while in the middle injection(second column)and bottom injection(third column)it is above the injection point.With time,the vortex area in the top injection decreases,and the vortex area in the middle and bottom injection tends to be constant.At the initial stage of the middle injection,there is a small swirling vortex under the slurry stream,and then it gradually disappears.However,the dispersion of the proppant particles within the fluid is also affected,resulting in the formation of relatively uniform proppant bank compared to bottom injection.In the bottom injection,due to the limitation of space below the injection point,the vortex is weak,and the proppant mainly moves along the dune surface by the erosion of the fluid.In addition,when the fluid reaches a certain velocity,a phenomenon of the proppant making a circular motion at the tail of the bank occurs.

    In conclusion,the swirling vortex generated below the injection point when injecting supercritical CO2slurry can increase the diffusion of the proppant,which in turn leads to a more uniform spreading of the proppant bank near the wellbore.After the vortex below the injection point disappears as the proppant bank grows,the proppants are fluid-driven and spread deeper into the fracture along the bank surface.Therefore,top injection is an effective way to improve the placement of the proppant in the fracture.

    3.3.Effect of multi-concentration injection order

    In this section,the placement of proppant under multiconcentration injection orders are investigated.Two different multi-concentration proppant injection orders are designed as shown Table 3.In order 1,the proppant concentration increases from stage 1 to stage 4.In order 2,the proppant concentration decreases from stage 1 to stage 4.The number of proppant particles injected in each stage of the two schemes is the same,so the total injection time of each stage is different.Other parameters are shown in Table 1,in which the injection temperature and proppant density are fixed at 338 K and 1250 kg/m3,respectively.

    Fig.5 shows the fluid velocity contours and proppant placement in the fracture after each injection stage in different injection orders.It can be observed that the shape of the proppant bank in the fracture is basically similar in each stage of the two orders.Compared to order 1,proppant transport in stage 1 of order 2 show longer distance under the same injection velocity conditions,which induce more smoothly proppant distribution in the fracture.This is because the proppant particles in the slurry with high concentration have higher velocity in the direction of flow and a lower settling velocity in the vertical direction in the initial injection stage(Blyton et al.,2015).Furthermore,from the fluid velocity contours in the fracture illustrated in Fig.5,we can observe that from stage 2 of the two injection orders,the proppant dunes gradually approached the high velocity flow region formed by supercritical CO2slurry injection,and the surface of the dunes began to be washed by the carrying fluid.This means that the shape of the proppant bank formed in the previous stage will also affect the proppant placement in the next stage by changing the distribution of the high-speed flow region.

    To determine the effect of multi-concentration injection order on the proppant placement in the fracture,we present the final proppant bank profile of each stage and the X-direction component of fluid velocity in the moving layer of the dune surface along the centerline,as shown in Fig.6.According to the coordinate system established in Fig.1,the velocity component in the X-direction is the average velocity of the fluid in the horizontal flow direction on the centerline.It can be observed that order 2 has a clear advantage in the spread distance of proppant banks,which is 11.21% higher than that in order 1.Furthermore,F(xiàn)ig.6 shows that there is a direct relationship between the proppant bank shape and the fluid velocity in the X-direction.For the fluid velocity in the X-direction,the flow to the fracture tip is positive,while flow to the wellbore is negative.The dashed line in the figure shows that there are two flow directions for the fluid near the surface of the proppant bank.In stage 1,the cutoff point for these two opposite directions of flow is approximately in the middle of the dune,and as the dune grows,the cutoff point in subsequent stages begins to move toward the wellbore.

    At the same time,we can find from Fig.6 that the fluid velocity in order 2 is higher than in order 1 in both stages 1 and 2,whether it is positive or negative.However,in stages 3 and 4,an inverse relationship displays:the fluid velocity in order 1 begins to be greater than in order 2.Combining the proppant concentration of the two injection orders in each stage,we can conclude that the velocity of the proppant bank surface fluid is proportional to the proppant injection concentration.Moreover,it can be also observed from Fig.6 that after a cone-shaped dune appears in the fracture(starting from stage 2),increasing the proppant concentration mainly affects the fluid velocity near the dune front.In the final stage 4 of the simulation,although the fluid velocity on the dune front in order 1 is about twice as high as that in order 2 due to the high injection concentration,the bank spreading distance is still shorter.In fracturing operation,proppant concentration is mostly injected in the way of step increase,while this study shows that a longer proppant spreading distance may be obtained by using a step decrease concentration injection method under the premise of effectively avoiding sand plugging at the fracture entrance.Therefore,the proppant concentration combination should be reasonably selected in the pumping design to obtain a longer dune placement distance.

    Fig.5.Fluid velocity contours and proppant placement in the fracture after each injection stage.

    Fig.6.Final proppant bank profile (solid line) of each stage and the X-direction component of fluid velocity in the moving layer of the dune surface along the centerline (dashed line).

    3.4.Effect of multi-density injection order

    For pumping slurry,the proppant with low density can reach the deeper fracture with the slurry,because its settling velocity in the fluid is decreased.However,the compressive strength of proppant also weaken with the decrease of density,which makes it easy to break when the fracture is closed,resulting in reduce of fracture conductivity.To investigate the effect of multi-density injection order on proppant placement,in this section we have designed two proppant pumping schedules with different injection orders as shown in Table 4.In order 1,the proppant density increases from stage 1 to stage 3.In order 2,the proppant density decreases from stage 1 to stage 3.Other parameters are shown in Table 1,in which the injection temperature and proppant concentration are fixed at 338 K and 0.04,respectively.

    Table 4Design of pumping schedules.

    To illustrate the placement process of the proppant under multidensity injection order,we present the proppant bank distributionafter each injection stage in Fig.7.It is evident that proppants of different densities show significantly different distribution structures in the fracture under the two designed injection orders.For the proppant placement in order 1,the low-density proppant spreads a relatively uniform proppant bed at the bottom of the fracture in the initial stage.Then,the medium-density proppant injected in stage 2 basically settles and accumulates on the proppant bed formed in stage 1 and does not over the peak of this bed.Subsequently,the high-density proppant injected in stage 3 cover the highest point of the dunes in the previous two stages and travel deep into the fracture,but still does not exceed the laying distance of the low-density proppant bed.In order 2,the high-density proppant first injected forms a distinct dune near the injection points.Then,the medium density proppant injected in the secondary stage completely cover the proppant dune and maintain the similar shape as before.When the low-density proppant was injected in the third stage,most of the proppant was washed to the back side of the dune except a small part of them accumulated in the front side of the dune,which promoted the proppant to be spread deep into the fracture.Comparing the placement process of proppant bank under the two injection schedules,it can be seen that there is obvious scouring phenomenon of the proppant dune in the stage 3 of the injection order 2,which makes the proppant easier to migrate to the depth of the fracture.

    Fig.7.The proppant bank distribution in the fracture after each injection stage.

    Additionally,we also extracted the profiles of the proppant bank along the fracture length direction after the end of each injection stage in the two schedules,to determine the effect of the injection order of different density proppant on the proppant bank spread,as shown in Fig.8.It can be clearly seen from Fig.8 that after the completion of the three injection stages in order 1,the placement distance of proppant in the fracture does not increase noticeably,which are 367,391,and 392 mm,respectively,which increased by 4.08%compared with the laying distance in order 1.In order 2,the proppant bank placement distance between different injection stages increased greatly,and the spread distance are 297,338,and 408 mm,respectively.The high-density proppant injected in order 1 can inhibit the scouring development to some extent by increasing the driving force of proppant particles needed,reduce the transport efficiency of the proppant.In contrast,the steep dune formed by the high-density proppant first injected in order 2 increased the velocity of the fluid flowing through it,enhancing the lifting effect on subsequent low-density proppant particles.Moreover,low-density proppant is more likely to be washed away when it is on top of the dune.Thereby,when pumping by multi-density injection,on the premise of avoiding premature accumulation of proppant bank to block the entrance,injecting high-density proppant first can,to a certain extent,facilitate the transport of the lowdensity proppant injected later in the fracture.

    Fig.8.The profiles of the proppant bank along the fracture length direction after the end of each stage in multi-density injection orders.

    3.5.Effect of injection temperature

    In a typical fracturing operation,the temperature of the injected fluid is lower than that of the reservoir.As the fluid flows in the fractures,heat transfer occurs between the fracturing fluid and the formation rock,resulting in variations in the properties of supercritical CO2.In this section,three sets of injection temperature,328,338,and 348 K are designed to investigate the effect of fluid temperature on proppant distribution under otherwise the same conditions.In each case,the proppant density and injection proppant concentration are 2650 kg/m3and 0.04,respectively.It is generally believed that viscosity is an important factor reflecting the proppant carrying capacity of the fluid.Fig.9 shows the viscosity distribution of supercritical CO2fluid and proppant placement in the fracture at different injection temperatures.It can be seen from the figure that for the same reservoir,the proppant slurry injection temperature will have a significant impact on the viscosity of supercritical CO2in the fracture.The viscosity of supercritical CO2depends on ambient temperature and pressure.With the same total injection time,outlet pressure and constant wall temperature,the temperature of supercritical CO2is directly related to the injection temperature when flowing in the fracture.Therefore,the viscosity of the fluid whose initial injection temperature is 328 K is higher than the other two temperatures in the fracture,which affects the placement of proppant.

    Fig.9.Viscosity distribution and proppant placement of supercritical CO2 fluid in fractures at different injection temperatures.

    Fig.10.Placement distance and the fluid viscosity of supercritical CO2 at different injection temperatures.

    To quantitatively compare the placement of the proppant in the three cases,the placement distance of proppant bank in the fracture and the fluid viscosity of supercritical CO2at different injection temperatures are given,as shown in Fig.10.The fluid viscosity is the average value of supercritical CO2fluid at the fluidization layer on the proppant bank surface.It can be seen from the figure that as the temperature of the injected slurry decreases from 348 to 338 K and 328 K,the viscosity of the fluid carrying proppant particles increases by 8.1% and 18.75%,respectively.At the same time,the placement distance of proppant bank in the fracture also increases,increases by 1.70%and 4.10%,respectively.This shows that reducing the injected slurry temperature can improve the viscosity of supercritical CO2fluid,increase its driving force on proppant particles,enhance the carrying capacity,and obtain a longer proppant placement distance.

    4.Conclusions

    In the present paper,a coupled CFD-DEM method was adopted to study the proppant transport and placement in a planar fracture with perforation tunnel under different pumping schedule during supercritical CO2fracturing.The effects of several groups of different pumping schemes on the placement of the proppant is analyzed,which has some guiding significance for the design of the proppant pumping schedule.Main conclusions are drawn as follows:

    (1) Different from conventional fracturing fluids,when supercritical CO2slurry is injected,it is the swirling vortex formed below the injection point that affects the diffusion of the proppant within the fracture,which in turn results in a more uniform distribution of proppant bank near the wellbore.The extent of the swirling vortex is positively correlated with the injection height,thus supercritical CO2slurry injection in the upper part of the reservoir is an effective way to improve the distribution of the proppant in the fracture.

    (2) For multi-concentration injection,the velocity of the fluid at the proppant bank surface is positively correlated with the concentration of the proppant injected.Under the conditions in this paper,a longer proppant bank placement distance is easily obtained with proppant injection concentration from high to low.

    (3) During multi-density injection schedule,the high-density proppant injected first creates a prominent dune that will increase the velocity of the fluid as it flows through the dune,and the low-density proppant injected later is more easily carried deeper into the fracture by the fluid washout,creating a longer proppant bank.

    (4) Proppant-carrying performance of supercritical CO2can be enhanced by decreasing the initial injection temperature.The lower initial temperature of supercritical CO2slurry contributes to obtain higher fluid viscosity,enhancing the driving effect of fluid on particles.

    Acknowledgements

    The authors would like to acknowledge the National Natural Science Foundation of China (Grant No.51874318,51922107,and 41961144026).Besides,support from the National Key Scientific Research Instrument Research Project of NSFC (Grant No.51827804) is appreciated.

    听说在线观看完整版免费高清| 亚洲人成网站高清观看| 国产成人一区二区三区免费视频网站| 19禁男女啪啪无遮挡网站| 欧美在线黄色| 国产成人欧美| 757午夜福利合集在线观看| 国产伦一二天堂av在线观看| 久久热在线av| 成熟少妇高潮喷水视频| av片东京热男人的天堂| 一级作爱视频免费观看| 麻豆av在线久日| 免费搜索国产男女视频| 狠狠狠狠99中文字幕| 国产亚洲精品久久久久5区| 色综合站精品国产| 中文字幕人妻丝袜一区二区| 亚洲 国产 在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人av激情在线播放| 99热这里只有精品一区 | 中文字幕高清在线视频| 色尼玛亚洲综合影院| 在线av久久热| 国产精品亚洲av一区麻豆| 免费在线观看黄色视频的| 男女那种视频在线观看| 免费看十八禁软件| 久久精品国产亚洲av高清一级| 欧美日韩乱码在线| 亚洲九九香蕉| 午夜福利免费观看在线| 美女高潮到喷水免费观看| 久久久国产精品麻豆| 黑丝袜美女国产一区| 久久久久久久午夜电影| 欧美性猛交黑人性爽| 免费高清视频大片| 免费在线观看完整版高清| 国产一卡二卡三卡精品| 岛国视频午夜一区免费看| 岛国视频午夜一区免费看| 亚洲av美国av| 免费电影在线观看免费观看| 亚洲成人国产一区在线观看| avwww免费| 国产主播在线观看一区二区| 欧美黑人精品巨大| 国产不卡一卡二| 美女免费视频网站| 国产av在哪里看| 亚洲精品在线观看二区| 欧美av亚洲av综合av国产av| 自线自在国产av| 麻豆久久精品国产亚洲av| 亚洲国产高清在线一区二区三 | 久热这里只有精品99| 女性生殖器流出的白浆| 亚洲专区中文字幕在线| 999久久久国产精品视频| 日韩精品青青久久久久久| 男人舔奶头视频| 国产又色又爽无遮挡免费看| 国产主播在线观看一区二区| 久久精品国产清高在天天线| 国内揄拍国产精品人妻在线 | 欧美激情极品国产一区二区三区| 国产精品一区二区三区四区久久 | 91国产中文字幕| 中亚洲国语对白在线视频| 精品免费久久久久久久清纯| 精品久久久久久久毛片微露脸| 黑丝袜美女国产一区| 搡老妇女老女人老熟妇| 亚洲精品中文字幕一二三四区| 51午夜福利影视在线观看| 国产精品久久久久久人妻精品电影| 国产高清有码在线观看视频 | 欧美日本亚洲视频在线播放| 亚洲中文字幕日韩| 亚洲 欧美 日韩 在线 免费| 国产午夜福利久久久久久| ponron亚洲| 自线自在国产av| 日韩欧美国产一区二区入口| 最好的美女福利视频网| 美女大奶头视频| avwww免费| 免费av毛片视频| 天堂动漫精品| 宅男免费午夜| 免费高清视频大片| 亚洲无线在线观看| 制服丝袜大香蕉在线| 免费高清在线观看日韩| 51午夜福利影视在线观看| 夜夜爽天天搞| 男人的好看免费观看在线视频 | 别揉我奶头~嗯~啊~动态视频| 亚洲色图 男人天堂 中文字幕| 欧美色视频一区免费| 可以在线观看的亚洲视频| 在线av久久热| 可以在线观看毛片的网站| 无限看片的www在线观看| 亚洲中文日韩欧美视频| 国产精品乱码一区二三区的特点| 可以免费在线观看a视频的电影网站| 非洲黑人性xxxx精品又粗又长| 桃色一区二区三区在线观看| 久久久久久亚洲精品国产蜜桃av| 欧美成人免费av一区二区三区| 91大片在线观看| 窝窝影院91人妻| 久久国产亚洲av麻豆专区| 久久久久久久久免费视频了| 夜夜躁狠狠躁天天躁| 国产精品 欧美亚洲| 亚洲男人的天堂狠狠| 久久人妻av系列| 波多野结衣巨乳人妻| 国产免费男女视频| 国产熟女xx| 免费在线观看视频国产中文字幕亚洲| 美国免费a级毛片| 成人国产综合亚洲| 在线观看www视频免费| 老熟妇仑乱视频hdxx| 久久久久亚洲av毛片大全| 亚洲第一av免费看| av欧美777| 亚洲激情在线av| 国产av不卡久久| 最近最新中文字幕大全免费视频| 久久精品91无色码中文字幕| 国产精品久久视频播放| 天天添夜夜摸| 又紧又爽又黄一区二区| 亚洲成人精品中文字幕电影| 亚洲国产欧洲综合997久久, | 777久久人妻少妇嫩草av网站| 国产精品一区二区免费欧美| 一级毛片高清免费大全| 成年版毛片免费区| 国产爱豆传媒在线观看 | 色综合亚洲欧美另类图片| 国产精品九九99| av视频在线观看入口| 91午夜精品亚洲一区二区三区| 国产免费一级a男人的天堂| 欧美日本亚洲视频在线播放| 国产成人91sexporn| 99九九线精品视频在线观看视频| 国产精品久久视频播放| 中文字幕精品亚洲无线码一区| 欧美成人免费av一区二区三区| 欧美潮喷喷水| 少妇高潮的动态图| 国内少妇人妻偷人精品xxx网站| 国产一区二区在线观看日韩| 久久精品影院6| 久久天躁狠狠躁夜夜2o2o| 少妇人妻精品综合一区二区 | 成人午夜高清在线视频| 成人性生交大片免费视频hd| 99久久成人亚洲精品观看| 国产成人一区二区在线| 亚洲国产高清在线一区二区三| 午夜福利18| 国产亚洲91精品色在线| 日日摸夜夜添夜夜爱| 51国产日韩欧美| 99久久成人亚洲精品观看| 欧美人与善性xxx| 人妻夜夜爽99麻豆av| 男女下面进入的视频免费午夜| 成人特级av手机在线观看| 亚洲精品国产av成人精品 | 日韩一本色道免费dvd| 菩萨蛮人人尽说江南好唐韦庄 | 女的被弄到高潮叫床怎么办| 看黄色毛片网站| 日本黄色视频三级网站网址| www日本黄色视频网| 国产精品一区二区三区四区久久| 日韩欧美免费精品| 老司机影院成人| 亚洲国产精品久久男人天堂| 日韩在线高清观看一区二区三区| 欧美日韩综合久久久久久| 深夜a级毛片| 国产黄a三级三级三级人| 嫩草影院入口| 国产精品一二三区在线看| 精品久久久噜噜| 女同久久另类99精品国产91| 美女大奶头视频| 成人特级av手机在线观看| 在线观看午夜福利视频| 国产视频内射| 精品国内亚洲2022精品成人| 国产精品久久久久久久电影| 精品人妻视频免费看| 国产精品女同一区二区软件| 最后的刺客免费高清国语| 一级av片app| a级一级毛片免费在线观看| 国产成人精品久久久久久| 精品国产三级普通话版| 激情 狠狠 欧美| 国产真实乱freesex| 欧美国产日韩亚洲一区| 成人永久免费在线观看视频| 日韩欧美 国产精品| 久久中文看片网| 亚洲七黄色美女视频| 国产视频内射| 国内精品久久久久精免费| 大香蕉久久网| 午夜福利在线在线| 亚洲国产欧洲综合997久久,| 在线观看午夜福利视频| 99久国产av精品| 97超碰精品成人国产| 久久久久久九九精品二区国产| 深夜精品福利| 精品久久久久久久末码| 国产91av在线免费观看| 看黄色毛片网站| 麻豆久久精品国产亚洲av| 久久精品国产清高在天天线| 日本a在线网址| 精品无人区乱码1区二区| 可以在线观看毛片的网站| 亚洲欧美精品自产自拍| 久久婷婷人人爽人人干人人爱| 欧美性感艳星| 能在线免费观看的黄片| 淫妇啪啪啪对白视频| 人人妻人人澡人人爽人人夜夜 | 一个人看的www免费观看视频| 神马国产精品三级电影在线观看| 日韩欧美三级三区| 久久6这里有精品| 麻豆精品久久久久久蜜桃| 99热精品在线国产| 国产熟女欧美一区二区| eeuss影院久久| 国产精品精品国产色婷婷| 麻豆成人午夜福利视频| 欧美潮喷喷水| 国产探花极品一区二区| 国产三级中文精品| 色哟哟·www| 亚洲国产精品成人久久小说 | 99热精品在线国产| 欧美日韩一区二区视频在线观看视频在线 | 丝袜美腿在线中文| 成年免费大片在线观看| 精品久久久久久久人妻蜜臀av| 日本-黄色视频高清免费观看| 精品一区二区免费观看| 国产高清有码在线观看视频| 欧美丝袜亚洲另类| 精品久久久噜噜| 亚洲中文字幕一区二区三区有码在线看| 男女做爰动态图高潮gif福利片| 搡老妇女老女人老熟妇| 国产69精品久久久久777片| 久久久欧美国产精品| 99国产极品粉嫩在线观看| 嫩草影院精品99| 国产黄色视频一区二区在线观看 | 亚洲国产欧洲综合997久久,| 最近手机中文字幕大全| 国产精品一区二区免费欧美| 日本-黄色视频高清免费观看| 久久精品国产自在天天线| 国产高清有码在线观看视频| 国产中年淑女户外野战色| 中文字幕av在线有码专区| 97碰自拍视频| 国产一区二区在线av高清观看| 麻豆成人午夜福利视频| 日韩大尺度精品在线看网址| 亚洲第一区二区三区不卡| 哪里可以看免费的av片| 国产精品嫩草影院av在线观看| 国产在线精品亚洲第一网站| 三级男女做爰猛烈吃奶摸视频| 国产伦精品一区二区三区视频9| 毛片一级片免费看久久久久| 日本在线视频免费播放| .国产精品久久| 久久久成人免费电影| 久久午夜福利片| 乱码一卡2卡4卡精品| 中文亚洲av片在线观看爽| 我的女老师完整版在线观看| 在线观看午夜福利视频| 日本一本二区三区精品| 精品熟女少妇av免费看| 91麻豆精品激情在线观看国产| 老司机影院成人| 亚洲三级黄色毛片| 久久久久久九九精品二区国产| 成人特级av手机在线观看| 久久人人爽人人片av| 波多野结衣巨乳人妻| 国产av在哪里看| 嫩草影院新地址| 国产91av在线免费观看| 全区人妻精品视频| 哪里可以看免费的av片| 性欧美人与动物交配| 免费看光身美女| 欧美日本亚洲视频在线播放| 搡老妇女老女人老熟妇| 亚洲人成网站在线播| 熟女电影av网| 国产男人的电影天堂91| 禁无遮挡网站| 欧美一级a爱片免费观看看| 国产熟女欧美一区二区| 少妇的逼好多水| 欧美色欧美亚洲另类二区| 欧美+日韩+精品| 一区福利在线观看| 一区二区三区免费毛片| 国产av一区在线观看免费| 小说图片视频综合网站| 国产精品,欧美在线| 一个人观看的视频www高清免费观看| 欧美高清性xxxxhd video| 在线天堂最新版资源| 欧美一区二区国产精品久久精品| 一级黄片播放器| 亚洲性久久影院| 日本 av在线| 男人的好看免费观看在线视频| 午夜精品一区二区三区免费看| 日本-黄色视频高清免费观看| 色综合站精品国产| 美女高潮的动态| 免费无遮挡裸体视频| 99久久中文字幕三级久久日本| 高清毛片免费观看视频网站| 精品一区二区三区av网在线观看| 国产女主播在线喷水免费视频网站 | 人妻制服诱惑在线中文字幕| 亚洲欧美日韩无卡精品| 日韩av在线大香蕉| 日韩欧美国产在线观看| 日本免费一区二区三区高清不卡| 免费看美女性在线毛片视频| 亚洲美女搞黄在线观看 | 美女内射精品一级片tv| 国产在线精品亚洲第一网站| 天天躁夜夜躁狠狠久久av| 我要搜黄色片| 校园春色视频在线观看| 寂寞人妻少妇视频99o| 欧美一区二区国产精品久久精品| 国产中年淑女户外野战色| 国产精品亚洲美女久久久| 男女下面进入的视频免费午夜| 成人无遮挡网站| 啦啦啦观看免费观看视频高清| 在线观看av片永久免费下载| 91在线精品国自产拍蜜月| 日韩欧美国产在线观看| 国产精品综合久久久久久久免费| 成年女人永久免费观看视频| 日本黄色片子视频| 欧美激情久久久久久爽电影| 在线免费观看的www视频| 日韩 亚洲 欧美在线| 天天一区二区日本电影三级| 久久中文看片网| 国产视频内射| .国产精品久久| 国产欧美日韩一区二区精品| aaaaa片日本免费| 欧美最新免费一区二区三区| 国产成年人精品一区二区| 夜夜爽天天搞| 国产高清三级在线| 又粗又爽又猛毛片免费看| 卡戴珊不雅视频在线播放| 成人亚洲精品av一区二区| 熟妇人妻久久中文字幕3abv| 日韩欧美一区二区三区在线观看| 久久午夜福利片| 亚洲av不卡在线观看| 欧美日韩一区二区视频在线观看视频在线 | 男女之事视频高清在线观看| 一本精品99久久精品77| 老熟妇仑乱视频hdxx| 亚洲欧美日韩卡通动漫| 又黄又爽又刺激的免费视频.| 久久欧美精品欧美久久欧美| 美女内射精品一级片tv| 成人欧美大片| 久久这里只有精品中国| 亚洲欧美日韩高清专用| 内射极品少妇av片p| 日本 av在线| 别揉我奶头~嗯~啊~动态视频| 欧美+日韩+精品| 亚洲精品色激情综合| 久久精品夜色国产| 国语自产精品视频在线第100页| 色综合站精品国产| 麻豆精品久久久久久蜜桃| 精品不卡国产一区二区三区| 亚洲内射少妇av| 亚洲精品456在线播放app| 欧美日韩综合久久久久久| 久久精品影院6| 国产v大片淫在线免费观看| 波多野结衣高清无吗| 性插视频无遮挡在线免费观看| 天堂av国产一区二区熟女人妻| 国产 一区精品| 日本欧美国产在线视频| 久久鲁丝午夜福利片| 日韩三级伦理在线观看| 精品人妻熟女av久视频| 亚洲av五月六月丁香网| 九色成人免费人妻av| 亚洲人成网站在线播| 久久久久久九九精品二区国产| 亚洲av熟女| 欧美一级a爱片免费观看看| 自拍偷自拍亚洲精品老妇| 成人无遮挡网站| 中文资源天堂在线| 丰满的人妻完整版| 桃色一区二区三区在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲国产欧洲综合997久久,| 一级av片app| 三级经典国产精品| 九九热线精品视视频播放| 久久精品国产鲁丝片午夜精品| 男人舔奶头视频| 好男人在线观看高清免费视频| 人妻制服诱惑在线中文字幕| 有码 亚洲区| 国产av一区在线观看免费| 亚洲精品在线观看二区| 99riav亚洲国产免费| 色5月婷婷丁香| 免费高清视频大片| 欧美丝袜亚洲另类| 亚洲经典国产精华液单| 听说在线观看完整版免费高清| 国产综合懂色| 欧美日本亚洲视频在线播放| 99久国产av精品国产电影| 欧美激情久久久久久爽电影| 久久久精品欧美日韩精品| 国产综合懂色| 97人妻精品一区二区三区麻豆| av女优亚洲男人天堂| 日韩,欧美,国产一区二区三区 | 真人做人爱边吃奶动态| 日韩制服骚丝袜av| 性欧美人与动物交配| 丰满的人妻完整版| 亚洲欧美清纯卡通| 老熟妇乱子伦视频在线观看| 狂野欧美激情性xxxx在线观看| 一本精品99久久精品77| www日本黄色视频网| 午夜激情福利司机影院| 亚洲国产高清在线一区二区三| 天天躁日日操中文字幕| 欧美绝顶高潮抽搐喷水| 久久久久久久午夜电影| 欧美一级a爱片免费观看看| 色综合亚洲欧美另类图片| 国产精品伦人一区二区| 国产精品久久久久久亚洲av鲁大| 一本久久中文字幕| 白带黄色成豆腐渣| 免费大片18禁| 午夜免费男女啪啪视频观看 | 久久久久精品国产欧美久久久| 国产 一区精品| 99热全是精品| 麻豆久久精品国产亚洲av| 18禁裸乳无遮挡免费网站照片| 免费无遮挡裸体视频| а√天堂www在线а√下载| 午夜视频国产福利| 国产精品人妻久久久影院| 亚洲天堂国产精品一区在线| 日本 av在线| 色av中文字幕| 欧美日韩在线观看h| 午夜激情欧美在线| ponron亚洲| 老司机影院成人| 夜夜看夜夜爽夜夜摸| 99国产精品一区二区蜜桃av| 欧美绝顶高潮抽搐喷水| 国产精品精品国产色婷婷| 小蜜桃在线观看免费完整版高清| 久久人人爽人人爽人人片va| 一本精品99久久精品77| 久久99热这里只有精品18| 亚洲成人av在线免费| 日本免费a在线| 久久久久免费精品人妻一区二区| 噜噜噜噜噜久久久久久91| 99久久精品一区二区三区| 日本精品一区二区三区蜜桃| 国产 一区 欧美 日韩| 免费av毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 最新中文字幕久久久久| 中出人妻视频一区二区| 亚洲欧美日韩无卡精品| aaaaa片日本免费| 久久欧美精品欧美久久欧美| 能在线免费观看的黄片| 国产爱豆传媒在线观看| 日韩亚洲欧美综合| 身体一侧抽搐| 男人和女人高潮做爰伦理| 久久精品91蜜桃| 国模一区二区三区四区视频| 亚洲成人av在线免费| 男女视频在线观看网站免费| 国产白丝娇喘喷水9色精品| 亚洲aⅴ乱码一区二区在线播放| 一个人看视频在线观看www免费| 国产一区二区三区在线臀色熟女| 天堂av国产一区二区熟女人妻| 99在线人妻在线中文字幕| 亚洲av成人精品一区久久| 国产精品一及| 男女那种视频在线观看| av免费在线看不卡| 在线播放国产精品三级| 久久精品国产亚洲av香蕉五月| 欧美极品一区二区三区四区| 在线观看一区二区三区| 欧美成人精品欧美一级黄| 搡老岳熟女国产| 精品午夜福利视频在线观看一区| 婷婷色综合大香蕉| 精品日产1卡2卡| 3wmmmm亚洲av在线观看| 91精品国产九色| 看非洲黑人一级黄片| 校园春色视频在线观看| 精品国产三级普通话版| 欧美国产日韩亚洲一区| 久久久成人免费电影| 乱码一卡2卡4卡精品| 免费不卡的大黄色大毛片视频在线观看 | 男女视频在线观看网站免费| 精品不卡国产一区二区三区| 禁无遮挡网站| 日韩欧美 国产精品| 国产人妻一区二区三区在| 久久欧美精品欧美久久欧美| 极品教师在线视频| 久久久久久伊人网av| 国产色婷婷99| 亚洲国产欧美人成| 国产亚洲精品综合一区在线观看| 久久综合国产亚洲精品| 国产伦一二天堂av在线观看| av在线播放精品| 黄色日韩在线| 99riav亚洲国产免费| av在线蜜桃| 欧美xxxx性猛交bbbb| av中文乱码字幕在线| 一级av片app| 少妇丰满av| 亚洲天堂国产精品一区在线| 国产精品人妻久久久影院| 可以在线观看毛片的网站| 欧美性猛交黑人性爽| 国产伦一二天堂av在线观看| 两个人视频免费观看高清| 三级经典国产精品| 日韩高清综合在线| 精品一区二区三区人妻视频| 亚洲人与动物交配视频| 日本与韩国留学比较| 一卡2卡三卡四卡精品乱码亚洲| 菩萨蛮人人尽说江南好唐韦庄 | 日本免费一区二区三区高清不卡| 亚洲人与动物交配视频| av在线观看视频网站免费| 国产视频内射| 黄片wwwwww| 中文在线观看免费www的网站| 国产成年人精品一区二区| 成人欧美大片| 免费大片18禁| 国产黄片美女视频| 美女cb高潮喷水在线观看| 欧美一区二区精品小视频在线| 少妇熟女欧美另类| 婷婷色综合大香蕉| 老女人水多毛片| 亚洲熟妇中文字幕五十中出| 成人欧美大片|