• <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.

    99re6热这里在线精品视频| 午夜久久久在线观看| 亚洲熟女精品中文字幕| 久久久久久久精品吃奶| 91精品三级在线观看| 免费黄频网站在线观看国产| 老司机靠b影院| 下体分泌物呈黄色| 亚洲 欧美一区二区三区| 91麻豆av在线| 女人被狂操c到高潮| 青草久久国产| 久久国产精品大桥未久av| 深夜精品福利| 欧美午夜高清在线| 久久国产精品男人的天堂亚洲| 黑人巨大精品欧美一区二区mp4| 黑人操中国人逼视频| 高潮久久久久久久久久久不卡| 免费久久久久久久精品成人欧美视频| 很黄的视频免费| 后天国语完整版免费观看| 国产亚洲精品久久久久5区| 国产三级黄色录像| 一进一出好大好爽视频| 亚洲国产看品久久| 丝瓜视频免费看黄片| 日韩欧美国产一区二区入口| 不卡一级毛片| 亚洲精品中文字幕在线视频| 精品国产亚洲在线| 一边摸一边抽搐一进一出视频| 搡老熟女国产l中国老女人| 亚洲人成伊人成综合网2020| 午夜福利免费观看在线| 亚洲熟妇熟女久久| 丰满的人妻完整版| 另类亚洲欧美激情| 精品视频人人做人人爽| tube8黄色片| 国产野战对白在线观看| 久久久久久人人人人人| 丰满饥渴人妻一区二区三| 建设人人有责人人尽责人人享有的| 麻豆成人av在线观看| 色老头精品视频在线观看| 女性被躁到高潮视频| 国产成人精品久久二区二区免费| 一级毛片高清免费大全| 免费在线观看日本一区| 91老司机精品| 高潮久久久久久久久久久不卡| 搡老乐熟女国产| 两个人看的免费小视频| 精品少妇久久久久久888优播| 国产日韩一区二区三区精品不卡| 正在播放国产对白刺激| 黄片播放在线免费| 人妻丰满熟妇av一区二区三区 | 两人在一起打扑克的视频| 国产精品美女特级片免费视频播放器 | 国产亚洲av高清不卡| 成年人午夜在线观看视频| 国产免费男女视频| 99国产精品一区二区三区| 国产亚洲精品一区二区www | 宅男免费午夜| 精品人妻在线不人妻| 午夜福利乱码中文字幕| 久久国产乱子伦精品免费另类| 婷婷精品国产亚洲av在线 | 国产成人影院久久av| 欧美日韩乱码在线| 俄罗斯特黄特色一大片| 日韩三级视频一区二区三区| 国产亚洲精品久久久久久毛片 | 国产一区二区三区综合在线观看| 午夜福利乱码中文字幕| 在线观看一区二区三区激情| 日韩精品免费视频一区二区三区| 国产主播在线观看一区二区| 欧美老熟妇乱子伦牲交| 欧美老熟妇乱子伦牲交| 久久久水蜜桃国产精品网| 五月开心婷婷网| 日韩欧美国产一区二区入口| 国产成人欧美| 国产欧美日韩一区二区三区在线| 欧美 日韩 精品 国产| 国产一区二区三区视频了| 国产无遮挡羞羞视频在线观看| 亚洲成人免费电影在线观看| 女同久久另类99精品国产91| 国产精品一区二区精品视频观看| 免费日韩欧美在线观看| 亚洲一区二区三区不卡视频| 高清欧美精品videossex| 国产成+人综合+亚洲专区| 手机成人av网站| 午夜激情av网站| 夜夜夜夜夜久久久久| 又大又爽又粗| 久久久久久久精品吃奶| 免费看十八禁软件| 亚洲五月婷婷丁香| 视频在线观看一区二区三区| 久久午夜亚洲精品久久| 国产av又大| 人人妻人人澡人人看| 亚洲免费av在线视频| 国产精品久久电影中文字幕 | 久久精品亚洲精品国产色婷小说| 美女国产高潮福利片在线看| 美女国产高潮福利片在线看| 精品少妇一区二区三区视频日本电影| 女同久久另类99精品国产91| 999精品在线视频| 亚洲美女黄片视频| 日韩精品免费视频一区二区三区| 国产av一区二区精品久久| 久久亚洲精品不卡| 国产不卡av网站在线观看| 精品乱码久久久久久99久播| 国产精品九九99| 99精国产麻豆久久婷婷| 欧美在线一区亚洲| 99riav亚洲国产免费| 免费观看a级毛片全部| 免费不卡黄色视频| 精品久久蜜臀av无| 欧洲精品卡2卡3卡4卡5卡区| 美女扒开内裤让男人捅视频| 黄色丝袜av网址大全| 黄色a级毛片大全视频| 91精品国产国语对白视频| 一本大道久久a久久精品| 国产精品一区二区精品视频观看| 欧美日韩黄片免| 亚洲精品在线美女| 日韩免费av在线播放| 精品亚洲成国产av| 国产一区在线观看成人免费| 国产精品二区激情视频| 色94色欧美一区二区| 69av精品久久久久久| 又大又爽又粗| 日韩有码中文字幕| 精品欧美一区二区三区在线| 日本撒尿小便嘘嘘汇集6| 最近最新免费中文字幕在线| 精品国产乱码久久久久久男人| 亚洲av日韩精品久久久久久密| 午夜亚洲福利在线播放| 最新的欧美精品一区二区| 一进一出抽搐动态| 亚洲综合色网址| 成人影院久久| 天天躁日日躁夜夜躁夜夜| 国产精品一区二区精品视频观看| 亚洲欧美激情在线| 亚洲国产欧美日韩在线播放| 美女午夜性视频免费| 午夜福利,免费看| 久久国产精品大桥未久av| 日韩人妻精品一区2区三区| 日韩一卡2卡3卡4卡2021年| 色尼玛亚洲综合影院| 国产激情欧美一区二区| 日韩免费高清中文字幕av| 精品国产亚洲在线| 国产一区有黄有色的免费视频| 在线免费观看的www视频| 成人亚洲精品一区在线观看| 满18在线观看网站| 一区二区三区精品91| 国产亚洲欧美98| 制服人妻中文乱码| 9191精品国产免费久久| 麻豆乱淫一区二区| 亚洲伊人色综图| 丰满的人妻完整版| 国产男女超爽视频在线观看| 人成视频在线观看免费观看| 国产国语露脸激情在线看| 欧美丝袜亚洲另类 | 久久久久精品人妻al黑| 中文字幕人妻丝袜制服| 久久性视频一级片| 女人被狂操c到高潮| 两个人看的免费小视频| 精品久久久精品久久久| 国产一区二区三区视频了| 91成人精品电影| 亚洲久久久国产精品| 高清黄色对白视频在线免费看| 欧美日韩视频精品一区| 老司机午夜福利在线观看视频| 脱女人内裤的视频| 激情在线观看视频在线高清 | 一夜夜www| 波多野结衣一区麻豆| 建设人人有责人人尽责人人享有的| 午夜91福利影院| 午夜两性在线视频| 国产男靠女视频免费网站| 在线免费观看的www视频| 午夜亚洲福利在线播放| 一区在线观看完整版| 欧美中文综合在线视频| 一级毛片精品| 老司机深夜福利视频在线观看| 黄网站色视频无遮挡免费观看| 久久国产乱子伦精品免费另类| 在线免费观看的www视频| 亚洲性夜色夜夜综合| 欧美精品一区二区免费开放| 成年女人毛片免费观看观看9 | 亚洲欧美日韩另类电影网站| av天堂在线播放| 老汉色av国产亚洲站长工具| 最近最新中文字幕大全免费视频| 亚洲人成伊人成综合网2020| 欧美久久黑人一区二区| 欧美乱色亚洲激情| 成人国产一区最新在线观看| 另类亚洲欧美激情| 亚洲av熟女| 国产成人免费观看mmmm| 国精品久久久久久国模美| 亚洲欧美激情在线| 欧美av亚洲av综合av国产av| 国产一区二区三区在线臀色熟女 | 欧美在线一区亚洲| 美女 人体艺术 gogo| 国产亚洲欧美精品永久| 亚洲av成人一区二区三| 欧美亚洲日本最大视频资源| 国产免费av片在线观看野外av| 精品国产乱码久久久久久男人| 亚洲人成电影免费在线| 久久影院123| 美女国产高潮福利片在线看| 亚洲少妇的诱惑av| 人人妻人人爽人人添夜夜欢视频| 国产精品久久久久久人妻精品电影| 国产区一区二久久| 免费观看精品视频网站| 女同久久另类99精品国产91| 新久久久久国产一级毛片| 国产一区二区激情短视频| 亚洲精华国产精华精| 国产单亲对白刺激| 女人爽到高潮嗷嗷叫在线视频| 婷婷成人精品国产| 国产不卡av网站在线观看| 狠狠狠狠99中文字幕| 日本vs欧美在线观看视频| xxxhd国产人妻xxx| 欧美成人午夜精品| 一级a爱片免费观看的视频| 久久中文字幕人妻熟女| 后天国语完整版免费观看| 91老司机精品| 女同久久另类99精品国产91| 天天影视国产精品| 免费高清在线观看日韩| 国产精华一区二区三区| 国产精品欧美亚洲77777| 国产在视频线精品| 中文字幕制服av| 黄色视频不卡| 午夜精品在线福利| 国产片内射在线| 国产精品99久久99久久久不卡| 国产在线观看jvid| 十八禁网站免费在线| 一区二区三区精品91| 久久久久精品国产欧美久久久| 9热在线视频观看99| 成年动漫av网址| 午夜福利在线观看吧| 亚洲国产中文字幕在线视频| 国产亚洲精品一区二区www | a在线观看视频网站| 国产精品免费大片| 女性被躁到高潮视频| 搡老岳熟女国产| 亚洲av第一区精品v没综合| 欧美丝袜亚洲另类 | 精品久久久久久,| 91精品三级在线观看| 欧美日韩成人在线一区二区| 国产男女内射视频| 国产在视频线精品| 国产男靠女视频免费网站| 欧美激情极品国产一区二区三区| 久久精品91无色码中文字幕| 亚洲精品国产一区二区精华液| 又黄又爽又免费观看的视频| 久久国产精品男人的天堂亚洲| 精品免费久久久久久久清纯 | 国产亚洲一区二区精品| 91麻豆精品激情在线观看国产 | 男女之事视频高清在线观看| 中文欧美无线码| 人人妻人人澡人人看| 又大又爽又粗| 叶爱在线成人免费视频播放| 黄片小视频在线播放| 每晚都被弄得嗷嗷叫到高潮| 亚洲片人在线观看| 久久香蕉精品热| 老汉色∧v一级毛片| 日韩 欧美 亚洲 中文字幕| 狠狠婷婷综合久久久久久88av| 久久久国产成人精品二区 | 91九色精品人成在线观看| 亚洲少妇的诱惑av| 日韩视频一区二区在线观看| 国产97色在线日韩免费| 不卡av一区二区三区| 亚洲国产毛片av蜜桃av| 99国产精品一区二区蜜桃av | 巨乳人妻的诱惑在线观看| 满18在线观看网站| 多毛熟女@视频| 欧美日韩av久久| 少妇被粗大的猛进出69影院| 亚洲五月天丁香| 亚洲欧洲精品一区二区精品久久久| 在线观看免费视频日本深夜| 亚洲,欧美精品.| tocl精华| 91精品国产国语对白视频| 国产深夜福利视频在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲人成伊人成综合网2020| 中文字幕高清在线视频| 国产在线一区二区三区精| 国产免费现黄频在线看| 亚洲一区高清亚洲精品| 色老头精品视频在线观看| 久久亚洲精品不卡| 久久久久久久久久久久大奶| 国产精品 欧美亚洲| 女同久久另类99精品国产91| 午夜免费观看网址| 成人手机av| xxxhd国产人妻xxx| 成人国产一区最新在线观看| 日韩欧美免费精品| 久久精品国产99精品国产亚洲性色 | 性少妇av在线| 亚洲成人手机| 大型黄色视频在线免费观看| 久久精品人人爽人人爽视色| 国产97色在线日韩免费| 欧美激情久久久久久爽电影 | 黄色a级毛片大全视频| 纯流量卡能插随身wifi吗| 亚洲av日韩精品久久久久久密| e午夜精品久久久久久久| 亚洲午夜理论影院| 99re在线观看精品视频| 亚洲国产欧美一区二区综合| √禁漫天堂资源中文www| 侵犯人妻中文字幕一二三四区| 人妻 亚洲 视频| 在线观看免费日韩欧美大片| 正在播放国产对白刺激| 午夜日韩欧美国产| 亚洲少妇的诱惑av| 女人久久www免费人成看片| 这个男人来自地球电影免费观看| av一本久久久久| 飞空精品影院首页| 国产男女超爽视频在线观看| 久久久久国产一级毛片高清牌| 麻豆成人av在线观看| 人人妻人人爽人人添夜夜欢视频| 一级作爱视频免费观看| 看片在线看免费视频| 天天影视国产精品| 久久久久久久久久久久大奶| 成人18禁高潮啪啪吃奶动态图| 丰满饥渴人妻一区二区三| 女人久久www免费人成看片| 久久精品国产99精品国产亚洲性色 | 亚洲成人手机| 国产成人系列免费观看| 精品福利永久在线观看| 亚洲情色 制服丝袜| 亚洲专区中文字幕在线| 男女之事视频高清在线观看| 久久国产精品人妻蜜桃| 日韩欧美国产一区二区入口| 男女之事视频高清在线观看| 亚洲精品国产色婷婷电影| 国产高清视频在线播放一区| 久久人妻熟女aⅴ| 最近最新中文字幕大全电影3 | 老司机深夜福利视频在线观看| 国产又色又爽无遮挡免费看| 99riav亚洲国产免费| 亚洲av片天天在线观看| 亚洲精品久久成人aⅴ小说| 国产精品成人在线| netflix在线观看网站| 在线永久观看黄色视频| av福利片在线| 欧美 亚洲 国产 日韩一| 精品人妻在线不人妻| 欧美日韩成人在线一区二区| 欧美乱码精品一区二区三区| av欧美777| 免费观看a级毛片全部| 免费在线观看视频国产中文字幕亚洲| 国产一卡二卡三卡精品| x7x7x7水蜜桃| 一二三四在线观看免费中文在| 免费高清在线观看日韩| 国产高清国产精品国产三级| 99国产极品粉嫩在线观看| 1024视频免费在线观看| 亚洲 欧美一区二区三区| 99re在线观看精品视频| 久99久视频精品免费| 欧美精品亚洲一区二区| 国产三级黄色录像| 一级片'在线观看视频| ponron亚洲| 99国产精品一区二区蜜桃av | 国产精品香港三级国产av潘金莲| 欧美日韩视频精品一区| 99国产精品一区二区三区| 黄色视频不卡| 久久久国产一区二区| 国产欧美日韩一区二区三区在线| 一本大道久久a久久精品| 国产不卡一卡二| 成人黄色视频免费在线看| 亚洲自偷自拍图片 自拍| 黄片小视频在线播放| 久久精品国产99精品国产亚洲性色 | 不卡一级毛片| 亚洲av日韩精品久久久久久密| 亚洲精品国产一区二区精华液| 丝袜在线中文字幕| 精品国产一区二区三区久久久樱花| 制服人妻中文乱码| 亚洲av成人不卡在线观看播放网| 中亚洲国语对白在线视频| 国产1区2区3区精品| 亚洲欧美色中文字幕在线| 精品高清国产在线一区| 老熟妇乱子伦视频在线观看| 男女床上黄色一级片免费看| 亚洲综合色网址| 18禁观看日本| 一区二区三区精品91| 国产免费男女视频| 久久香蕉国产精品| 多毛熟女@视频| 免费黄频网站在线观看国产| 久99久视频精品免费| 精品少妇一区二区三区视频日本电影| 少妇裸体淫交视频免费看高清 | 男女午夜视频在线观看| 亚洲五月婷婷丁香| 精品福利观看| 丝袜美腿诱惑在线| 亚洲欧美日韩高清在线视频| 日韩精品免费视频一区二区三区| 精品无人区乱码1区二区| 亚洲五月色婷婷综合| 成在线人永久免费视频| 国产日韩一区二区三区精品不卡| 亚洲美女黄片视频| 欧美人与性动交α欧美精品济南到| 女性被躁到高潮视频| videosex国产| 久久国产精品人妻蜜桃| 村上凉子中文字幕在线| 后天国语完整版免费观看| 国产精品国产高清国产av | 制服人妻中文乱码| 国产一区在线观看成人免费| 日本一区二区免费在线视频| 日韩中文字幕欧美一区二区| 久久精品国产亚洲av高清一级| 亚洲精品av麻豆狂野| 一进一出好大好爽视频| 亚洲第一欧美日韩一区二区三区| 日韩制服丝袜自拍偷拍| 香蕉国产在线看| 激情在线观看视频在线高清 | av国产精品久久久久影院| 久热这里只有精品99| 久久久久久亚洲精品国产蜜桃av| 久久久久久人人人人人| 18在线观看网站| 在线观看66精品国产| 亚洲一码二码三码区别大吗| 色综合婷婷激情| 国产一区二区激情短视频| 久久国产精品人妻蜜桃| 国产精品香港三级国产av潘金莲| 亚洲一区高清亚洲精品| 丝袜在线中文字幕| 久久久久国内视频| 激情在线观看视频在线高清 | 啦啦啦在线免费观看视频4| 久久久精品免费免费高清| 午夜老司机福利片| 欧美亚洲日本最大视频资源| 久久精品亚洲av国产电影网| 国产精品98久久久久久宅男小说| 一夜夜www| 亚洲一区二区三区不卡视频| 女人久久www免费人成看片| 91老司机精品| 午夜免费观看网址| 亚洲avbb在线观看| 精品亚洲成a人片在线观看| 久9热在线精品视频| 动漫黄色视频在线观看| 亚洲五月天丁香| 亚洲欧美精品综合一区二区三区| 亚洲国产精品一区二区三区在线| 黑人操中国人逼视频| 亚洲av日韩精品久久久久久密| 精品少妇久久久久久888优播| 精品国产美女av久久久久小说| 国产精品免费大片| av欧美777| 免费久久久久久久精品成人欧美视频| 啦啦啦免费观看视频1| 夜夜夜夜夜久久久久| 日本五十路高清| 日本wwww免费看| 好男人电影高清在线观看| 别揉我奶头~嗯~啊~动态视频| 悠悠久久av| 香蕉国产在线看| 妹子高潮喷水视频| 久久久久国内视频| 侵犯人妻中文字幕一二三四区| 一级片免费观看大全| 色老头精品视频在线观看| 久久精品国产亚洲av高清一级| 亚洲自偷自拍图片 自拍| 国产亚洲精品久久久久5区| 天天躁日日躁夜夜躁夜夜| 99国产综合亚洲精品| 免费在线观看视频国产中文字幕亚洲| 狠狠婷婷综合久久久久久88av| 欧美人与性动交α欧美软件| 亚洲性夜色夜夜综合| 美国免费a级毛片| 中文字幕人妻丝袜制服| 99国产精品一区二区三区| 亚洲国产精品一区二区三区在线| 色婷婷久久久亚洲欧美| cao死你这个sao货| 国产成+人综合+亚洲专区| 欧美日韩乱码在线| av不卡在线播放| 亚洲国产精品一区二区三区在线| 欧美激情久久久久久爽电影 | 亚洲欧美日韩另类电影网站| 中亚洲国语对白在线视频| 侵犯人妻中文字幕一二三四区| 国产精品久久视频播放| 一级a爱视频在线免费观看| 日韩人妻精品一区2区三区| 在线观看免费视频日本深夜| 国产成人影院久久av| 在线观看免费视频日本深夜| 国产精华一区二区三区| 欧美亚洲日本最大视频资源| 波多野结衣一区麻豆| 国产欧美日韩一区二区三区在线| 美女 人体艺术 gogo| 国产精品影院久久| 18禁裸乳无遮挡动漫免费视频| 麻豆乱淫一区二区| 大型av网站在线播放| 在线观看www视频免费| 精品无人区乱码1区二区| 亚洲美女黄片视频| 精品欧美一区二区三区在线| avwww免费| 成人影院久久| 亚洲精品粉嫩美女一区| 色综合婷婷激情| 午夜亚洲福利在线播放| 欧美性长视频在线观看| 夜夜爽天天搞| 天堂俺去俺来也www色官网| 69精品国产乱码久久久| 欧美 日韩 精品 国产| 亚洲人成电影免费在线| 久久精品国产a三级三级三级| 操美女的视频在线观看| 国产一区有黄有色的免费视频| 美女扒开内裤让男人捅视频| 国产91精品成人一区二区三区| 亚洲精品久久午夜乱码| 成人免费观看视频高清| 国产成+人综合+亚洲专区| 亚洲av日韩在线播放| 激情在线观看视频在线高清 | 欧美人与性动交α欧美软件| 国产精品av久久久久免费| 欧美日韩瑟瑟在线播放|