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

    A modified soil water content measurement technique using actively heated fiber optic sensor

    2020-07-12 12:36:24MengWngXuLiLihongChenSenqunHouGuiynWuZhilinDeng

    Meng Wng, Xu Li,b, Lihong Chen,*, Senqun Hou, Guiyn Wu, Zhilin Deng

    a Key Laboratory of Urban Underground Engineering of Ministry of Education, Beijing Jiaotong University, Beijing, China

    b Qinghai Research and Observation Base, Key Laboratory of Highway Construction & Maintenance Technology in Permafrost Regions, Ministry of Transport,Xining, Qinghai, China

    Abstract Soil water content measurement is critical in practical engineering. The actively heated fiber Bragg grating optic sensor (FBGS) has great potential of multi-point measurement for soil water content measurement in field.In this study,the effect of heating time on the measurement accuracy is discussed,and modifications are made for actively heated fiber optic(AHFO) sensors.The results demonstrate that if an integration data analysis method is used, the accuracy and reliability of soil water content measurement with AHFO sensors will be improved. Both a short fiber length and a short-term heating pattern are effective and can help to reduce soil disturbance. With the proposed integration method, a short heating time is guaranteed for measuring the soil water content. Such improvements will reduce the thermal disturbance to soil sample and improve the reliability of measurement.

    2020 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    Keywords:Fiber bragg grating Carbon fiber heated sensor Unsaturated soil Field monitoring Soil water content measurement

    1. Introduction

    In recent years, fiber optic sensing technology has been developed rapidly and has achieved great success in both field monitoring system and laboratory model test. Compared to other sensing technologies, the advantages of fiber optic technology are obvious, e.g. it is capable of both distributed long-distance measurements and indoor tests,with anti-corrosion and anti-jamming characteristics (Freifeld et al., 2008; Selker et al., 2006; Westhoff et al., 2007;Tyler et al., 2008, 2009; Kim et al., 2010). Grattan and Sun, (2000) reviewed the types of fiber optic sensors and the scope of application of various types in 1960e2000. In a wide variety of fiber optic sensors,Bragg gratings have some advantages in the measurement of strain and temperature due to its ease for use.For indoor tests, Lee et al. (2010) had assembled the fiber optic testing system into the triaxial testing equipment and successfully measured the pressure/force and linear displacement.In such tests,conventional displacement sensors are easily short-circuited when immersed underwater or subjected to electromagnetic interference, but fiber optic sensors will not hold.

    It is accepted that the specific heat capacity and thermal conductivity of soil are significantly affected by soil water content(e.g.Mitchell and Soga, 2005; Lu et al., 2007). On one hand, thermal properties are essential in the thermal related problem of soils,e.g.frost deformation (Liu et al., 2019), expansion of saline soil (Fang et al., 2018), and thermal consolidation (Bai, 2006). On the other hand,the thermal properties of unsaturated soils can be measured and further used to infer soil water content.Many efforts have been made to estimate the water content by measuring the thermal property values of the soil(Kersten,1949;Johansen,1977;Brandon and Mitchell,1989;Ochsner et al.,2001;Abu-Hamdeh,2003;C?té and Konrad, 2005; Lu et al., 2007; Usowicz et al., 2017; Oyeyemi et al., 2018). Halloran et al. (2016) also reviewed the thermal methods in unsaturated soils. In addition, the thermal methods applicable in unsaturated soils were reviewed by Halloran et al.(2016) and the methods can be divided into two categories: the passive heat monitoring methods and the active heat tracing methods(including heating).

    Fiber Bragg grating optic sensor (FBGS) is a fast developing method for temperature and moisture measurement, from simply monitoring the temperature(passive)to involving a heater(active).Recently, the actively heated fiber optic (AHFO) and the passive fiber optic (PFO) sensors have been used to measure soil water content based on thermal effect. Steele et al. (2010) used PFO sensors to measure soil temperature changes caused by diurnal temperature cycle difference,in order to estimate water content in the soil, which could be exactly the passive heat tracing methods.However,this method could obtain good measurement results only when the temperature difference in the diurnal cycle is large.

    Many researchers have designed AHFO sensors with different structures and verified their suitability. Sayde et al. (2010)improved Steele et al. (2010)’s method by replacing the original fiber with a kind of AHFO sensor. Moreover, they proposed an empirical relation between the temperature variation characteristic values and soil water content. A large number of embedded AHFO sensors were used to predict the occurrence of seepage by monitoring the temperature change of the dam (C?té et al., 2007; Xiao and Huang, 2013). Later, Cao et al. (2015, 2018) developed an AHFO sensor to detect the volumetric water content,which shows much higher accuracy compared to that used by Sayde et al.(2010).Because of their efforts,measurement of soil water content by fiber sensors has become practical both in laboratory and field. This technique will be a good alternative to the conventional measurement methods such as time domain reflectometer (TDR) method(Dalton et al.,1984;Malicki et al.,1992,1996;Roth et al.,1992;Chen et al., 2007, 2009, 2013) or frequency domain reflectometry (FDR)method(Heimovaara et al.,1996;Veldkamp and O’Brien,2000;Lin,2003; Jones and Or, 2004; Tsheko and Savage, 2005), which has been stocked for many years and cannot overcome the difficulties of measurement accuracy and sensitivity to the contact condition between sensor and soil matrix.

    As a young technique with great economic advantage for multipoint measurement, AHFO sensor has a great potential for soil water content measurement and is continuously improved. However, it still has some limitations, for example:

    (1) The previous researches mainly focused on measuring the thermal conductivity () to calculate the water content.However, it is a rather complicated process (Kersten,1949;Johansen, 1977; Sepaskhah and Boersma, 1979; Lu et al.,2007; Ciocca et al., 2012; Striegl and Loheide, 2012). Sayde et al. (2010) proposed a simple response variable (Tcum),which monotonically varies with soil water content, to calibrate the moisture content, but they did not propose a specific calibration equation forTcumand water content. Later,Cao et al. (2015, 2018) proposed a modified method for sensors.However,the size of the sample used to calibrate the sensor is too large, making it difficult to control the dry density (400 mm in diameter and 460 mm in height (Cao et al., 2015); 1 m in height and 0.15 m in diameter (Cao et al., 2018)).

    (2) Soil disturbance may be still high due to the long-term heating.If the heating process is conducted with a high-power and long time, migration of water content in the soil sample will definitely occur(heating time and power are 120 s and 20 W/m, respectively (Sayde et al., 2010); and 300 s and 36 W/m,respectively (Ciocca et al., 2012)). Attempts to reduce soil disturbance need to be studied,such as lowering the heating power or reducing the heating time. It is a balance between remarkable signals and limiting soil disturbance.

    (3) The data interpretation technique should be further improved.In soil water content measurement, a well-accepted form for characterizing the relation between electric constant and soil water content has been established. However, the relations between temperature parameters and soil water content in Sayde et al. (2010) and Cao et al. (2018) cannot satisfy the users’ need. For example, Cao et al. (2015, 2018)’s relation depended on human judgment of a rising stage and a stable stage and has a complex expression. Sayde et al. (2010)’s relation was lack of data verification.

    In other words, the measurement accuracy, data analysis method, sensor design and heating pattern can be further improved.In this study,some improvements on the AHFO sensors are discussed.To achieve this,a series of laboratory tests is carried out, and four modifications are examined, including data analysis method, heating pattern, calibration method, and sensor length.

    2. Soil water content measurement by AHFO sensors

    Change in water content will lead to variations in the specific heat capacity and thermal conductivity of the soil, which is the basis for measuring the water content by the AHFO sensors(Mitchell and Soga, 2005; Cao et al., 2015, 2018). Thus, the AHFO sensors are essentially the combination of optical fiber temperature measurement technology and small heaters.

    2.1. Temperature measurement by FBGS

    Hartog (1983) firstly proposed the distributed temperature sensing (DTS) technology based on Raman optical time-domain reflectometry. Spontaneous Raman scattering yields by a collision occurring between fiber molecules and photons after a pumped laser pulse when a certain amount of energy is launched into a sensing optical fiber. As a result, two light photons with different wavelengths are generated, i.e. Stokes light and anti-Stokes light(Grattan and Sun,2000).According to the relationship between the anti-Stokes and Stokes intensity ratio, the temperature measurement of the optical fiber can be realized (Tyler et al., 2009).

    The FBGSs used by Hartog(1983)and C?té et al.(2007)have an accuracy of only 1C. However, with the innovation of the interrogator, the measurement accuracy of the FBGSs can reach 0.1C(Cao et al.,2018),and higher accuracy of temperature measurement will reduce the heating energy required for recording valuable data during water content measurement.

    2.2. Temperature duration curve during active heating process

    The unsaturated soil at room temperature is generally composed of the solid phase, liquid phase, and gas phase.Furthermore,the thermal conductivity and specific heat capacity of each component are different. When soil is heated by a constant linear heat source in a semi-infinite space,the temperature change curve of the soil is shown in Fig.1 (Campbell,1985; Shiozawa and Campbell, 1990). There are obvious segmented characteristics of the heating and cooling curves, thus it is necessary to explore the data of different segments separately.

    Fig.1. Temperature change curve of the soil heated by a constant linear heat source.The expression of the rapid rising stage and tapering stage comes from Cao et al.(2015,2018).

    2.3. Calibrated relation between soil water content and temperature characteristic value

    The typical relationship between thermal conductivity and water content for unsaturated soil undergoing drying process is shown in Fig.2.If the water content falls below a certain value(or a certain range),the soil thermal conductivity decreases significantly(Johansen, 1977; C?té and Konrad, 2005; Lu et al., 2007). In the authors’opinion,the critical water content in some models is close to the residual one of soil. As soil water content goes below the residual water content, water phase changes from a continuous distribution to a discontinuous distribution, subsequently the conductivity decreases significantly.

    The factors affecting soil thermal conductivity mainly include volumetric water content, soil density, and grain size distribution.Among them, the volumetric water content is the most important(Kersten,1949; Sepaskhah and Boersma,1979). Some empirical or theoretical relations have been proposed. For example, Kersten(1949) proposed the following empirical equation to calculate the thermal conductivity using soil water content for an unsaturated soil:

    Johansen (1977) proposed a predictive model for estimation of soil thermal conductivity () based on the degree of saturation,porosity,and bulk density of the soil.Later,C?té and Konrad(2005)simplified Johansen (1977)’s model by introducing a soil texture parameter that depends on soil type and a concise equation for estimating the thermal conductivity of dry soil. Lu et al. (2007)expanded C?té and Konrad (2005)’s work and accomplished the model in the entire range of soil water content.

    The original equations proposed by Johansen (1977)were

    Fig. 2. Typical relationship between thermal conductivity and water content for an unsaturated soil.

    The modified equations for calculatingKeanddryproposed by Lu et al. (2007)were

    whereKeis the normalized thermal conductivity proposed by Johansen (1977);satanddryare the thermal conductivities of saturated and dry soils(W m1K1),respectively;sis a parameter determined by the quartz content (q);qis the thermal conductivity of quartz, andq7.7 W m1K1;ois the thermal conductivity of other minerals, which is taken as 2 W m1K1whenq> 0.2, and 3 W m1K1whenq0.2;wis the conductivity of water, andw0:594 W m1K1at 20C;Sris the degree of saturation;eis the base of natural logarithm; a is a soil texture dependent parameter;1.33 is a shape parameter value;nis the soil porosity;andaandbare the empirical parameters.

    Definitely, the thermal conductivity highly depends on soil water content. However, it is difficult to accurately measure thermal conductivity and use it to infer soil water content. Hence, a substituting strategy is used to infer soil water content from soil thermal phenomena, e.g. the heating or decaying process of soil.Sayde et al.(2010)tried to infer soil water content from a response variable (Tcum) that monotonically varies with soil water content and can be measured by the DTS measurement method:

    whereTcumis the cumulative temperature increase(C s);t0is the total time of integration; andTis the DTS reported temperature change from the prepulse temperature (C). This is the first time that the water content is calibrated by using a variable that does not contain thermal conductivity ().

    Later, based on the test data of the tapering stage (Fig. 1), a temperature characteristic valueTtis proposed by Cao et al.(2015,2018) as

    wherenis the total record number of the tapering stage;iis the record number; andTiis theith record of temperature. The soil water content is linked toTtby a piece-wise function (Cao et al.,2018) as

    In fact,a linear relation was used in Cao et al.(2015)and it had a lower accuracy compared to the piece-wise function expressed by Eqs.(6)and(7).However,the drawbacks of this calibration method are also obvious:

    Fig. 3. Relation between soil water content and temperature characteristic value.

    (1) It is hard to accurately identify the two stages shown in Fig.1,that is to say, the calculation ofTtmay depend on human judgments and have some errors.

    (2) Only the data at the tapering stage are used.

    (3) In Eq.(7),Ttis expressed as a function ofw.In the application of this technique,wshould be inferred fromTt. An inverse function of Eq.(7)is to be solved and some other errors may be induced.

    3. Soil water content measurement technique modifications by AHFO sensors

    3.1. An integration data analysis method

    The integration data analysis method was first proposed by Sayde et al. (2010); however, due to the limited resolution and precision of the fiber interrogation, this method can only qualitatively measure the water content. The integration data analysis method used by the authors can be described as

    whereFtis an integrated temperature characteristic value,tis the total heating time, andfis a function to be determined.

    AfterFtis obtained, the temperature characteristic valueFtandwhave the following relations (Fig. 4):

    wherea,b,k,andcare the fitting parameters,and these parameters could be determined by indoor calibration test.Similar to0in Cao et al. (2018),tis a turning point onFt-wcurve. The test results demonstrate that above and belowt,w-Ftcurve shows quite different trend.tcan be determined by curve fitting method,such as the least square error method.Interestingly,tseems to be close to the water content at the air entry value on soil water characteristic curve (SWCC) (Fig. 4b) and corresponds to the degree of saturation for the soil around 0.8e0.9. Hence, if a small water content interval is used when the degree of saturation for the soil ranges between 0.8 and 0.95, a more accuratetcan be obtained.

    Fig.4. Comparison of w-Ft curve and SWCC:(a)Relation between soil water content and integrated temperature characteristic value; and (b) Typical curve of volumetric water content versus matric suction.

    Ciocca et al. (2012) calibrated the volumetric water content using the temperature data from the cooling process in Fig.1, and proved that the data from the cooling part could be used to calibrate the soil water content. In this paper, an integration data analysis method was used to analyze the data of the heating and cooling processes, and its effectiveness was compared.

    3.2. Modification of AHFO sensor

    The size of the sensor directly determines its application scope.Literature review indicated that the existing AHFO sensors were too bulky and long (Sayde et al., 2010; Cao et al., 2015, 2018), which caused great inconvenience to the calibration and use of the sensor.Hence, there is an urgent need for a kind of AHFO sensor which is suitable for both in situ and indoor calibrations.

    The length and diameter of the AHFO sensors modified by the authors are 40 mm and 4 mm,respectively,as shown in Fig.5a and b.Furthermore,the shell of the sensor is made of corundum,which is well-conductive and corrosion-resistant. Meanwhile, the ordinary heating wire (made by resistance heater) is replaced by a carbon fiber rod because of its marked thermal stability, i.e. its resistance value will not change much with rising temperature.

    3.3. Sensitive study of heating time

    Fig.5. Experimental equipment:(a)FBGS and soil sample;and(b)Cross-section of the FBGS (not to scale) and fiber interrogation.

    In the study of Sayde et al.(2010),it is necessary to heat the soil for 120 s with 20 W/m power to achieve an appropriate temperature rise for water content measurement.The heating time for the test conducted by Cao et al. (2015) even reached 30 min. Such a long heating time is likely to cause a considerable disturbance to water distribution in the soil sample(Sayde et al.,2010).Then,Cao et al. (2018) shortened the heating time to 300 s, but the effect of heating time on the test results was not discussed in their study.

    It will be much more helpful for application of AHFO sensors in practical engineering if the heating time could be shortened without reducing accuracy.Therefore,the influence of heating time on the accuracy of water content measurement will be discussed.

    4. Test technique and plan

    4.1. Soil samples and test equipment

    In this study,some laboratory tests were carried out to check the performance of AHFO sensors during soil water content measurement.The structure of sensors used in this study was similar to that in Cao et al. (2018) but with a shorter length. The cross-section of the sensor is shown in Fig.5b,and the sensor itself has a length of 40 mm and a diameter of 4 mm. The complete measurement system consists of: (1) a rheostat; (2) 12 V DC power supply; (3) an ammeter;(4)automatic clock switch;and(5)a fiber interrogation,type NZS-FBG-A01.

    The soil used has a plastic limit and a liquid limit of 15.2% and 30.7%, respectively. It can be classified as silty clay. The basic properties of soil samples are listed in Table 1.

    4.2. Test plan for laboratory calibration

    The soil samples used in this study were prepared with the following procedure:

    (1) First, soil samples with various water contents were prepared by a statistic compaction method. Each sample had a diameter of 10 cm and a height of 10 cm.

    (2) Then, soil samples were put into a thermostat where the environmental temperature was maintained at 25C to reduce the influence of heat convection on soil water content measurement.

    (3) Third, a short-term heating action was implemented by a short carbon fiber rod,which was around 4 cm in length and coated with a corundum shell. At the same time, the temperature increasing process of surrounding soil was monitored by the core of FBGS (Fig. 5). In this study, the heating power was maintained at 2.7 W/m and the total heating duration was 10e450 s.

    The influence range of AHFO sensors depends on the heating time, power, and soil thermal conductivity (Johansen,1977; Abu-Hamdeh, 2003; Mitchell and Soga, 2005). For a given condition,the influence range is supposed to only depend on the heating time(C?té and Konrad, 2005; Sayde et al.,2010). According to Fourier’s law, for a transient conduction phenomenon, the influence range should be linear to the root of time or the logarithm of time.

    The heating process of AHFO sensors in a large sample was monitored to check the influence range. As the soil thermal conductivity increased with decrease of water content, a soil sample with 10 cm in both diameter and height, and 3.4% of volumetric water content was used to explore the maximum influence range of the AHFO sensors.AHFO sensors were inserted into the sample,as shown in Fig. 6a, one in the center and the others in the radial direction with an interval of 1 cm. After 20 min of heating with a power of 2.7 W/m, only inner 4 sensors (within 4.5 cm) demonstrated some increment of temperature, as shown in Fig. 6bed.Referring to Fig.6e,if the heating time is within 30 s,the influence range of AHFO sensors should be within 1 cm. In this study, the heating time was shorter than 20 min, and the diameter of the sample used was 10 cm, thus the influence of the outer boundary could be ignored.

    5. Results

    Eleven samples with different water contents were tested. The recorded temperature duration curves are illustrated in Fig. 7. As shown in Fig.7,the temperature duration curve recorded by AHFO sensors during the heating process can be divided into two stages:rapid rising stage and tapering stage.This trend is consistent with the report in Figs. 4 and 9 in Cao et al. (2015).

    Table 1Physical properties of the tested silty clay.

    Fig.6. Determination of the maximum influence range of heating process:(a)Sensor installation plan;(b,c,d)Contour maps measured after 5 s,30 s and 20 min of heating time;and(e)Radial temperature distribution after heating for 20 min.The spacing of the measuring points is 1 cm,and the temperature of the position between the measuring points is obtained by interpolation of adjacent measuring points. The point where the temperature increases by 0.5 C is the point affected.

    Fig. 7. Temperature duration curves of silty clay measured by FBGS.

    5.1. Data analysis using temperature characteristic value method in Cao et al. (2018)

    Based on the test data in Fig.7,the relation betweenTtandwfor the silty clay tested can be determined by Eqs.(5)e(7),as shown in Fig. 8. Such relation can be applied to further laboratory/field measurements to infer the soil water content from the measured temperature characteristic valueTt, which was the average temperature during the period of 150e400 s in the heating process at the tapering stage in Fig. 7.

    5.2. Data analysis using integration data analysis method

    Referring to Fig. 9, the relation between temperature characteristic value (Ft) and volumetric water content (w) can be divided into two segments, i.e. a linear segment (Eq. (11)) and a curve segment (Eq. (10)). A possible reason can be that whenwis larger thant,the bulky water, which is mainly stored in the macropores,controls the thermal conductivity of soil;whereas whenwis lower thant,the bulky water has been drained. Therefore, only the capillary water, which is mainly stored in the micropores or intraaggregate pores,exists and controls the thermal conductivity of soil.

    Fig. 8. Calibrated relation between soil water content and temperature characteristic value.

    Fig. 9. Results by the integration data analysis method: (a) Heating time of 10 s, 20 s,30 s and 40 s;(b) Heating time of 50 s,100 s,150 s and 200 s;and (c)Heating time of 250 s, 300 s, 350 s and 400 s.

    Fig.10. Calibrated relation using the data recorded from the cooling process.

    As shown in Fig. 9, Eqs. 8e11 can characterize the relation betweenFtandwvery well, i.e. the fitted curves haveR2values larger than 0.97 except for curve at heating time of 10 s and 20 s in Fig.9a.That is to say,even if the heating time is 30 s,Eqs.(10)and(11)can still be well fitted to calculate the water content.Moreover,Eq. (10) with only two fitting parameters will have a better fit convergence than Eq. (7).

    5.3. Data analysis in cooling process

    As shown in Fig. 7, a sharp drop begins immediately after the heating stops,and the curve gradually tends to be stable after 350 s.In order to avoid the error caused by the sharp drop of the curve,the paper used the integral value of the entire cooling process(450e950 s) to calibrate the water content.

    Fig.10 reveals that the integration data analysis method is also valid for the data recorded from the cooling process.The turn point between the linear part and the nonlinear part occurs at water content of around 30%, which is similar to Fig. 9a. Nevertheless, it seems feasible to use the data recorded from the cooling process to calibrate theFt-wcurve,but this is not convenient compared with using the data recorded from the heating process. Moreover, thecorrelation coefficient (R2) from the cooling process is slightly less than that in the heating process (Fig. 9b and c).

    Table 2Comparison of the accuracy of original and integration methods.

    5.4. Comparison of the accuracy of two data analysis methods

    As listed in Table 2,the root mean square errors(RMSE)values of the integration method are no more than 2%.Meanwhile,RMSE of the original method(Eqs.(5)e(7)and Fig.7)is 3.4%.As lower RMSE indicates higher accuracy, the integration method has a better performance compared to the original method.

    As for calibrating the FBGS using the data from the cooling process, the RMSE value is around 1.7%, and the correlation coefficient is higher than 95%. Even though it takes a relatively long time for soil cooling-down,the method can still calibrate the water content effectively.

    Figs.11 and 12 show the results of the error of the data obtained after heating of 50e400 s. The RMSE values of the integration method are no larger than 1.5%compared with those of the original method (Table 2). Of all the data in Figs.11 and 12, the maximum error value does not exceed 2.9%. However, there are many inevitable factors that might affect test accuracy,e.g.soil density,particle shape,and pore distribution(C?té and Konrad,2005).Furthermore,with increase of heating time, the maximum error is stable with value within 2.7%, and the fluctuation trend of error is gradually reduced (Fig.12).

    Fig. 13 revels the comparison of the integration data analysis method (in the heating and cooling processes) and Cao et al.(2018)’s method. The maximum error of calibration result by Cao et al. (2018)’s method is around 8.98%, which is nearly three times the error by the integration data analysis method.As shown in Fig.13,whether in the heating or cooling process,the calibration value from the integration method is closely around the 1:1 line,suggesting that this method is a good way to calibrate the AHFO sensors.

    In fact,since the integration data analysis method was proposed by Sayde et al. (2010), the factor controlling the accuracy of the method was the accuracy of the fiber itself. Hence, they did not finalize the integration analysis method and offer good test results(Note:the minimum scan interval was 5 s at that time,and it is now less than 1 s).

    Fig. 11. Results of error analysis (comparison of data obtained after heating of 50e 200s).

    Fig.12. Results of error analysis (comparison of data obtained after heating of 250e 400 s).

    5.5. Verification of repeatability

    Parallel tests were used to check the repeatability of this technique.Three groups of soil samples with different volumetric water contents(3.6%,19.8%and 40.5%)were tested.The results are illustrated in Fig.14. The temperature errors are within0.15C.With the same data analysis technique, the maximum errors of measured water content for each sample are 3.4%,1.9% and 1.4%,respectively.

    One should be noted that the start time of heating process is found to have some influence on the test results. If the wavelength does not reach a stable value (the fluctuation of the wavelength is within 0.01 nm), a thermal equilibrium between FBGS and the surrounding soil may be not reached and the results may involve some error. Hence, in the integration data analysis method, only after the wavelength does not change, the heating process could be started and the equilibrium temperature is recorded asT0.During the integration,onlyTTT0should be integrated. Hence, the start time for the integration does not affect the result.

    Fig. 13. Results of error analysis (comparison between integration data analysis method (in heating and cooling processes) and Cao et al. (2018)’s method).

    Fig.14. Verification of the repeatability of the test:(a)w 3.6%, (b)w 19.8%, and(c) w 40.5%. Each soil sample was heated by AHFO sensor at a constant power of 2.7 W/m for 450 s and then stopped heating.Each sample was tested three times with an interval of 30 min. The explanation of “400 s -Test 1” is as follows: “400 s” means the integration in the heating curve from 0 to 400 s,and the “Test 1”indicates the first test.

    Fig.15. Integration analysis results of the data in Fig. 8 from Cao et al. (2015).

    6. Discussion

    6.1. Revisit of data reported in Cao et al. (2015)

    The data in Fig. 8 of Cao et al. (2015) were analyzed by the integration analysis method. As shown in Fig. 15, the integration data analysis method has good performance and can avoid the trouble for the rough division of the rapid rising stage and the tapering stage.

    6.2. Influence of fiber length

    The actively heated FBGSs used by Cao et al. (2015) had the length and diameter of 15 cm and 5 cm, respectively. The sensor length used in Cao et al.(2018)even reached 1 m.A long fiber can be easily broken upon bending or shearing, especially when it is buried in coarse soil or upon compaction. Besides, a long fiber sensor is hard to calibrate and lack of flexibility in actual engineering applications.

    In this study, the length and diameter of fiber sensors were reduced to 4 cm and 4 mm,respectively,and high-quality data were obtained as shown in Fig.7.That is to say,a short-length sensor can still be effective for this technique.If a multi-point measurement is in demand, the short sensor can be easily connected into long clusters by hot pyro-condensation pipes.Such clusters can achieve the same function as the long fiber, but avoid the broken problem and the difficulty of calibration.

    6.3. Influence of heating time

    Sayde et al.(2010)pointed out that a heat pulse with a long time might cause significant disturbance to soil samples, e.g. the soil water content migration. At the same time, the measurement results in Fig.14b and c also confirm this.As the number of repeated heating increases,the measured water content decreases.That is to say, a shorter heating time during the measurement can help to reduce the disturbance to the soil water content field and improve the accuracy of soil water content measurement.

    Using the integration data analysis method,w-Ftrelation with a heating time of 10e400 s is calculated,as compared in Fig.9.The results demonstrate that the relations betweenFtandwhave similarity when heating time varies. Even at the heating time of 20 s, an effective relation ofw-Ft can be obtained.

    Fig.16. Variations of fitting parameters involved in Eqs. (6)and (7) with heating time t: (a) a and c, and (b) k and b.

    Referring to Fig.9,the fitting parameters seem to depend on the heating time. As illustrated in Fig.16, the fitting parametersa,b,kandcvary with heating time.Obviously,the parametersa,bandktend to stabilize after the heating time exceeds 100 s except for the parameterb.

    Fig.17. A recommended technique for the application of AHFO sensor in field.

    That is to say, if the heating time is larger than 75 s, the fitting parameters will not be sensitive to the heating time and the reliability of calibration will be high.With this consideration,a heating time of 75 s is recommended but not always, because when rapid measurements are needed in practical engineering applications,the results of 30 s-heating are also acceptable in view ofR2and RMSE values from Fig. 9a and Table 2. Such a recommendation of heating time is much lower than that used by Cao et al.(2015)(25e 30 min)and Cao et al.(2018)(300 s).In fact,shortening the heating time can also benefit continuous measurement.

    6.4. Suggestions on engineering application of AHFO sensors

    The major limitation of AHFO sensor is that the relation between volumetric water content and thermal phenomenon varies with soil types. Hence, prior to AHFO sensor being used practically,calibration tests should be carried out to obtainw-Ftcurve.Another issue is that the test results can only reflect the average thermal properties of soil. The homogeneity of soil cannot be considered in this technique. Therefore, a representative soil sample in field should be used as a benchmark. A recommended technique for the application of AHFO sensor in field is illustrated in Fig.17.The corundum is hard and has poor ductility.Therefore,it is not suitable for measurement in gravel soil. If a soil with a higher content of gravel is encountered in actual engineering applications,some supplementary measures such as pre-drilling should be adopted to avoid the damage of corundum.

    For multi-point measurement outdoors, individual AHFO sensors should be connected into clusters, and the spacing between the sensors can be easily adjusted,as shown in Fig.17.The weatherresistant adhesive at the end of the flexible hot pyrocondensation pipe is impervious to water and can protect the heating wire of the sensor from rusting. However, the embedding of sensors is a technical problem that needs further exploration. At the construction site,any excavation and backfilling will disturb the intact soil. The dry density of the soil after backfilling must be different from the undisturbed soil,which will inevitably bring errors to the measurement.

    7. Summary

    By reviewing the effect of heating time on measuring the volumetric water content using the FBGS, the following improvements were made to the measurement method of water content by FBGSs:

    (1) An integration data analysis method and new fitting functions can be applied, instead of the original temperature characteristic value method proposed by Cao et al. (2015,2018).

    (2) A short-term heating pattern can be used instead of the long term heating pattern.

    (3) A short fiber can be used instead of the long fiber.

    Declaration of Competing Interest

    The authors wish to confirm that there are no known conflicts of interests associated with this publication and there has been no significant financial support for this work that could have influenced its outcome.

    Data availability statement

    Some or all data, models or codes that support the findings of this study are available from the corresponding author upon reasonable request.

    Acknowledgments

    This study is supported by the National Natural Science Foundation of China (Grant No.51979002).

    国产v大片淫在线免费观看| 日本一区二区免费在线视频| 午夜免费成人在线视频| 亚洲片人在线观看| 久久99热这里只有精品18| 亚洲激情在线av| 国产亚洲欧美在线一区二区| 国产高清有码在线观看视频 | 超碰成人久久| 夜夜躁狠狠躁天天躁| 99国产精品99久久久久| 国产精品亚洲美女久久久| 搡老熟女国产l中国老女人| 日韩中文字幕欧美一区二区| 女同久久另类99精品国产91| 身体一侧抽搐| 美女午夜性视频免费| 91麻豆精品激情在线观看国产| 国产免费av片在线观看野外av| 免费在线观看完整版高清| 亚洲国产精品sss在线观看| 韩国精品一区二区三区| 最好的美女福利视频网| 巨乳人妻的诱惑在线观看| 少妇被粗大的猛进出69影院| 日韩欧美三级三区| 亚洲精品粉嫩美女一区| 极品教师在线免费播放| 国产成人精品无人区| 成人手机av| 久久这里只有精品19| www国产在线视频色| 高清在线国产一区| 啦啦啦韩国在线观看视频| 国产伦在线观看视频一区| 在线观看免费视频日本深夜| 午夜福利高清视频| 一进一出抽搐gif免费好疼| 精品国产亚洲在线| 精品国产一区二区三区四区第35| 88av欧美| 久久久久久大精品| 久久久国产精品麻豆| 欧美激情高清一区二区三区| 黄色视频,在线免费观看| 欧美成人性av电影在线观看| 久久久久九九精品影院| 别揉我奶头~嗯~啊~动态视频| xxxwww97欧美| 亚洲国产精品久久男人天堂| 精品第一国产精品| 亚洲av中文字字幕乱码综合 | 精品久久久久久久毛片微露脸| 亚洲aⅴ乱码一区二区在线播放 | x7x7x7水蜜桃| 黑人欧美特级aaaaaa片| 每晚都被弄得嗷嗷叫到高潮| 91成年电影在线观看| 国产成人av教育| 看黄色毛片网站| 精品电影一区二区在线| av有码第一页| 一进一出抽搐动态| 又黄又粗又硬又大视频| 日韩欧美一区视频在线观看| 日韩av在线大香蕉| 最新美女视频免费是黄的| 女同久久另类99精品国产91| 变态另类丝袜制服| 两人在一起打扑克的视频| 久久久久免费精品人妻一区二区 | 1024手机看黄色片| 一区二区三区国产精品乱码| 亚洲成人精品中文字幕电影| 18禁美女被吸乳视频| 亚洲欧美精品综合久久99| 亚洲一区中文字幕在线| 国产欧美日韩一区二区三| 这个男人来自地球电影免费观看| 中文在线观看免费www的网站 | 欧美黄色片欧美黄色片| tocl精华| 丁香六月欧美| 成人三级黄色视频| 99riav亚洲国产免费| 亚洲欧洲精品一区二区精品久久久| 少妇 在线观看| 国产91精品成人一区二区三区| 日韩高清综合在线| 国产激情偷乱视频一区二区| 男女午夜视频在线观看| 日本 欧美在线| 制服人妻中文乱码| 老熟妇仑乱视频hdxx| 久久久国产精品麻豆| xxxwww97欧美| 欧美日韩瑟瑟在线播放| 一区福利在线观看| 国产欧美日韩一区二区精品| 亚洲一区二区三区不卡视频| 日韩精品免费视频一区二区三区| 1024视频免费在线观看| 亚洲精品av麻豆狂野| 久久欧美精品欧美久久欧美| 12—13女人毛片做爰片一| 女性生殖器流出的白浆| 色av中文字幕| 一级a爱视频在线免费观看| 亚洲激情在线av| 国产色视频综合| 久久青草综合色| 天天躁夜夜躁狠狠躁躁| 欧美激情高清一区二区三区| 黄色丝袜av网址大全| 亚洲精品在线美女| 女生性感内裤真人,穿戴方法视频| 非洲黑人性xxxx精品又粗又长| 久久精品国产综合久久久| 天天躁狠狠躁夜夜躁狠狠躁| 最近最新中文字幕大全免费视频| 91老司机精品| 国产在线观看jvid| 成年人黄色毛片网站| 免费观看精品视频网站| 亚洲成人精品中文字幕电影| 国产精华一区二区三区| 成人一区二区视频在线观看| 成人三级黄色视频| 18禁裸乳无遮挡免费网站照片 | 一区二区日韩欧美中文字幕| av福利片在线| 曰老女人黄片| 亚洲熟妇熟女久久| 91麻豆av在线| 12—13女人毛片做爰片一| 90打野战视频偷拍视频| 成人18禁在线播放| 此物有八面人人有两片| 日本a在线网址| 午夜成年电影在线免费观看| 欧美成人午夜精品| 国产在线观看jvid| 亚洲国产高清在线一区二区三 | ponron亚洲| 日韩精品免费视频一区二区三区| 极品教师在线免费播放| 亚洲五月婷婷丁香| 老司机午夜十八禁免费视频| 日韩欧美在线二视频| 热re99久久国产66热| 中文亚洲av片在线观看爽| 人人妻人人澡欧美一区二区| 国产一区二区在线av高清观看| 国产精品久久视频播放| 久久久国产成人精品二区| 国产激情欧美一区二区| 国产精品免费视频内射| 在线视频色国产色| 久久久国产精品麻豆| 久久午夜亚洲精品久久| 美女扒开内裤让男人捅视频| 国产精品亚洲美女久久久| 国产男靠女视频免费网站| 亚洲中文av在线| 黄片大片在线免费观看| 麻豆国产av国片精品| 国产又黄又爽又无遮挡在线| 国产在线精品亚洲第一网站| 亚洲最大成人中文| 亚洲美女黄片视频| 免费高清视频大片| 午夜福利在线观看吧| 热99re8久久精品国产| 麻豆一二三区av精品| 日本成人三级电影网站| 男女下面进入的视频免费午夜 | 波多野结衣高清作品| 免费高清在线观看日韩| 老鸭窝网址在线观看| 亚洲国产精品久久男人天堂| 一二三四在线观看免费中文在| 日韩 欧美 亚洲 中文字幕| 欧美日韩中文字幕国产精品一区二区三区| 在线观看免费日韩欧美大片| 99国产精品一区二区三区| 在线永久观看黄色视频| 18禁黄网站禁片午夜丰满| 欧美激情久久久久久爽电影| 亚洲熟妇熟女久久| svipshipincom国产片| 成人av一区二区三区在线看| 日本 欧美在线| 成人欧美大片| 欧美在线黄色| x7x7x7水蜜桃| 成人特级黄色片久久久久久久| 亚洲狠狠婷婷综合久久图片| 精品一区二区三区av网在线观看| 欧美在线黄色| 精品一区二区三区四区五区乱码| 亚洲人成网站在线播放欧美日韩| 亚洲av第一区精品v没综合| 日本一本二区三区精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产av又大| 日韩国内少妇激情av| 国产精品,欧美在线| 99热6这里只有精品| 手机成人av网站| 亚洲国产欧洲综合997久久, | 久久99热这里只有精品18| 国产极品粉嫩免费观看在线| 黄网站色视频无遮挡免费观看| 国产精品精品国产色婷婷| 色尼玛亚洲综合影院| 亚洲 国产 在线| 久久亚洲真实| av电影中文网址| 宅男免费午夜| 久久久久国产一级毛片高清牌| 国产99白浆流出| 免费女性裸体啪啪无遮挡网站| 最新美女视频免费是黄的| 午夜影院日韩av| 黄片大片在线免费观看| 精品福利观看| 日韩欧美免费精品| 看片在线看免费视频| 日韩精品免费视频一区二区三区| 国产国语露脸激情在线看| 国产乱人伦免费视频| 国产一区二区激情短视频| 国产亚洲欧美精品永久| 国产av在哪里看| 国产精华一区二区三区| 欧美精品亚洲一区二区| 少妇裸体淫交视频免费看高清 | 精品国产国语对白av| 国产91精品成人一区二区三区| 国产区一区二久久| 久久久久久九九精品二区国产 | 成在线人永久免费视频| 国产激情欧美一区二区| 岛国视频午夜一区免费看| 精品久久蜜臀av无| 国产精品久久久久久亚洲av鲁大| 男人操女人黄网站| 欧美成人午夜精品| av超薄肉色丝袜交足视频| 亚洲国产日韩欧美精品在线观看 | 亚洲中文日韩欧美视频| 亚洲成av片中文字幕在线观看| 精品午夜福利视频在线观看一区| 久久久久久人人人人人| 亚洲精华国产精华精| 久久精品国产综合久久久| 99热只有精品国产| 精品久久久久久久久久久久久 | 一级作爱视频免费观看| 女同久久另类99精品国产91| 一个人观看的视频www高清免费观看 | 国产高清激情床上av| 国内精品久久久久久久电影| 国产99白浆流出| 色播亚洲综合网| 国产亚洲精品综合一区在线观看 | 最近最新免费中文字幕在线| 欧美日韩瑟瑟在线播放| 国产伦人伦偷精品视频| 亚洲国产欧美网| 在线免费观看的www视频| 免费看日本二区| 日本免费一区二区三区高清不卡| 大香蕉久久成人网| 精品少妇一区二区三区视频日本电影| 免费女性裸体啪啪无遮挡网站| 人人澡人人妻人| 日韩欧美国产一区二区入口| 国内毛片毛片毛片毛片毛片| 99久久精品国产亚洲精品| 久久国产亚洲av麻豆专区| 国产成人系列免费观看| 国产区一区二久久| 天天躁夜夜躁狠狠躁躁| 黄片大片在线免费观看| 精品少妇一区二区三区视频日本电影| 亚洲国产欧美一区二区综合| 久久久久免费精品人妻一区二区 | 国产成人啪精品午夜网站| 女人高潮潮喷娇喘18禁视频| x7x7x7水蜜桃| 每晚都被弄得嗷嗷叫到高潮| 桃红色精品国产亚洲av| 国产视频一区二区在线看| 亚洲av熟女| 精品国产国语对白av| АⅤ资源中文在线天堂| 国产精品一区二区免费欧美| 色综合亚洲欧美另类图片| 中出人妻视频一区二区| 亚洲欧美日韩无卡精品| 亚洲成人久久爱视频| 亚洲精品在线美女| 色播亚洲综合网| 欧美一级a爱片免费观看看 | 亚洲九九香蕉| 欧美乱码精品一区二区三区| 国产一卡二卡三卡精品| 国产片内射在线| www日本黄色视频网| 两人在一起打扑克的视频| 欧美一区二区精品小视频在线| 亚洲五月婷婷丁香| 欧美三级亚洲精品| 久久精品91无色码中文字幕| 一二三四在线观看免费中文在| 亚洲精品中文字幕一二三四区| 91在线观看av| 99国产综合亚洲精品| 精品久久久久久久人妻蜜臀av| 亚洲,欧美精品.| 国产久久久一区二区三区| 精品人妻1区二区| 欧美黑人巨大hd| bbb黄色大片| 亚洲avbb在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产一卡二卡三卡精品| 麻豆成人午夜福利视频| 欧美日本视频| 中文字幕人妻丝袜一区二区| 我的亚洲天堂| 日本精品一区二区三区蜜桃| 50天的宝宝边吃奶边哭怎么回事| 中文字幕人妻丝袜一区二区| 国产aⅴ精品一区二区三区波| av欧美777| 亚洲欧美精品综合久久99| 久久人人精品亚洲av| 人人妻人人澡人人看| 少妇粗大呻吟视频| 后天国语完整版免费观看| 欧美成人午夜精品| 久久九九热精品免费| 一级黄色大片毛片| 桃红色精品国产亚洲av| 久久精品91蜜桃| 十八禁人妻一区二区| 日本成人三级电影网站| 在线观看午夜福利视频| 熟女电影av网| 国产亚洲av嫩草精品影院| 女同久久另类99精品国产91| 男人的好看免费观看在线视频 | 婷婷精品国产亚洲av在线| 免费在线观看黄色视频的| 不卡一级毛片| 99在线视频只有这里精品首页| 不卡一级毛片| 日韩有码中文字幕| 免费在线观看影片大全网站| 精品熟女少妇八av免费久了| 久久精品人妻少妇| 老司机午夜十八禁免费视频| 午夜福利18| 首页视频小说图片口味搜索| 精华霜和精华液先用哪个| 国产精品亚洲av一区麻豆| 午夜福利18| 亚洲成国产人片在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 色av中文字幕| 在线永久观看黄色视频| 亚洲人成网站高清观看| 19禁男女啪啪无遮挡网站| 窝窝影院91人妻| 欧美激情极品国产一区二区三区| 午夜日韩欧美国产| 欧美日韩乱码在线| 欧美不卡视频在线免费观看 | svipshipincom国产片| 午夜福利一区二区在线看| 亚洲第一电影网av| 午夜福利一区二区在线看| 丰满人妻熟妇乱又伦精品不卡| 国产高清激情床上av| 欧美国产日韩亚洲一区| 99热只有精品国产| 成人欧美大片| 欧美性猛交╳xxx乱大交人| 国产亚洲精品一区二区www| 丝袜美腿诱惑在线| 欧美性长视频在线观看| a级毛片在线看网站| 麻豆成人av在线观看| 国产成人影院久久av| 久久精品91蜜桃| 好男人在线观看高清免费视频 | 亚洲第一欧美日韩一区二区三区| avwww免费| √禁漫天堂资源中文www| 日本 av在线| 亚洲精品国产一区二区精华液| 一区二区日韩欧美中文字幕| 成人18禁高潮啪啪吃奶动态图| 久久中文看片网| 国产成人av激情在线播放| 国产av不卡久久| 99国产精品99久久久久| 中文字幕另类日韩欧美亚洲嫩草| 成人国产综合亚洲| 俺也久久电影网| e午夜精品久久久久久久| 午夜影院日韩av| 欧美一级毛片孕妇| 国产蜜桃级精品一区二区三区| 一本精品99久久精品77| 啦啦啦 在线观看视频| 一级毛片高清免费大全| 精品久久久久久久毛片微露脸| 成人永久免费在线观看视频| 男女视频在线观看网站免费 | 久久99热这里只有精品18| 制服人妻中文乱码| 国产精品九九99| 亚洲午夜理论影院| а√天堂www在线а√下载| 亚洲午夜精品一区,二区,三区| 啪啪无遮挡十八禁网站| 久久精品91蜜桃| 欧美日韩福利视频一区二区| 亚洲精品一区av在线观看| 成年女人毛片免费观看观看9| 国产不卡一卡二| 色尼玛亚洲综合影院| 亚洲av美国av| 久久国产亚洲av麻豆专区| av在线天堂中文字幕| 超碰成人久久| 久久久久国内视频| 久热爱精品视频在线9| 亚洲七黄色美女视频| 美女 人体艺术 gogo| 久久久久久久久中文| 色综合站精品国产| 日日爽夜夜爽网站| 哪里可以看免费的av片| 一级毛片女人18水好多| 久久婷婷人人爽人人干人人爱| 国产精品亚洲av一区麻豆| 最好的美女福利视频网| 免费在线观看黄色视频的| 国产精品免费视频内射| 18禁美女被吸乳视频| 中文字幕精品免费在线观看视频| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 亚洲精品国产区一区二| 在线观看免费日韩欧美大片| 国产精品1区2区在线观看.| 亚洲精品中文字幕一二三四区| 午夜日韩欧美国产| 久久中文看片网| 99热这里只有精品一区 | 黄片大片在线免费观看| 国产精品一区二区精品视频观看| 狠狠狠狠99中文字幕| 中文资源天堂在线| 国产亚洲欧美在线一区二区| 精品国产乱子伦一区二区三区| 在线永久观看黄色视频| 亚洲精品国产一区二区精华液| 美女大奶头视频| 十分钟在线观看高清视频www| 男人舔奶头视频| 一进一出抽搐动态| 亚洲国产欧美日韩在线播放| 高清毛片免费观看视频网站| 国产精品1区2区在线观看.| 欧美zozozo另类| 国产精品乱码一区二三区的特点| 精品国产乱码久久久久久男人| 9191精品国产免费久久| 国产精品日韩av在线免费观看| 亚洲 欧美一区二区三区| 亚洲av电影在线进入| 又黄又粗又硬又大视频| 国产成人影院久久av| 91字幕亚洲| 久久草成人影院| 国产一卡二卡三卡精品| 色播在线永久视频| a级毛片在线看网站| 午夜视频精品福利| 色综合欧美亚洲国产小说| 国产私拍福利视频在线观看| 色播在线永久视频| 精品久久久久久久末码| 国产一级毛片七仙女欲春2 | 12—13女人毛片做爰片一| 欧美人与性动交α欧美精品济南到| 一区二区三区高清视频在线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产在线观看jvid| 人人澡人人妻人| 极品教师在线免费播放| 国产亚洲av高清不卡| 午夜老司机福利片| 黄色女人牲交| 人人妻人人看人人澡| 99久久无色码亚洲精品果冻| 午夜福利高清视频| 精品久久久久久成人av| 亚洲精品国产精品久久久不卡| 超碰成人久久| 一本一本综合久久| 97人妻精品一区二区三区麻豆 | 中文字幕人成人乱码亚洲影| 50天的宝宝边吃奶边哭怎么回事| 午夜亚洲福利在线播放| 国产99白浆流出| 久久久国产成人精品二区| 亚洲第一青青草原| 国产熟女xx| 一边摸一边抽搐一进一小说| 午夜福利在线在线| 国产人伦9x9x在线观看| 成在线人永久免费视频| 国产亚洲精品久久久久久毛片| 很黄的视频免费| 在线永久观看黄色视频| 99久久99久久久精品蜜桃| 午夜福利在线观看吧| www.自偷自拍.com| 黄色丝袜av网址大全| 国产一区二区三区在线臀色熟女| 麻豆成人av在线观看| 免费在线观看成人毛片| 色综合亚洲欧美另类图片| 国产亚洲精品综合一区在线观看 | 久久午夜综合久久蜜桃| 久久中文字幕一级| 妹子高潮喷水视频| 成年女人毛片免费观看观看9| 国产精品二区激情视频| 亚洲最大成人中文| 中文字幕人妻熟女乱码| 亚洲片人在线观看| 男男h啪啪无遮挡| 麻豆av在线久日| av中文乱码字幕在线| 成人18禁高潮啪啪吃奶动态图| 叶爱在线成人免费视频播放| 精品国产超薄肉色丝袜足j| 少妇的丰满在线观看| 亚洲精品美女久久久久99蜜臀| 欧美绝顶高潮抽搐喷水| 午夜精品在线福利| 搡老岳熟女国产| 啦啦啦韩国在线观看视频| 人成视频在线观看免费观看| 欧美一区二区精品小视频在线| 欧美中文综合在线视频| 精品久久久久久久毛片微露脸| 国产视频一区二区在线看| 桃色一区二区三区在线观看| 男女床上黄色一级片免费看| 女性被躁到高潮视频| 亚洲成人久久爱视频| 国产精品久久久人人做人人爽| 9191精品国产免费久久| 久久性视频一级片| 美女免费视频网站| 黄色毛片三级朝国网站| 中国美女看黄片| 亚洲美女黄片视频| 中文字幕久久专区| 母亲3免费完整高清在线观看| 久久欧美精品欧美久久欧美| 最近最新中文字幕大全电影3 | 热re99久久国产66热| 99riav亚洲国产免费| 波多野结衣av一区二区av| 黄色a级毛片大全视频| 成人午夜高清在线视频 | 美女扒开内裤让男人捅视频| 90打野战视频偷拍视频| 亚洲 欧美 日韩 在线 免费| 精品欧美一区二区三区在线| 国产三级黄色录像| 欧美不卡视频在线免费观看 | 好看av亚洲va欧美ⅴa在| tocl精华| 午夜福利欧美成人| 国产亚洲精品综合一区在线观看 | 精品国产乱子伦一区二区三区| cao死你这个sao货| 国产精品久久久久久亚洲av鲁大| 九色国产91popny在线| 最近最新中文字幕大全免费视频| 自线自在国产av| 男女那种视频在线观看| cao死你这个sao货| 亚洲真实伦在线观看| 啦啦啦观看免费观看视频高清| 日本 欧美在线| 久久国产乱子伦精品免费另类| 精品国产乱子伦一区二区三区| av超薄肉色丝袜交足视频| 午夜福利18| 黄色毛片三级朝国网站| 十八禁人妻一区二区| 亚洲国产看品久久| netflix在线观看网站| 国产欧美日韩一区二区精品|