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

    Empirical Equations of P-Wave Velocity in the Shallow and Semi-Deep Sea Sediments from the South China Sea

    2021-06-25 11:13:08LIGuanbaoHOUZhengyuWANGJingqiangKANGuangmingandLIUBaohua
    Journal of Ocean University of China 2021年3期

    LI Guanbao, HOU Zhengyu, WANG Jingqiang, KAN Guangming, and LIU Baohua, 7)

    Empirical Equations of P-Wave Velocity in the Shallow and Semi-Deep Sea Sediments from the South China Sea

    LI Guanbao1), 6), HOU Zhengyu2), 3), 4), 5), *, WANG Jingqiang1), 6), KAN Guangming1), 6), and LIU Baohua1), 6), 7)

    1),,,266061,2),,,,510301,3),,511458,4),,100190,5),511458,6),,266237,7),,266237,

    Three-hundred and thirty-one sediment cores, including six sediment types (clayey- and sandy-silt, silty- and sandy-clay, clayey- and silty-sand) were collected from the shallow and semi-deep seas of the South China Sea, and the P-wave velocities and physical properties of core sediments were measured under standard laboratory conditions. To eliminate the influence of environmental factors, the empirical sound speed ratio equations were established. Compared with several equations from literature, the porosity and wet bulk density empirical equations established in this paper agree well with Richardson and Briggs (2004)’sequations, which implies that our empirical equations can be used in the similar region of world’s oceans under certain conditions and will be useful in areas lacking first-hand P-wave speed data. However, the mean grain size equations established in this study, similar to the previous studies, have low accuracy, which may be due to the different particle arrangements and degrees of compaction in sediments. The results also show that for different sediment types, the equation based on all sediment data is in good agreement with the measured data in the study area, as there are both siliciclastic and carbonate sediments on the studied seabed. It is suggested that appropriate empirical equations should be selected according to sediment types and sedimentary environment in future works, and the empirical equation of porosity or the two-parameter equation of porosity and grain size should be preferred.

    seafloor sediment; geo-acoustic property; sound speed; empirical equation; South China Sea

    1 Introduction

    In the waveguide of the shallow and semi-deep sea (SSD),the seabed is a dominant factor controlling underwater acoustic propagation by bottom reflection and scattering losses (Zhou and Zhang, 2010). Sediment geo-acoustic pa- rameters, such as P-wave speeds, and their relationships with physical properties are important in many marine- geophysical applications (Liu., 2019), such as seismic exploration, seismic inversion and marine engineering. The geo-acoustic properties of SSD sediments are affected by many factors, such as sediment source, ocean current, water column and benthic organisms. Successful estimation of geo-acoustic properties in seabed sediments, such as the P-wave speed (hereafter called sound speed or SS), has been the goal of many studies involving inverse techniques (Zou and Kan, 2011; Suzuki., 2015; Mi- chalopoulou and Gerstoft, 2019). Core sediments are the primary samples for obtaining geo-acoustic parameters, which provide valuable insight into physical properties of the seabed (Konno., 2015; Ballard., 2019).

    Both empirical and physical models have been used to predict the relationships between seafloor geo-acoustic parameters and physical properties; however, the former has significant advantages since some of parameters required by physical models are either very difficult or nearlyimpossible to be obtained (Jackson and Richardson, 2007). To acquire the accurate sound speed data, many researchers have established several empirical equations based on the sediment cores from different regions. In the last century, Hamilton and co-workers (Hamilton, 1970; Hamilton and Bachman, 1982; Bachman, 1989) studied the acoustic characteristics of marine sediment worldwide and established empirical relationships between sound speed and physical properties for various topographical regions (con- tinental terrace, abyssal hill, and abyssal plain). Consi- derable efforts have been made to develop empirical equations (Richardson., 2002; Goff., 2004; Endler., 2016) that can be used to predict the acoustic behavior of marine sediments. Few studies concerning the geo-acoustic and physical properties of sediments in the South China Sea (SCS) have been published. Basic characteristics and general distribution of sediments in the southeast offshore sea area off the Hainan Island were analyzed by Lu. (2008), and the regional empirical equations for sediments were established based on the data. Wang. (2018) used thesystem to measure seafloor sediment acoustic properties in SCS, and com- pared the results with several empirical equations, and found that the measured data were closed to the results calculated by the equation of Richardson and Briggs (2004). Hou. (2015, 2018a, 2018b) investigated the vertical and planar distribution of acoustic parameters of sediments in the southern and central areas of SCS and studied the effects of ocean currents and sediment source on the sound speed of sediments. Geo-acoustic parameters are influenced by multiple factors, and different seabed environments have different sediment acoustic properties. Due to the limited number of sediment samples and no standard measurement methods, the empirical equations for SCS given in many of these studies, which are based on variety of techniques (different measurement methods or frequencies), are limited in practical application. In this study, the same measurement method and frequency were used to measure the geo-acoustic parameters for more than 300 different types of core sediment samples collected in the SSD of the SCS at the same batch. Empirical equations between sediment geo-acoustic and physical parameters were established and compared with several equations from the literature.

    2 Materials and Methods

    Sediment cores were collected by a gravity corer in the SSD area in the north and west of SCS (water depth less than 1500m), among which the shortest sediment core was only 25cm long, and the longest core was 3m long. There was a plexiglass hard tube in the gravity corer and the two ends of sediment cores were sealed, which was conducive to the preservation of moisture and transport of cores. The sediment acoustic parameters were measured by the coaxial differential distance method under the standard laboratory conditions (23℃, salinity 35, and atmosphe- ric pressure). Details about the measurement method are described in Hou. (2018a, 2018b). The transmitting waveform was produced from a signal generator, and the waveform was received by the receiving transducer and recorded by a digital recorder at a sample rate of 10MHz. The long cores are cut into sections every 25cm and each section is acoustically measured. The transducer frequency was 100kHz, and the compressional sound speed was cal- culated by using the sediment length and the propagation time of the acoustic signals. The speed in pore-water was reported under the standard laboratory conditions favored by Hamilton (23℃, salinity 35, and atmospheric pressure). These conditions correspond to a water sound speed of 1529.4ms?1and a water density of 1024kgm?3, as calculated by the Del Grosso equation (1974). The sound speed ratio (SSR) is defined as the ratio of the measured sediment sound speed to the pore-water sound speed at same temperature, salinity, and pressure.

    All sections were sub-sampled at their tops and ends, and additional samples were taken at depths with special interests. The porosity and water content were determined from the weight loss of sediments. Grain density was determined using a pycnometer. Sediment bulk density is calculated from the porosity and densities of pore water and sediment particles (Jackson and Richardson, 2007). The sediment particle size parameters were analyzed by the British Mastersizer 2000 laser diffraction particle size analyzer, and according to the percentage of each fraction in the sediment, the sediment type can be determined.

    3 Results

    The sediments collected in this study consisted of clay- ey-silt, sandy-silt, silty-clay, sandy-clay, clayey-sand and silty-sand (based on Shepard ternary diagrams, 1954). The geo-acoustic and physical parameters measured in this study included sound speed (SS), wet bulk density (ρ), meangrain size (Mz), and porosity (n). Three-hundred and thirty- one samples were examined in the study area, with 135 silty-clay, 131 clayey-silt, 37 silty-sand, 15 clayey-sand, 8 sandy-silt, and 5 sandy-clay. The statistical acoustic and physical parameters for different sediment types are listed in Table 1. The results show that the SS of sandy sediments (clayey-sand and silty-sand) was the highest, with the average larger than 1600ms?1, and the sandy sediments had higher wet bulk densities, particle sizes, and lower porosity. The lowest SS occurred in the clayey sediments (silty-clay and sandy-clay), with the average of about 1500ms?1, and the clayey sediments had lower wet bulk density, particle size, and higher porosity. The SS of silty sediment (sandy-silt and clayey-silt) was between clayey and sandy sediments. The results also show that the SS was closely related with physical properties, and the SS was proportional to the wet bulk density and particle size but was inversely proportional to porosity, which is consistent with previous studies (Hamilton, 1972; Hamilton and Bachman, 1982). Sound speed in sediments varies with seabed environment (temperature, pressure, and pore water salinity), and the laboratory environment is usually different from theenvironment. In order to eliminate the influence of the environment on the speed of sediments, Hamilton proposed the term ‘sound speed ratio (SSR)’, which was defined as the ratio of SS of sediments to the SS of water at the same temperature, salinity, and pressure (Jackson and Richardson, 2007). We used a least squares method to perform regression analysis on all data and established empirical equations for the SSR and some physical parameters (,and). The empirical equations are listed in Table 2.

    Table 1 The acoustic and physical parameters of six types of sediments

    4 Discussion

    As previous studied, SS exhibits strong inverse relation- ship with(Fig.1). Also shown in Fig.2, SSR decreased with the increase of. Whenwas approximately 75%, an inflection was encountered, and SSR began to increase with. The differences between this paper and other stu- dies (Bachman, 1989; Richardson and Briggs, 2004; Jack- son and Richardson, 2007) are plotted in Fig.2, and the equations used in this paper are listed in Table 2. The curve of this study fit well with the datameasured by using ISSAMS (Jackson and Richardson, 2007). The curve of Bachman (1989) is much higher than our curve, especially for high SS sediments. The curves of Richardson and Briggs (2004) are slightly higher than ours, while the one determined by all sediment equation was closer to the measured data. These results agree well with the work done by Wang. (2018) in the same sea area. The results obtained from the carbonate equation were only a bit less than those of all sediment equations when porosity wassmaller than 75%. However, when the porosity was greater than 75%, the prediction accuracy of carbonate equation was less than all sediment equations, which demonstrates that sediment type is an important factor influencing the SS.

    Table 2 Regression equations used in this paper

    Fig.3 illustrates the relationship between SSR andfor sediments from the SSD of SCS. It can be seen that the SSR has a positive relationship with. Whenwas larger than 1.40gcm?3, SSR began to increase with. As in the case of Fig.3, the curve of Bachman equation is higher than other studies, while the curve of this study is basically consistent with theone. After comparing Fig.2 with Fig.3, we can see that the trends ofandwere basically inverse becauseandhad a negative linear relationship (Fig.1). Just as the relationship ofand, thealso had a linear relationship with, which was more dispersed (Fig.1). The relationship between SSR andis display- ed in Fig.4. Generally, SSR increased with increasing. Kim. (2012) suggested that grain size affects sound speed through its influence on porosity, density, and other factors. The relationship between SSR andis sometimes important when grain size analyses can be made on dried sediments, on which density, porosity and sound speed measurements cannot be made (Hamilton and Bachman, 1982). In Fig.4, the curve of this study was not close to the curve of the Bachman equation or the curve of themeasurement (ISSAMS). Whenwas between 3 and 5.5f, the measured SS were relatively scattered. In this case, the Bachman equation had the highest predicted SS, and the equation established in this paper and ISSAMS equations give the results in the middle of the scattered points. When thewas between 5.5 and 9f, the measured SS was relatively concentrated. Our equation was con- sistent with the measured data; however, the predicted SS of other equations were higher than the measured data, and the siliciclastic equation has the highest predicted SS. Due to the uncertain relationship between the SSR and, the equation of the SSR based onof the sediment often has a low prediction accuracy for the data in SCS, especially whenis small. This may be due to the different particle arrangement and compaction degree in se- diments. The mean grain size (Mz) alone cannot express the characteristics of the particles (Wang., 2019). The different measuring methods of mean grain size used in these studies, which bring different errors, may be one of the reasons for not being able to use thealone.

    Fig.1 Scatter plots of acoustic and physical parameters of different type sediments.

    Fig.2 Sound speed ratio vs. porosity for sediments from the SSD of SCS. The empirical equation in this paper is: VSSR=1.587?1.591n+1.042n2, with n ranging from 0.43 to 0.82. The root mean square error (RMSE) for this equation is 0.0135, the sum of squared estimate of errors (SSE) is 0.05982 and R-square is 0.8232.

    Fig.3 Sound speed ratio vs. wet bulk density for sediments from the SSD of SCS. The empirical equation of this paper is: VSSR=1.692?1.011ρ+0.3591ρ2, in which ρ ranges 1.30 to 1.97gcm?3. The RMSE for this equation is 0.0144, the SSE is 0.06798 and R-square is 0.7991.

    Fig.4 Sound speed ratio vs. mean grain size for sediments from the SW of SCS. The empirical equation of this paper is: VSSR=1.307?0.07763Mz+0.004578Mz2, in which Mzranges 3.04 to 8.60f. The RMSE for this equation is 0.01377, the SSE is 0.06215 and R-square is 0.8163.

    Richardson and Briggs (2004) collected nearly 800 cores from 67 shallow-water sites worldwide (12 carbonate and 55 siliciclastic sediment sites), and they obtained nearly 4500 data points resulting from measurements or calculations (3922 siliciclastic and 621 carbonate) to yield the SSR regression equations (Table 2). From Fig.2, we can see that the SSR curve calculated by the equation for sili- ciclastic sediments was nearly parallel to the curve of this paper, and the maximum difference was about 0.0216. In Fig.3, the equations for siliciclastic and all sediments almost coincided, and both had a maximum difference of about 0.0212. We found in Fig.4 that the equations for siliciclastic and all sediments had a large difference in our study area with the maximum up to 0.0405. However, the discrepancy between the equations for all sediments and those in this paper is smaller, and for the low SS sediments, their curves almost coincided. The results calculated by the regression equations between SSR and physical properties are similar for carbonate and siliciclastic sediments, and Richardson and Briggs (2004) suggest a universal applicability of the SSR regression equations for all sediments. The sediments in the shallow and semi-deep areasof the SCS are mostly siliciclastic sediments and often con- tain more than 30% CaCO3(Li., 2004). Therefore, all sediment equations fit better with the data for the SSD in SCS than the siliciclastic equation, especially for porosity. All sediment equations were based on data for both siliciclastic and carbonate sediments. They have a wide range of applications and are suitable for all sediment types.

    Another probable reason for the differences among equations may be the differences in instrumentation and computation. Data from Bachman (1989) were measured at a frequency of 200kHz, whereas data from Richardson and Briggs (2004) were measured at a single high frequency (400kHz) and may not be compared with our lower frequency data (100kHz). According to dispersion characteristics, the sound speed decreases with decreasing fre- quency. It is reasonable that the values we measured could have been smaller than their measurement values. The data for theequations was measured by ISSAMS at frequency of 38 or 58kHz, and therefore have higher accuracy. Until now, no systematic analysis has been con- ducted to yield standardized sediment sound speed data. This process is out of the scope of this study, but a standard measurement procedure or anmeasuring me- thod should be developed for geo-acoustic databases.

    Hamilton and Bachman (1982) believed that marine sediments consist of only a few common mineral species, and there are few variations in the physical properties of bottom water. So, for a given sediment type (in a given environment), the bulk elastic properties should be almost the same worldwide. The common geo-acoustic relationship could be used to estimate the sound speed for unknown regions with similar porosity. We also confirmed that for given the same sediment types, the porosity-SSR equation can be used to estimate SSR for different regions. However, when the SSR of sediments was relatively low (SSR<1), there was an apparent envelope for sediments within the range from 1.30 to 1.76gcm?3,of 0.57 to 0.82, andof 5.78 to 8.60f(Phi), indicating that the SSR was relatively insensitive to changes of sediment phy- sical parameters. When the SSR of sediments was relatively high (SSR>1), there is an apparent envelope for se- diments within the range from 1.68 to 1.97gcm?3,of 0.43 to 0.58, andof 3.04 to 5.93f, which indicates that the SSR was relatively sensitive to changes of sediment physical parameters. The sound speed can be affected by multiple physical parameters although porosity is the primary one, given that the porosity may not produce an accurate speed. With the same porosity, the sediments may still have different grain sizes, resulting in different sound speeds. For example, in Fig.5, when the porosity was about 0.55 to 0.58, the different mean grain sizes (or different sediment types) account for the different sound speeds of sediments.

    Fig.5 Sound speed ratio vs. mean grain size and porosity. The red cross is for samples with Mz ranging from 3 to 5f, the yellow diamond is for Mz ranging from 5 to 7f, and the blue hollow circle is for Mz ranging from 7 to 9f.

    In order to yield a creditable equation that takes into account both porosity and sediment types, we recommend adding a parameter, mean grain size, in the above formula because the grain size can not only represent the sediment type, but also is invariant for measurements both in the seabed and in the laboratory. The two-parameter equation (SSR--equation) in our study was that:

    in whichranged from 0.43 to 0.82, andranged from 3.04 to 8.60f. The RMSE for this equation was 0.01024, the SSE was 0.03529 and-square was 0.8957. This equa- tion has a higher-square and a smaller RMSE than the single parameter equations (SSR-, SSR-equation). How- ever, it can generate an accurate estimate of speed for dif- ferent regions and sediment types when porosity and mean grain size data are available.

    5 Conclusions

    Seafloor sediments from the SSD of SCS were collected to study the relationship between geo-acoustic parameters and physical properties. The results show that the empirical equations established in this paper accord well with theequations. This study confirms the sound speed is more correlated with wet density and porosity than with mean grain size. The discrepancy among different equations may be caused by the different carbo- nate contents, the various measured frequency for sound speeds, and different measurement methods for physical properties.

    Due to the limitations of samples, the empirical equations given in this study are more suitable for the fine sediments in shallow and semi-deep sea. This study suggests that moredata should be measured, and it is necessary to strengthen the comparative study between sampling measurements andmeasurements since the sound speed of sediments is affected by many factors, such as sediment source, ocean current, and sediment types. In the next stage, we will compare and study the geo- acoustic properties of seafloor sediments in other sea areas of China, such as the South China, Yellow, and East China Seas. The applicability of empirical equations established in different areas, and the influence mechanism of sedimentary environment on geo-acoustic properties will be investigated.

    Acknowledgements

    This study was funded by the State Key Laboratory of Acoustics, Chinese Academy of Sciences (No. SKLA202 007), the National Natural Science Foundation of China (Nos. 41706045, 42076082, 41706062), the Director Fund of Qingdao National Laboratory for Marine Science and Technology (No. QNLM201713), the Guangdong Natural Science Foundation (No. 2017A030313237), the Taishan Scholar Project Funding (No. tspd20161007), the Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDA13010102), and the Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (No. GML2019ZD0104).

    Bachman, R. T., 1989. Estimating velocity ratio in marine sediment.,89: 2029-2032.

    Ballard, M. S., Lee, K. M., McNeese, A. R., Wilson, P. S., Chay- tor, J. D., Goff, J. A., and Reed, A. H., 2019.measurements of compressional wave speed during gravity coring operations in the New England mud patch., 45 (1): 26-38.

    Del Grosso, V. A., 1974. New equation for the speed of sound in natural waters (with comparisons to other equations)., 56 (4): 1084-1091.

    Endler, M., Endler, R., Wunderlich, J., Bobertz, B., Leipe, T., Moros, M., Jensen, J., and Arz, H., 2016. Geo-acoustic mo- delling of late and postglacial sedimentary units in the Baltic Sea and their acoustic visibility., 376: 86- 101.

    Goff, J. A., Kraft, B. J., Mayer, L. A., Schock, S. G., Sommerfield, C. K., Olson, H. C., and Nordfjord, S., 2004. Seabed cha- racterization on the New Jersey middle and outer shelf: Correlatability and spatial variability of seafloor sediment pro- perties., 209 (1-4): 147-172.

    Hamilton, E. L., 1970. Sound velocity and related properties of marine sediments, North Pacific., 75 (23): 4423-4446.

    Hamilton, E. L., 1972. Compressional-wave attenuation in marine sediments., 37 (4): 620.

    Hamilton, E. L., and Bachman, R. T., 1982. Sound velocity and related properties of marine sediments., 72 (6): 1891-1904.

    Hou, Z., Chen, Z., Wang, J., Zheng, X., Yan, W., Tian, Y., and Luo, Y., 2018a. Acoustic characteristics of seafloor sediments in the abyssal areas of the South China Sea., 156: 93-100.

    Hou, Z., Chen, Z., Wang, J., Zheng, X., Yan, W., Tian, Y., and Luo, Y., 2018b. Acoustic impedance properties of seafloor se- diments off the coast of southeastern Hainan, South China Sea., 154: 1-7.

    Hou, Z., Guo, C., Wang, J., Chen, W., Fu, Y., and Li, T., 2015. Seafloor sediment study from South China Sea: Acoustic & physical property relationship., 7 (9): 11570- 11585.

    Jackson, D. R., and Richardson, M. D., 2007.. Springer Science & Business Media, New York, 1-605.

    Kim, D. C., Kim, G. Y., Yi, H. I., Seo, Y. K., Lee, G. S., Jung, J. H., and Kim, J. C., 2012. Geoacoustic provinces of the South Sea shelf off Korea., 263: 139-147.

    Konno, Y., Jin, Y., Yoneda, J., Kida, M., Egawa, K., Ito, T., and Nagao, J., 2015. Effect of methane hydrate morphology on compressional wave velocity of sandy sediments: Analysis of pressure cores obtained in the eastern Nankai Trough., 66: 425-433.

    Li, X., Chen, F., Liu, J., and Huang, X., 2004. Distribution and its dissolution of carbonate in seafloor surface sediment in the western South China Sea., 33 (3): 35-41 (in Chinese with English abstract).

    Liu, J., He, Z. L., Liu, X. W., Huo, Z. Z., and Guo, P., 2019. Usingfrequency-dependent AVO inversion to predict the ‘sweet spots’ of shale gas reservoirs., 102: 283-291.

    Lu, B., Li, G., Liu, Q., Huang, S., and Zhang, F., 2008. Sea floor sediment and its acouso-physical properties in the southeast open sea area of Hainan Island in China., 26 (3): 129-144.

    Michalopoulou, Z. H., and Gerstoft, P., 2019. Multipath broadband localization, bathymetry, and sediment inversion., 45 (2): 92-102.

    Richardson, M. D., and Briggs, K. B., 2004. Empirical predictions of seafloor properties based on remotely measured se- diment impedance. In:. Porter, M. B., and Siderius, M., eds., AIP, Melville, NY, 12-21.

    Richardson, M. D., Briggs, K. B., Bentley, S. J., Walter, D. J., and Orsi, T. H., 2002. The effects of biological and hydrodynamic processes on physical and acoustic properties of sediments off the Eel River, California., 182 (1): 121-139.

    Shepard, F. P., 1954. Nomenclature based on sand-silt-clay ratios., 24: 151-158.

    Suzuki, K., Schultheiss, P., Nakatsuka, Y., Ito, T., Egawa, K., and Holland, M., 2015. Physical properties and sedimento- logical features of hydrate-bearing samples recovered from the first gas hydrate production test site on Daini-Atsumi Knoll around eastern Nankai Trough., 66: 346-357.

    Wang, J., Hou, Z., Li, G., Kan, G., and Liu, B., 2019. A new compressional wave speed inversion method based on granularity parameters., 7: 185849-185856

    Wang, J., Li, G., Liu, B., Kan, G., Sun, Z., Meng, X., and Hua, Q., 2018. Experimental study of the ballastsediment acoustic measurement system in South China Sea., 36 (5): 515-521.

    Zhou, J. X., and Zhang, X. Z., 2010. Physical parameters for four seabed geoacoustic models from low-frequency measure- ments., 1272 (1): 163-172.

    Zou, D., and Kan, G., 2011. Application study onacoustic measurement system of seafloor sediments., 9 (4): 1507-1510.

    February 7, 2020;

    May 6, 2020;

    July 24, 2020

    ? Ocean University of China, Science Press and Springer-Verlag GmbH Germany 2021

    . E-mail: zyhou@scsio.ac.cn

    (Edited by Chen Wenwen)

    免费在线观看成人毛片| 一本综合久久免费| 国产精品伦人一区二区| 日本黄色片子视频| 国产高清有码在线观看视频| 亚洲电影在线观看av| 嫁个100分男人电影在线观看| 久久人人爽人人爽人人片va | 又粗又爽又猛毛片免费看| 在线观看免费视频日本深夜| av黄色大香蕉| 色综合欧美亚洲国产小说| 日韩高清综合在线| 综合色av麻豆| 天堂√8在线中文| 精品一区二区免费观看| 精品不卡国产一区二区三区| 欧美精品啪啪一区二区三区| 久久伊人香网站| 国产亚洲精品久久久com| 精品久久久久久久久av| 欧美又色又爽又黄视频| 白带黄色成豆腐渣| 人妻制服诱惑在线中文字幕| 窝窝影院91人妻| 久99久视频精品免费| 亚洲av熟女| 久久久久久久久大av| 国产精品亚洲一级av第二区| 美女被艹到高潮喷水动态| 久久久国产成人免费| 国产一区二区亚洲精品在线观看| 亚洲国产精品999在线| 欧美不卡视频在线免费观看| 免费黄网站久久成人精品 | 男插女下体视频免费在线播放| 精品人妻视频免费看| 一本精品99久久精品77| 欧美黑人巨大hd| 男女下面进入的视频免费午夜| 美女被艹到高潮喷水动态| 午夜免费激情av| 狠狠狠狠99中文字幕| 日韩欧美国产一区二区入口| 精品久久久久久久人妻蜜臀av| 一区二区三区免费毛片| 少妇被粗大猛烈的视频| 亚洲 国产 在线| 两个人的视频大全免费| 国产探花在线观看一区二区| 99热这里只有是精品在线观看 | 极品教师在线免费播放| 此物有八面人人有两片| 成人欧美大片| 女人被狂操c到高潮| 宅男免费午夜| 亚洲专区国产一区二区| 国产精品美女特级片免费视频播放器| 亚洲无线在线观看| 一进一出好大好爽视频| 永久网站在线| 亚洲中文日韩欧美视频| 一级av片app| 欧美xxxx黑人xx丫x性爽| 国产91精品成人一区二区三区| 日韩精品青青久久久久久| 人妻丰满熟妇av一区二区三区| xxxwww97欧美| 国产美女午夜福利| 国产av不卡久久| 国产乱人视频| 一本综合久久免费| 老司机福利观看| 夜夜爽天天搞| 国产精品自产拍在线观看55亚洲| 激情在线观看视频在线高清| 熟女人妻精品中文字幕| 美女黄网站色视频| 伊人久久精品亚洲午夜| 中文字幕久久专区| 老司机福利观看| 老司机午夜十八禁免费视频| 免费看美女性在线毛片视频| 国产高潮美女av| 国产精品国产高清国产av| 一进一出抽搐动态| 偷拍熟女少妇极品色| av在线蜜桃| 亚洲国产日韩欧美精品在线观看| 亚洲精品456在线播放app | 国产黄片美女视频| 亚洲美女搞黄在线观看 | 国产精品人妻久久久久久| 国产欧美日韩一区二区精品| 99在线人妻在线中文字幕| 成人三级黄色视频| 美女免费视频网站| 中文亚洲av片在线观看爽| 国产成人a区在线观看| 麻豆av噜噜一区二区三区| 99riav亚洲国产免费| 哪里可以看免费的av片| 国产精品久久久久久久电影| 成人一区二区视频在线观看| av女优亚洲男人天堂| 久久久久久久精品吃奶| 日韩欧美 国产精品| 看片在线看免费视频| 欧美高清性xxxxhd video| 一级作爱视频免费观看| 久久午夜福利片| 国产一区二区三区在线臀色熟女| 国产精品,欧美在线| 国内少妇人妻偷人精品xxx网站| 亚洲色图av天堂| 国产精品电影一区二区三区| 成年免费大片在线观看| 亚洲国产色片| 给我免费播放毛片高清在线观看| 淫秽高清视频在线观看| 热99在线观看视频| 亚洲电影在线观看av| 看十八女毛片水多多多| 日本一本二区三区精品| 婷婷丁香在线五月| 午夜激情欧美在线| 午夜福利视频1000在线观看| 亚洲五月天丁香| 国产亚洲av嫩草精品影院| 午夜亚洲福利在线播放| 色综合站精品国产| 国产欧美日韩精品亚洲av| 日韩中字成人| 国产男靠女视频免费网站| www日本黄色视频网| 又粗又爽又猛毛片免费看| 欧美激情久久久久久爽电影| 午夜视频国产福利| 又黄又爽又刺激的免费视频.| 夜夜夜夜夜久久久久| av女优亚洲男人天堂| 日本免费a在线| 美女黄网站色视频| 精品人妻偷拍中文字幕| 日本五十路高清| 波多野结衣高清无吗| 深夜精品福利| 日本a在线网址| 亚洲av第一区精品v没综合| 高清毛片免费观看视频网站| 一级黄片播放器| www.熟女人妻精品国产| 国产精品电影一区二区三区| 欧美一级a爱片免费观看看| 99久久精品一区二区三区| 成人美女网站在线观看视频| 窝窝影院91人妻| 国产高清三级在线| 成人欧美大片| 波野结衣二区三区在线| 久久精品国产自在天天线| 国产欧美日韩精品亚洲av| 男女之事视频高清在线观看| 听说在线观看完整版免费高清| 国产精品,欧美在线| 久久精品夜夜夜夜夜久久蜜豆| 午夜福利欧美成人| 日韩欧美三级三区| 我的女老师完整版在线观看| 91麻豆精品激情在线观看国产| 国产蜜桃级精品一区二区三区| 国产乱人视频| 中文字幕av成人在线电影| 国产熟女xx| 级片在线观看| 久久久精品大字幕| 怎么达到女性高潮| 免费观看精品视频网站| 日本一本二区三区精品| 99精品久久久久人妻精品| 简卡轻食公司| 老司机午夜十八禁免费视频| 精品日产1卡2卡| 级片在线观看| 国产激情偷乱视频一区二区| 99久久99久久久精品蜜桃| 国产探花在线观看一区二区| 免费观看精品视频网站| 窝窝影院91人妻| 久久久久久久亚洲中文字幕 | 天堂影院成人在线观看| 日韩欧美精品免费久久 | 亚洲成人久久性| 免费一级毛片在线播放高清视频| 日韩欧美在线二视频| bbb黄色大片| 欧美xxxx性猛交bbbb| 午夜激情福利司机影院| 欧美黑人巨大hd| 中亚洲国语对白在线视频| 欧美成人一区二区免费高清观看| 久9热在线精品视频| 国产美女午夜福利| 国产成人aa在线观看| 国产人妻一区二区三区在| 国产成人aa在线观看| 少妇被粗大猛烈的视频| 日韩成人在线观看一区二区三区| 成人毛片a级毛片在线播放| 在线观看av片永久免费下载| 成人三级黄色视频| 免费av毛片视频| 久久人人爽人人爽人人片va | 午夜福利欧美成人| 国产高清激情床上av| av在线观看视频网站免费| 日韩成人在线观看一区二区三区| 男女下面进入的视频免费午夜| 国产一区二区三区视频了| 午夜福利在线观看免费完整高清在 | 麻豆久久精品国产亚洲av| 全区人妻精品视频| 999久久久精品免费观看国产| 亚洲色图av天堂| 免费看光身美女| 91字幕亚洲| 在线观看一区二区三区| 欧美午夜高清在线| 国产一区二区亚洲精品在线观看| 毛片一级片免费看久久久久 | 欧美不卡视频在线免费观看| 日本成人三级电影网站| 日韩免费av在线播放| 2021天堂中文幕一二区在线观| 一a级毛片在线观看| 久久中文看片网| 在线免费观看的www视频| 床上黄色一级片| 女人十人毛片免费观看3o分钟| www.www免费av| 毛片女人毛片| 亚洲欧美激情综合另类| 尤物成人国产欧美一区二区三区| 赤兔流量卡办理| 亚洲国产精品合色在线| 老司机深夜福利视频在线观看| 中文资源天堂在线| 国产爱豆传媒在线观看| 久久精品久久久久久噜噜老黄 | 69av精品久久久久久| 偷拍熟女少妇极品色| 亚洲精品影视一区二区三区av| 美女cb高潮喷水在线观看| 国产日本99.免费观看| 三级毛片av免费| 精品日产1卡2卡| 午夜视频国产福利| 久久精品91蜜桃| 97人妻精品一区二区三区麻豆| 国产黄a三级三级三级人| 久久久久国内视频| 两人在一起打扑克的视频| 校园春色视频在线观看| 真实男女啪啪啪动态图| 午夜日韩欧美国产| 亚洲熟妇熟女久久| 亚洲熟妇中文字幕五十中出| 亚洲七黄色美女视频| 国产精品美女特级片免费视频播放器| 淫妇啪啪啪对白视频| .国产精品久久| av欧美777| 嫩草影院新地址| 男插女下体视频免费在线播放| 久久久色成人| 少妇人妻精品综合一区二区 | 精品国产三级普通话版| 免费电影在线观看免费观看| 国产高清有码在线观看视频| 亚洲成人久久爱视频| 国产又黄又爽又无遮挡在线| 免费无遮挡裸体视频| 亚洲一区二区三区色噜噜| 国产精品影院久久| 五月伊人婷婷丁香| 噜噜噜噜噜久久久久久91| 午夜福利在线观看吧| 国产视频内射| 窝窝影院91人妻| 国产极品精品免费视频能看的| av国产免费在线观看| 露出奶头的视频| 1024手机看黄色片| 国内精品久久久久精免费| 变态另类成人亚洲欧美熟女| 亚洲自拍偷在线| 男女床上黄色一级片免费看| 欧美在线一区亚洲| 午夜激情欧美在线| 免费在线观看日本一区| 亚洲av二区三区四区| 91久久精品电影网| 97人妻精品一区二区三区麻豆| 91av网一区二区| or卡值多少钱| 亚洲成av人片免费观看| 高清在线国产一区| 少妇熟女aⅴ在线视频| 日本黄色视频三级网站网址| 五月伊人婷婷丁香| 成年女人永久免费观看视频| 久久久久九九精品影院| 国产视频内射| 特级一级黄色大片| 国产三级在线视频| or卡值多少钱| 日韩有码中文字幕| 亚洲av中文字字幕乱码综合| 午夜福利高清视频| av天堂中文字幕网| 亚洲av美国av| 男人的好看免费观看在线视频| 国产av一区在线观看免费| 免费在线观看影片大全网站| 欧美精品国产亚洲| 很黄的视频免费| 国内精品一区二区在线观看| 亚洲美女黄片视频| 夜夜爽天天搞| 夜夜躁狠狠躁天天躁| 久久精品人妻少妇| 黄色视频,在线免费观看| 亚洲人成伊人成综合网2020| 美女被艹到高潮喷水动态| 在线天堂最新版资源| 婷婷六月久久综合丁香| 啪啪无遮挡十八禁网站| 国产精品亚洲av一区麻豆| 观看免费一级毛片| 少妇的逼水好多| 国产爱豆传媒在线观看| 欧美xxxx黑人xx丫x性爽| 可以在线观看毛片的网站| 日韩免费av在线播放| 国产精品女同一区二区软件 | 丁香欧美五月| 综合色av麻豆| 一进一出抽搐动态| 99热精品在线国产| 亚洲成人久久性| 久久精品国产清高在天天线| 高清毛片免费观看视频网站| 又爽又黄a免费视频| 波野结衣二区三区在线| 成人三级黄色视频| 露出奶头的视频| aaaaa片日本免费| 国产激情偷乱视频一区二区| 麻豆久久精品国产亚洲av| 90打野战视频偷拍视频| 亚洲五月婷婷丁香| 成人精品一区二区免费| 少妇人妻一区二区三区视频| 51午夜福利影视在线观看| 久久久色成人| 免费观看精品视频网站| 99国产极品粉嫩在线观看| 两个人视频免费观看高清| 亚洲精品影视一区二区三区av| 欧美日韩亚洲国产一区二区在线观看| 网址你懂的国产日韩在线| 男人舔女人下体高潮全视频| 亚洲狠狠婷婷综合久久图片| 在线看三级毛片| 国产熟女xx| 国产乱人视频| 亚洲av日韩精品久久久久久密| 日韩高清综合在线| 国产在线精品亚洲第一网站| 变态另类成人亚洲欧美熟女| 高清日韩中文字幕在线| 青草久久国产| 脱女人内裤的视频| 97人妻精品一区二区三区麻豆| 中亚洲国语对白在线视频| 成人美女网站在线观看视频| 国产极品精品免费视频能看的| 欧美黑人欧美精品刺激| 亚洲最大成人手机在线| 啦啦啦观看免费观看视频高清| 欧美三级亚洲精品| 色吧在线观看| 在现免费观看毛片| 亚洲av中文字字幕乱码综合| 国内精品美女久久久久久| 波多野结衣高清作品| 一级黄片播放器| 91在线观看av| 中文字幕人成人乱码亚洲影| 一区二区三区免费毛片| 精品国内亚洲2022精品成人| 精品人妻偷拍中文字幕| 超碰av人人做人人爽久久| 欧美另类亚洲清纯唯美| 国产伦一二天堂av在线观看| 久久精品久久久久久噜噜老黄 | 欧美日韩亚洲国产一区二区在线观看| 少妇熟女aⅴ在线视频| 91麻豆精品激情在线观看国产| 午夜激情欧美在线| 人人妻人人看人人澡| 草草在线视频免费看| 欧美成人一区二区免费高清观看| 中出人妻视频一区二区| 免费av不卡在线播放| 久久久久久久久大av| 男插女下体视频免费在线播放| 日本撒尿小便嘘嘘汇集6| 91九色精品人成在线观看| 亚洲精品在线美女| 亚洲综合色惰| 成人性生交大片免费视频hd| 亚洲av中文字字幕乱码综合| 在线观看66精品国产| 香蕉av资源在线| 国产成人aa在线观看| 搞女人的毛片| 美女高潮喷水抽搐中文字幕| 国产亚洲欧美98| 国产精品98久久久久久宅男小说| 美女高潮喷水抽搐中文字幕| 在线免费观看不下载黄p国产 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 两个人的视频大全免费| 日本免费a在线| 久久精品国产亚洲av天美| 亚洲中文字幕日韩| 十八禁网站免费在线| 一级av片app| 不卡一级毛片| 亚洲性夜色夜夜综合| 中文在线观看免费www的网站| 波多野结衣高清无吗| 亚洲精品影视一区二区三区av| 啦啦啦观看免费观看视频高清| 国产午夜精品久久久久久一区二区三区 | 999久久久精品免费观看国产| 久久精品综合一区二区三区| 色av中文字幕| 伊人久久精品亚洲午夜| 他把我摸到了高潮在线观看| 亚洲人与动物交配视频| 亚洲精品一卡2卡三卡4卡5卡| 久久人人爽人人爽人人片va | 日本一本二区三区精品| 国产三级中文精品| 色播亚洲综合网| 日日摸夜夜添夜夜添小说| 国产一区二区在线观看日韩| 麻豆成人午夜福利视频| 成人国产综合亚洲| 色吧在线观看| 三级毛片av免费| 在线观看66精品国产| 国内揄拍国产精品人妻在线| 久久伊人香网站| 在线播放国产精品三级| 亚洲 国产 在线| 午夜免费激情av| 欧美激情久久久久久爽电影| АⅤ资源中文在线天堂| 内地一区二区视频在线| 亚洲乱码一区二区免费版| 一本一本综合久久| 国产伦在线观看视频一区| 亚洲av免费在线观看| 精品日产1卡2卡| 国产高潮美女av| 精品久久久久久久末码| 亚洲国产日韩欧美精品在线观看| 日韩欧美一区二区三区在线观看| 国产69精品久久久久777片| 国产一级毛片七仙女欲春2| 国产高清有码在线观看视频| 美女被艹到高潮喷水动态| 免费av毛片视频| 夜夜夜夜夜久久久久| 日本黄色视频三级网站网址| 我要搜黄色片| 天堂av国产一区二区熟女人妻| 精品乱码久久久久久99久播| 中出人妻视频一区二区| 日韩免费av在线播放| 男女视频在线观看网站免费| 看片在线看免费视频| 亚洲综合色惰| 国产69精品久久久久777片| 午夜福利高清视频| 天天躁日日操中文字幕| 在线播放无遮挡| 亚洲精品色激情综合| 色尼玛亚洲综合影院| 此物有八面人人有两片| 午夜福利在线观看吧| 亚洲成av人片在线播放无| 高清日韩中文字幕在线| 麻豆av噜噜一区二区三区| 日本黄大片高清| 色精品久久人妻99蜜桃| 在线十欧美十亚洲十日本专区| 亚洲国产欧美人成| 久久伊人香网站| 国产亚洲精品久久久com| 成人亚洲精品av一区二区| 日本免费一区二区三区高清不卡| 色播亚洲综合网| 天堂av国产一区二区熟女人妻| 国产私拍福利视频在线观看| 老司机午夜福利在线观看视频| ponron亚洲| 99riav亚洲国产免费| 亚洲一区二区三区不卡视频| 床上黄色一级片| 国产综合懂色| 久久国产乱子伦精品免费另类| 国产高清有码在线观看视频| 乱码一卡2卡4卡精品| 欧美在线黄色| 国产日本99.免费观看| 国产免费一级a男人的天堂| 久久久久国内视频| 老司机深夜福利视频在线观看| 国产精品爽爽va在线观看网站| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品国产亚洲在线| 成人无遮挡网站| 一个人看的www免费观看视频| 欧美性猛交╳xxx乱大交人| 亚洲欧美日韩高清在线视频| 久久精品国产亚洲av香蕉五月| 国产精品亚洲av一区麻豆| 久久人人精品亚洲av| 99精品久久久久人妻精品| 成人特级黄色片久久久久久久| 一区二区三区免费毛片| 国产白丝娇喘喷水9色精品| 搡老妇女老女人老熟妇| 亚洲av成人精品一区久久| 久久99热这里只有精品18| 色av中文字幕| 久久国产乱子免费精品| 成人鲁丝片一二三区免费| 欧美日本亚洲视频在线播放| 久99久视频精品免费| 日本免费一区二区三区高清不卡| 最近在线观看免费完整版| 国产精品1区2区在线观看.| 嫩草影院新地址| 欧美激情国产日韩精品一区| 嫩草影院精品99| 国产精品野战在线观看| АⅤ资源中文在线天堂| 怎么达到女性高潮| 九色成人免费人妻av| 看黄色毛片网站| 成年女人永久免费观看视频| 国产欧美日韩精品一区二区| 长腿黑丝高跟| 免费观看精品视频网站| 88av欧美| 午夜久久久久精精品| avwww免费| 亚洲第一欧美日韩一区二区三区| 久久欧美精品欧美久久欧美| 一进一出抽搐动态| 中文字幕av成人在线电影| bbb黄色大片| 91字幕亚洲| 中文字幕精品亚洲无线码一区| 国产三级中文精品| 18禁黄网站禁片午夜丰满| 成人国产一区最新在线观看| 中文字幕精品亚洲无线码一区| 国产极品精品免费视频能看的| 一本一本综合久久| 最好的美女福利视频网| 亚洲av不卡在线观看| 久久99热这里只有精品18| 丰满人妻熟妇乱又伦精品不卡| 亚洲18禁久久av| 9191精品国产免费久久| 日韩av在线大香蕉| 国产精品亚洲一级av第二区| 午夜福利视频1000在线观看| 在线天堂最新版资源| 亚洲欧美日韩高清专用| 麻豆久久精品国产亚洲av| 麻豆成人午夜福利视频| 动漫黄色视频在线观看| 十八禁人妻一区二区| 91狼人影院| 亚洲七黄色美女视频| netflix在线观看网站| 日韩欧美免费精品| 黄色配什么色好看| 国产高清有码在线观看视频| 黄色配什么色好看| 亚洲精华国产精华精| 精品一区二区三区视频在线观看免费| 精品一区二区三区视频在线| 嫩草影视91久久| 成年女人毛片免费观看观看9| 少妇人妻一区二区三区视频| 亚洲av电影不卡..在线观看| 国产主播在线观看一区二区|