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

    Field Observations of Near-Surface Wind Flow Across Expressway Embankment on the Qinghai-Tibet Plateau

    2022-10-18 12:29:22YanhuMuWiMaZhaohuiJoyYangXiaolinLiKunZhangYunchngMao
    Engineering 2022年7期

    Yanhu Mu, Wi Ma,*, Zhaohui (Joy) Yang*, Xiaolin Li, Kun Zhang, Yunchng Mao

    a State Key Laboratory of Frozen Ground Engineering, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China

    b University of Chinese Academy of Sciences, Beijing 100049, China

    c College of Engineering, University of Alaska Anchorage, Anchorage, AK 99508, USA

    d College of Civil Engineering and Architecture, Jiaxing University, Jiaxing 314001, China

    e School of Civil Engineering, Northwest Minzu University, Lanzhou 730024, China

    Keywords:Near-surface wind flow Field observation Air-cooling structures Linear transportation infrastructure Qinghai-Tibet Plateau

    ABSTRACT Crushed rock layers(CRLs),ventilation ducts(VDs)and thermosyphons are air-cooling structures(ACSs)widely used for maintaining the long-term stability of engineered infrastructures in permafrost environments.These ACSs can effectively cool and maintain the permafrost subgrade’s frozen state under climate warming by facilitating heat exchange with ambient air in cold seasons.As convection is a crucial working mechanism of these ACSs,it is imperative to understand the near-surface wind flow(NSWF)across a constructed infrastructure,such as an embankment.This article describes a yearlong field observation of the NSWF across an experimental expressway embankment, the first of its kind on the Qinghai-Tibet Plateau(QTP).The wind speed and direction along a transect perpendicular to the embankment on both the windward and leeward sides and at four different heights above the ground surface were collected and analyzed. The results showed that the embankment has a considerable impact on the NSWF speed within a distance of up to ten times its height, and in the direction on the leeward side. A power law can well describe the speed profiles of NSWF across the embankment, with the power-law indices(PLIs)varying from 0.14 to 0.40.On an annual basis,the fitted NSWF PLI far away from the embankment was 0.19,which differs substantially from the values widely used in previous thermal performance evaluations of ACSs on the QTP. Finally,the significance of the NSWF to the thermal performance of the ACSs,particularly the CRLs and VDs, in linear transportation infrastructure is discussed. It is concluded that underestimating the PLI and neglecting wind direction variations may lead to unconservative designs of the ACSs.The results reported in this study can provide valuable guidance for infrastructure engineering on the QTP and other similar permafrost regions.

    1. Introduction

    One-quarter of the Northern Hemisphere and 17% of Earth’s exposed land surface are underlain by permafrost [1]. These vast permafrost regions play an important role in the world economy due to their rich natural resources, including an estimated 25% of the world’s undiscovered oil and gas reserves [2,3]. In recent decades, economic development has resulted in an extensive expansion of hydrocarbon extraction, transportation networks,communication lines, industrial projects, civil facilities, and engineering maintenance systems in permafrost regions [4].However, ongoing climate warming and consequent permafrost degradation have severely damaged thousands of infrastructure components and led to significant economic losses [5-10]. The most recent assessments show that by 2050,33%of the infrastructure built in the Northern Hemisphere’s permafrost regions will experience damage from ground subsidence and loss of structural bearing capacity induced by permafrost degradation [11]. Thus,how to economically maintain the existing infrastructure and scientifically construct future infrastructures on degrading permafrost is a great engineering challenge. It is critical to regulate heat exchange between a permafrost foundation and the ambient environment because the permafrost thermal regime is a key variable that determines the bearing capacity and deformation behaviors of foundations [12-14].

    Construction activity affects the natural ground surface conditions and the ground thermal regime by changing the energy exchange between the ground and the atmosphere.For permafrost roadway embankments,such an impact can be substantial as great changes in material properties, surface albedo, super-permafrost hydrology and vegetation covers can occur during embankment construction and operation. High-porosity crushed rock layers(CRLs),ventilation ducts(VDs)and thermosyphons have been proven to be effective air-cooling structures(ACSs)for maintaining the long-term thermal stability of permafrost subgrades beneath various engineered infrastructures [15]. By facilitating heat exchange with ambient air in cold seasons, these ACSs can cool down and maintain permafrost foundations to protect engineered infrastructures from damage induced by permafrost warming and thawing during climate warming. In cold seasons, CRLs and VDs can enhance the natural ventilation driven by wind pressure and temperature differences. Because convection is a much more efficient heat-transfer mechanism,heat in the underlying permafrost could be transferred into the ambient environment more rapidly [15-18].A thermosyphon is a heat exchanger that can effectively transfer a heat source to a heat sink through the evaporation and condensation of a working fluid [19]. Because of their superior heattransfer ability, thermosyphons have been extensively used in the thermal management of a wide range of engineering systems[20-27].When used for permafrost cooling,thermosyphons generally function in winter and transfer the heat from the evaporation section buried in permafrost to the condensation section exposed in the air [20,22,25]. In the past few decades, these three types of ACSs have been widely used for roadway embankments, airfield runways, oil and gas pipeline systems, buildings and power transmission lines in northern Canada,Alaska,Siberia,and the Qinghai-Tibet Plateau (QTP) [12,13,28-30] (Fig. 1).

    The near-surface wind flow(NSWF)is a factor that controls the thermal performance of these ACSs when used in linear transportation infrastructures. The speed and direction of the NSWF determine not only the heat convection type within the CRLs, namely,wind-driven convection or buoyancy-driven convection, but also the intensity of heat convection [16,17,31-34]. Wu et al. [33]studied the cooling mechanisms of the CRL in embankments through distributed temperature measurements within an experimental railway section. The results showed that, in cold seasons, wind-driven convection occurred within the CRL when the ambient wind speed was large, whereas weak buoyancydriven convection occurred when the ambient wind speed was small. They concluded that the cooling mechanisms and performance of the CRL were controlled by the ambient wind speed and direction. For VDs used for roadway and railway embankments, with a given diameter and a burial height within the embankments, the airflow speed within the duct is closely related to the ambient wind speed and direction[17,35,36].The NSWF also represents a crucial operating condition for thermosyphons used for permafrost foundation cooling[37,38].In coupled air-thermosyphon-soil heat transfer systems, the NSWF speed determines the convective heat transfer of the condenser section and hence its heat loss rate during the operation of thermosyphons.

    However, to date, field-observed data of NSWFs across engineered infrastructures equipped with ACSs are very scarce. Thus,the wind flow characteristics required to assess the long-term thermal performance of the infrastructure with ACSs were generally assumed based on a prevailing wind direction and a mean or a seasonally varying wind speed [22,25,30,35-43]. These assumptions can result in either over or underestimation of the heat exchange efficiency of the ACSs and hence, unreliable prediction of the permafrost subgrade’s long-term thermal regime, which may further lead to an unconservative or uneconomical foundation design.

    This article describes systemic field observations of the NSWF across an experimental expressway embankment in the Beiluhe Basin,the central area of the QTP.A total of six meteorological towers were established at both sides of the embankment to collect the NSWF data, including wind speeds at four levels (0.5, 1.5, 3.0, and 4.5 m) above the ground surface and wind directions. Data collection at these observation sites and an adjacent meteorological station serving as a reference was carried out for an entire year. The characteristics of the NSWF across the embankment, including its distribution and seasonal variations, were analyzed based on the field observed data.Based on field-observed results,the impacts of the wind directions and speed profiles on the thermal performances of the ACSs were discussed. This study’s observation data and results will provide an essential scientific basis for the design and maintenance of the ACSs used for linear transportation infrastructure built on the QTP, for example, the existing Qinghai-Tibet Railway and Qinghai-Tibet Highway and the proposed Qinghai-Tibet Expressway.

    Fig.1. The individual or combined use of ACSs in a wide range of engineered infrastructures on the QTP:(a)CRLs,(b)CRLs and VDs,and(c)thermosyphons used in railway embankments; (d) CRLs, (e) CRLs and VDs, and (f) thermosyphons used in roadway embankments; and thermosyphons in foundations of (g) power transmission lines,(h) railway bridges, and (i) buildings. The building base was also elevated above the ground surface to promote natural ventilation.

    2. Field site and methods

    2.1. Field site

    The field site(92°55′36′′E,34°49′21′′N)is located in the Beiluhe Basin,the interior of the QTP,as shown in Fig.2(a).At an elevation of 4618 m above sea level,the local terrain is slightly sloped,with 70%-90%of the ground surface covered by an alpine meadow.The permafrost at the site is predominately continuous, with mean annual ground temperatures at depths of 10-15 m ranging from-0.5 to-1.0°C.The active layer thickness(the maximum seasonal thaw depth) varies from 2 to 3 m. Ice-rich permafrost layers are encountered immediately beneath the permafrost table based on borehole exploration data. A meteorological station is located approximately 2 km in the northeast direction from the field site.The data from this station show that,from 2003 to 2015,the mean annual air temperature ranged from-2.8 to-4.1°C and the mean annual precipitation ranged from 230 to 460 mm. Data, including air temperature, wind speed, and wind direction at 2 m above the ground, were collected from July 1, 2018 to June 30, 2019(referred to as the study period hereafter),at this station to provide a reference in this study.

    2.2. Experimental expressway embankment

    An experimental expressway embankment with asphalt pavement was built nearly parallel to the Qinghai-Tibet Railway at the field site (refer to Fig. 2(b)) in August 2009 to collect data for constructing expressways in a warm permafrost environment on the QTP.Various ACSs were included in the embankment to investigate their long-term performances in cooling the underlying permafrost. The experimental embankment has a total length of 300 m, an azimuth angle of 197°, a height of 3 m, a crest width of 13 m, and 1:1.5 side slopes. The embankment was equally divided into ten segments, each 30 m in length. Among them, a single ACS, or a combination of ACSs, including CRLs, VDs, thermosyphons,and hollow concrete brick layers,was installed in nine segments, and the segment without any cooling structures was used as a reference [44].

    2.3. Measurement of NSWF

    Six meteorological towers were established along a transect perpendicular to the embankment route on both sides of the reference section, as illustrated in Fig. 3. On each side, the distances between the three towers and the embankment slope toe were 0.5,30,and 60 m,respectively,and as shown in Fig.3(a),the three towers on the right side (from Golmud to Lhasa, China) were labeled TR-0.5,TR-30,and TR-60,respectively,and those on the left side were TL-0.5, TL-30, and TL-60, respectively. The two towers adjacent to the embankment slope toe (i.e., TR-0.5 and TL-0.5)were 4.5 m tall, and the remaining four (i.e., TR-30, TR-60, TL-30,and TL-60)were 3.0 m tall.The site is relatively flat,and the differences among elevations at the locations of the six towers could be neglected. Considering the prevailing wind direction at the study area and the embankment route, towers TR-60, TR-30, and TR-0.5 were regarded as being on the windward side, whereas TL-60, TL-30, and TL-0.5 were on the leeward side.

    Fig.2. (a)The location of the field site;(b)an aerial view of the experimental expressway embankment.The distance between the railway and the experimental expressway ranges from 168 to 186 m. Mt.: mountains.

    Fig.3. Deployment of the NSWF observation towers:(a)schematic diagram of the six meteorological tower locations across the embankment and views of the towers on the(b) leeward and (c) windward sides.

    Fig. 4. Wind speed and direction sensors mounted on (a) 3.0 m tall and (b) 4.5 m tall towers.

    Wind speeds and wind directions at six towers were measured by 010C wind-speed sensors and 020C wind-direction sensors,respectively. All sensors were manufactured and calibrated by Met One Instruments,Inc.(USA).The 010C wind-speed sensor uses a lightweight three-cup anemometer to measure the horizontal wind speed. This sensor has a calibrated operating range of 0-50 m·s-1, a resolution of <0.1 m·s-1, and an accuracy of ±1%. On the 3.0 m tall towers,the three wind-speed sensors were mounted at heights of 0.5,1.5,and 3.0 m above the ground surface(Fig.4(a)).On the 4.5 m tall towers near the embankment toe, four windspeed sensors were mounted at heights of 0.5, 1.5, 3.0, and 4.5 m(Fig. 4(b)). On each tower, the 020C wind-direction sensor, with a resolution of <0.1° and an accuracy of ±3°, was paired with the uppermost wind speed sensor to provide azimuth data (0°-360°)by coupling a lightweight airfoil vane and a single-precision potentiometer. Note that such a sensor is visible on the right end of a horizontal beam in Fig. 4(b). On each side of the embankment, all the sensors were connected to a CR3000 datalogger manufactured and calibrated by Campbell Scientific Inc. (USA) to automatically collect,compute,and store the data at 10 min intervals throughout the study period.

    3. Results and analysis

    3.1. Characteristics of natural wind flow in the study area

    3.1.1. Wind direction

    Fig.5(a)shows the rose diagram of the wind frequency obtained at the meteorological station in the study period.In the study area,west (W), west-southwest (WSW), and west-northwest (WNW)were the prevailing wind directions,and their occurrence frequencies accounted for 18.5%,14.9%,and 13.0%,respectively,during the study period. The frequencies of the other directions were relatively low, ranging from 0.9% to 8.5%.

    The characteristics of the wind directions in cold and warm seasons differed significantly.The cold and warm seasons were determined based on the monthly average air temperature. In the cold season, which was October through April of the following year in this study, the composition of the wind direction was simple(Fig.5(b)).The frequencies of the three prevailing wind directions,namely, W, WSW, and WNW, were all above 18.5%, with a sum of 65.5%. Following these three directions, the SW winds accounted for 11.5%, and the others were all less than 4.3%. However, in the warm season from May to September,the wind direction was more dispersed(Fig.5(c)).The peak frequency was less than 12.4%,with four wind directions,namely,east(E),east-northeast(ENE),northeast (NE), and north-northeast (NNE), accounting for 9.6%-12.4%.The frequencies of the other wind directions varied from 1.6% to 8.1%, with a total of 55.4%.

    3.1.2. Wind speed and air temperature

    As mentioned above, wind speed, air temperature, and their seasonal variations are crucial for the heat transfer processes within the ACSs and the cooling effect on the underlying permafrost subgrade.

    Fig. 6 shows the variations in wind speed and air temperature obtained at the meteorological station in the study period.Throughout this period,wind speed varied from 0 to 18 m·s-1,with a mean annual value of 4.4 m·s-1. Air temperature ranged from-26.2 to 17.5 °C, with a mean annual value of -3.6 °C. Both wind speed and air temperature varied seasonally,exhibiting sinusoidal patterns with opposite trends. The wind was stronger in the cold season and weaker in the warm season, which is crucial for applications of ACSs, particularly VDs, for permafrost subgrade cooling on the QTP [13,16,17].

    Two indices, mean monthly wind speed (MMWS) and mean monthly air temperature (MMAT), were calculated and are listed in Table 1. From November to March of the following year,MMWSs reached or exceeded 4.8 m·s-1, with a maximum of 7.1 m·s-1in February, while MMATs varied between -14.6 and-9.2 °C. From June to September, MMWSs fell below 3.6 m·s-1,with a minimum of 3.1 m·s-1in September, while MMATs ranged from 3.6 to 7.9 °C.

    Fig. 5. Rose diagrams of natural wind frequency in the study area during (a) the entire study period, (b) the cold season, and (c) the warm season. The red dotted line represents the embankment route.

    Fig. 6. Variations in wind speed and air temperature in the study period.

    Table 1 The mean monthly values of wind speed and air temperature in the study period.

    3.2. Variation in wind direction across the embankment

    3.2.1. Observation towers located far from the embankment

    Fig. 7 compares the wind rose diagrams for the towers located farthest away from the embankment, namely, TR-60 and TL-60,in the study period.Note that the two towers were located at a distance of 20 times the height of the embankment.

    The wind roses in both locations were very close to each other.Specifically, the prevailing winds at TR-60, namely, W, WSW, and WNW, accounted for 17.8%, 15.2%, and 14.1% of the total, respectively, with a sum of 47.1%. In comparison, these three winds at TL-60 represented 18.5%, 13.9%, and 14.6% of the total, respectively,with an aggregate of 47%.Moreover,there is no considerable difference among the three rose charts in Figs. 5(a) and 7, indicating that the existence of an embankment had no meaningful impact on the NSWF direction at a distance of 20 times the embankment height away.

    3.2.2. Observation towers located close to the embankment

    Figs. 8(a)-(d) show wind rose diagrams at the four towers located relatively close to the embankment (at a distance of ten times of the embankment height or less), namely, TR-30, TL-30,TR-0.5, and TL-0.5, respectively, during the study period.

    On the windward side, the prevailing winds were the same as those at TR-60 and TL-60. For instance, Fig. 8(a), the wind rose for TR-30 (ten times the embankment height away), reveals that the three prevailing winds accounted for 18.2%, 14.9%, and 14.5%of the total.In comparison,Fig.8(c),the wind rose for TR-0.5(near the embankment toe), shows that the prevailing winds comprised 17.5%,14.3%,and 14.1%of the total.These values were very close to those at TR-60 and TL-60,suggesting that the embankment’s existence did not significantly influence the NSWF directions on the windward side.

    Fig. 7. The wind rose diagrams for (a) TR-60 (windward) and (b) TL-60 (leeward) during the study period. The red dotted line represents the embankment route.

    Fig. 8. Wind rose graphs for (a) TR-30, (b) TL-30, (c) TR-0.5, and (d) TL-0.5 during the study period. The red dotted line represents the embankment route.

    However, on the leeward side, the wind directions were significantly altered at both TL-30 and TL-0.5. Fig. 8(b) displays that the frequencies of the WSW, W, and WNW winds for TL-30 were 17.2%, 16.8%, and 11.9%, respectively, with the prevailing wind switching from W to WSW. Compared with the wind rose at TL-60 (Fig. 7(b)), the frequency of the WSW wind increased by 3.3%, while the frequency of the WNW wind decreased by 2.7%. The impact on wind directions intensifies as the observation towers are located closer to the embankment. Fig. 8(d)shows that the WSW wind frequency at TL-0.5 reached 18.9%,whereas the WNW wind frequency fell to 11.7%. The results reveal that the embankment itself could significantly alter the NSWF directions within a distance of ten times the embankment height on the leeward side. Apparently, the embankment obstructed the winds perpendicular to its course and facilitated those parallel to its course.

    3.3. Variation in wind speed across the embankment

    Fig.9 exhibits the variations in MMWS at heights of 0.5,1.5,and 3.0 m along the six towers in the study period. The seasonal variations in all MMWSs were synchronous with their values,generally reaching their respective maximum values in February and falling to their minimum in September. However, the embankment clearly impacted the wind speeds, especially for the two towers close to the slope toe, namely, TR-0.5 and TL-0.5. The closer to the ground surface it is, the greater the embankment influence on the NSWF speed is.

    Fig.9. Variations in the MMWS at heights of(a)0.5 m,(b)1.5 m,and(c)3.0 m along the six towers in the study period.

    At the 0.5 m height, the MMWS values at TR-60, TL-60, TR-30,and TL-30 shown in Fig. 9(a) were close. Among the four towers,the MMWS values at TR-30 in the study period were approximately 0.2-0.3 m·s-1lower than those of the rest. At the observation location near its slope toe, the embankment significantly slowed the wind, particularly on the leeward side. On the windward side, the maximum MMWS at TR-0.5 was 1.8 m·s-1or 40%smaller than that of TR-60.However,on the leeward side,the peak MMWS at TL-0.5 was 2.2 m·s-1or 50% lower than that of TL-60.

    The embankment impact on the wind speeds appeared to decrease with an increasing observation height. At 1.5 m above the ground surface, Fig. 9(b) shows that the MMWS data at TR-60, TR-30, TL-60, and TL-30 almost overlapped throughout the entire study period. However, the differences between the two towers near the embankment slope toe (i.e., TR-0.5 and TL-0.5)and the rest(i.e.,TR-60,TR-30,TL-60,and TL-30)were still significant.For example,the MMWS values at TR-0.5 were 0.8-1.5 m·s-1smaller than those at the other four locations. When the observation height reached 3.0 m (Fig. 9(c)), the MMWS data differences at all six towers in the same month were within 0.7 m·s-1.

    Notably,at 0.5 and 1.5 m,the MMWS values at TR-0.5 were generally greater than those at TL-0.5, except for July and August. As seen from Fig. 5(c), the prevailing wind during the warm season switched from W to E, rendering TL-0.5 on the windward side and TR-0.5 on the leeward side, and hence, the MMWSs values were higher at TL-0.5 in these two summer months. However,the MMWS values at 3.0 m at TL-0.5 were slightly but consistently greater than those at TR-0.5 throughout the study period.The considerable differences in the MMWS at 0.5 and 1.5 m above the ground surface along the six towers clearly illustrated the embankment’s obstruction effect on the NSWF.

    Fig. 10 compares MMWS at a height of 4.5 m along TR-0.5 and TL-0.5 in the study period.At this height,the MMWS values at the two locations were very close throughout the study period. In the summertime, the MMWS values at the two locations were almost identical. In contrast, in the cold season, the MMWS values at TL-0.5 were marginally greater than those at TR-0.5. This result can be explained by the acceleration of wind flow when it passes over the embankment. The seasonal variations in the MMWS differences between the two locations can be attributed to the more consistent prevailing winds during the cold season and the more dispersed winds during the warm season in the study area.

    3.4. Wind speed profiles across the embankment

    As mentioned above,the NSWF speed is one of the main factors affecting the thermal performance of ACSs used in permafrost subgrade cooling. However, the conventional practice uses wind data measured at a given height,for example,10 m,from a weather station.A power law is typically applied to interpolate wind speeds at various heights above the ground surface for engineering applications because of its simplicity[45-47].The power law is described as follows:

    where y and v10are the wind speeds at heights of x and 10 m above the ground surface,respectively,and a is the power-law index(PLI).The value of the PLI is considered to be a constant that depends on ground roughness.

    Fig.10. Variations in the MMWS at a height of 4.5 m along TR-0.5 and TL-0.5 in the study period.

    The MMWS data in February and September from the six towers were analyzed to examine the impact of existing structures,such as an embankment,on the wind speed profiles.The v10values for February and September obtained from the nearby weather station were 8.6 and 3.4 m·s-1,respectively.Figs.11(a)-(d)depict the wind speed profiles at TR-60, TL-60, TR-30, and TL-30. At the four locations, the wind speed variations with increasing height could be well fitted using the power law.However,the PLI values varied considerably with the location and the season. Among the four observation towers, the PLI values in September varied between 0.14 and 0.18, whereas the PLI values in February were greater,ranging from 0.21 to 0.24.

    Fig. 11. Wind speed profiles at (a) TR-60, (b) TL-60, (c) TR-30, and (d) TL-30 in September and February and fitted curves using the power law. R2: correlation coefficient.

    Fig. 12. Wind speed profiles at (a) TR-0.5 and (b) TL-0.5 in September and February and fitted curves using the power law.

    Fig. 12 compares the wind speed profiles at TR-0.5 and TL-0.5,which are closest to the embankment and therefore directly affect the ACSs’ performance. The results show that even though the wind speeds near the embankment slope toe were considerably affected (refer to Figs. 9 and 10), the wind speed profiles can still be described by the power law. However, the PLI values at the two locations were significantly greater than those at the four towers farther away from the embankment as discussed above.The PLI values at TR-0.5 and TL-0.5 reached 0.28 and 0.23 in September and jumped to 0.40 and 0.38 in February, respectively.These substantially greater PLI values at the two towers clearly reflect the embankment’s accelerating effect on the ambient wind flow.

    To quantify the spatial and temporal variations in the PLI values of wind speed profiles across the embankment, Fig. 13 shows the fitted PLI values at TR-60,TR-30,TL-60,and TL-30 during the entire study period. The figure shows the seasonal variations in the PLIs are synchronous and significant. The PLI values are smaller for lower wind speeds in the warm season and larger for higher wind speeds in the cold season,consistent with observations from urban terrains[45-47].Additionally,the PLI values on the windward side towers are slightly greater than those on the leeward side at the same distance from the embankment in the same month.

    The PLI values at 60 m away from the embankment in the cold season ranged from 0.18 to 0.23,with an average of 0.21,whereas the values in the warm season varied from 0.13 to 0.22, with an average of 0.16. At 30 m away from the embankment, the PLI values increase slightly,with an average PLI rising to 0.25 in the cold season and to 0.18 in the warm season. Among the four locations,through the year, only a few PLI values were being smaller than 0.16, which was typically used in previous numerical evaluations of the long-term thermal performances of ACSs used in engineered infrastructures on the QTP [16,22,25,27,30,35-43].

    Fig. 13. NSWF PLI at TR-60, TR-30, TL-60, and TL-30 in the study period.

    4. Discussion

    4.1. Significance of the NSWF speed profile to the long-term thermal performances of ACSs

    A reliable evaluation of the long-term thermal performances of ACSs is essential for the design and maintenance of engineered infrastructures built in a permafrost environment. In numerical evaluations, the NSWF speed is a crucial boundary condition for the heat exchange between the ACSs and ambient air. As mentioned previously,the wind speed observed at a 10 m height is generally used to determine the NSWF speed according to a power law,and the two PLI values(i.e.,0.12 and 0.16)were typically used in previous numerical evaluations of the long-term thermal performances of ACSs used in engineered infrastructures on the QTP[16,22,25,27,30,35-43]. However, the field-observed data in this study showed that the PLI not only varied considerably with wind speed but also substantially exceeded the assumed values(Fig.11).

    The mean annual wind speed (MAWS) values at three heights(0.5,1.5,and 3.0 m)along TR-60 and TL-60 were analyzed to quantify the PLI values on an annual basis in the study period,as shown in Fig. 14. The fitted results showed that the annual average PLI value was 0.19.In the figure,the power-law predicted wind speed profiles with v10from the reference meteorological station and PLI values of 0.12 and 0.16 were also plotted for comparison. The NSWF speed values estimated using previous PLI values are considerably greater than the field-observed values. For instance, the wind speed is overestimated by 12%-25% and 5%-13% at 0.5 and 1.5 m above the ground surface, respectively. At a height of 3.0 m, the predicted wind speed with a PLI of 0.12 would still be overvalued by 10%.

    In permafrost regions on the QTP,most of the railway and roadway embankments are generally under 5 m in height [48,49],meaning that most numerical evaluations of the long-term thermal performances of embankments constructed with ACSs, including CRLs, VDs, and thermosyphons, applied overrated NSWF speeds;these overrated speeds led to an overestimated ACS heatexchange efficiency and hence an unconservative assessment of their long-term permafrost cooling effects.For example,an overestimated NSWF speed would result in an overrated airflow speed within VDs and an inflated heat extraction rate for thermosyphons.More importantly, for CRLs used in railway and roadway embankments, an overestimated NSWF speed may not only erroneously determine the convection type, namely, wind-driven versus buoyancy-driven, but also, the convection intensity within the CRLs.The overvalued heat-exchange efficiency of these ACSs would lead to insufficient designs of these ACSs, which may jeopardize the long-term thermal stability of the permafrost subgrade.

    Fig. 14. Observed mean annual wind speed profiles along (a) TR-60 and (b) TL-60 and the power-law predicted profiles with PLI values of 0.12 and 0.16.

    4.2. Significance of the NSWF direction on the long-term thermal performances of ACSs

    The NSWF direction is also a vital factor affecting the long-term thermal performances of the CRLs and VDs used in linear transportation infrastructures. Taking a roadway embankment with VDs as an example,the NSWF direction,together with the embankment route azimuth,determines the angle between the NSWF and the duct axis, and hence, determines the relationship between the environmental wind speed and the duct airflow speed.

    The field-observed data on two days, namely, August 15 and December 15 of 2018, were selected to represent the duct airflow conditions in warm and cold seasons, respectively. Fig. 15 shows the half-hour average environmental wind speed and airflow speed within a duct in the experimental expressway embankment. The duct has an inner diameter of 0.4 m and is buried at a height of 1.0 m above the ground surface. The environmental wind speed data were collected from the weather station at a height of 2 m,and the airflow speed data were measured by an ultrasonic anemometer installed along the duct axis. Fig. 16 illustrates the half-hour average environmental wind direction on the two selected days. On December 15, the environmental wind direction was relatively persistent, ranging from 242° to 332°, with a daily average of 282°, which is roughly perpendicular to the embankment route (197°). On this day, there was a clear positive correlation between the environmental wind speed and the duct airflow speed, with the latter reaching a maximum of 5.3 m·s-1. On the other hand,the environmental wind direction on August 15 varied greatly, ranging from 4° to 360°. For instance, the environmental wind direction at 18:00 was nearly perpendicular to the embankment route or almost parallel to the duct axis, resulting in a duct airflow speed of 2.1 m·s-1, which was about half of the environmental wind speed (4.3 m·s-1). However,the environmental wind direction at 21:30 was approximately parallel to the embankment route or almost perpendicular to the duct axis, generating a duct airflow speed of only 0.2 m·s-1,which accounts for approximately 1/30 of the environmental wind speed (6.8 m·s-1).

    These results demonstrate that both the environmental wind direction and speed need to be considered in evaluating the long-term thermal performance of linear transportation infrastructure with ACSs. An overestimation of the ACS heat-exchange efficiency would likely occur otherwise. However, most existing numerical simulations[16,22,25,27,30,35-43]ignored wind direction variations due to their complexities. Fortunately, the wind direction on the QTP in wintertime is relatively persistent (Fig.5(b)), which instills confidence in assuming a predominant wind direction in numerical simulations.Regarding engineering applications, the angle between the actual environmental wind direction and the route of linear transportation infrastructure should be considered in the design of VDs. For example, if the prevailing wind direction is close to the route of roadway embankment, thermosyphons or CRLs would be a better choice than VDs for permafrost subgrade cooling. Otherwise, the VDs should be aligned with the prevailing wind direction in the cold season as much as possible.

    5. Conclusions

    The NSWF is crucial in determining the heat-exchange efficiency of the ACSs used for permafrost subgrade cooling.This study focuses on NSWF across linear transportation infrastructures in a high-altitude permafrost environment. A yearlong field observation of the NSWF across an embankment located in the interior of the QTP,including the observation system layout,measurement methods, and collected data, was described. Based on the fieldobserved data, the environmental air and wind characteristics,wind direction and speed from both the windward and leeward sides of the embankment, and their variations with varying distances from the embankment were analyzed. Furthermore, the significance of the NSWF to the thermal performance of ACSs,including CRLs, VDs, and thermosyphons, was discussed in the context of linear transportation infrastructure. The following conclusions were obtained from the field observations, analyses,and discussion:

    (1) The wind direction in the cold season at the study site was relatively persistent. The occurrence frequencies of W, WNW,and WSW winds each accounted for more than 18.5%, with a sum of 65.5%. In contrast, the wind direction in the warm season was considerably more dispersed, with no single wind direction exceeding 15% in its occurrence frequency. The wind speed in the cold winter was substantially higher than that in the warm summer. The ACSs can significantly benefit from cold winds in the cold season to preserve the permafrost subgrade.

    (2)The embankment’s existence did not impact the wind direction on the windward side. However, the wind perpendicular to the embankment course on the leeward side was slowed,and that parallel to the embankment course was accelerated. The embankment obstructed the wind and resulted in a considerable decrease in wind speed at elevations lower than its height on both the windward and leeward sides. Above the embankment crest, the wind speed on the leeward side was enhanced. The NSWF speed could recover at approximately 15-20 times the embankment height.

    Fig. 15. Environmental wind speed versus duct airflow speed on (a) December 15 and (b) August 15 of 2018.

    Fig. 16. Environmental wind direction on (a) December 15 and (b) August 15 of 2018.

    (3) The wind speed profiles across the embankment can be fitted well with a power law. The PLI values changed considerably with the seasonally varied wind speed and distance from the embankment.At locations 20 times the embankment height away,the average PLI was 0.21 in wintertime and 0.16 in summertime.On an annual basis,the PLI value was 0.19,which was significantly larger than the values (e.g., 0.12 and 0.16) used in previous numerical simulations. The underestimation of the PLI in these numerical simulations would lead to an overestimated heatexchange efficiency of the ACSs and, therefore, an unconservative design for permafrost subgrade cooling.

    (4)The wind direction is also a crucial factor affecting the thermal performances of the CRLs and VDs commonly used in linear transportation infrastructures. Ignoring the wind direction variations would also yield an overrated heat-exchange efficiency and an unconservative design.It is highly recommended that the angle between the prevailing environment wind direction and the linear transportation infrastructure route be considered in the design of VDs.

    The presented results are obtained from one case, and more field observations on different settings are needed.Factors including ground surface conditions, landforms, the geometry of embankment have a significant impact on the characteristics of the NSWF across an expressway embankment and should be considered when applying the results or selecting future field observation sites.

    Acknowledgment

    This study is supported by the National Natural Science Foundation of China (41630636 and 41772325) and China’s Second Tibetan Plateau Scientific Expedition and Research(2019QZKK0905).

    Compliance with ethics guidelines

    Yanhu Mu,Wei Ma,Zhaohui(Joey)Yang,Xiaolin Li,Kun Zhang,and Yuncheng Mao declare that they have no conflict of interest or financial conflicts to disclose.

    国产精品蜜桃在线观看| 国产精品一及| 日韩欧美在线乱码| 亚洲激情五月婷婷啪啪| 日韩高清综合在线| 久久人人爽人人片av| 床上黄色一级片| 午夜福利高清视频| 国产亚洲最大av| 久久久久久久久久久丰满| 午夜免费激情av| 国产精品综合久久久久久久免费| 欧美潮喷喷水| 麻豆乱淫一区二区| 亚洲成人精品中文字幕电影| 老司机福利观看| 国产av在哪里看| 久久久a久久爽久久v久久| 国内揄拍国产精品人妻在线| 国产精品一二三区在线看| 少妇的逼好多水| 国产不卡一卡二| 18+在线观看网站| 极品教师在线视频| 婷婷六月久久综合丁香| 欧美成人免费av一区二区三区| 久久婷婷人人爽人人干人人爱| 亚洲在久久综合| 国产中年淑女户外野战色| av又黄又爽大尺度在线免费看 | 亚洲婷婷狠狠爱综合网| 少妇丰满av| 久久久久免费精品人妻一区二区| 午夜福利视频1000在线观看| 精品熟女少妇av免费看| 久久草成人影院| 欧美性猛交黑人性爽| 国产午夜精品一二区理论片| 婷婷六月久久综合丁香| 亚洲,欧美,日韩| 国产成人a区在线观看| 国产精品久久久久久久久免| 欧美丝袜亚洲另类| 一本久久精品| 国产精品国产高清国产av| 99在线视频只有这里精品首页| 久久韩国三级中文字幕| 一区二区三区四区激情视频| 国产一区二区在线av高清观看| 只有这里有精品99| 国产精品综合久久久久久久免费| 久久精品国产亚洲av天美| 精品人妻视频免费看| 少妇的逼水好多| 简卡轻食公司| 国产成人91sexporn| 欧美日韩国产亚洲二区| 久热久热在线精品观看| 亚洲最大成人手机在线| 97人妻精品一区二区三区麻豆| 亚洲在久久综合| 久久亚洲国产成人精品v| 成年女人永久免费观看视频| 成人综合一区亚洲| 久久久久性生活片| 日韩成人伦理影院| 中文在线观看免费www的网站| 国产精品一区二区三区四区免费观看| 国产精品伦人一区二区| 欧美激情在线99| 精品国产三级普通话版| 国产高潮美女av| 大香蕉久久网| 久久久国产成人精品二区| 超碰av人人做人人爽久久| 久久久久久久久大av| 三级国产精品欧美在线观看| 特级一级黄色大片| 午夜精品一区二区三区免费看| 免费不卡的大黄色大毛片视频在线观看 | 国产一区有黄有色的免费视频 | 三级男女做爰猛烈吃奶摸视频| 一级黄色大片毛片| 国产精品电影一区二区三区| 成人一区二区视频在线观看| 高清在线视频一区二区三区 | 亚洲av电影在线观看一区二区三区 | 视频中文字幕在线观看| 国产亚洲5aaaaa淫片| 男的添女的下面高潮视频| 精品酒店卫生间| 欧美+日韩+精品| 国产一区有黄有色的免费视频 | 亚洲国产欧美在线一区| 国产伦精品一区二区三区视频9| 老司机影院毛片| 久久精品国产亚洲av涩爱| 国产精品99久久久久久久久| 国产午夜精品论理片| 成人鲁丝片一二三区免费| 亚洲四区av| 99久久无色码亚洲精品果冻| 爱豆传媒免费全集在线观看| eeuss影院久久| 一区二区三区高清视频在线| 国语自产精品视频在线第100页| 久久亚洲国产成人精品v| 午夜激情欧美在线| 最后的刺客免费高清国语| 伦理电影大哥的女人| 日本爱情动作片www.在线观看| 色吧在线观看| 亚洲国产高清在线一区二区三| 久久久a久久爽久久v久久| 欧美高清性xxxxhd video| 久久国内精品自在自线图片| 一区二区三区四区激情视频| 噜噜噜噜噜久久久久久91| 久久人人爽人人爽人人片va| 国产日韩欧美在线精品| 日韩av在线大香蕉| 性插视频无遮挡在线免费观看| 久久国内精品自在自线图片| 国产黄色视频一区二区在线观看 | 国产亚洲精品久久久com| 久久久久久久久久成人| 日本黄色片子视频| 国产亚洲一区二区精品| 69av精品久久久久久| 国产精品,欧美在线| 免费av观看视频| 国产亚洲一区二区精品| 国产一区二区三区av在线| 狂野欧美白嫩少妇大欣赏| 国产一级毛片在线| 在现免费观看毛片| 在现免费观看毛片| 好男人在线观看高清免费视频| 成人av在线播放网站| 久久人人爽人人片av| 少妇人妻精品综合一区二区| 午夜老司机福利剧场| 最近中文字幕高清免费大全6| 欧美性猛交╳xxx乱大交人| 一级毛片电影观看 | 久久精品久久久久久噜噜老黄 | 99热6这里只有精品| 国产人妻一区二区三区在| 日韩欧美精品免费久久| 青青草视频在线视频观看| 国产国拍精品亚洲av在线观看| 99热网站在线观看| 久久亚洲精品不卡| 国产av在哪里看| 自拍偷自拍亚洲精品老妇| 亚洲激情五月婷婷啪啪| 91午夜精品亚洲一区二区三区| 亚洲第一区二区三区不卡| 美女脱内裤让男人舔精品视频| 成人三级黄色视频| 免费观看a级毛片全部| 亚洲国产精品专区欧美| 一个人观看的视频www高清免费观看| 韩国av在线不卡| 国产精品野战在线观看| 高清在线视频一区二区三区 | 99视频精品全部免费 在线| 亚洲真实伦在线观看| 国产成人福利小说| 亚洲va在线va天堂va国产| 免费黄网站久久成人精品| 一个人看视频在线观看www免费| 午夜福利网站1000一区二区三区| 国产精品人妻久久久影院| 久久99蜜桃精品久久| 国产亚洲av片在线观看秒播厂 | 亚洲欧美成人精品一区二区| 青春草视频在线免费观看| 啦啦啦韩国在线观看视频| 天堂av国产一区二区熟女人妻| 免费人成在线观看视频色| 久久久久久久久久成人| or卡值多少钱| 婷婷色麻豆天堂久久 | 久久久久久久久久黄片| 精品午夜福利在线看| 亚洲精品国产成人久久av| 七月丁香在线播放| 亚洲精品国产成人久久av| 热99在线观看视频| 日本免费一区二区三区高清不卡| 亚洲国产高清在线一区二区三| h日本视频在线播放| 大香蕉久久网| 美女脱内裤让男人舔精品视频| 天堂中文最新版在线下载 | 非洲黑人性xxxx精品又粗又长| 成人性生交大片免费视频hd| 久久人妻av系列| 精品国内亚洲2022精品成人| 日韩欧美在线乱码| 国产激情偷乱视频一区二区| 久久鲁丝午夜福利片| 亚洲第一区二区三区不卡| 91久久精品电影网| 午夜福利在线观看免费完整高清在| 2021天堂中文幕一二区在线观| 最新中文字幕久久久久| 欧美一区二区国产精品久久精品| 色哟哟·www| 菩萨蛮人人尽说江南好唐韦庄 | 99久久精品热视频| 亚洲伊人久久精品综合 | 亚洲五月天丁香| 99久久精品国产国产毛片| 国产精品爽爽va在线观看网站| 国产一区二区在线观看日韩| 精华霜和精华液先用哪个| 国产精品av视频在线免费观看| 日韩欧美三级三区| 亚洲国产欧美人成| 亚洲精品国产成人久久av| 午夜亚洲福利在线播放| 亚洲高清免费不卡视频| 99久久精品国产国产毛片| 国产黄色视频一区二区在线观看 | 纵有疾风起免费观看全集完整版 | 欧美潮喷喷水| 麻豆成人av视频| 黄色日韩在线| 永久免费av网站大全| 日本午夜av视频| 精品久久久久久电影网 | 久久久久网色| 亚洲人成网站在线播| 赤兔流量卡办理| 在线天堂最新版资源| 一个人看的www免费观看视频| 亚洲天堂国产精品一区在线| 国产精品一区www在线观看| 乱系列少妇在线播放| 午夜激情福利司机影院| 精品久久久噜噜| 亚洲图色成人| 亚洲丝袜综合中文字幕| 国产伦一二天堂av在线观看| 特级一级黄色大片| 欧美精品一区二区大全| 人人妻人人澡欧美一区二区| 男女下面进入的视频免费午夜| 村上凉子中文字幕在线| 国产一区亚洲一区在线观看| 国国产精品蜜臀av免费| 中国国产av一级| 好男人在线观看高清免费视频| 国产成人aa在线观看| 高清毛片免费看| 老女人水多毛片| 国产亚洲精品av在线| 国内精品一区二区在线观看| 你懂的网址亚洲精品在线观看 | 一个人免费在线观看电影| 国产精品久久视频播放| 国产真实乱freesex| 欧美成人免费av一区二区三区| 男女那种视频在线观看| 久久久久免费精品人妻一区二区| 一级黄色大片毛片| 神马国产精品三级电影在线观看| 最新中文字幕久久久久| 久久精品影院6| 一级二级三级毛片免费看| 丰满少妇做爰视频| 水蜜桃什么品种好| 欧美成人午夜免费资源| 精品久久久噜噜| 欧美激情国产日韩精品一区| 欧美三级亚洲精品| 久久热精品热| 搡女人真爽免费视频火全软件| 色哟哟·www| 亚洲五月天丁香| 久久人人爽人人片av| 日韩欧美国产在线观看| 午夜免费男女啪啪视频观看| 国产精品99久久久久久久久| 全区人妻精品视频| 永久网站在线| 大话2 男鬼变身卡| 少妇丰满av| 国产精品日韩av在线免费观看| 最近视频中文字幕2019在线8| 我要看日韩黄色一级片| 国产精品女同一区二区软件| 国产亚洲精品久久久com| 亚洲天堂国产精品一区在线| 亚洲成色77777| 色综合亚洲欧美另类图片| 欧美激情在线99| 亚洲国产成人一精品久久久| 九九爱精品视频在线观看| 特大巨黑吊av在线直播| 国产成人午夜福利电影在线观看| 亚洲av日韩在线播放| 亚洲精品色激情综合| 午夜福利在线观看吧| 久久精品国产亚洲网站| 中文字幕亚洲精品专区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品亚洲一区二区| 日韩大片免费观看网站 | 国产老妇伦熟女老妇高清| 精品无人区乱码1区二区| 色尼玛亚洲综合影院| 成年版毛片免费区| 亚洲精品,欧美精品| 在线免费观看的www视频| 性色avwww在线观看| 亚洲在线观看片| av.在线天堂| 国产av在哪里看| 欧美又色又爽又黄视频| 97超视频在线观看视频| 爱豆传媒免费全集在线观看| 少妇猛男粗大的猛烈进出视频 | 夫妻性生交免费视频一级片| 亚洲乱码一区二区免费版| 亚洲内射少妇av| 午夜福利视频1000在线观看| 免费黄网站久久成人精品| 久久久久免费精品人妻一区二区| 1000部很黄的大片| 日韩制服骚丝袜av| 国产高清有码在线观看视频| 亚洲欧美中文字幕日韩二区| 99久久人妻综合| 欧美人与善性xxx| 神马国产精品三级电影在线观看| 欧美精品国产亚洲| 免费观看性生交大片5| 日韩一区二区三区影片| 啦啦啦韩国在线观看视频| 韩国av在线不卡| 男女视频在线观看网站免费| 99热全是精品| 伦理电影大哥的女人| 男女下面进入的视频免费午夜| 精品免费久久久久久久清纯| 欧美成人午夜免费资源| 成人亚洲欧美一区二区av| 一边亲一边摸免费视频| 简卡轻食公司| 亚洲久久久久久中文字幕| 欧美高清成人免费视频www| 国产午夜精品久久久久久一区二区三区| 成人漫画全彩无遮挡| 日本一二三区视频观看| 欧美极品一区二区三区四区| 国产又色又爽无遮挡免| 亚洲精品aⅴ在线观看| av免费观看日本| 亚洲精品国产成人久久av| 免费黄色在线免费观看| 精品久久久久久久人妻蜜臀av| 亚洲国产精品合色在线| 国产精品av视频在线免费观看| 高清日韩中文字幕在线| 美女内射精品一级片tv| av黄色大香蕉| 亚洲经典国产精华液单| 熟女电影av网| 亚洲无线观看免费| 免费无遮挡裸体视频| 99视频精品全部免费 在线| 免费播放大片免费观看视频在线观看 | 高清视频免费观看一区二区 | 久久精品91蜜桃| 我的女老师完整版在线观看| 国产成人a∨麻豆精品| 成年av动漫网址| 色综合色国产| eeuss影院久久| 国产免费视频播放在线视频 | 麻豆av噜噜一区二区三区| 日日摸夜夜添夜夜添av毛片| 中文字幕av成人在线电影| 97在线视频观看| 99热6这里只有精品| 国产精品不卡视频一区二区| 亚洲av中文av极速乱| 久久久精品欧美日韩精品| 日韩一区二区三区影片| 中文资源天堂在线| 亚洲自拍偷在线| 狠狠狠狠99中文字幕| 欧美极品一区二区三区四区| 国产精品.久久久| 99久久精品热视频| 18禁裸乳无遮挡免费网站照片| 乱码一卡2卡4卡精品| 亚洲激情五月婷婷啪啪| 国产视频首页在线观看| 一级av片app| www.av在线官网国产| 国产私拍福利视频在线观看| 最近2019中文字幕mv第一页| 黄色配什么色好看| 蜜臀久久99精品久久宅男| av视频在线观看入口| 高清毛片免费看| 午夜福利网站1000一区二区三区| 日本免费a在线| 在线播放国产精品三级| 国语自产精品视频在线第100页| 国产精品一区二区三区四区久久| 性插视频无遮挡在线免费观看| 非洲黑人性xxxx精品又粗又长| 伊人久久精品亚洲午夜| 亚洲乱码一区二区免费版| 欧美xxxx黑人xx丫x性爽| 国产欧美日韩精品一区二区| 哪个播放器可以免费观看大片| 中国国产av一级| 欧美另类亚洲清纯唯美| 亚洲国产高清在线一区二区三| 激情 狠狠 欧美| 成年av动漫网址| 老司机影院毛片| 色吧在线观看| 国产亚洲一区二区精品| 在线观看66精品国产| 国产亚洲午夜精品一区二区久久 | 久久精品国产亚洲av涩爱| 国产欧美另类精品又又久久亚洲欧美| 国产精品美女特级片免费视频播放器| 国产美女午夜福利| 久久久久免费精品人妻一区二区| 只有这里有精品99| 日本免费a在线| 久久精品国产鲁丝片午夜精品| av福利片在线观看| 国产成年人精品一区二区| 免费黄色在线免费观看| 日韩,欧美,国产一区二区三区 | 人妻系列 视频| 波野结衣二区三区在线| 有码 亚洲区| 婷婷色av中文字幕| 久久久久国产网址| 菩萨蛮人人尽说江南好唐韦庄 | 日本一本二区三区精品| 2022亚洲国产成人精品| 午夜福利在线观看吧| 国产在线男女| 国产精品av视频在线免费观看| 久久精品国产亚洲网站| 美女脱内裤让男人舔精品视频| 精品久久久久久久久av| 国产老妇伦熟女老妇高清| 国产午夜精品一二区理论片| 久久久久久久久大av| 韩国高清视频一区二区三区| 精品不卡国产一区二区三区| av在线天堂中文字幕| 久久久久久久久久久丰满| 国产精品1区2区在线观看.| av又黄又爽大尺度在线免费看 | 久久久a久久爽久久v久久| 天天躁夜夜躁狠狠久久av| 在线观看美女被高潮喷水网站| 内地一区二区视频在线| 国产一级毛片七仙女欲春2| 日韩制服骚丝袜av| 99热6这里只有精品| 小说图片视频综合网站| 最近视频中文字幕2019在线8| 午夜久久久久精精品| 午夜福利在线在线| 一边摸一边抽搐一进一小说| 精品久久久久久久人妻蜜臀av| 精华霜和精华液先用哪个| 国产精品美女特级片免费视频播放器| 国产精品1区2区在线观看.| 久久综合国产亚洲精品| 免费人成在线观看视频色| 国产一区二区在线观看日韩| 麻豆成人午夜福利视频| 伊人久久精品亚洲午夜| 中文精品一卡2卡3卡4更新| 国产一级毛片在线| 国产免费一级a男人的天堂| 狠狠狠狠99中文字幕| 啦啦啦韩国在线观看视频| 男人舔女人下体高潮全视频| 在线免费十八禁| 天天一区二区日本电影三级| 久久人人爽人人爽人人片va| 欧美极品一区二区三区四区| 国产爱豆传媒在线观看| 嫩草影院入口| 日韩成人伦理影院| 亚洲18禁久久av| 国产一级毛片在线| 91久久精品电影网| 六月丁香七月| 国产激情偷乱视频一区二区| 自拍偷自拍亚洲精品老妇| 国产精品女同一区二区软件| 欧美极品一区二区三区四区| 边亲边吃奶的免费视频| videossex国产| 色5月婷婷丁香| 国产精品三级大全| 亚洲美女搞黄在线观看| 国产极品天堂在线| av国产久精品久网站免费入址| 一级爰片在线观看| 午夜福利在线在线| 中文乱码字字幕精品一区二区三区 | 我要搜黄色片| 欧美另类亚洲清纯唯美| 寂寞人妻少妇视频99o| 色综合站精品国产| 啦啦啦观看免费观看视频高清| 久久久成人免费电影| 免费看日本二区| 亚洲乱码一区二区免费版| 精品人妻熟女av久视频| www.av在线官网国产| 成人漫画全彩无遮挡| 久久精品熟女亚洲av麻豆精品 | 国产免费又黄又爽又色| 亚洲精品久久久久久婷婷小说 | 蜜桃亚洲精品一区二区三区| 亚洲精品aⅴ在线观看| 欧美zozozo另类| 一个人看视频在线观看www免费| 九色成人免费人妻av| 久久精品国产亚洲av天美| 亚洲av日韩在线播放| 观看美女的网站| 成人美女网站在线观看视频| 国产av在哪里看| 国产精品国产三级国产专区5o | 毛片一级片免费看久久久久| 日韩欧美在线乱码| 只有这里有精品99| 精品欧美国产一区二区三| 99久久九九国产精品国产免费| 男女国产视频网站| 亚洲在线观看片| 国产私拍福利视频在线观看| 精品99又大又爽又粗少妇毛片| 国产极品天堂在线| 1024手机看黄色片| 夜夜看夜夜爽夜夜摸| 亚洲五月天丁香| 两性午夜刺激爽爽歪歪视频在线观看| 综合色av麻豆| 青春草视频在线免费观看| 亚洲一级一片aⅴ在线观看| 欧美成人精品欧美一级黄| 亚洲精品,欧美精品| 亚洲欧美一区二区三区国产| 精品一区二区三区人妻视频| 亚洲婷婷狠狠爱综合网| 伦理电影大哥的女人| 天美传媒精品一区二区| 尾随美女入室| 丰满乱子伦码专区| 日韩一区二区视频免费看| 亚洲av免费在线观看| 亚洲天堂国产精品一区在线| 国产一区二区在线观看日韩| 一级av片app| 只有这里有精品99| 久久精品国产亚洲网站| 日日啪夜夜撸| 欧美xxxx黑人xx丫x性爽| 波野结衣二区三区在线| 欧美激情在线99| 看免费成人av毛片| 欧美性感艳星| 亚洲精华国产精华液的使用体验| 久久久久久国产a免费观看| 啦啦啦观看免费观看视频高清| 国产高清视频在线观看网站| 十八禁国产超污无遮挡网站| 欧美又色又爽又黄视频| 久久99热6这里只有精品| 日本五十路高清| 又爽又黄无遮挡网站| 免费观看a级毛片全部| 如何舔出高潮| 97人妻精品一区二区三区麻豆| 免费看日本二区| 最新中文字幕久久久久| 免费黄网站久久成人精品| www.色视频.com| 精品人妻一区二区三区麻豆| 国产精品av视频在线免费观看| 久久久久久久国产电影| 午夜a级毛片| 高清在线视频一区二区三区 | 国产视频首页在线观看| 国产一区有黄有色的免费视频 | av.在线天堂| 大香蕉97超碰在线| 两个人的视频大全免费| 成人性生交大片免费视频hd| 久久热精品热| 国产人妻一区二区三区在| 日韩av不卡免费在线播放| 国产黄色视频一区二区在线观看 | 麻豆成人午夜福利视频|