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

    A distributed measurement method for in-situ soil moisture content by using carbon-fiber heated cable

    2015-10-09 07:10:03DingfengCoBinShiHonghuZhuGungqingWeiShenEnChenJunfnYn
    關(guān)鍵詞:機(jī)加成株吡蚜

    Dingfeng Co,Bin Shi,*,Honghu Zhu,*,Gungqing Wei,Shen-En Chen,Junfn Yn

    aSchool of Earth Sciences and Engineering,Nanjing University,Nanjing,210023,China

    bSuzhou NanZee Sensing Technology Co.,Ltd.,Suzhou,215123,China

    cDepartment of Civil and Environmental Engineering,University of North Carolina at Charlotte,NC 28223,USA

    A distributed measurement method for in-situ soil moisture content by using carbon-fiber heated cable

    Dingfeng Caoa,Bin Shia,*,Honghu Zhua,*,Guangqing Weib,Shen-En Chenc,Junfan Yana

    aSchool of Earth Sciences and Engineering,Nanjing University,Nanjing,210023,China

    bSuzhou NanZee Sensing Technology Co.,Ltd.,Suzhou,215123,China

    cDepartment of Civil and Environmental Engineering,University of North Carolina at Charlotte,NC 28223,USA

    A R T I C L E I N F O

    Article history:

    in revised form

    20 August 2015

    Accepted 24 August 2015

    Available online xxx

    In-situ soil moisture content

    Distributed measurement

    Carbon-fiber heated cable(CFHC)

    Fiber-optic sensing

    A B S T R A C T

    Moisture content is a fundamental physical index that quantifies soil property and is closely associated with the hydrological,ecological and engineering behaviors of soil.To measure in-situ soil moisture contents,a distributed measurement system for in-situ soil moisture content(SM-DTS)is introduced. The system is based on carbon-fiber heated cable(CFHC)technology that has been developed to enhance the measuring accuracy of in-situ soil moisture content.Using CFHC technique,a temperature characteristic value(Tt)can be defined from temperature-time curves.A relationship among Tt,soil thermal impedance coefficient and soil moisture content is then established in laboratory.The feasibility of the SM-DTS technology to provide distributed measurements of in-situ soil moisture content is verified through field tests.The research reported herein indicates that the proposed SM-DTS is capable of measuring in-situ soil moisture content over long distances and large areas.

    ?2015 Institute of Rock and Soil Mechanics,Chinese Academy of Sciences.Production and hosting by

    Elsevier B.V.All rights reserved.

    1.Introduction

    Soil moisture content is an important and fundamental physical index that correlates to the hydrological,ecological and engineering properties of soil(Maleki and Bayat,2012;Mittelbach et al., 2012;Cho and Choi,2014).To obtain accurate and reliable in-situ soil moisture content is imperative for solving many engineering problems.In laboratory,the most direct and common method for measuring soil moisture content is the drying method.Obviously, this method is not directly applicable to field conditions,where the capability of measuring over large areas and long distances is an essential requirement(Schmugge et al.,1980;Lunt et al.,2005;Yin et al.,2013).

    There are also several in-situ methods for measuring soil moisture content that have been suggested,including the time domain reflectometry(TDR)method,the capacitance and frequency domain reflectometry(FDR)method,the ground penetratingradar(GPR)method,andremotesensingmethods(Dobriyal et al.,2012).Each proposed method has its own advantages,disadvantages and limiting scope.All the abovementioned methods,however,do share a common weakness, that they provide only point-mode measurements and cannot satisfy the real-time measuring requirement of in-situ soil moisture content(Dobriyal et al.,2012;Susha Lekshmi et al.,2014;Sun et al.,2015).Although GPR and remote sensing methods can obtain soil moisture content distribution on ground surface over an area,their measuring accuracy is poor due to environmental factors,such as vegetation,soil surface temperature,surface reflectivity,and thus the measuring depth is limited to the soil surface(Doolittle and Collins,1995;Weihermüller et al.,2007;Schrott and Sass,2008).

    Distributed fiber-optic sensing(DFOS)technology for condition monitoring of engineering systems has been rapidly developed in recent years.Compared with other sensing technologies, DFOS technology is superior in that it is capable of distributed and long-distance measurements,characterized by anti-corrosion and anti-jamming.One of the DFOS technologies is the Raman optical time domain reflectometry(ROTDR)technique which has extensive applications to integrity monitoring of electric power structures,tunnels and oil wells.ROTDR has also been applied as fire warning systems and for seepage monitoring in dams.For example,C?té et al.(2007)used optical fibers to detect water leakage at the Peribonka dam,and Vogt et al.(2010)estimated the seepage rates in a losing stream by means of a fiber-optic,highresolution vertical temperature profiling technique.In applicationof moisture monitoring,some distributed technology and equipment were proposed.These technologies include two types:passive fiber optics(PFO)method and actively heated fiber optics(AHFO)method.For the PFO,there is no external power to heat cable,and soil moisture content is determined by sensing the natural temperature of soil.While for AHFO,an external power is needed to heat cable and produce thermal pulses on soil,and the moisture content is inferred by the spreading characteristics of these thermal pulses.

    Because the PFO does not need the heat source,it is an energysaving and intrinsic-safe method.Steele-Dunne et al.(2010)demonstrated the feasibility of using buried fiber-optic cables as a distributed measurement technique to determine the thermal response of soil moisture to diurnal temperature cycles.Although the PFO is a low-cost method,its precision can be affected bya large number of factors,e.g.temperature difference,plants cover, weather and seasons.

    By using an actively heated fiber-optic cable,Weiss(2003)demonstrated the potential for collecting distributed soil moisture content data but with errors>6%.A laboratory test conducted by Sayde et al.(2010)showed that it is possible to use ROTDR system for measurements of active soil moisture content with high precision and short measurement time.They fitted an exponential relation between soil moisture and cumulative temperature increase.Striegl and Loheide(2012)used the similar AHFO method in the field to measure soil moisture content.Ciocca et al.(2012)systematically introduced the basic principles of the AHFO method in heating and cooling phases.Gil-Rodríguez et al.(2012)used the AHFO method to estimate wetting bulbs and water distribution under drip emitters in the laboratory.Benítez-Buelga et al.(2014)obtained volumetric heat capacity and water content by using dual-probe approach.Read et al.(2014)quantified vertical flow in boreholes by the AHFO technology.Sayde et al.(2014)firstly used the AHFO method in a field across 1-1000 m scales.

    Although these previous experimental results have demonstrated the feasibility of applying optical fibers for distributed measurement of soil moisture content,further studies are needed to improve the technique including possibly integrating laboratory and field measurements into more efficient methods or developing hybrid cable configurations to optimize heating schemes,etc.(Sayde et al.,2010;Striegl and Loheide,2012).This is because the thermal property of hybrid cable directly influences the measurement accuracy of soil moisture content,and the heating mode or power supply determines the measurement distance.A perfect heated cable should have uniformity in its axial direction,i.e.the electric resistance should be identical for each section with the same length.Furthermore,the calibration procedure proposed by Sayde et al.(2010)is very complicated and the relationship between soil thermal conductivity and moisture content has not been established.A large number of functions were used to fit the relation between soil thermal conductivity and moisture content,and the most one appears in a form of exponential function,however, all of these are empirical or semi-empirical equations which cannot be obtained from theoretical analytical solution(Chung and Horton,1987;C?té and Konrad,2005;Lu et al.,2007;Jougnot and Revil,2010).

    A distributed measurement system for in-situ soil moisture content(SM-DTS)is introduced in this paper.This work aims at promoting the application of optical fiber sensing to soil moisture content in the field,simplifying the calibration procedure,and obtaining a more maneuverable and effective relation between thermal conductivity and moisture content.The following sections will describe the fundamental principles of SM-DTS,the structure of the carbon-fiber heated cable(CFHC),the definition of the temperature characteristic value(Tt)and the relationship among Tt, soil thermal impedance coefficient and soil moisture.

    2.Fundamental principle of SM-DTS

    Distributed temperature sensing(DTS)technology based on the ROTDR was first proposed by Hartog(1983)that when a pumped laser pulse with a certain amount of energy is launched into a sensing optical fiber,a collision occurs between fiber molecules and photons,which yields spontaneous Raman scattering. 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/ Stokes intensity ratio and the temperature of the optical fiber, the distributed measurement of temperature along the sensing fiber can be realized(Tyler et al.,2009).The most prominent advantage of the Raman scattering spectrum is that it is only sensitive to temperature and single-end measurement,which is very fitting for geotechnical engineering monitoring where measurements are usually deep in the earth.Modern DTS systems have relatively high temperature measurement accuracies at ±0.1°C,and are capable of long measuring distances at dozens of kilometers.In this paper,an attempt is made to adopt DTS technology to measure in-situ soil moisture content by using the SMDTS technology.The underlying principle of the SM-DTS technology is based on the assumption that a straightforward relationship between the temperature characteristic value(Tt)and soil moisture content can be established,which can be used to quantify in-situ soil moisture content using a heated cable buried in the tested soil and a DTS sensing unit.

    2.1.The carbon-fiber heated cable

    In order to enhance the measurement accuracy of the soil moisture content by a DTS unit,a CFHC is suggested and is developed for field applications.Fig.1 shows a cross-sectional view of the CFHC with 4 mm in diameter,which consists of a multimode optical fiber,a carbon-fiber and a carbon-fiber jacket. The optical fiber contains its own core,cladding,coating and fiberoptic jacket.The optical fiber is then encased by 24,000 carbonfibers.

    Fig.1.Cross-sectional view of the CFHC(not to scale).

    The CFHC is heated by the conduction of the carbon-fibers,and the heating principle is shown in Fig.2.The carbon-fiber is connected to several close-loop electric circuits consisting of two copper wires that are electrically collocated and connected to the carbon-fibers through annular binding posts.The annular binding post is installed with a distance of D.One end of the copper wire is connected with the annular binding post,and the other end is connected to a power supply.The heating power of each section of the carbon-fiber between two adjacent annular binding posts should be identical.It is assumed that the carbon-fiber with thesame length will have an equivalent resistance of R.The equivalent circuit diagram is shown in Fig.2a.

    Fig.2.Diagram for the heating principle of the CFHC.(a)Equivalent circuit diagram of CFHC,and(b)Circuit connection style of CFHC.

    In both the laboratory and field tests,the heating power(Q)is set as 2.74 W/m,and the electric current flowing through every equivalent resistance(R)is determined to be 0.378 A.Heat spreading characteristic of the CFHC afterconducting heat is shown in Fig.3 where Q1is the energy to heat CFHC itself,and Q2is the energy that has been spread into surrounding soil.The sum of Q1and Q2is the heating power(Q).Since soils with different moisture contents will have disparate thermal impedance coefficients, consequently,the heating temperature-time curves of soils with different moisture contents around the CFCH can also be different.

    A typical heating temperature-time curve for soil is shown in Fig.4.According to the heating rate,the curve can be divided into two stages.The first stage involves temperature rapidly rising within the soil when the conducting carbon-fiber starts to heat up;and the second stage involves the temperature rising rate slowing down(after an initial heating period),eventually becoming plateaued(C?té et al.,2007).

    Fig.3.Diagram for the principle of the soil moisture content measurement technique using SM-DTS.

    2.2.Temperature characteristic value

    Sayde et al.(2010)have already proposed a temperature variable,cumulative temperature increase,Tcum,which is defined as the temperature integration during the total time.They have established an empirical and complicated relation between Tcumand soil moisture content.Even though the Tcumhas many advantages,it is not defined scrupulously,because the temperature is recorded by DTS discretely,and we cannot calculate the integration of these discrete data directly.In order to establish a simple but effective relationship between the temperature distribution along the CFHC,and to overcome the weakness of Tcum,a temperature characteristic value(Tt)is proposed,which is defined as the arithmetic mean value of the temperatures measured by the DTS within a defined characteristic time interval of[ta,tb].The temperatures at taand tbare defined as Taand Tb,respectively in Fig.4.The time interval is selected based on the temperaturetime curve.Ttis calculated as

    where Tiis the temperature value measured by DTS within a characteristic time interval[ta,tb]and n is the number of measurements during the time interval.On the temperature-time curve,the characteristic time interval is selected based on the tapering stage of the time-temperature curve.Ttis then calculated as shown in Fig.4.

    The factors that affect Ttmay come from two aspects.The first one is the intrinsic property of soil itself,such as mineral composition,grain size and dry density;and the second one comes from external elements,such as moisture and environmental temperature(Wierenga et al.,1969;Abu-Hamdeh and Reeder,2000;C?té and Konrad,2005;Markle et al.,2006;Abbasy et al.,2014).For a specific in-situ soil,the mineral composition and structural performance typically remain the same.The effect of ambient temperature can be easily compensated by sensor calibration.So the main factor influencing the soil thermal impedance coefficient and Ttis only the soil moisture content.Therefore,once the relationshipbetween soil moisture content and Ttis established,the soil moisture content can be easily measured by Tt.

    Fig.4.Temperature-time curve of CFHC during the heating process.

    2.3.Relationship between Ttand soil moisture content

    Assuming that the soil is uniform and isotropic,for a line thermal source such as CFHC buried in a boundless soil medium,the temperature(T)obtained byDTS satisfies the following relationship(C?té et al.,2007):

    where qLis the heat source strength per unit length(W/m);ρis the thermal impedance coefficient((m K)/W);T0is the initial temperature of the line thermal source;D is a constant for a specific soil,which is independent of time but a function of the medium thermal diffusivity.Eq.(2)can be combined into Eq.(1)and result in:

    Once the characteristic time interval is selected,ln tbe

    i comes a constant.Then Eq.(3)can be transformed into the following form:

    Tt=Aρ+B(4)

    where

    Both A and B are values to be determined in laboratory or during field tests.The soil thermal impedance coefficient depends on the thermal impedance of solid particles,dry density and moisture content,but for a specific soil,the solid particles and dry density change slightly with time.Therefore,the factor affecting Ttis only the moisture content.According to the results obtained by Sayde et al.(2010),the relation between Ttand moisture content is nonlinear,especially when soil is very dry,but for a large range of soil moisture content,this relation can be approximately regarded as linear:

    ρ=Cω+D(ω1<ω<ω2)(6)

    whereωis the soil moisture content,ω1is the lower limit of moisture content for the application scope of Eq.(6),andω2is its upper limit.After Eq.(6)is simplified,the relationship between Ttandωcan be written as

    Tt=Eω+F(7)

    where E=AC,F=AD+B.

    According to Eq.(7),it can be obtained that:

    ω=kTt+b(8)

    where k=1/E and b=-F/E.Both k and b are constants that can be determined through laboratory tests.According to Eq.(8),the soil moisture content can be determined through Ttmeasurements.

    3.SM-DTS design

    Based on the description of the SM-DTS principle in previous section,a distributed measurement system for soil moisture content is devised which consists of three subsystems:a heatingsensing unit,a demodulating unit and a data processing unit,as shown in Fig.5.

    The heating-sensing subsystem consists of the carbon-fiber,the sensing optical fiber,the copper wires,the power supply,the annular binding posts and the cable coating,all bundled together. The function of the subsystem is both sensing and transmitting of temperature along the CFHC.The demodulating subsystem includes the DTS demodulation instrument and associated software. Its main function is to transform the light signal in the optical fiber sensor into temperature signal.The data processing subsystem includes both the data processing instrument and associated software,which is used to analyze the data measured by the DTS,to generate the temperature-time curves,and to determine Tt.The soil moisture content can then be calculated from the computed Tt.

    As shown in Eq.(8),in the field measuring process,it is not necessary to investigate the basic physical indices such as composition,dry density and grain size as long as the coefficients k and b are calibrated before measurement.When SM-DTS is applied in field for a long time,the temperature sensed by CFHC can be affected by environmental factors,such as climate,season and weather.To avoid this adverse influence,as a rule of thumb,the buried depth of CFHC should be more than 0.2 m.

    4.Laboratory test

    In order to verify the feasibility of SM-DTS in measuring soil moisture content,calibration tests have been performed in the laboratory.

    鉤藤蚜蟲可用0.3%苦參堿水劑200~300倍,或0.6%乙基多殺菌素懸浮劑1 000~1 500倍,或35%吡蟲啉懸浮劑10 000~13 000倍,或25%吡蚜酮懸浮劑2 000~3 000倍在苗期噴霧,或在成株期密集田塊用煙霧機(jī)加化學(xué)藥劑產(chǎn)煙防控。

    4.1.Testing scheme

    Because the spatial resolution of DTS is only 1 m,so,to improve the spatial resolution,a CFHC with 5 m in length is coiled onto a polyvinyl chloride(PVC)tube of 5 cm in diameter,forming a measuring tube,as shown in Fig.6.

    The spatial resolution of the measuring tube can then be calculated as

    Fig.5.Schematic diagram of the SM-DTS including heating-sensing subsystem,demodulating subsystem and data processing subsystem.

    Fig.6.Forming of the measuring tube.

    where M is the spatial resolution of the measuring tube,φis the diameter of CFHC(4 mm in this study),P is the diameter of the measuring tube(5 cm in this study),and S is the spatial resolution of the DTS(1 m in this case).From Eq.(9),M is computed to be 2.5 cm,showing a significant improvement.

    The laboratory test involved the measurements of clayey soil with different moisture contents.Before measuring,the clayey soil was filled layerwise into an iron drum of 40 cm in diameter and 46 cm in height,and the measuring tube was buried into the center of the iron drum,as shown in Fig.7.The chemical composition of soil is described in Table 1.The soil was prepared into six groups with moisture contents of 10%,15%,25%,30%and 35%,respectively, and Ttwas measured by the SM-DTS.

    Fig.7.Basic configuration of laboratory tests.

    4.2.Results and discussion

    Fig.8 shows the temperature-time curves for the soil samples with different moisture contents.During moisture content measurement,a nonheating,temperature sensing optical fiber was pasted on the inner wall of the iron drum to measure the boundary temperaturewithin the iron drum.Notemperatureperturbation on the inner wall was found;hence,the boundary effect of the iron drum can be neglected.The temperature rising time is in the range of 0-23 min,and the temperature will taper at around the time range of 23-30 min.On the temperature tapering stage,the interval of 26-30 min is selected as the characteristic time interval[ta,tb]to calculate Tt.The time of the selected interval is 4 min and n is 5.

    Table 1 Chemical composition of the test soil.

    As can be seen in Fig.8,the relationship between T and time t appears to be a logarithmic function,which is consistent with the theoretical analysis(C?té et al.,2007).It is found that the deviations between any two neighboring curves are approximately consistent except for that between the 10%and 15%curves.This phenomenon can be established as an optimal proportion of water absorption of 10%moisture content within the soil.When the absorbed water is heated,it did not move freely like free water through molecular movement to spread energy.As a result,the thermal impedance ishigher,so the corresponding T is higher(Luo and Li,1990;Barry-Macaulay et al.,2013).

    Fig.8.Temperature-time curves of the measuring tube buried in soils with different moisture contents for 30 min heating under 2.74 W/m power.

    Fig.9.Relation between Ttand moisture content,ω,with 2.74 W/m heat power.Solid line represents the linear fitting result for soil samples except the sample with 10% moisture content,and the R2of fitting function is 0.99.Dashed line represents the linear fitting result for all soil samples,and the R2of fitting function is 0.95.

    It can be seen from Fig.9 that there is a linear relationship betweenωand Ttobtained from the soils with moisture contents of 15%,20%,25%,30%and 35%.The coefficient of determination(R2)is 0.99,and root mean squared error(RMSE)is 0.25°C,shown as the solid line in Fig.9.If we consider the Ttobtained from 10%moisture soil into the fitting process,R2will decrease to 0.95,and RMSE will increase to 0.38°C,shown as the dashed line in Fig.9.Thus,the relationship betweenωand Ttcan be express as

    Tt=-0.13ω+13.33(10)

    Eq.(10)can be expressed as Eq.(11)by some linear transformation steps:

    ω=-7.69Tt+103.53(11)

    On the basis of Eq.(11),the soil moisture content can be determined with a determined Tt.Since actual soil composition and dry density have influence on Tt,it is necessary to calibrate the coefficients k and b.

    5.Field test

    The following field test is conducted to verify the effectiveness of SM-DTS in the field.

    5.1.Schematic design of field test

    The field test was conducted in a 19 m long trench with silty clay,and its mineral composition is the same as that used in the laboratory tests.The CFHC was buried in the soil directly instead of the measuring tube.Because the field is 19 m in length,1 m spatial resolution DTS can meet the measurement requests.The heating power in the field was kept the same as that in the laboratory tests. Fig.10 shows the soil trench with a cross-section of 20 cm×20 cm. The CFHC was paved at the bottom of the trench,then backfilled and compacted.The trench from A to B was lengthwise divided into three sections where 1-7 m and 13-19 m were with the same moisture content and 7-13 m with less moisture content.After the measurement was made using the SM-DTS,the soil moisture contents of each sectionweremeasured using drying method in the laboratory.

    5.2.Results and discussion

    Comparison of soil moisture content measured by drying method(ωD)and SM-DTS(ωS)is shown in Fig.11.It can be seen from Fig.11 that the variation ofωS(by SM-DTS)agrees with that of ωD(by drying method)and both reflect the soil moisture content distribution from A to B.IfωDis regarded as the precise value,the absolute errors(AE)ofωSfrom 19 samples collected in every meter are shown in Table 2.

    Fig.10.Diagram of the soil trench with a 19 m length and a 20 cm×20 cm crosssection in the field test.

    As shown inTable 2,the absolute errors of the 16 samples are all less than 3%,and the minimum is 0.01%.The other absolute errors are slightly higher,ranging from 2.97%to 6.57%.These relativelyhigh errors may be caused by variations in soil composition and properties including root of plants,non-uniformity of dry densities, etc.Fig.12 shows the error analysis of the results measured by SMDTS as compared toωD.

    Fig.11.Comparison of soil moisture content measured by drying method(ωD)and SMDTS(ωS)in a 19 m long field test.

    Table 2 Absolute errors between SM-DTS and drying method.

    Fig.12.Error analysis of the results measured by SM-DTS.

    As can be seen from Fig.12,the data measured by the two methods have a relevance with R2=0.72 and RMSE=2.85%. Comparing with the measurement accuracy of other methods proposedbyotherresearchers(e.g.Calamitaetal.,2012;Rowlandson et al.,2013;Yin et al.,2013),current study has higher relevance.The measuring result indicates that the SM-DTS has good measurement accuracy for the in-situ soil moisture content.

    6.Conclusions

    In conclusion,the following observations can be made about the SM-DTS method suggested in this paper:

    (1)The CFHC is made with embedded sensing optical fiber and is shown to enhance measurement accuracy of soil moisture content. Temperature-time curves for soil moisture content correlations have been obtained in the laboratory.

    (2)A temperature characteristic value(Tt)is proposed,and the relationship among Tt,soil thermal impedance and moisture content is established.Based on this relationship,the soil moisture content can be measured by the proposed DTS system.

    (3)A distributed soil moisture measurement system(SM-DTS)is designed and made,which includes three subsystems:heatingsensing,demodulating and data processing.

    (4)The parameters k and b of the SM-DTS are essential in soil moisture content measurement and can be determined in the laboratory.

    The field results show that the SM-DTS method is viable for distributed measurements of in-situ soil moisture content,which has better results than previously proposed methods.For future studies,it is important to theoretically analyze the influences of intrinsic factors such as soil density and soil composition,and to establish relationship between temperature and soil moisture content by Richard's equation of soil moisture content and Fick's law.

    Conflict of interest

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

    Acknowledgments

    The authors wouldlike tothank all participants in the laboratory and field tests.The financial supports provided by the National Natural Science Foundation of China(Grant Nos.41230636, 41372265,41427801),and National Basic Research Program of China(973 Project)(Grant No.2011CB710605)are gratefully acknowledged.

    References

    Abbasy F,Hassani FP,Madiseh SAG,C?té J,Nokken MR.An experimental study on the effective parameters of thermal conductivity of mine backfill.Heat Transfer Engineering 2014;35(13):1209-24.

    Abu-Hamdeh N,Reeder RC.Soil thermal conductivity:effects of density,moisture, salt concentration,and organic matter.Soil Science Society of America Journal 2000;64(4):1285-90.

    Barry-Macaulay D,Bouazza A,Singh RM,Wang B,Ranjith PG.Thermal conductivity of soils and rocks from the Melbourne(Australia)region.Engineering Geology 2013;164:131-8.

    Benítez-Buelga J,Sayde C,Rodríguez-Sinobas L,Selker JS.Heated fiber optic distributed temperature sensing:a dual-probe heat-pulse approach.Vadose Zone Journal 2014;13(11).http://dx.doi.org/10.2136/vzj2014.02.0014.

    Calamita G,Brocca L,Perrone A,Piscitelli S,Lapenna V,Melone F,Moramarco T. Electrical resistivity and TDR methods for soil moisture estimation in central Italy test-sites.Journal of Hydrology 2012;454-455:101-12.

    Cho E,Choi M.Regional scale spatio-temporal variability of soil moisture and its relationship with meteorological factors over the Korean peninsula.Journal of Hydrology 2014;516:317-29.

    Chung S,Horton R.Soil heat and water flow with a partial surface mulch.Water Resources Research 1987;23(12):2175-86.

    Ciocca F,Lunati I,van de Giesen N,Parlange MB.Heated optical fiber for distributed soil-moisture measurements:a lysimeter experiment.Vadose Zone Journal 2012;11(4):1-10.

    C?té A,Carrier B,Leduc J,No?l P,Beauchemin R,Soares M,Garneau C,Gervais R. Water leakage detection using optical fiber at the Peribonka dam.In:Proceedings of the 7th International Symposium on Field Measurements in Geomechanics.American Society of Civil Engineers;2007.p.1-12.

    C?té J,Konrad JM.A generalized thermal conductivity model for soils and construction materials.Canadian Geotechnical Journal 2005;42(2):443-58.

    Dobriyal P,Qureshi A,Badola R,Hussain SA.A review of the methods available for estimating soil moisture and its implications for water resource management. Journal of Hydrology 2012;458-459:110-7.

    Doolittle JA,Collins ME.Use of soil information to determine application of ground penetrating radar.Journal of Applied Geophysics 1995;33(1-3):101-8.

    Gil-Rodríguez M,Rodríguez-Sinobas L,Benítez-Buelga J,Sánchez-Calvo R.Application of active heat pulse method with fiber optic temperature sensing for estimation of wetting bulbs and water distribution in drip emitters.Agricultural Water Management 2012;120:72-8.

    Grattan KTV,Sun T.Fiber optic sensor technology:an overview.Sensors Actuators A:Physical 2000;82(1-3):40-61.

    Hartog AH.A distributed temperature sensor based on liquid core optical fibers. Journal of Lightwave Technology 1983;1(3):498-509.

    Jougnot D,Revil A.Thermal conductivity of unsaturated clay-rocks.Hydrology and Earth System Sciences 2010;14(1):91-8.

    Lu S,Ren T,Gong Y,Horton R.An improved model for predicting soil thermal conductivity from water content at room temperature.Soil Science Society of America Journal 2007;71(1):8-14.

    Lunt IA,Hubbard SS,Rubin Y.Soil moisture content estimation using groundpenetrating radar reflection data.Journal of Hydrology 2005;307(1-4):254-69.

    Luo GY,Li SL.Principles of engineering geology.Nanjing:Nanjing University Press;1990.

    Maleki M,Bayat M.Experimental evaluation of mechanical behavior of unsaturated silty sand under constant water content condition.Engineering Geology 2012;141-142:45-56.

    Markle JM,Schincariol RA,Sass JH,Molson JW.Characterizing the two-dimensional thermal conductivity distribution in a sand and gravel aquifer.Soil Science Society of America Journal 2006;70(4):1281-94.

    Mittelbach H,Lehner I,Seneviratne SI.Comparison of four soil moisture sensor types under field conditions in Switzerland.Journal of Hydrology 2012;430-431:39-49.

    Read TO,Bour JS,Selker VF,Bense T,Borgne R,Hochreutener LN.Active-distributed temperature sensing to continuously quantify vertical flow in boreholes.Water Resources Research 2014;50(5):3706-13.

    Rowlandson TL,Berg AA,Bullock PR,Ojo ER,McNairn H,Wiseman G,Cosh MH. Evaluation of several calibration procedures for a portable soil moisture sensor. Journal of Hydrology 2013;498:335-44.

    Sayde C,Buelga JB,Rodriguez-Sinobas L,Khoury LE,Engish M,van de Giesen N, Selker JS.Mapping variability of soil water content and flux across 1-1000 m scales using the actively heated fiber optic method.Water Resources Research 2014;50(9):7302-17.

    Sayde C,Gregory C,Gil-Rodriguez M,Tufillaro N,Tyler S,van de Giesen N,Selker JS. Feasibility of soil moisture monitoring with heated fiber optics.Water Resources Research 2010;46(6):1-8.

    Schmugge TJ,Jackson TJ,McKim HL.Survey of methods for soil moisture determination.Water Resources Research 1980;16(6):961-79.

    Schrott L,Sass O.Application of field geophysics in geomorphology:advances and limitations exemplified by case studies.Geomorphology 2008;93(1-2):55-73. Steele-Dunne SC,Rutten MM,Krzeminska DM,Hausner M,Tyler SW,Selker J, Bogaard TA,van de Giesen N.Feasibility of soil moisture estimation using passive distributed temperature sensing.Water Resources Research 2010;46(3):1-12.

    Striegl AM,Loheide SP.Heated distributed temperature sensing for field scale soil moisture monitoring.Ground Water 2012;50(3):340-7.

    Sun YK,Li Q,Li XY,Yang DX.Progress of real-time monitoring technology in oil and gas industry based on fiber Bragg grating sensing.Science&Technology Review 2015;33(12).http://dx.doi.org/10.3981/j.issn.1000-7857.

    Susha Lekshmi SU,Singh DN,Baghini MS.A critical review of soil moisture measurement.Measurement 2014;54:92-105.

    Tyler SW,Selker JS,Hausner MB,Hatch CE,Torgersen T,Thodal CE,Geoffrey Schladow S.Environmental temperature sensing using Raman spectra DTS fiber-optic methods.Water Resources Research 2009;45(4):1-11.

    Vogt T,Schneider P,Hahn-Woernle L,Cirpka OA.Estimation of seepage rates in a losing stream by means of fiber-optic high-resolution vertical temperature profiling.Journal of Hydrology 2010;380(1-2):154-64.

    Weihermüller L,Huisman JA,Lambot S,Herbst M,Vereecken H.Mapping the spatial variation of soil water content at the field scale with different ground penetrating radar techniques.Journal of Hydrology 2007;340(3-4):205-16.

    Weiss JD.Using fiber optics to detect moisture intrusion into a landfill cap consisting of a vegetative soil barrier.Journal of the Air&Waste Management Association 2003;53(9):1130-48.

    Wierenga PJ,Nielsen DR,Hagan RM.Thermal properties of a soil based upon field andlaboratorymeasurements.SoilScienceSocietyofAmericaJournal 1969;33(3):354-60.

    Yin Z,Lei T,Yan Q,Chen Z,Dong Y.A near-infrared reflectance sensor for soil surface moisture measurement.Computers and Electronics in Agriculture 2013;99:101-7.

    Dingfeng Cao is currently pursuing a Ph.D.degree in geological engineering from Nanjing University,China.He is mainly engaged in engineering geological research and the application of distributed fiber optic sensors in monitoring soil moisture.Contact information:School of Earth Sciences and Engineering,Nanjing University,163 Xianlin Avenu,Nanjing,210023,China.Tel.:+86-18205187501. E-mail:caodingfeng2014@gmail.com

    15 April 2015

    *Corresponding authors.Tel.:+86 25 85386640.

    E-mail addresses:shibin@nju.edu.cn(B.Shi),zhh@nju.edu.cn(H.Zhu).

    Peer review under responsibility of Institute of Rock and Soil Mechanics, Chinese Academy of Sciences.

    1674-7755?2015 Institute of Rock and Soil Mechanics,Chinese Academy of Sciences.Production and hosting by Elsevier B.V.All rights reserved.

    http://dx.doi.org/10.1016/j.jrmge.2015.08.003

    猜你喜歡
    機(jī)加成株吡蚜
    柔性智能化檢測技術(shù)助力連桿機(jī)加質(zhì)量提升
    汽車工藝師(2024年6期)2024-06-24 06:15:32
    研究發(fā)現(xiàn)黃瓜成株期耐熱性新基因
    Cr12Mo1V1鋼制塔盤沖孔模具制造
    2018年貴州省玉米區(qū)域試驗L組六枝點(diǎn)總結(jié)
    壓力容器機(jī)加裝夾方式改進(jìn)實踐
    機(jī)加零件精益單元構(gòu)建與應(yīng)用
    控制釋放農(nóng)藥顆粒劑吡蚜酮防治稻飛虱田間藥效試驗總結(jié)
    更 正
    黃瓜耐弱光性鑒定指標(biāo)研究及其苗期與成株期的相關(guān)分析
    中國蔬菜(2016年8期)2017-01-15 14:23:36
    50%吡蚜酮水分散粒劑防治水稻白背飛虱效果研究
    耕作與栽培(2015年1期)2015-07-22 08:27:00
    欧美亚洲日本最大视频资源| 热re99久久国产66热| 男女无遮挡免费网站观看| 我的亚洲天堂| 亚洲av成人精品一二三区| 国产毛片在线视频| 亚洲欧洲日产国产| 视频在线观看一区二区三区| 国产毛片在线视频| 国产成人欧美| 日本vs欧美在线观看视频| 中文字幕色久视频| 亚洲在久久综合| 成人二区视频| 2021少妇久久久久久久久久久| 亚洲精品国产色婷婷电影| 在线观看免费日韩欧美大片| 国产男女超爽视频在线观看| 在线观看一区二区三区激情| 少妇人妻精品综合一区二区| 少妇人妻久久综合中文| 国产精品香港三级国产av潘金莲 | 国产又色又爽无遮挡免| 三上悠亚av全集在线观看| 亚洲精品av麻豆狂野| 人人妻人人澡人人看| a级毛片在线看网站| 欧美日本中文国产一区发布| 亚洲图色成人| 少妇 在线观看| 我要看黄色一级片免费的| 2018国产大陆天天弄谢| 免费高清在线观看视频在线观看| 欧美日本中文国产一区发布| 亚洲,欧美,日韩| 亚洲图色成人| a 毛片基地| 少妇人妻 视频| 免费高清在线观看日韩| 国产有黄有色有爽视频| 考比视频在线观看| 蜜桃国产av成人99| 中文字幕最新亚洲高清| 哪个播放器可以免费观看大片| 亚洲一码二码三码区别大吗| 麻豆精品久久久久久蜜桃| 欧美中文综合在线视频| 国产熟女午夜一区二区三区| 在线天堂最新版资源| videos熟女内射| 丰满迷人的少妇在线观看| 日本wwww免费看| 久久精品国产亚洲av高清一级| 欧美激情极品国产一区二区三区| 一级黄片播放器| 国产av国产精品国产| 亚洲成国产人片在线观看| 日韩制服丝袜自拍偷拍| 哪个播放器可以免费观看大片| 国产成人精品在线电影| 日韩电影二区| 午夜免费观看性视频| 热99久久久久精品小说推荐| 欧美国产精品va在线观看不卡| 国产精品女同一区二区软件| 亚洲国产av影院在线观看| 在线观看免费日韩欧美大片| 男女下面插进去视频免费观看| 免费看不卡的av| 久久人人97超碰香蕉20202| 最近的中文字幕免费完整| 在线观看美女被高潮喷水网站| 美女主播在线视频| 少妇熟女欧美另类| 色播在线永久视频| 精品少妇内射三级| 国产极品粉嫩免费观看在线| 天堂中文最新版在线下载| 国产伦理片在线播放av一区| 久久精品aⅴ一区二区三区四区 | 王馨瑶露胸无遮挡在线观看| 日本wwww免费看| 久久久久久久久久久久大奶| 永久免费av网站大全| 欧美av亚洲av综合av国产av | 女人精品久久久久毛片| 高清黄色对白视频在线免费看| 99精国产麻豆久久婷婷| 日本91视频免费播放| 亚洲成色77777| 一个人免费看片子| 久久久久久人妻| 777久久人妻少妇嫩草av网站| 91精品国产国语对白视频| 国产成人午夜福利电影在线观看| videos熟女内射| 天天影视国产精品| 亚洲一区中文字幕在线| www日本在线高清视频| 亚洲 欧美一区二区三区| 老司机影院成人| 国产成人精品久久久久久| av在线老鸭窝| 国产在视频线精品| 18在线观看网站| 国产成人精品在线电影| av在线老鸭窝| 亚洲国产欧美日韩在线播放| 国产精品无大码| 熟女av电影| 老汉色∧v一级毛片| 欧美日韩精品网址| 亚洲成人手机| 亚洲情色 制服丝袜| xxx大片免费视频| 这个男人来自地球电影免费观看 | 在现免费观看毛片| 久久久久久久久久人人人人人人| 国产成人aa在线观看| 9色porny在线观看| 一级片'在线观看视频| 亚洲 欧美一区二区三区| 亚洲精品国产一区二区精华液| 日韩一本色道免费dvd| 亚洲欧美日韩另类电影网站| 老司机亚洲免费影院| 另类精品久久| 中文字幕亚洲精品专区| 在线免费观看不下载黄p国产| 亚洲,一卡二卡三卡| 日本猛色少妇xxxxx猛交久久| 少妇人妻久久综合中文| 美女大奶头黄色视频| 国产女主播在线喷水免费视频网站| 欧美激情高清一区二区三区 | 国产亚洲精品第一综合不卡| 中文字幕亚洲精品专区| 人妻系列 视频| a级毛片黄视频| 少妇 在线观看| 久久久久网色| 中文字幕精品免费在线观看视频| 色婷婷久久久亚洲欧美| 精品卡一卡二卡四卡免费| 最近中文字幕高清免费大全6| 久久久久久久久免费视频了| 91久久精品国产一区二区三区| 91久久精品国产一区二区三区| 男女高潮啪啪啪动态图| 中文字幕最新亚洲高清| 中国国产av一级| 搡老乐熟女国产| 久久久久国产一级毛片高清牌| 亚洲成人av在线免费| 国产欧美亚洲国产| 黄片小视频在线播放| 97精品久久久久久久久久精品| 蜜桃国产av成人99| 激情视频va一区二区三区| 啦啦啦啦在线视频资源| 国产精品嫩草影院av在线观看| 国产片特级美女逼逼视频| 国产精品熟女久久久久浪| 国产片特级美女逼逼视频| 久久精品国产综合久久久| 少妇的丰满在线观看| 日本午夜av视频| av福利片在线| 亚洲av成人精品一二三区| 欧美精品一区二区大全| 大话2 男鬼变身卡| 国产福利在线免费观看视频| 亚洲国产日韩一区二区| 黄色毛片三级朝国网站| 国产片内射在线| 亚洲国产欧美在线一区| 亚洲一区中文字幕在线| 精品国产一区二区三区久久久樱花| 捣出白浆h1v1| 国产av精品麻豆| 精品少妇一区二区三区视频日本电影 | 在线观看免费日韩欧美大片| 欧美人与善性xxx| 欧美在线黄色| 日韩欧美一区视频在线观看| 国产精品二区激情视频| 国产精品二区激情视频| 国产又爽黄色视频| 国产亚洲精品第一综合不卡| 天天躁夜夜躁狠狠躁躁| 777久久人妻少妇嫩草av网站| 9热在线视频观看99| 亚洲av免费高清在线观看| 热re99久久精品国产66热6| 久久精品国产鲁丝片午夜精品| 涩涩av久久男人的天堂| 欧美av亚洲av综合av国产av | 欧美日韩一级在线毛片| 亚洲国产日韩一区二区| 久久久国产欧美日韩av| 国产成人aa在线观看| 国产精品女同一区二区软件| 国产日韩一区二区三区精品不卡| 国产成人精品婷婷| 免费观看性生交大片5| 啦啦啦在线免费观看视频4| 欧美日韩视频高清一区二区三区二| 国精品久久久久久国模美| 久久ye,这里只有精品| 亚洲,一卡二卡三卡| 久久久久久久精品精品| 王馨瑶露胸无遮挡在线观看| 天堂8中文在线网| 国产精品免费视频内射| 曰老女人黄片| 日韩一区二区三区影片| 国产成人精品一,二区| 天天影视国产精品| 国产不卡av网站在线观看| 看免费成人av毛片| 青草久久国产| 久久青草综合色| 国产精品国产三级国产专区5o| 美女午夜性视频免费| 欧美日本中文国产一区发布| 国产福利在线免费观看视频| 久久久精品国产亚洲av高清涩受| 少妇人妻久久综合中文| 久久韩国三级中文字幕| 免费黄色在线免费观看| 国产在线一区二区三区精| 亚洲第一青青草原| 亚洲av综合色区一区| 999精品在线视频| 妹子高潮喷水视频| 男人舔女人的私密视频| 国产免费又黄又爽又色| 大陆偷拍与自拍| 午夜福利网站1000一区二区三区| 麻豆av在线久日| 久久99蜜桃精品久久| 99热网站在线观看| 桃花免费在线播放| 免费播放大片免费观看视频在线观看| 久久精品久久久久久噜噜老黄| 各种免费的搞黄视频| 成年人免费黄色播放视频| 大香蕉久久网| 香蕉丝袜av| 国产精品国产三级专区第一集| 久久精品国产自在天天线| 亚洲精品乱久久久久久| 黄色毛片三级朝国网站| 免费日韩欧美在线观看| 男人舔女人的私密视频| 免费大片黄手机在线观看| 激情五月婷婷亚洲| 高清av免费在线| 综合色丁香网| 国产免费一区二区三区四区乱码| 久久国内精品自在自线图片| 久久久久久久大尺度免费视频| 免费黄网站久久成人精品| 精品人妻一区二区三区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 成年女人在线观看亚洲视频| 黄片小视频在线播放| 九草在线视频观看| 叶爱在线成人免费视频播放| 亚洲国产av新网站| 国产精品秋霞免费鲁丝片| 一本—道久久a久久精品蜜桃钙片| 亚洲国产av新网站| 成年女人毛片免费观看观看9 | av电影中文网址| 免费在线观看视频国产中文字幕亚洲 | 熟女少妇亚洲综合色aaa.| 在线观看国产h片| 99九九在线精品视频| 亚洲一区中文字幕在线| 午夜精品国产一区二区电影| 精品人妻一区二区三区麻豆| 秋霞伦理黄片| 九草在线视频观看| 欧美精品人与动牲交sv欧美| av福利片在线| 男女边吃奶边做爰视频| 国产男人的电影天堂91| 久久久久视频综合| 亚洲欧美一区二区三区国产| 999久久久国产精品视频| 久久久久久久亚洲中文字幕| 亚洲,欧美精品.| 国产高清不卡午夜福利| 成人毛片60女人毛片免费| 欧美日韩视频精品一区| 啦啦啦在线观看免费高清www| 成人18禁高潮啪啪吃奶动态图| av国产精品久久久久影院| 老司机亚洲免费影院| 熟妇人妻不卡中文字幕| 两个人免费观看高清视频| 波野结衣二区三区在线| 交换朋友夫妻互换小说| 看非洲黑人一级黄片| 综合色丁香网| 妹子高潮喷水视频| 看十八女毛片水多多多| 一级,二级,三级黄色视频| 一本—道久久a久久精品蜜桃钙片| 国产日韩一区二区三区精品不卡| 午夜福利影视在线免费观看| 熟女少妇亚洲综合色aaa.| 热re99久久国产66热| 婷婷色av中文字幕| 少妇的逼水好多| av.在线天堂| 日韩精品有码人妻一区| 日韩在线高清观看一区二区三区| 狠狠婷婷综合久久久久久88av| 免费高清在线观看视频在线观看| 日本猛色少妇xxxxx猛交久久| 一级a爱视频在线免费观看| 亚洲三区欧美一区| 亚洲国产欧美日韩在线播放| 免费少妇av软件| 午夜精品国产一区二区电影| 少妇的逼水好多| 夜夜骑夜夜射夜夜干| 亚洲图色成人| 久久精品国产a三级三级三级| 亚洲,欧美,日韩| 久久久久网色| 亚洲图色成人| 黑丝袜美女国产一区| 国产熟女午夜一区二区三区| 亚洲欧洲日产国产| 97精品久久久久久久久久精品| 国产精品麻豆人妻色哟哟久久| 欧美激情极品国产一区二区三区| 18在线观看网站| 9色porny在线观看| 97人妻天天添夜夜摸| 日韩视频在线欧美| 丝袜美腿诱惑在线| 欧美bdsm另类| 我的亚洲天堂| 性色avwww在线观看| av国产精品久久久久影院| 午夜激情av网站| 一区福利在线观看| 国产男女超爽视频在线观看| 爱豆传媒免费全集在线观看| 在现免费观看毛片| 王馨瑶露胸无遮挡在线观看| 一区二区三区激情视频| 久久婷婷青草| 国产精品国产三级专区第一集| 男女高潮啪啪啪动态图| 国产精品99久久99久久久不卡 | 久久99精品国语久久久| 精品久久蜜臀av无| 春色校园在线视频观看| 日本91视频免费播放| 99精国产麻豆久久婷婷| 亚洲综合色网址| 国产乱来视频区| 久久久欧美国产精品| 日韩一本色道免费dvd| 一级毛片 在线播放| 亚洲婷婷狠狠爱综合网| 亚洲美女搞黄在线观看| 在线精品无人区一区二区三| 91午夜精品亚洲一区二区三区| 亚洲精品第二区| 人妻一区二区av| 久久久国产一区二区| 大片免费播放器 马上看| 丝瓜视频免费看黄片| 女人久久www免费人成看片| 精品酒店卫生间| 亚洲人成电影观看| 一边亲一边摸免费视频| 曰老女人黄片| 久热这里只有精品99| 天天躁夜夜躁狠狠久久av| 国产精品免费大片| 超碰97精品在线观看| 亚洲精华国产精华液的使用体验| 啦啦啦在线观看免费高清www| 亚洲国产最新在线播放| 在线观看人妻少妇| 国产精品二区激情视频| 大香蕉久久成人网| 丰满迷人的少妇在线观看| 性色av一级| 人妻一区二区av| 亚洲一区二区三区欧美精品| 午夜老司机福利剧场| 亚洲美女搞黄在线观看| 久久人人97超碰香蕉20202| a级毛片黄视频| 亚洲欧美成人综合另类久久久| 一本久久精品| 中文字幕精品免费在线观看视频| 欧美在线黄色| 老鸭窝网址在线观看| 成人影院久久| 一区二区三区激情视频| 男女边吃奶边做爰视频| 国产精品嫩草影院av在线观看| 成人午夜精彩视频在线观看| 久久免费观看电影| 午夜免费鲁丝| 自线自在国产av| 高清黄色对白视频在线免费看| 成年动漫av网址| 男女午夜视频在线观看| 纵有疾风起免费观看全集完整版| 色视频在线一区二区三区| 三上悠亚av全集在线观看| 久久久国产精品麻豆| 日韩大片免费观看网站| 国产黄色视频一区二区在线观看| 精品少妇黑人巨大在线播放| 亚洲美女视频黄频| 热99国产精品久久久久久7| 少妇被粗大的猛进出69影院| 欧美日韩一级在线毛片| av不卡在线播放| 欧美老熟妇乱子伦牲交| 五月伊人婷婷丁香| 午夜91福利影院| 成年女人在线观看亚洲视频| 亚洲美女黄色视频免费看| 亚洲精品av麻豆狂野| 少妇精品久久久久久久| 国产女主播在线喷水免费视频网站| 校园人妻丝袜中文字幕| 美女午夜性视频免费| 亚洲美女黄色视频免费看| 卡戴珊不雅视频在线播放| 黄色 视频免费看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av中文av极速乱| 亚洲精品国产色婷婷电影| 飞空精品影院首页| 精品一区二区三卡| 视频在线观看一区二区三区| 国产综合精华液| 久久这里有精品视频免费| 国产极品粉嫩免费观看在线| 国产精品麻豆人妻色哟哟久久| 精品少妇一区二区三区视频日本电影 | 亚洲中文av在线| 久久这里只有精品19| 日韩免费高清中文字幕av| 两个人免费观看高清视频| 男女无遮挡免费网站观看| 午夜福利网站1000一区二区三区| 精品国产一区二区久久| 在线 av 中文字幕| 欧美激情 高清一区二区三区| 91aial.com中文字幕在线观看| 热re99久久国产66热| 亚洲国产看品久久| 汤姆久久久久久久影院中文字幕| 日韩成人av中文字幕在线观看| 天天躁夜夜躁狠狠躁躁| 久久精品国产亚洲av高清一级| 大码成人一级视频| 国产精品国产三级国产专区5o| 国产在线免费精品| 自拍欧美九色日韩亚洲蝌蚪91| 国产无遮挡羞羞视频在线观看| 天天躁日日躁夜夜躁夜夜| 人人妻人人添人人爽欧美一区卜| 国产亚洲一区二区精品| 一个人免费看片子| 这个男人来自地球电影免费观看 | 看免费成人av毛片| 国产成人91sexporn| 在线免费观看不下载黄p国产| 在线天堂中文资源库| 午夜激情久久久久久久| 久热久热在线精品观看| 高清黄色对白视频在线免费看| 校园人妻丝袜中文字幕| av视频免费观看在线观看| 久久这里只有精品19| 日韩 亚洲 欧美在线| 亚洲精品视频女| 一本大道久久a久久精品| 国产在线免费精品| 欧美国产精品一级二级三级| 久久精品国产亚洲av天美| 不卡视频在线观看欧美| 精品一区二区三卡| 国产精品一国产av| 欧美 亚洲 国产 日韩一| 九草在线视频观看| av网站在线播放免费| av线在线观看网站| 亚洲精品久久久久久婷婷小说| 在现免费观看毛片| 日韩精品免费视频一区二区三区| 免费人妻精品一区二区三区视频| 春色校园在线视频观看| 国产在线免费精品| 男女高潮啪啪啪动态图| 亚洲国产毛片av蜜桃av| 大香蕉久久成人网| 亚洲欧美一区二区三区国产| 欧美人与善性xxx| 欧美日韩视频精品一区| 亚洲精品成人av观看孕妇| 国产成人精品无人区| 黄色怎么调成土黄色| 波多野结衣av一区二区av| 久久久久精品久久久久真实原创| 久久精品夜色国产| av.在线天堂| 少妇猛男粗大的猛烈进出视频| a级片在线免费高清观看视频| 亚洲av男天堂| 国产精品蜜桃在线观看| 亚洲精品日韩在线中文字幕| 99久久精品国产国产毛片| 日韩,欧美,国产一区二区三区| 高清av免费在线| 亚洲图色成人| 亚洲美女黄色视频免费看| 亚洲精品av麻豆狂野| 日韩成人av中文字幕在线观看| 伊人久久国产一区二区| 交换朋友夫妻互换小说| 久久精品国产综合久久久| 亚洲av国产av综合av卡| 亚洲国产最新在线播放| 欧美日韩综合久久久久久| 日韩成人av中文字幕在线观看| 国产不卡av网站在线观看| av网站在线播放免费| 天天影视国产精品| 男女边摸边吃奶| 久久精品夜色国产| 国产免费福利视频在线观看| 91精品伊人久久大香线蕉| 一区福利在线观看| 中国三级夫妇交换| 伊人久久大香线蕉亚洲五| 亚洲天堂av无毛| 久久精品夜色国产| 男女下面插进去视频免费观看| 久久韩国三级中文字幕| 午夜福利,免费看| 欧美中文综合在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 最近最新中文字幕免费大全7| 熟妇人妻不卡中文字幕| 久久99蜜桃精品久久| 久久精品夜色国产| 午夜福利网站1000一区二区三区| 免费日韩欧美在线观看| 最近最新中文字幕大全免费视频 | 欧美激情高清一区二区三区 | av天堂久久9| 青春草视频在线免费观看| 国产国语露脸激情在线看| 亚洲精品在线美女| 黑人猛操日本美女一级片| 天天影视国产精品| 99久久中文字幕三级久久日本| 欧美成人午夜免费资源| 69精品国产乱码久久久| 欧美亚洲日本最大视频资源| 国产精品成人在线| 一区二区三区激情视频| 汤姆久久久久久久影院中文字幕| 国产黄色视频一区二区在线观看| 男女免费视频国产| 美国免费a级毛片| 美女脱内裤让男人舔精品视频| 国产1区2区3区精品| 成年人免费黄色播放视频| 黄色怎么调成土黄色| 人人妻人人爽人人添夜夜欢视频| 国产精品一区二区在线不卡| 热99久久久久精品小说推荐| 2018国产大陆天天弄谢| 免费少妇av软件| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产精品国产精品| 亚洲国产av影院在线观看| 亚洲国产成人一精品久久久| 啦啦啦中文免费视频观看日本| 伦精品一区二区三区| xxxhd国产人妻xxx| 亚洲在久久综合| xxxhd国产人妻xxx| 伦理电影大哥的女人| 91精品国产国语对白视频| 亚洲成色77777| 亚洲精品国产av成人精品| 中文字幕另类日韩欧美亚洲嫩草| 一本久久精品| 亚洲第一区二区三区不卡| 热99国产精品久久久久久7| 欧美xxⅹ黑人| 老司机亚洲免费影院| 国产成人精品久久久久久| 亚洲精品国产av成人精品| 国产成人精品久久久久久| 欧美日韩av久久| 亚洲精品日韩在线中文字幕|