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

    Assessment of recoverable oil and gas resources by in-situ conversion of shale-Case study of extracting the Chang 73 shale in the Ordos Basin

    2022-06-02 04:59:32LinHuHouXiLuoSenHuLinYongXinLiLiJunZhngWeiJio
    Petroleum Science 2022年2期

    Lin-Hu Hou ,Xi Luo ,Sen-Hu Lin ,Yong-Xin Li ,**,Li-Jun Zhng ,Wei-Jio M

    a Research Institute of Petroleum Exploration &Development,PetroChina,Beijing,100083,China

    b National Energy Tight Oil &Gas R&D Center,Beijing,100083,China

    c School of Geosciences,China University of Petroleum,Qingdao,Shandong,266580,China

    Keywords:Shale oil In-situ conversion processing Thermal simulation experiment Hydrocarbon generation mechanism Hydrocarbon resource Ordos basin

    ABSTRACT The purpose of this study is to investigate the entire evolution process of shales with various total organic contents (TOC) in order to build models for quantitative evaluation of oil and gas yields and establish methods for assessing recoverable oil and gas resources from in-situ conversion of organic matters in shale.Thermal simulation experiments under in-situ conversion conditions were conducted on Chang 73 shales from the Ordos Basin in a semi-open system with large capacity.The results showed that TOC and Ro were the key factors affecting the in-situ transformation potential of shale.The remaining oil and gas yields increased linearly with TOC but inconsistently with Ro.Ro ranged 0.75%-1.25%and 1.05%-2.3%,respectively,corresponding to the main oil generation stage and gas generation stage of shale in-situ transformation.Thus a model to evaluate the remaining oil/gas yield with TOC and Ro was obtained.The TOC of shale suitable for in-situ conversion should be greater than 6%,whereas its Ro should be less than 1.0%.Shales with 0.75% (Ro) could obtain the best economic benefit.The results provided a theoretical basis and evaluation methodology for predicting the hydrocarbon resources from in-situ conversion of shale and for the identification of the optimum“sweet spots”.The assessment of the Chang 73 shale in the Ordos Basin indicated that the recoverable oil and gas resources from in-situ conversion of organic matters in shale are substantial,with oil and gas resources reaching approximately 450×108 t and 30×1012 m3,respectively,from an area of 4.27×104 km2.

    1.Introduction

    With the future increase in the global population and economy,the energy demand will grow substantially in the next few decades(Ma et al.,2020a).The global oil consumption has already risen from 41.43×108t in 2008 to 46.62×108t in 2018 (BP,2019).Energy security has drawn the attention of many countries to their oil and gas supply.The successful development of shale oil and gas in the US has reshaped the energy landscape in the world (Jia et al.,2012;Zendehboudi and Bahadori et al.,2017;Yang and Zou et al.,2019;Zhao et al.,2020).In 2019,the crude oil production reached 44.65 billion barrels in the US,of which,29.14 billion barrels were from shale oil,accounting for 65.12%of the total oil production that year(EIA,2020-01).Shale oil is playing an increasingly important role in the dramatic growth of American oil production(EIA,2018,2019a,2019b;Zhou et al.,2019;Zou et al.,2019) and exerts a significant influence on the global energy and geopolitical layout.

    The technologies of horizontal wells and multi-stage fracturing have driven the rapid growth of oil and gas production from shales with medium-high maturity in the US(IHS,2016;Jin et al.,2019a;Perrin,2019-06-06).Consequently,China has attached great importance to shale oil and gas(Zou et al.,2010;Lu et al.,2012;Xue et al.,2015;Du et al.,2019;Ju et al.,2020;Kuang et al.,2020;Song et al.,2020).Surveys indicate that two types of shale,i.e.,lacustrine shale and marine shale,are developed in China.Marine shales are mainly distributed in South China with relatively high maturity,in which the vitrinite reflectance (Ro) is generally greater than 2.0%(Zou et al.,2011,2015;Zhao et al.,2012,2016;Nie et al.,2016).These areas contain abundant shale gas resources and some have been developed successfully with a production of 150×108m3in 2019.However,lacustrine shales are primarily distributed in North China,including the Ordos,Songliao,Bohai Bay and Junggar Basins,covering an area of more than 20×104km2(Dang et al.,2015;Zhao et al.,2018a,b;Yang and Jin,2019;Zou et al.,2020) and characterized by relatively low maturity for which the Rois normally less than 1.0%.China has been attempting to exploit the shale oil in its lacustrine formations for more than a decade by drilling nearly 500 horizontal wells.However,the economic development of shale oil with low and medium maturity has not yet been realized,and the undeveloped organic pores,poor connectivity,low gas-oil ratio(GOR) and high viscosity of oil are the main reasons,all of which lead to an extremely low recovery factor of only 1%-3%(Saif et al.,2017;Jin et al.,2019b;Hou et al.,2020a) by simply relying on volume fracturing in horizontal wells.Prospecting history in the US also shows that conventional volume fracturing in horizontal wells is unable to realize the economic development of shale oil with low to medium maturity.Instead,in-situ conversion must be used to achieve economic development.Shell,ExxonMobil,Total,and many other giant oil companies have conducted relevant research and established four types of in-situ heating technologies,including electric heating,convection heating,radiant heating,and combustion heating(Wang et al.,2013).

    In-situ conversion of shale oil,including physical and chemical process,enables the organic matter(including solid organic matter and asphalt)in shale or oil shale to be rapidly converted into light oil and natural gas by in-situ heating.It is then exploited,while coke and other residues are left in the subsurface.This process can be regarded as an“underground refinery”(Zhao et al.,2018a,b).Insitu transformation of shale has four technical advantages:(1) the oil obtained through the in-situ heating and conversion is of good quality(35-49°API);(2)the overpressure fluid and micro-fracture network generated by the in-situ conversion process increase the driving force,the permeable channels,and the drainage system,resulting in an ultimate recovery surpassing 60%;(3) there are no tailings or waste produced,with little groundwater pollution and no hydraulic fracturing needed,which attains actual clean mining by minimizing the damage to the ecosystem caused by harmful byproducts during conventional production;and (4) the occupied area is very small,which helps to protect farmland.The in-situ conversion technologies can successfully tackle the challenges in exploiting the shale with low-medium maturity,and lead to efficient and environmentally friendly utilization (hereafter we use shale to represent both the shale and oil shale).Nonetheless,some key issues must be addressed to achieve commercial development of the shale oil with low-medium maturity,including but not limited to criteria for the selection and evaluation of the recovery targets.To address these issues,we took the shale of the 3rd section in the Chang 7 member of the Tertiary in the Ordos Basin(hereafter referred to as the Chang 73shale) as an example to simulate the hydrocarbon generation and expulsion process under in-situ conditions with a semi-open system.The purpose was to establish an approach with which to evaluate the recoverable oil and gas resources from in-situ conversion of shale,to provide a theoretical basis and evaluation methodology for predicting the hydrocarbon resources from in-situ conversion of shale,and to identify the optimum “sweet spots.”

    2.Samples and experimental conditions

    2.1.Samples

    The Chang 73shale developed in the Late Triassic sediment in semi-deep to deep lacustrine facies.It is characterized by continuous sedimentation with large thickness,high organic matter abundance(average total organic content(TOC)is 14.32 wt%),and a good organic matter type (TypeII kerogen (Ma et al.,2020b;Hou et al.,2020b)).A series of hydrocarbon generation and expulsion thermal simulations were conducted by analyzing various Chang 73shale outcrop samples with different TOC values.Altogether,nine samples were acquired from 5 m below the surface (to avoid the effects of epigenesis,such as weathering and leaching) in the southeast of the Ordos Basin(Fig.1).The Roof the samples was less than 0.5%.The samples were crushed into 40-60 mesh size,mixed,and divided evenly into several parts.Then,some samples were randomly chosen for TOC,Rock-Eval,and vitrinite reflectance (Ro)tests to obtain fundamental geochemical information(see Table 1).

    Table 1Basic geochemical data of the original unheated shale samples.

    Table 2Characteristics of the residual samples after pyrolysis experiments.

    2.2.Experimental device and process

    Pyrolysis experiment is one of the most direct and efficient ways to investigate the oil and gas generation process.Previous researchers have conducted a large number of simulation studies by changing experimental conditions,including the sample feed(from 10 mg (You et al.,2019) to 200 g (Tang et al.,2015)),heating rate(from 0.5 to 25°C/s to 12°C/h (Doan et al.,2013;Lan et al.,2015;Siramard et al.,2017),pyrolysis temperature (max temperature from 420°C (Siramard et al.,2017;Ma et al.,2020c) to 1000°C(Shao et al.,2019)),and pressure (0.06-68 MPa (Lan et al.,2015;Sun et al.,2019;Yan et al.,2019)).Based on these studies,kinetic models have been built of hydrocarbon generation of different types of organic matter (Han et al.,2014),impacts of temperature and pressure on hydrocarbon generation have been discussed,and it was concluded that a rapid heating rate reduces oil generation,increase gas generation(Siramard et al.,2017),and exerts very little influence on the kinetic parameters (Bai et al.,2015).Thermal simulation experiments deepen our understanding of the hydrocarbon generation process.However,they also have disadvantages like other technologies.For instance,the error in assessing oil and gas generation can increase if the sample feed is small,and the heating rate is rapid.Additionally,there has not been much discussion regarding the oil and gas expulsion process constrained by experimental instruments.In this paper,thermal simulations with large sample feed and a slow heating rate were designed in a semiopen system using newly developed reactors with a large capacity.A total of 99 experiments were completed with nine samples at 11 temperatures to study the hydrocarbon generation features of the low maturity shale with various TOC values.

    The experimental devices are diagramed in Fig.2.The reactor was made of a special alloy,resistant to corrosion by H2S,CO2,and H2.A special gasket and coating were used between the reactor body and cap to ensure long-term operation without leakage.The highest temperature that the reactor could bear was 700°C,the maximum pressure was 40 MPa,the container capacity was 1.36 L,and samples weighing approximately 2 kg could be accommodated.The instrument meets the experimental requirements of thermal simulation of large-volume shale under different temperatures and pressures.The reactor has two heating systems,one at the bottom and the other on the side.They can operate simultaneously.Additionally,multiple thermocouples and pressure sensors are equipped to precisely control the temperature and pressure inside the reactor,such that the temperature error is below 1°C and the fluid pressure error is less than 0.1 MPa during experiments.

    Fig.1.Shale distribution in the Ordos Basin and the sampling point.

    Fig.2.Schematic diagram of the pyrolysis device.(1) autoclave,(2) sealing cover,(3) outer cover,(4) electric heater,(5) thermo couple,(6) pressure meter,(7) piston valve controlled by three-way solenoid valve,(8) needle valve,(9) heating tape,(10) synchronous collector,(11) floating piston,(12) liquid collector,(13) condenser,(14) gas bag,(15)measuring cylinder.

    The experimental steps are discussed below.

    (1) Sample loading.The crushed and evenly mixed shale samples were loaded to the reactor.

    (2) Heating simulation.①Leak test.After sealing the reactor,vacuumed it,filled it with 20 MPa helium,and placed it for leak test.When there was no leakage,released the gas.Repeated the process for 3-4 times,and vacuumized the reactor for the last time.②Heating.5 MPa and 7 MPa were set as the fluid pressure and hydrocarbon expulsion pressure respectively.A total of 11 thermal simulation temperature points were designed (Table 2),covering the whole process of oil and gas generation.The heating procedure was as follows:the sample was heated from room temperature to the temperature before each set point at the heating rate of 20°C/d (the temperature before the first temperature point is 200°C),and then the sample was heated to the set temperature at the heating rate of 5°C/d,and kept at this temperature for 10 h.

    (3) Simulation,product collection and quantification.①hydrocarbon expulsion.The process of hydrocarbon expulsion was controlled by the three-way solenoid valve.At the beginning,the piston valve was closed and the whole reactor was sealed.When the generated products accumulated continuously in the container and the fluid pressure exceeded the set value,the piston valve would open automatically,and the products were released from the top of the reactor until the pressure dropped to the preset value.After that,the piston valve would close automatically and the reactor would be closed again.②Collection and quantification of pyrolysis products.When products were discharged,their release rate was controlled by the needle valve,so that the pressure in the container was not reduced rapidly,and the synchronous collector could also play a buffer role,so as to ensure that the oil and gas could be separated in time when they were generated and discharged in large quantities.The hydrocarbon expulsion pipeline in the apparatus was wrapped with a heating belt,and the temperature of the heating belt was synchronized with the temperature in the reactor to prevent the condensation of products from blocking the pipeline.When the product passed through the condenser (water circulation at 20°C),water and hydrocarbon liquid(including C5) were collected in the liquid collector.Hydrocarbon gas(C1-C4) and non-hydrocarbon gas were collected in the air bag immersed in the water tank.The volume of the gas bag was determined by the volume of water discharged from the expansion of the air bag.Comparing the mass or volume of produced oil or gas before thermal simulation,the oil and gas yields of shale sample could be obtained.③Collection and quantification of retained oil and solid residual sample.After the autoclave cooled down to room temperature,the simulated source rock sample was taken out,weighed and extracted with dichloromethane,which was regarded as retained oil.The gas obtained by purging was the retained gas(the retained oil and gas were not discussed here,as they were not involved in our product modelling).The residual samples were tested for TOC and Roafter dichloromethane extraction (Table 2).

    (4) Reliability analysis of simulation results.Theoretically,the mass of the reactants before and after the thermal simulation was conserved.Nevertheless,measurement errors and equipment leakage may cause product loss and result in disequilibrium.The mass balance rate,the ratio of the product mass after simulation(including oil,gas,water,and residuals)to that of the shale sample before simulation,was used to evaluate the reliability of the experimental results(Equation (1)).

    where mass balance is the mass balance rate,residual rock is the mass of the residual sample after simulation,oil is the mass of the produced and remaining oil after simulation,water is the mass of the produced and remaining water after simulation,gas is the mass of the produced and remaining gas after simulation,and unheated rock is the mass of the sample in the reactor before simulation.

    The mass balance rate should be no less than 99.5% if there is only a measurement error.Otherwise,leakage occurred during the simulation and the results are not reliable.Under such circumstances,the thermal simulation must be conducted again.In this study,the mass balance rates of all 99 simulations exceeded 99.75%(Fig.3),indicating the results of the pyrolysis experiments were reliable.

    3.Experimental results

    3.1.Relationship between pyrolysis and Ro

    Thermal maturity is an important parameter with which hydrocarbon generation and expulsion (Tissot and Welte,1984) can be assessed.Of the various evaluation indexes,Rois the only indicator put forward as an international measurement standard(Stach et al.,1982) and is widely used.The Roof the residual samples at each simulation temperature was analyzed,and the results indicated that there were very highly positive correlations between Roand temperature was analyzed,and the results indicated that there were very highly positive correlations between Roand the thermal simulation temperature (Fig.4).Furthermore,the measured Rovalues of different samples at the same temperature were very close.Therefore,the average Roof the nine samples at each temperature was used to establish the model to express the relationship between thermal simulation temperature and Ro(Equation(2)).

    where Rois vitrinite reflectance,T is thermal simulation temperature,and a1and a2are empirical coefficients,equal to 0.1380 and 0.0057,respectively.

    3.2.Characteristics of hydrocarbon yield from in-situ conversion of shale

    Fig.5 shows the yields of oil and gas with increasing Ro.Generally,the oil and gas yields of shale increased linearly with TOC,that was,the higher the abundance of organic matter(TOC)in source rock,the higher the yields.For example,the oil and gas yields of the sample with TOC of 25.988 wt% are 66.85 mg/g·rock and 33.30 mL/g·rock at Ro=3.65%.However,the variation trend of oil and gas yields were different from that of Ro.For oil yield,it increased rapidly at the beginning of thermal evolution and reached the peak of oil generation at Ro?1.2%,which was similar with that of Lucaogou shale (Hou et al.,2021),another typical lacustrine type-II shale in China.After that,it almost remained the same.For gas yield,the value first increased rapidly with Roand then increased slowly,showing that it was affected by the cracking of crude oil and wet gas during high evolution stage.However,for samples with TOC ≤13.34 wt%,the gas production increased rapidly before Ro?1.67%,and then increased very slowly.While for the two samples with TOC>20 wt%,the trend of increasing gas production was very obvious before Ro≈2.4%,and then it showed a gentle increasing trend.

    4.Discussions

    4.1.Characteristics of the remaining hydrocarbon yield

    Fig.3.Evolution of mass fractions of major products with increasing thermal maturity.

    Fig.4.Correlation between pyrolysis temperature and Ro.

    Fig.5.Yields of oil and gas with increasing Ro.(a) yield of oil,(b) yield of gas.

    Fig.6.Yield of the remaining hydrocarbons with increasing Ro.(a) remaining yield of oil,(b) remaining yield of gas.

    Fig.7.Yield of remaining hydrocarbons with original TOC.(a) remaining yield of oil,(b) remaining yield of gas.

    Fig.8.Yield of remaining hydrocarbons with residual TOC.(a) remaining yield of oil,(b) remaining yield of gas.

    The hydrocarbons generated but retained in shale,together with the organic matters that were not yet degraded constituted the major source of oil and gas during the in-situ conversion of shale.Subsurface shales have generally undergone thermal evolution.Some oil and gas may have been generated and discharged already.Therefore,in this paper,the hydrocarbon yield from in-situ conversion of shale refers to the remaining hydrocarbon yield.With this in mind,we have estimated the hydrocarbon yield from in-situ conversion of shale by calculating the obtained oil and gas at each experimental temperature (yields of the remaining oil and gas are plotted in Fig.6).It shows that the remaining oil and gas yields exhibit consistent trends as Roincreases and both were negatively correlated with Ro,whereas regarding specific performance,these data reduced slowly at first and then decreased dramatically when Roreached a certain value.However,oil and gas have different evolution stages.Regarding the remaining oil yield,it decreased slowly when Roranged from 0.58%to 0.75%,then it dropped almost linearly when Rowas between 0.75% and 1.25% and decreased slowly again until it reached zero when Rowas greater than 1.25%.In contrast,the remaining gas yield declined slowly when Rowas from 0.58% to 1.05% and diminished rapidly in a near power function pattern when Rowas between 1.05%and 2.3%.Additionally,the major oil and gas generation stages were different,and it can be concluded from the yield rates that the stage between 0.75% and 1.25%was the major oil generation period,which accounted for 75%of the total oil yield.However,the amount of oil generated during the high evolution stage(Ro1.7%-2.3%)was small and approached to zero when Roexceeded 2.3%.On the contrary,the major gas generation period was 1.05%-2.3%.The remaining gas yield declined slowly to a constant level when Ro>2.3%.At this stage,the gas was mainly generated from thermal expansion because of increasing temperature.

    TOC was another key factor influencing the potential of shale transformation,as well as Ro.The organic matter provided a material basis for hydrocarbons,and the original TOC in shale determined its ultimate yield of oil and gas.Fig.7 shows that both the remaining oil yield and gas yield increased with the original TOC for the samples with the same Ro,and there were strong linear correlations between the yield and the original TOC (Fig.7).Furthermore,the remaining oil yield and gas yield decreased as Roincreased under the same original TOC,which was consistent with previous findings.

    When we evaluated the recoverable petroleum resources from in-situ conversion of shale under geologic conditions,we had to understand that the evaluated objects had generally experienced thermal evolution.Because of this,the organic matter in shale had already been transformed making it difficult to obtain the original TOC value.Therefore,the relationship between the remaining oil/gas yield and the residual TOC was a focus in this paper(Fig.8).The remaining oil yield from shale samples with different TOC values increased as the residual TOC increased when Rowas below 1.7%and changed consistently with TOC when Rowas below 1.1%(Fig.8a).The remaining gas yield increased as the residual TOC increased when Roremained the same.However,when TOC was unchanged,the remaining gas yield increased at first and then decreased as Roincreased,and the turning point was at Ro=1.1%(Fig.8b).This feature could be attributed to the evolution process because first,the organic matter generated oil and gas as the temperature rised,then some oil was cracked into coke at high maturity,and these actions caused TOC to decline rapidly at first,then increased slightly,and ultimately remained stable.

    Based on the correlations between the remaining oil/gas yield(per mass unit of shale and mass unit of TOC) versus TOC and Ro(Fig.9a),it could be concluded that,as for the same shale sample,the remaining oil yield increased slightly at first,then decreased steadily and finally diminished rapidly as Roincreased.However,for the remaining gas yield,Roincreased (Fig.9b).

    4.2.Models for evaluating the remaining oil/gas yield from in-situ shale transformation

    As previously described,the experimental results demonstrated that there was a strong linear correlation between the remaining oil/gas yield and TOC.Accordingly,a model could be built to express this relationship.Some key parameters were extracted and then used to establish links with Ro,such that a model to evaluate the remaining oil/gas yield with TOC and Rocould be obtained (Equation (3)).

    where Qpo is the remaining oil yield per mass unit of shale,Rois the vitrinite reflectance,and TOC is the total organic carbon.a1,a2,a3,a4,a5,and a6are empirical parameters,where a1and a6equal 0.99892 and 0.01538,respectively.When Ro≤0.76%,a2and a3equal 0.4265 and 0.7516,when 0.76%<Ro≤0.95%,a2and a3equal 0.4593 and 1.41,when Ro>0.95%,a2and a3equal 4.164 and 5.3161,when Ro≤0.77%,a4and a5equal 0.068 and 1.1297,when 0.77% <Ro≤1.06%,a4and a5equal 2.6881 and 3.2629,and when Ro>1.06%,a4and a5equal 3.5488 and 4.1449,respectively.

    where Qpg is the remaining gas yield per mass unit of shale,Rois the vitrinite reflectance,and TOC is the total organic carbon.b1to b15are empirical parameters and equal 1.0062,0.9478,0.5744,-0.0997,-1.1745,3.4118,2.1756,1.5235,-2.3651,-0.2334,2.9012,-2.9174,-0.0967,0.5035,and -0.4776,respectively.

    To test the confidence of the models,the remaining oil and gas yields were estimated by Equations (3) and (4) using TOC and Rodata obtained from this study.Through comparisons with experimental results,the correlation coefficients of both the oil and gas yields were greater than 0.997 (Fig.10).Specifically,the average absolute error between the calculated remaining oil yield and the experimental yield was 0.0001 mg/g·rock,and the absolute values of the absolute errors averaged 0.4796 mg/g·rock.The average absolute error between the calculated remaining gas yield and the experimental one was-0.0002 m3/t·rock,and the absolute values of the absolute errors averaged 0.3707 m3/t·rock.Such figures showed that the models built in this paper were reliable.

    Fig.9.Yield of remaining hydrocarbons with TOC and Ro.(a) remaining yield of oil,(b) remaining yield of gas.

    Fig.10.Relationship between calculated yield and experimental yield.(a) oil yield,(b) gas yield.

    Fig.11.Relationship of the cutoffs of the remaining oil yield versus Ro and TOC.

    Fig.12.Relations of the remaining oil yield cutoff versus minimum TOC cutoff and Ro

    4.3.Cutoffs of TOC and remaining oil yield

    Shale to be developed by in-situ transformation into oil have normally experienced some thermal evolution,and thus,have different remaining oil yields.Furthermore,shale with commercial development value also differs in TOC.Taking the shale samples used in our experiments as an example,the cutoffs of the remaining oil yield of shale were 50 mg/g·rock,45 mg/g·rock,40 mg/g·rock,35 mg/g·rock,30 mg/g·rock,25 mg/g·rock,20 mg/g·rock,and 15 mg/g·rock,respectively.The TOC cutoffs under different Rovalues could be obtained by using Equation (1).The results are shown in Fig.11.When Roremained unchanged,the TOC cutoff declined as the yield increased.When the cutoff of the remaining oil yield remained the same,the TOC cutoff declined slightly at first and then increased as Roincreased.The minimum value of the TOC cutoff occurred when Rowas approximately 0.75%.Therefore,the shale with Roof approximately 0.75% offered the best commerciality for development by in-situ transformation,and it was untrue that lower Rowas better.The reasons might be that after oil and/or gas are generated,they are absorbed first by organic matter(Baker,2016;Tissot et al.,1971;Stainforth and Reinders,1990;Larsen and Li,1997;Ertas et al.,2006;Kelemen et al.,2006;Han et al.,2017)and minerals (Schettler and Parmely,1991;Li et al.,2020).Only when the amount of the generated hydrocarbons is large enough,can microfractures be formed within the source rocks because of increased pressure (Jarvie,2012),and the oil and/or gas generated from kerogens can then move and be discharged from the source rocks.In summary,the source rock needs to reach a certain maturity to generate enough oil and gas and form an effective hydrocarbon expulsion channel.At this time,the oil and gas generated during in-situ transformation can be directly exploited from shale.As a result,the time to produce oil and/or gas is advanced,which helps to shorten the investment recovery period and improve project economics.

    Fig.13.Relationship between the maximum recoverable oil volume and the maximum recoverable gas volume of the shales.

    Fig.12 plots the relationship of the remaining oil yield cutoffs versus minimum TOC cutoffs and Ro,where it can be observed that the minimum TOC cutoff increased as the remaining oil yield cutoff increased,and there was a strong linear correlation between them.In contrast,the Rocorresponding to the minimum TOC cutoff decreased as the cutoff of the remaining oil yield increased.For instance,when the cutoffs of the remaining oil yield were 15 and 50 mg/g·rock,respectively,the Rovalues corresponding to the minimum TOC cutoff were 0.76% and 0.74%,respectively.

    5.Assessment of recoverable oil and gas resources in Chang 73 shale,Ordos Basin

    5.1.Geological settings and exploration

    Fig.14.Scheme of TOC low limit vs Ro.

    The Ordos Basin is located in the central part of North China(Fig.1),covering an area of approximately 37×104km2and contains 146.5×108t and 15.68×1012m3of conventional oil and gas resources,respectively.The Triassic Yanchang Formation.and the Jurassic Yan’an Formation are the major oil-bearing series in the basin.The oil source rocks belong to the Yanchang Formation deposited in semi-deep to deep lacustrine facies (Zhang et al.,2006;Li et al.,2019).Specifically,organic-rich shale is mainly developed in the Chang 73section.During the sedimentary period of this section,the maximum water depth of the ancient lake basin could reach 60-120 m,and the distribution area of organic-rich shale is approximately 5×104km2.By the end of 2019,the total resources in the Ordos Basin were more than 5.6 trillion tons,and the annual oil production was more than 36 million tons.The reservoirs of the Yanchang Formation are tight with low permeability.The average porosity and permeability were 9.9% and 1.23 mD,respectively.Approximately 81.6% of the resources were in reservoirs with a permeability lower than 5 mD,among which 54.3%of the reservoirs featured permeability lower than 1 mD.The estimated ultimate recovery factor of the resources was 15.8%.

    The Chang 73shale has low maturity with a Rovalue average of 0.82%.Tests indicate that it is impossible to realize commercial development using the existing volume fracturing of horizontal wells.The initial oil rates of the seven vertical wells were only 0.60-3.23 t/d after volume fracturing and declined rather rapidly.Taking Well G295 with the slowest decline rate as an example,the daily oil production was 1.09 t/d after producing for 1 year,and the annual oil production was 415.03 t.After producing 912 d,its cumulative oil production was only 623.6 t.All other wells were shut after producing 1-8 months.Estimation of GOR showed that the peak GOR of the Chang 73shale ranged from 60 m3/m3to 98 m3/m3,and averaged 83 m3/m3.There was liquid petroleum in the shale.In-situ conversion was feasible,although the existing technologies could not realize commercial development.

    5.2.Determination of key parameters

    The methods for determining the key parameters,such as TOC,Ro,and net pay of shale for assessing the recoverable oil/gas resources from in-situ conversion of shale are introduced in this paper using the Chang 73shale in the Ordos Basin as an example.

    5.2.1.Determination of TOC cutoff

    To evaluate whether or not the target is economical,it is necessary to determine the cutoffs of recoverable oil and gas volumes of the shale.Because there was a strong linear correlation between the maximum recoverable oil volume and the maximum recoverable gas volume (Fig.13),it was enough to determine only the cutoff of the recoverable oil volume.

    During the commercial development via in-situ shale transformation,the cutoff of recoverable oil volume could be determined with a model (Equation (5)) based on the oil production cutoff of a well group and the mass of the shale in the effective heating area controlled by this group.

    Where Qpo_limitis the cutoff (or lower limit) of recoverable oil volume per mass unit of shale,Qoil_limitis the lower limit of the cumulative oil produced by a well group,and Wtrockis the mass of shale in the effective heating area controlled by the well group.

    Horizontal wells were used in this paper.Specifically,there were 10 heating wells and one production well in the well group.The spacing between two heating wells was 15 m,and the lateral length of a horizontal well was 1,200 m.Based on these data,the rock mass in the effective heating area was 7,080,000 tons,the cutoff of the oil production was 100,000 tons,and the cutoff of recoverable oil volume of shale was 14 mg/g·rock as estimated with Equation (5).

    Fig.15.Distribution of TOC in Chang73 of Triassic in Ordos basin.

    The TOC cutoff corresponding to the cutoff of recoverable oil resources per mass unit of shale could be determined with Equation (3) by using the cutoff of recoverable oil volume determined above (Equation (6)).

    where TOClimitis the TOC lower limit (i.e.,cutoff) of the effective shale,Qpo_limitis the lower limit of recoverable oil volume per mass unit of shale,Rois the shale to be tested,and c1,c2,c3,c4,c5,and c6are empirical parameters.c1and c4are 1.0011 and -0.0154,whereas c2and c3are 0.068 and 1.1297 when Ro≤0.77%,and c2and c3equal 2.6881 and 3.2629 when 0.775%<Ro≤1.06%,and c2and c3are 3.5488 and 4.1449 when Ro>1.06%.

    Regarding c5 and c6,they were 0.4265 and 0.7516 when Ro≤0.76%,equaled 0.4593 and 1.41 when 0.76% <Ro≤1.0%,and they are 4.164 and 5.3161 when Ro>1.0%.

    Assuming the cutoff of recoverable oil volume per mass unit of rock is 14 mg/g·rock,then the plot of Roversus the TOClimitof the effective shale(Fig.14)could be obtained using Equation(5),from which it can be seen that the cutoff of TOC is approximately 6%when Rois approximately 0.8%.

    5.2.2.Determination of other key parameters

    5.2.2.1.TOC.The analysis of 8690 cores taken from the Chang 73 shale from 271 wells indicated that their TOC ranged from 6 to 39 wt%,and averaged 14.32 wt%.In this paper,TOC and Ro obtained from the laboratory tests were used to calibrate the log data using the ΔlogR method(Passey et al.,1990).Then,the TOC values of the target intervals in 791 wells were obtained by the interpretation model built by zoning.It should be noted that the TOC used here refers to the average TOC of the shale intervals satisfying the TOC cutoff criteria.To test the reliability of the calculation results,referring to the standard of effective shale section,the area and thickness weighted statistics of TOC values in the control area of each well point showed that the variation range was from 6 to 32.4 wt%,and averaged 14.12 wt%,which was very close to the value obtained from laboratory cores.This result indicated that the calculations were reliable.As Fig.15 shows,the Chang 73 shale with TOC greater than 6 wt% covered an area of approximately 5.1×104 km2,and the organic-rich region (TOC>18 wt%) was mainly distributed in Zhengning-Miaowan,Ningxian-Jingchuan,Heshui-Qingcheng,Maling-Huachi-Changguanmiao,and Maling-Gengwan-Jiyuan.Overall,organic-rich shale sections were distributed in succession,but there was still strong heterogeneity,which reflected the effects of the paleoenvironment on the development of high-quality source rocks.Nevertheless,the Chang 73 shale had a good material basis for in-situ transformation.

    Fig.16.Distribution of Ro in Chang73 of Triassic in Ordos Basin

    5.2.2.2.Ro.Altogether,2200 shale samples were taken from 321 wells to test and analyze the vitrinite reflectance of the Chang 73 shale.Based on the analysis results,the Ro distribution map was drawn for the Chang 73 shale.According to the results,the Chang 73 shale experienced uneven thermal evolution.Its Ro was mainly distributed in 0.52%-1.25%,with an average value of 0.82%.Among these,93% of the areas had a Ro less than 1.0% and 54.3% had a Ro less than 0.82%;these were mainly distributed in the south of Qincheng,Yanwu-Wuqi,and Gengwan-Yanchi.According to the previous analysis,shales in the above areas have good hydrocarbon generation capacity and should be used as primary targets for insitu transformation (Fig.16).

    Fig.17.Distribution of effective shale thickness in the Chang 73 of the Triassic in Ordos Basin.

    5.2.2.3.Net paythickness.The effective shale sections suitable for in-situ transformation were divided using the TOC cutoff of effective shale and the TOC value of logging interpretation to obtain the effective shale thickness and determine the favorable shale distribution area.The principles were as follows.①the TOC of shale was higher than its lower limit value (6 wt%),and the continuous thickness was greater than 2 m;or ②the cumulative thickness of shale section was greater than 2 m,and the accumulated thickness of TOC less than 6 wt%was less than 20%of the thickness of this section,and the single-layer thickness with average TOC less than 6 wt% was less than 1 m.According to the above principles,the plane distribution of effective shale thickness of the Chang 73member in the study area was drawn using the logging interpretation results of 791 wells(Fig.17).The results showed that the effective thickness of the Chang 73shale ranged from 2 m to 49 m,with an average of 17.5 m.Among these,the distribution area with an effective thickness greater than 20 m exceeded 1.5×104km2was mainly located in Jiyuan-Maling,the east of Wuqi-Nanliang,and Zhengning-Miaowan,where the average thickness of effective shale is as high as 27.5 m.Additionally,from the profile,the spatial distribution of the effective shale of Chang 73was also very stable,especially in the long axis direction of the paleo deep lake semi-deep lake.The effective shale distribution was continuous and stable,and the thickness change was small (Fig.18).In the short axis direction,affected by the ancient landform,the effective thickness of the shale changed substantially,but the main body was still relatively continuous with a large thickness (Fig.19).

    5.3.Evaluation method of recoverable oil and gas resources

    The recoverable oil and gas resources from in-situ shale transformation were assessed using the following steps:(1)Estimate the remaining oil and gas yield per mass unit of shale through Equations (3) and (4) according to TOC and Rodistribution of effective shale distribution area.(2) Estimate the abundance (or density) of recoverable oil and gas resources in the assessment area by Equations (7) and (8) using the effective shale thickness and shale density (obtained through log data).(3) Estimate the recoverable oil and gas resources based on the control area of evaluation points.

    Fig.18.Organic rich shale section of the Chang7 of Triassic in Ordos Basin (NW-SE)(Location is shown in Fig.1).

    Fig.19.Organic rich shale section of the Chang7 of Triassic in Ordos Basin (SW-NE)(Location is shown in Fig.1).

    Fig.20.Distribution of recoverable oil resource abundance in Chang73 shale of Triassic in Ordos Basin.

    where AOR is the abundance of recoverable oil resources in the assessment area,AGR is the abundance of recoverable gas resources in the assessment area,Qpois the oil recovery volume per mass unit of shale within the effective interval,Qpgis the gas recovery volume per mass unit of shale in the effective interval,and ρ is the shale density in the effective interval.

    The evaluation results showed that the recoverable oil resources of the Chang 73shale by in-situ conversion were approximately 450×108t,natural gas was approximately 30×1012m3,and the distribution area was approximately 4.27×104km2.The abundance of recoverable oil resources in the effective shale distribution area was 50×104t/km2-595×104t/km2,with an average value of 146×104t/km2.Among these,the area with technically recoverable oil resource abundance greater than 160×104t/km2was approximately 1.55×104km2,mainly distributed in Jiyuan-Gengwan-Maling,Wuqi and Qingcheng-Miaowan (Fig.20),which were the key areas of the Chang 73shale in-situ conversion.The distribution of recoverable natural gas resource abundance was consistent with this finding.The abundance of recoverable natural gas resources was 3×108m3/km2-33×108m3/km2,with an average value of 8.1×108m3/km2.The area with the resource abundance greater than 8×108m3/km2was approximately 1.66×104km2,mainly distributed in Jiyuan-Gengwan-Maling,Huachi and Qingcheng-Miaowan (Fig.21).It was demonstrated that the Chang 73shale in Ordos Basin has a large amount of recoverable oil and gas resources.If it could be effectively utilized,the crude oil production of in-situ transformation could have a resource base of 100 t/year and stable production for more than 100 years.

    6.Conclusions

    (1) In this paper,the hydrocarbon generation and expulsion process of shale under in-situ conversion were simulated using a large-capacity reactor,slow heating rate,and a semiopen system.It was revealed that TOC and Rowere the key factors affecting the in-situ transformation potential of shale.The remaining oil and gas yield of shale increased linearly with TOC,and the remaining oil yield increased slightly at first and then declined steadily and finally decreased rapidly as Roincreased;however,the remaining gas yield increased as Roincreased.

    Fig.21.Distribution of recoverable gas resource abundance in Chang73 shale of Triassic in Ordos Basin.

    (2) Roranged 0.75%-1.25% and 1.05%-2.3%,respectively,corresponding to the main oil generation stage and gas generation stage of shale in-situ transformation.The TOC of shale suitable for in-situ conversion should be greater than 6%,whereas its Roshould be less than 1.0%.Shales with 0.75%(Ro) could obtain the best economic benefit.An evaluation model of the remaining oil/gas yield was established.

    (3) The organic-rich intervals in the Triassic Chang 73shale of the Ordos Basin were characterized as being continuous and having a large thickness and high TOC.They were the optimal targets for in-situ conversion.The resource assessment results of Chang 73shale showed that the recoverable oil resources of in-situ transformation were approximately 450×108t,whereas the recoverable gas resources were 30×1012m3,and the distribution area was approximately 4.27×104km2.If it could be used effectively,the crude oil output of in-situ transformation in the Ordos Basin has a resource base of 100 t/year and stable production for more than 100 years.

    Acknowledgments

    This work was supported by PetroChina Co Ltd.(Grant number:2015D-4810-02;2018YCQ03;2021DJ52) and National Natural Science Foundation of China(Grant number:42172170).We thank the two anonymous reviewers for their valuable suggestions.

    黄色毛片三级朝国网站| 国产精品一区二区在线观看99| 国产免费一级a男人的天堂| av福利片在线| 91久久精品国产一区二区三区| 亚洲高清免费不卡视频| 国产精品久久久久久久久免| 91精品三级在线观看| 男女边摸边吃奶| 日韩大片免费观看网站| 美国免费a级毛片| 香蕉精品网在线| 蜜桃在线观看..| 寂寞人妻少妇视频99o| 26uuu在线亚洲综合色| 青春草亚洲视频在线观看| 高清欧美精品videossex| 精品99又大又爽又粗少妇毛片| kizo精华| 亚洲精品乱码久久久久久按摩| 精品久久蜜臀av无| 国内精品宾馆在线| 国产乱人偷精品视频| 国产又色又爽无遮挡免| 亚洲成av片中文字幕在线观看 | 十八禁网站网址无遮挡| 日本欧美视频一区| 波野结衣二区三区在线| 视频区图区小说| 久久久久久久精品精品| 久热久热在线精品观看| 国产在线视频一区二区| 日韩伦理黄色片| 国产成人av激情在线播放| 曰老女人黄片| 亚洲欧美色中文字幕在线| 国产成人精品在线电影| 99热网站在线观看| 大香蕉97超碰在线| 精品国产一区二区三区四区第35| 纵有疾风起免费观看全集完整版| 夫妻午夜视频| 99热全是精品| 精品人妻熟女毛片av久久网站| 欧美亚洲 丝袜 人妻 在线| 国产精品欧美亚洲77777| 日本欧美国产在线视频| h视频一区二区三区| 日本黄大片高清| 一本—道久久a久久精品蜜桃钙片| 99九九在线精品视频| 中文字幕免费在线视频6| 日韩成人伦理影院| 69精品国产乱码久久久| 国产成人免费无遮挡视频| 亚洲成人一二三区av| 色5月婷婷丁香| 亚洲精品久久久久久婷婷小说| 欧美国产精品一级二级三级| 青春草国产在线视频| 欧美老熟妇乱子伦牲交| 男女边摸边吃奶| 丰满乱子伦码专区| 另类精品久久| 女人久久www免费人成看片| 欧美人与性动交α欧美精品济南到 | av播播在线观看一区| av一本久久久久| 狂野欧美激情性bbbbbb| 尾随美女入室| 免费播放大片免费观看视频在线观看| 亚洲精品av麻豆狂野| 午夜精品国产一区二区电影| 大陆偷拍与自拍| 一区二区av电影网| 欧美人与性动交α欧美软件 | 丰满乱子伦码专区| 在线免费观看不下载黄p国产| 成年动漫av网址| 一级,二级,三级黄色视频| √禁漫天堂资源中文www| 亚洲国产精品成人久久小说| 国产伦理片在线播放av一区| 最黄视频免费看| av国产精品久久久久影院| 亚洲精品456在线播放app| 十八禁网站网址无遮挡| 菩萨蛮人人尽说江南好唐韦庄| 99九九在线精品视频| 在线观看人妻少妇| 永久免费av网站大全| 九色成人免费人妻av| 日韩 亚洲 欧美在线| 赤兔流量卡办理| 日产精品乱码卡一卡2卡三| 中国国产av一级| 国产精品国产三级国产专区5o| a级片在线免费高清观看视频| 国产极品天堂在线| 一二三四中文在线观看免费高清| 老女人水多毛片| 寂寞人妻少妇视频99o| 啦啦啦在线观看免费高清www| 欧美国产精品va在线观看不卡| 久久国产精品男人的天堂亚洲 | 热re99久久精品国产66热6| 亚洲精品美女久久久久99蜜臀 | 一本久久精品| 国产成人av激情在线播放| 高清不卡的av网站| 欧美国产精品va在线观看不卡| 中文乱码字字幕精品一区二区三区| 国产深夜福利视频在线观看| 寂寞人妻少妇视频99o| 国产成人精品婷婷| 一边摸一边做爽爽视频免费| 91精品三级在线观看| 最黄视频免费看| 91精品国产国语对白视频| 亚洲av在线观看美女高潮| 亚洲成人av在线免费| 亚洲av.av天堂| 国产精品久久久久成人av| 欧美精品高潮呻吟av久久| 亚洲国产精品一区二区三区在线| 久久热在线av| 九九爱精品视频在线观看| 国产男女超爽视频在线观看| 九色成人免费人妻av| 少妇熟女欧美另类| 色94色欧美一区二区| 桃花免费在线播放| 少妇 在线观看| 深夜精品福利| av片东京热男人的天堂| 99热6这里只有精品| 人人妻人人爽人人添夜夜欢视频| 精品国产一区二区久久| 中文乱码字字幕精品一区二区三区| 成人午夜精彩视频在线观看| 如日韩欧美国产精品一区二区三区| 一区二区三区四区激情视频| 夜夜爽夜夜爽视频| 成人手机av| 91在线精品国自产拍蜜月| 国产在线一区二区三区精| 你懂的网址亚洲精品在线观看| 久久久国产一区二区| 日韩 亚洲 欧美在线| 国产毛片在线视频| 最近最新中文字幕免费大全7| 亚洲四区av| 男女无遮挡免费网站观看| 国产极品天堂在线| 夜夜骑夜夜射夜夜干| 日韩一区二区视频免费看| 国产精品成人在线| 中文字幕精品免费在线观看视频 | 久久精品久久精品一区二区三区| 国产精品久久久av美女十八| 五月伊人婷婷丁香| 中文乱码字字幕精品一区二区三区| 亚洲,欧美精品.| 26uuu在线亚洲综合色| 99久国产av精品国产电影| 久久久欧美国产精品| 99香蕉大伊视频| 美国免费a级毛片| 80岁老熟妇乱子伦牲交| 99热网站在线观看| 亚洲精品成人av观看孕妇| 国产在线免费精品| 国国产精品蜜臀av免费| 丝瓜视频免费看黄片| 欧美精品一区二区大全| 久久99热这里只频精品6学生| 精品人妻一区二区三区麻豆| 欧美亚洲 丝袜 人妻 在线| 亚洲成人手机| 18禁国产床啪视频网站| 五月玫瑰六月丁香| 免费看av在线观看网站| 亚洲精品aⅴ在线观看| 国产成人av激情在线播放| 精品亚洲成a人片在线观看| 久久久亚洲精品成人影院| 亚洲色图综合在线观看| 精品人妻在线不人妻| 如日韩欧美国产精品一区二区三区| av国产精品久久久久影院| 国产精品蜜桃在线观看| 国产午夜精品一二区理论片| 精品国产一区二区久久| 国产精品嫩草影院av在线观看| 亚洲精品久久成人aⅴ小说| 欧美精品人与动牲交sv欧美| 久久久久久久精品精品| 亚洲国产精品国产精品| 日韩av在线免费看完整版不卡| 色哟哟·www| 人妻系列 视频| 久久av网站| 91aial.com中文字幕在线观看| 两个人看的免费小视频| 国产一级毛片在线| 成年美女黄网站色视频大全免费| 在线观看国产h片| 久久99热6这里只有精品| 精品少妇黑人巨大在线播放| 成人国产av品久久久| 大片免费播放器 马上看| 91成人精品电影| 成人毛片60女人毛片免费| www.熟女人妻精品国产 | 午夜av观看不卡| 日韩成人av中文字幕在线观看| 999精品在线视频| 亚洲精品美女久久av网站| 亚洲av成人精品一二三区| 97在线视频观看| 2021少妇久久久久久久久久久| 久久婷婷青草| 中文欧美无线码| 久久国产亚洲av麻豆专区| 日韩制服丝袜自拍偷拍| 免费观看av网站的网址| 久久97久久精品| 99久久人妻综合| 日本免费在线观看一区| www.熟女人妻精品国产 | 国产女主播在线喷水免费视频网站| 成人国产av品久久久| 亚洲成av片中文字幕在线观看 | 亚洲国产精品一区二区三区在线| av在线观看视频网站免费| 18禁动态无遮挡网站| 巨乳人妻的诱惑在线观看| 香蕉丝袜av| 精品久久蜜臀av无| 一二三四中文在线观看免费高清| 99久久精品国产国产毛片| 婷婷成人精品国产| 中国美白少妇内射xxxbb| 亚洲一区二区三区欧美精品| 国产一区亚洲一区在线观看| 美女福利国产在线| 香蕉精品网在线| 成人免费观看视频高清| 制服诱惑二区| 中国国产av一级| 欧美xxxx性猛交bbbb| 国产又爽黄色视频| 亚洲欧美一区二区三区黑人 | 久久久久久久久久久免费av| 成年美女黄网站色视频大全免费| 视频在线观看一区二区三区| 天天操日日干夜夜撸| 精品久久久久久电影网| 一二三四中文在线观看免费高清| 亚洲五月色婷婷综合| 午夜福利视频在线观看免费| 99精国产麻豆久久婷婷| 少妇的逼好多水| 国产片特级美女逼逼视频| 欧美3d第一页| 色94色欧美一区二区| 国产精品无大码| 美女内射精品一级片tv| 三上悠亚av全集在线观看| 国产一区二区三区综合在线观看 | 国产一区二区三区av在线| 嫩草影院入口| 丁香六月天网| 熟女电影av网| 国产麻豆69| 精品久久久精品久久久| 亚洲人成网站在线观看播放| 亚洲精品视频女| 久久ye,这里只有精品| 久久久久国产精品人妻一区二区| 亚洲成av片中文字幕在线观看 | 宅男免费午夜| 午夜av观看不卡| 天堂8中文在线网| 九色亚洲精品在线播放| 中国美白少妇内射xxxbb| 亚洲精品乱码久久久久久按摩| 多毛熟女@视频| 丝袜脚勾引网站| 久久婷婷青草| 欧美精品一区二区免费开放| 丰满迷人的少妇在线观看| 狂野欧美激情性bbbbbb| 日韩一区二区视频免费看| 少妇高潮的动态图| 国产成人91sexporn| 又黄又粗又硬又大视频| 纯流量卡能插随身wifi吗| 国产一区二区激情短视频 | 丰满乱子伦码专区| 国产在视频线精品| a级毛色黄片| 国产日韩欧美在线精品| 精品国产一区二区三区久久久樱花| 亚洲经典国产精华液单| av免费观看日本| 久久久精品94久久精品| 一级毛片我不卡| 在线观看www视频免费| 极品少妇高潮喷水抽搐| 国产成人精品福利久久| 国产黄频视频在线观看| 亚洲情色 制服丝袜| videos熟女内射| 国产极品天堂在线| 欧美精品亚洲一区二区| 国产有黄有色有爽视频| 精品久久久久久电影网| 三上悠亚av全集在线观看| 中文天堂在线官网| 好男人视频免费观看在线| 国产精品人妻久久久久久| 亚洲激情五月婷婷啪啪| 久久人妻熟女aⅴ| 黑丝袜美女国产一区| 亚洲美女视频黄频| 91精品国产国语对白视频| 免费少妇av软件| 欧美精品一区二区免费开放| 国产成人免费无遮挡视频| 精品人妻在线不人妻| 国产精品偷伦视频观看了| 又黄又粗又硬又大视频| 午夜影院在线不卡| 久久久久久久久久人人人人人人| 成人手机av| 久久久久视频综合| 欧美xxxx性猛交bbbb| 久久精品aⅴ一区二区三区四区 | 国产免费一区二区三区四区乱码| 乱码一卡2卡4卡精品| 久久久久久久久久人人人人人人| 精品一区二区三区视频在线| 亚洲欧美一区二区三区国产| 99国产综合亚洲精品| 国产日韩一区二区三区精品不卡| 深夜精品福利| 国产激情久久老熟女| 女的被弄到高潮叫床怎么办| 欧美日韩精品成人综合77777| 欧美日韩av久久| 久久99一区二区三区| 亚洲高清免费不卡视频| 免费大片黄手机在线观看| 看免费av毛片| 国产精品不卡视频一区二区| 97在线视频观看| 久久久国产欧美日韩av| 最新中文字幕久久久久| 日日啪夜夜爽| 免费观看性生交大片5| 免费在线观看完整版高清| 中文字幕最新亚洲高清| 欧美97在线视频| 日韩中字成人| 国产欧美日韩一区二区三区在线| 久久99精品国语久久久| 亚洲精品第二区| 永久网站在线| 啦啦啦视频在线资源免费观看| 日韩,欧美,国产一区二区三区| 一边亲一边摸免费视频| 久久久久久久亚洲中文字幕| 男人操女人黄网站| 建设人人有责人人尽责人人享有的| 最近中文字幕2019免费版| 欧美日韩国产mv在线观看视频| 免费在线观看完整版高清| 欧美性感艳星| 免费看不卡的av| av又黄又爽大尺度在线免费看| 丝瓜视频免费看黄片| 菩萨蛮人人尽说江南好唐韦庄| 国产乱来视频区| 亚洲第一区二区三区不卡| 视频中文字幕在线观看| av有码第一页| 91在线精品国自产拍蜜月| 国产精品人妻久久久久久| 各种免费的搞黄视频| 欧美日韩亚洲高清精品| 亚洲伊人久久精品综合| 国产一级毛片在线| 2022亚洲国产成人精品| av女优亚洲男人天堂| 天天躁夜夜躁狠狠久久av| 在现免费观看毛片| 亚洲成人av在线免费| 高清不卡的av网站| 91精品三级在线观看| 久久午夜综合久久蜜桃| 亚洲精品久久成人aⅴ小说| 亚洲av成人精品一二三区| 国产乱人偷精品视频| 制服诱惑二区| 国产69精品久久久久777片| 精品福利永久在线观看| 最近的中文字幕免费完整| av.在线天堂| 中文字幕av电影在线播放| 亚洲精华国产精华液的使用体验| 2018国产大陆天天弄谢| 国产精品国产av在线观看| 少妇的逼好多水| 日日啪夜夜爽| 亚洲欧洲国产日韩| 超碰97精品在线观看| 人成视频在线观看免费观看| 日韩欧美精品免费久久| 夫妻性生交免费视频一级片| 久久97久久精品| 欧美最新免费一区二区三区| 日本av手机在线免费观看| 两个人免费观看高清视频| 亚洲国产欧美在线一区| 亚洲图色成人| 黑人高潮一二区| 亚洲美女搞黄在线观看| 人人妻人人澡人人看| 欧美3d第一页| 最新中文字幕久久久久| 中文字幕制服av| 日本与韩国留学比较| 亚洲av免费高清在线观看| 久久免费观看电影| 午夜激情av网站| 肉色欧美久久久久久久蜜桃| 国产精品女同一区二区软件| 国产一区二区激情短视频 | 色网站视频免费| 一二三四中文在线观看免费高清| 不卡视频在线观看欧美| 毛片一级片免费看久久久久| 免费看光身美女| 搡老乐熟女国产| 成年动漫av网址| 中文字幕最新亚洲高清| 一区二区日韩欧美中文字幕 | 99热全是精品| 国产精品 国内视频| 亚洲成人一二三区av| 午夜日本视频在线| 日韩一区二区三区影片| 亚洲av成人精品一二三区| 精品99又大又爽又粗少妇毛片| 最近最新中文字幕免费大全7| videos熟女内射| 香蕉国产在线看| 色婷婷久久久亚洲欧美| 天堂中文最新版在线下载| 美女内射精品一级片tv| 男的添女的下面高潮视频| 女性被躁到高潮视频| 亚洲av国产av综合av卡| 欧美xxxx性猛交bbbb| 天天操日日干夜夜撸| a 毛片基地| 国产女主播在线喷水免费视频网站| 麻豆精品久久久久久蜜桃| 欧美日韩成人在线一区二区| 欧美 日韩 精品 国产| 美女主播在线视频| 超碰97精品在线观看| 中文欧美无线码| av电影中文网址| 91午夜精品亚洲一区二区三区| 亚洲欧美一区二区三区国产| 欧美日韩视频精品一区| 欧美亚洲日本最大视频资源| 女人精品久久久久毛片| 丁香六月天网| 一本—道久久a久久精品蜜桃钙片| 少妇被粗大的猛进出69影院 | 欧美日韩av久久| 街头女战士在线观看网站| 国产乱来视频区| 超色免费av| 欧美xxxx性猛交bbbb| 99香蕉大伊视频| 午夜精品国产一区二区电影| 国产欧美日韩综合在线一区二区| 国产精品久久久av美女十八| 亚洲精品乱久久久久久| 青春草视频在线免费观看| 欧美日韩视频高清一区二区三区二| 免费在线观看完整版高清| xxxhd国产人妻xxx| 久久久久久久久久人人人人人人| 最后的刺客免费高清国语| 大陆偷拍与自拍| 最近2019中文字幕mv第一页| 黑丝袜美女国产一区| 搡老乐熟女国产| 午夜免费男女啪啪视频观看| 看免费av毛片| 国产精品人妻久久久久久| 夫妻午夜视频| 久热这里只有精品99| 91久久精品国产一区二区三区| 午夜视频国产福利| 高清不卡的av网站| 久久精品人人爽人人爽视色| 大片免费播放器 马上看| 国产精品 国内视频| 国产精品99久久99久久久不卡 | 久久精品国产鲁丝片午夜精品| 国产精品免费大片| 免费看不卡的av| 中国国产av一级| 国产男女超爽视频在线观看| 国内精品宾馆在线| 欧美精品一区二区大全| 国产精品一区www在线观看| 99国产综合亚洲精品| www.av在线官网国产| 亚洲国产av新网站| 黄网站色视频无遮挡免费观看| 国产高清不卡午夜福利| 国产成人av激情在线播放| 一区二区日韩欧美中文字幕 | 最近手机中文字幕大全| 亚洲人成网站在线观看播放| 免费在线观看黄色视频的| 亚洲精品日本国产第一区| 看十八女毛片水多多多| 欧美日韩视频高清一区二区三区二| 免费看不卡的av| 国产高清国产精品国产三级| 中文字幕亚洲精品专区| 欧美老熟妇乱子伦牲交| 成人国语在线视频| 国产精品久久久久久久久免| 欧美精品一区二区大全| 国产成人a∨麻豆精品| 午夜视频国产福利| 久久久国产一区二区| 另类精品久久| 欧美日韩一区二区视频在线观看视频在线| 日韩成人av中文字幕在线观看| 国产日韩欧美在线精品| 国产熟女午夜一区二区三区| 中文字幕精品免费在线观看视频 | 久久99一区二区三区| 久久韩国三级中文字幕| 日韩在线高清观看一区二区三区| 七月丁香在线播放| 国产成人欧美| 午夜福利在线观看免费完整高清在| 国产av国产精品国产| 国产午夜精品一二区理论片| 桃花免费在线播放| 999精品在线视频| 日韩一本色道免费dvd| av电影中文网址| 国产免费现黄频在线看| 青青草视频在线视频观看| 免费人成在线观看视频色| 大片电影免费在线观看免费| 日韩欧美一区视频在线观看| av卡一久久| 全区人妻精品视频| 日本欧美国产在线视频| 99精国产麻豆久久婷婷| 久久99热6这里只有精品| 伊人亚洲综合成人网| 在线免费观看不下载黄p国产| 免费高清在线观看视频在线观看| 久久久国产欧美日韩av| 多毛熟女@视频| 精品一品国产午夜福利视频| 黄色怎么调成土黄色| 久久久久久久久久人人人人人人| 国产精品女同一区二区软件| 欧美变态另类bdsm刘玥| 交换朋友夫妻互换小说| 有码 亚洲区| 久久久国产欧美日韩av| 777米奇影视久久| av线在线观看网站| 中文欧美无线码| 好男人视频免费观看在线| 99视频精品全部免费 在线| 男男h啪啪无遮挡| 看免费成人av毛片| 女性生殖器流出的白浆| 日韩精品有码人妻一区| 久久99热6这里只有精品| 满18在线观看网站| 国内精品宾馆在线| 久久久久久久久久成人| 一级a做视频免费观看| 成人综合一区亚洲| 有码 亚洲区| 夜夜爽夜夜爽视频| 色婷婷av一区二区三区视频| 青春草国产在线视频| 中文欧美无线码| 国产69精品久久久久777片| 亚洲精品,欧美精品| 午夜精品国产一区二区电影| 看非洲黑人一级黄片| 蜜桃国产av成人99| 99热网站在线观看| 丝袜喷水一区| 波多野结衣一区麻豆|