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

    Forecasting method of induced passenger flow for Sichuan-Tibet railway based on improved MD model

    2020-05-10 09:11:54ZhangJinXuJunxiangGuoJingni

    Zhang Jin Xu Junxiang Guo Jingni

    (School of Transportation and Logistics, Southwest Jiaotong University, Chengdu 611756, China)

    Abstract:Aiming at the limitation of forecasting induced passenger flow with the modal demand(MD) model, which assumes that the travel time value and travel utility obey the fixed probability distribution, a forecasting method of induced passenger flow is proposed by using an improved MD model based on the prospect theory. Combined with the actual traffic flow survey data of the Sichuan Tibet passenger transport channel, the improved MD model is used to forecast the total passenger flow and induced passenger flow along the Sichuan Tibet railway. Finally, the improved forecasting model is compared with the traditional MD model and the Logit model, and the effects of the parameters in the foreground theory on the forecasting results are analyzed. The results show that the Logit model forecasting result is 8.5% higher than the improved MD model forecasting result, the traditional MD model forecasting result is 3.7% higher than the improved MD model forecasting result, but the difference between the forecasting results of each model is within 10%.The parameter value in the prospect theory has a significant effect on the passenger flow forecasting results, especially in the two types of travelers, the variable and the conservative, where the forecasted difference reaches 23%.

    Key words:induced passenger flow; passenger flow forecasting; prospect theory; modal demand (MD) model

    In recent years, an increasing number of experts and scholars have paid attention to the research of induced passenger flow forecasting. Forecasting of induced passenger flow for new-built passenger dedicated lines can be constructed based on the RBF neural network model for short-term traffic flow forecasting according to the characteristics of high forecasting accuracy of the RBF neural network in complex network systems[1]. By analyzing the mechanism of induced passenger flow, a passenger flow forecasting model based on the TOD model is proposed, and a growth curve forecasting model based on the TOD model is established. This model mainly forecasts induced passenger flow through land use change[2], but this method is not suitable for forecasting the induced passenger flow in the initial stage of new-built high-speed railways due to the explosive manifestation of potential demand. According to the autoregressive moving average model and utility theory, the diversion passenger flow and induced passenger flow can be distinguished from each other[3]. Some scholars have quantitatively analyzed the spatial and temporal distribution characteristics of induced passenger flow, and proposed a method of using the fractal dimension theory to describe the complex characteristics and variation characteristics of induced passenger flow in urban rail transit. This method is based on the variation and distribution of passenger flow data to depict the characteristics of induced passenger flow, without considering the influencing factors of passenger travel behavior, so that it can be applied to the intersection. The research of traffic flow distribution is not suitable for the forecasting of induced passenger flow[4-5], and the limitation of this method is that the forecasting model can be applied to Beijing, Shanghai and other economically developed areas as well as areas with an obvious polarization effect in the growth pole theory. If it is a poor area, this model is no longer applicable. There are also many methods for forecasting induced passenger flows, for example using the trip distribution-mode choice (TDMC) model to forecast the increased passenger flow by public transport[6].Aiming at the train delay condition, the neural network model is used to forecast the passenger flow[7-8]. The traditional forecast method of induced passenger flow adopts the four-stage method[9-13].

    The MD forecasting model is a new passenger flow forecasting method developed by the Japan Transport Administration in the 1970s. It was first applied to the passenger flow forecasting of the Shinkansen Line. Its principle is to adopt the utility theory of economics and take the individual traveling as the object to analyze the potential passenger demand, so it is very suitable for the induced increase of new-built railway lines[14-17]. The prospect theory can describe people’s risk attitudes in the face of gains or losses and explain the impact of travel utility considering travelers’ psychological expectations on passenger flow forecasting results. Therefore, this paper considers using the improved MD model based on the prospect theory to forecast the induced passenger flow after the completion of the Sichuan-Tibet railway. The main innovations are as follows:

    1) The travel time value and travel utility in the traditional MD model are recharacterized by the prospect theory value function and weight function, and the travel psychological behavior of travelers is taken into consideration in the study.

    2) The improved MD model is used to forecast the induced passenger flow after the completion of the Sichuan Tibet railway, and it is compared with the traditional forecasting model to verify the rationality of this method.

    3) This research shows the influence of different parameters in the value function and weight function of the prospect theory on the forecasting results, and analyzes the influence of travelers’ different travel psychological preferences for the forecasting results.

    1 Theory and Method

    The prospect theory was proposed by Kahneman and Tversky in 1979[18]. It is a descriptive theory, which assumes that the decision maker is limited in rationality, that the individual’s ability to calculate and think is limited, and that the judgment information of the options is not fully determined. The prospect theory holds that the decision maker will receive the first satisfactory result compared with the expected reference point, but not necessarily the result of the greatest benefit[19].In the prospect theory, value functionv(x) and decision weight functionw(p) are usually used to describe individual decision-making behavior.The value function is a subjective function, which represents the difference between the current value and the reference point, in which the reference point is determined by the subjective judgment of the decision maker. The decision weight function is a subjective probability function, which is the subjective judgment made by the decision maker according to the probability in the expected utility theory.The value function and the decision weight function together constitute the prospect function, which affects the final decision of the decision maker. The specific form of the prospect function is

    Vi=w+(p)v+(x)+w-(p)v-(x)

    (1)

    whereVirepresents the comprehensive prospect value;v+(x) represents the profit value function;w+(p) represents the profit weight function;v-(x) represents the loss value function;w-(p) represents the loss weight function.

    Profit value function(v+(x)) and loss value function (v-(x)) are expressed as

    (2)

    where Δxrepresents the difference between the actual decision result and the reference point;σrepresents the profit pursuit coefficient;λrepresents the loss avoidance factor. In the actual decision-making process, there is a natural endowment effect. When people own an item, their evaluation value of the item is much higher than before. Therefore, when people lose the item, they expect to sell it at a high price. Therefore, the acquiescence of profit pursuitσis 1, so the value function can be expressed as

    (3)

    Whenλ>1, it means that decision makers are more sensitive to losses than to gains.αrepresents a diminishing marginal sensitivity parameter for revenue, andβrepresents a marginally sensitive decreasing parameter for loss, 0<α<1, 0<β<1.

    The prospect theory holds that people tend to overestimate probability events and underestimate probability events in decision-making, but they are relatively insensitive to the change of intermediate probability. The expression of decision weight function is

    (4)

    whereωis the risk-return attitude coefficient;μis the attitudinal coefficient of risk loss.

    At present,the prospect theory has been widely applied in many fields, such as route selection[20], mode selection[21], dynamic traffic flow allocation[22], and emergency plan decision-making[23], and many valuable research results have been obtained.

    Several basic concepts in the MD model include the potential passenger demand, travel utility, travel sacrifice, time value and potential demand manifestation rate, which can describe the mechanism of induced passenger flow as described above[24].Therefore, it is suitable to use the MD model to research the induced passenger flow forecast caused by the explosive growth of demand in the near future after the completion of high-speed rail lines[25].The variation of the induced passenger flow is shown in Fig.1.

    Fig.1 The variation of induced passenger flow

    The interpretations of the potential passenger demand, travel utility, travel sacrifice, time value and potential demand manifestation rate in the MD model are as follows:

    1) Potential passenger demand.In the MD model, the potential passenger demand is defined as the total number of potential passengers who may displace from the origin to the destination. In the MD model, potential passenger demand refers to the number of people who do not consider the ability to pay between OD as long as they have the intention to travel.

    2) Travel utility.Travel utility in the MD model refers to the value or benefit that travelers gain through this trip, such as travel, trade, school, etc.The MD model holds that people are willing to travel because they want to obtain travel utility. The level of social and economic development of the region determines the general level of travel utility. The travel utility is also related to the characteristics of travelers. The traditional MD model assumes that the evaluations of travel utility obeys the lognormal distribution, and all OD pairs present this kind of distribution.However, it should be noted that the travel utility in the MD model is different from that in the traditional four-stage forecasting model. The value of travel utility in the MD model is only related to the traveler’s own factors. It is a subjective psychological evaluations of whether the traveler can satisfy his/her inner desire to travel. Price has no direct impact on the specific mode of transportation and its service characteristics, while the travel utility in the traditional four-stage forecasting model is related to both[26].

    3) Travel sacrifice. In the MD model, the traveler’s cost to obtain the travel value is called the travel sacrifice. This cost is a generalized cost including travel time consumption, travel cost, travel fatigue, etc.

    4) Time value.Time value refers to the monetary manifestation of the added value of benefits and the loss of benefits caused by the consumption of time. Different travelers have different evaluations of time value, which is due to the different evaluations of time value by different travelers for different purposes and time-consumption[27]. In the MD model, it is also assumed that the evaluation of time value also obeys lognormal distribution.

    5) Potential demand manifestation rate.When the traveler’s travel utility is greater than the traveling sacrifice of the mode, that is to say, the benefit is more than the cost, the potential transport demand will be realized, which is called the manifestation of the potential transport demand. This part of passenger flow formed after the manifestation of the potential demand is called the passenger demand volume, and its proportion in the total potential passenger demand volume is called the latent demand manifestation rate.

    2 Forecasting Process of Induced Passenger Flow

    2.1 Determination of time value

    In the specific transportation process, if travel time is not consumed in this way, then it can be used for other labor productions and create economic benefits and value. Therefore, the authors assume that the economic value created in unit time is equivalent to the time value of travelers. The expression of the decrease in income caused by the time consumed in travel which leads to the incapacity to work is

    (5)

    (6)

    2.2 Determination of travel sacrifice quantity

    1) Travel time consumingT. Travel time is the total time spent by the chosen means of transport in the whole OD, including the travel time during the journey and the connection time between the start and end of the journey and the station. As the research object of this paper is the newly built railway, without considering other modes of transport, it can be expressed as

    (7)

    f(T)=viT

    (8)

    2) Travel expensesC. Travel expenses are the total expenses that travelers spend in the process of traveling. They can be expressed as

    (9)

    whereCirepresents the total cost in the travel process. Assuming thatf(Ci) represents the contribution value of passenger travel expenses to travel sacrifices, the monetary unit dimension can be expressed as

    (10)

    3) Travel comfortG(t). Different modes of transportation will bring different feelings of comfort to passengers. There are many factors affecting the comfort degree of passengers, such as stability, seat comfort, the internal environment, service quality, etc. The authors consider using the travel fatigue recovery time of travelers to quantify the comfort degree, which can be expressed as

    (11)

    whereRmaxrepresents the longest time it takes to recover from travel fatigue;γis a dimensionless parameter, representing the minimum travel fatigue recovery time;τrepresents the fatigue recovery time intensity factor per unit travel time. According to Ref.[28], whenRmax=15, the value of parameterτis 0.28. Assume thatf(G) represents the contribution value of the comfort degree to travel sacrifice, which is the same as that of travel time quantification. The comfort degree is converted into a monetary unit dimension by unit time value, which can be expressed as

    f(G)=viG(t)

    (12)

    4) Travel safety.Travelers need to bear the risk of travel when choosing the mode of transportation. The authors quantify the safety of travel by using the rate of travel accidents. The rate of travel accidents is defined as the number of passenger casualties that occur within 1.0×108km of passenger turnover in a certain period of time. The following expression can be used to express the number of passenger casualties that occur within 1.0×108km of passenger turnover.

    (13)

    whereεrepresents the accident rate;Arepresents the passenger traffic volume;Lrepresents mileage. Assume thatf(ε) represents the contribution value of safety loss to travel sacrifice. The amount of compensation after the accident is converted into a monetary unit dimension, which can be expressed as

    f(ε)=θε

    (14)

    According to the national policy in China, the limit of liability for personal injury and death of railway passengers caused by accidents is 150 000 yuan, so the value is 150 000 yuan.

    In the MD model, the amount of travel sacrifice is affected by travel time, travel cost, travel comfort and travel safety. The formula for calculating the amount of travel sacrifice can be obtained by adding the contribution values:

    (15)

    2.3 Determination of travel utility

    The traditional MD model considers that the travel utility of travelers is the same as the time value and obeys lognormal distribution. It has been analyzed above that the travel utility of the MD model is only related to the traveler himself, and different travelers have different travel utilities. Therefore, the traveler’s travel benefit can be described by the value function and decision weight function in the prospect theory. The authors believe that the traveler’s expectation of travel benefit can be reflected by the comprehensive prospect value function of travel, which can be expressed as

    (16)

    (17)

    (18)

    f(Vi)=(f(lnu))αw+(f(lnu))+

    (f(lnu))βw-(f(lnu))

    (19)

    Travel utility can be described by Eqs.(16) to (19). Eq. (16) represents the expression of travel utility in the traditional MD model. Assume thatηlnuis the logarithmic mean of travel utility andφl(shuí)nuis the logarithmic variance of travel utility.μandσare the traditional mean variance symbols, which are used in Eqs.(2) and (4), so the symbols are replaced here.Eq.(17) denotes the value function in the prospect theory, Eq.(18) denotes the decision weight function in the prospect theory, and Eq.(19) denotes the comprehensive prospect function, which is used here to denote travel utility.

    2.4 Determination of potential demand and potential demand explicit rate

    According to the forecasting principle of the MD model, the combination formula can be expressed as

    (20)

    whereRn,ijrepresents the potential passenger manifestation rate for selecting modenfrom nodeitoj;qn,ijdenotes the passenger volume from nodeitojto choose the mode of transportn;Qijrepresents the potential passenger traffic volume of modenchosen fromitoj;vn+1,nandvn,n-1represent the border substitution rates between modes of transportation. Potential passenger demand can be expressed as follows and potential passenger demand changes are shown in Fig.2.

    (21)

    Fig.2 Potential passenger demand change chart

    The forecast steps of the induced passenger flow based on the improved MD model are shown in Fig.3.

    Fig.3 Flow chart of induced passenger flow forecasting

    3 Forecasting of Induced Passenger Flow

    3.1 Introduction for the Sichuan-Tibet Railway

    Located in Sichuan Province and the Tibet Autonomous Region of China, the Sichuan-Tibet railway, with a total length of 1 838 km, is a fast railway connecting Sichuan Province and Tibet Autonomous Region. It is the second railway to enter the Tibet in China and one of the trunk railways in southwest China[29]. The traffic zones along the Yalin section of Sichuan-Tibet Railway can be divided into the Ganzi Tibetan Autonomous Prefecture, Changdu area, Linzhi area and Shannan area.

    3.2 Current situation of traffic passenger volume

    3.2.1 Highway passenger flow

    The Highway along the Sichuan-Tibet Railway is dominated by the G317 and G318 national highways. Combining the traffic investigations in G317 and G318 and collating the OD survey data of motor vehicles in relevant station statistical reports, the passenger flow tables of highway passage sections along the Sichuan-Tibet Railway in 2018 (base year) are shown in Tab.1.

    Tab.1 OD table for highway passengers in alpine areas along the Sichuan-Tibet railway in the base year 104person/a

    ODGanzi TibetanChangdu areaLinzhi areaShannan areaLhasaGanzi Tibetan0845642102Changdu area780413882Linzhi area665002680Shannan area211116019Lhasa785221140

    3.2.2 Air passenger flow

    Since there are almost no direct flights on the air airport routes along the Sichuan-Tibet Railway in alpine areas, the Kangding-Lhasa flights on the only route are very few, and the annual passenger flow is very low. It is not useful to count the air traffic in the alpine areas along the Sichuan-Tibet Railway. Therefore, the research considers that for the main statistics of air transport corridors, the passenger flow of the whole line of the Chuan-Tibet Railway has practical significance for the future research of the passenger flow transfer caused by the new construction of the Sichuan-Tibet Railway.The specific air passenger flow statistics are shown in Tab.2(The statistical passenger flow is direct passenger flow). The actual total air passenger flow of the Sichuan-Tibet Railway in the frigid areas studied in this paper only includes Kangding, Changdu, Linzhi and Lhasa.

    Tab.2 Air OD table of Sichuan-Tibet area in the base year 104person/a

    ODChengduKangdingChangduLinzhiLhasaXigazeChengdu0220124246396282Kangding114000560Changdu68000240Linzhi14200000Lhasa212840000Xigaze12400000

    3.2.3 Passenger flow in corridor

    Combining the passenger flow density of the above two modes of transportation, the passenger flow in the attraction range of the transportation modes in the Sichuan-Tibet area in 2018 is shown in Tab.3.

    Tab.3 Passenger flow and market share of two transportation modes in Sichuan-Tibet region

    TransportationPassenger flow/(104 person·a-1)Market share/%Highway transport97758.3Trans70041.7Total1 677100

    3.3 Improvement of MD model forecasting based on the prospect theory

    According to the data of the current investigation, the passenger volume operation of the Sichuan-Tibet Railway in the alpine regions after its completion is studied and forecast. The forecast base year is 2018, and the total OD of the current passenger flow is 16.77 million. Since the Sichuan-Tibet Railway is expected to open to traffic in 2030, the forecast year is set to be 2035 in the near future and 2040 in the mid-term.

    The calculation steps of the improved MD model based on the prospect theory are as follows:

    1) Calculation of time value

    After completing the preparatory work of the data collection, questionnaire survey and traffic district division, the expected reference points of time value for people in this area are determined by Eq.(5) according to the regional economy and population, and the values ofλandβin Eq.(6) are obtained by a questionnaire analysis, so as to obtain the the information of travelers in this area. This paper includes a questionnaire designed to help understand the travel characteristics of people in Sichuan and Tibet. Through data analysis, we obtain the values of parametersλandβof the prospect theory in Eq.(5), which are 1.38 and 0.59, respectively.

    2) Calculation of traffic sacrifice

    The travel time and fare of transportation corridors in Sichuan and Tibet were collected, and the travel comfort of people was investigated by designing questionnaires for different travel modes. The contribution values of different factors affecting travel sacrifice are calculated by combining Eqs.(7) to (14). Eq.(15) is used to sum the contribution values of each influencing factor, and the base and target annual values of the line sacrifice are calculated.

    3) Calculation of traffic utility

    Through the analysis of the questionnaire and the calculation results of passenger flow and time value of different modes of transportation between OD pairs, the values of parametersα,β,ω,μandλbased on the foreground theory model are obtained. The travel utility can be obtained by substituting Eq.(16) into Eq.(19). In this paper, the values of theoretical parameters of foreground in Eq.(16) and Eq.(19) are obtained by data analysis through questionnaires,α=0.34,β=0.46,ω=1.12,μ=0.82,λ=1.58.

    4) Calculation of the potential demand manifestation rate

    According to the calculation of the apparent rate of potential demand, the relevant parameters obtained above are substituted into Eq.(20), and the upper and lower bounds of the integral are substituted for the substitution rate of the boundaries of different modes of transportation. Then the basic annual valueR0n,ijand the target annual valueRtn,ijof the apparent rate of different OD pairs of transportation modes can be obtained.

    5) Estimation of potential passenger demand in the base year

    By dividing the passenger demand between base year intervaliandjby the apparent rate of potential demand in the base year, the potential passenger demand in the base year, i.e. Eq. (21) is obtained.

    6) Calculation of potential passenger transport demand in the future year

    Based on the socio-economic and demographic data of the base year of a traffic district, the potential passenger demand of the base year calculated in the previous step is used to determine the potential passenger demand function of the target year. Then, the forecasting value of the socio-economic population index of the target year is substituted into the potential passenger transport demand function, and the potential passenger transport demand of the target year is obtained.

    7) Calculation of passenger transport demand in the coming year

    Passenger demand in the coming year is calculated according to the potential demand dominance rateRtn,ijand obtained potential passenger demandQt,ij.

    qtn,ij=Qt,ijRtn,ij

    (22)

    8) Transfer of passenger traffic

    The transfer of passenger traffic is due to the change of traffic mode between OD and OD after the opening of the railway line. This part of passenger flow is divided into a new railway, that is, the transfer of passenger traffic. The calculation formula is

    M=Qt-Q0

    (23)

    whereMdenotes the transfer of passengers after the construction of a new railway;Qtindicates the passenger flow of existing modes of transport when there are no new modes of transport in the coming year;Q0indicates the passenger flow of existing modes of transport when new modes of transport are added in the coming year.

    9) Induced passenger flow

    The induced passenger flow of the newly built railway can be expressed by the following formula according to the description above:

    Qt,ij=qt,ij-Mt,ij

    (24)

    whereQt,ijdenotes the induced passenger flow of the newly built railways in intervalsitojin the coming year;qt,ijrepresents passenger flow betweeniandjin the next year;Mt,ijrepresents the sum of passenger flow transfers from existing modes of transport in the coming year.

    From the forecasting steps above, we can obtain the forecast of passenger flow in the alpine area along the Yalin section of Sichuan-Tibet Railway in the future. The results are shown in Tab.4.

    Tab.4 Forecast results of passenger flow along Sichuan-Tibet Railway in the alpine region 104person/a

    Forecast year20352040Sichuan-Tibet Highway9471 120Transport Aviation8621 094Sichuan-Tibet Railway1 5802 350

    From this, we can obtain the forecast results of the induced passenger flow in the alpine areas along the Sichuan-Tibet Railway as shown in Tab.5.

    Tab.5 Forecast results of induced passenger flow in the alpine areas along the Sichuan-Tibet Railway 104person/a

    Forecast year20352040Induced passenger flow 1 7122 887

    3.4 Comparative analysis of forecasting model

    In order to verify the accuracy and reliability of the MD model for the forecasting results of induced passenger flow of the newly-built Sichuan-Tibet Railway, this paper uses the improved MD model based on the prospect theory, the Logit model in the four-step forecasting method, and the traditional MD model to compare the forecasting of total induced passenger flow in 2035 in the areas along the line. The results are shown in Fig.4.

    It can be seen from Fig.4 that the results of the improved MD model are lower than those of the traditional MD model and Logit model. For induced passenger flow forecasting, the Logit model forecasting results are 8.5% higher than the results in this paper, and the traditional MD model forecasting results are 3.7% higher than the results in this paper, but the difference between the forecasting results of each method is within 10%. This shows that the forecasting method based on the improved MD model proposed in this paper shows reliability in the research of induced passenger flow forecasting.

    3.5 Analysis of the influence of parameters α and β on the results

    The values of parametersαandβin the value function of the prospect theory reflect the travelers’ psychological behavior of travel decision-making. This paper discusses the effect of parameters on the forecasting results by comparing the fluctuation of passenger flow in different modes of transportation and the changes of induced total passenger flows under different values. The specific parameter values are shown in Tab.6, and the impact of the forecasting results is shown in Fig.5.

    (a)

    (b)

    From Fig.5, it can be seen that the effect of different parameters and values on the forecasting results is quite different, especially for the two types of traveler pairs, which are variable and conservative. Therefore, it can be seen that the improved MD model based on the prospect theory is important for the analysis of the effect of parameter values on the forecasting results in the prospect theory. Using this method, we need to fully investigate the

    Tab.6 Parameter value characteristics of travelers between traffic OD and accurately describe the parameter value of the prospect theory.

    Traveler categoryParameter value rangeParameter valueVariable type0<α,β<1α=0.34,β=0.46Conservative typeα>1,β>1α=2.25,β=1.62Intermediate typeα=β=1α=β=1

    (a)

    (b)

    4 Conclusions

    1) Based on the prospect theory, it is feasible and applicable to improve the forecasting method of the MD model to forecast the traffic flow of the railway under construction.

    2) The forecasting model based on the prospect theory describes travelers’ psychological behavior, so the forecasting result will be more accurate than those of the traditional methods.

    3) The parameter value in the prospect theory reflects travelers’ travel decision preference, which has a great impact on the passenger flow forecasting results. Therefore, it is necessary to accurately depict the characteristics of regional travelers in the future research.

    久久人人爽人人爽人人片va| 亚洲成人久久爱视频| 亚洲国产成人一精品久久久| 欧美高清成人免费视频www| 日本免费在线观看一区| 97超视频在线观看视频| 一二三四中文在线观看免费高清| 久久精品久久久久久久性| 精品国内亚洲2022精品成人| 成人毛片60女人毛片免费| 国产精品av视频在线免费观看| 国产精品爽爽va在线观看网站| 久久久色成人| 国产日韩欧美在线精品| 免费观看的影片在线观看| 蜜桃久久精品国产亚洲av| 亚洲在久久综合| 神马国产精品三级电影在线观看| 久久精品综合一区二区三区| 久久人人爽人人片av| 极品少妇高潮喷水抽搐| 国产高清不卡午夜福利| 亚洲欧美成人综合另类久久久| 视频中文字幕在线观看| 欧美成人精品欧美一级黄| 可以在线观看毛片的网站| 日本-黄色视频高清免费观看| 精品一区二区免费观看| 国产在视频线精品| 日韩视频在线欧美| 久久久亚洲精品成人影院| 免费av毛片视频| 色播亚洲综合网| 99久国产av精品| 亚洲无线观看免费| 国产一区二区三区av在线| 我要看日韩黄色一级片| 日韩一区二区三区影片| 熟妇人妻久久中文字幕3abv| 欧美bdsm另类| 国产精品久久视频播放| 别揉我奶头 嗯啊视频| 国产成人精品福利久久| 亚洲伊人久久精品综合| 久久国内精品自在自线图片| 成人一区二区视频在线观看| 国产精品综合久久久久久久免费| 搡老乐熟女国产| 熟女人妻精品中文字幕| 日韩成人伦理影院| 亚洲精品成人久久久久久| 日本免费a在线| 麻豆精品久久久久久蜜桃| 80岁老熟妇乱子伦牲交| 22中文网久久字幕| 搡老乐熟女国产| 亚洲精品成人av观看孕妇| 一级毛片aaaaaa免费看小| 成年女人看的毛片在线观看| 亚洲精品乱久久久久久| 精品一区二区三区视频在线| 大又大粗又爽又黄少妇毛片口| 国产精品久久久久久久电影| 国产成人精品婷婷| 欧美xxxx性猛交bbbb| 联通29元200g的流量卡| 亚洲激情五月婷婷啪啪| av在线天堂中文字幕| 少妇人妻精品综合一区二区| 美女被艹到高潮喷水动态| 日本黄大片高清| 国产精品国产三级国产av玫瑰| 国产精品久久视频播放| 国产在线一区二区三区精| 精品酒店卫生间| 3wmmmm亚洲av在线观看| 日日摸夜夜添夜夜爱| 国产精品久久久久久久久免| 久久久久久久久久人人人人人人| 国产成年人精品一区二区| 久久久久九九精品影院| 国产日韩欧美在线精品| 国产高潮美女av| 在线免费观看的www视频| 91精品伊人久久大香线蕉| 久久草成人影院| 伦精品一区二区三区| 日日啪夜夜撸| 一个人免费在线观看电影| 日本黄色片子视频| 熟妇人妻久久中文字幕3abv| 日韩中字成人| 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| 人妻制服诱惑在线中文字幕| 美女国产视频在线观看| 丝袜美腿在线中文| 久99久视频精品免费| 女的被弄到高潮叫床怎么办| 国产精品99久久久久久久久| 亚洲成人中文字幕在线播放| 春色校园在线视频观看| 18+在线观看网站| 免费看a级黄色片| 亚洲av一区综合| 久久草成人影院| 成人亚洲精品一区在线观看 | 丝瓜视频免费看黄片| 亚洲成色77777| 一级毛片 在线播放| 2021少妇久久久久久久久久久| 日韩视频在线欧美| 免费观看a级毛片全部| 国精品久久久久久国模美| 精品久久久久久成人av| 亚洲精品aⅴ在线观看| 亚洲最大成人av| 成人鲁丝片一二三区免费| 欧美激情久久久久久爽电影| 国产精品人妻久久久影院| 成人av在线播放网站| 免费av观看视频| 久久99精品国语久久久| 国产黄色小视频在线观看| 久热久热在线精品观看| 黄色配什么色好看| 免费av毛片视频| 午夜精品在线福利| 国产色婷婷99| 国产亚洲5aaaaa淫片| 青青草视频在线视频观看| 国产69精品久久久久777片| 不卡视频在线观看欧美| 99热网站在线观看| 国产乱来视频区| 国产淫语在线视频| av在线观看视频网站免费| 青春草亚洲视频在线观看| 中文乱码字字幕精品一区二区三区 | 麻豆成人av视频| 看非洲黑人一级黄片| 嫩草影院入口| 国产伦精品一区二区三区四那| 国产精品福利在线免费观看| 午夜免费激情av| 美女被艹到高潮喷水动态| 亚洲自拍偷在线| 黑人高潮一二区| 久久久色成人| 黄色配什么色好看| 日韩欧美三级三区| 国产亚洲一区二区精品| 欧美日韩综合久久久久久| 中文精品一卡2卡3卡4更新| 久久精品夜色国产| 亚洲色图av天堂| 国产精品人妻久久久影院| 禁无遮挡网站| 伦精品一区二区三区| 亚洲国产精品专区欧美| 午夜爱爱视频在线播放| 波野结衣二区三区在线| 有码 亚洲区| 欧美3d第一页| 熟女电影av网| 久久综合国产亚洲精品| 真实男女啪啪啪动态图| 黄片无遮挡物在线观看| 国产老妇伦熟女老妇高清| 亚洲欧美成人精品一区二区| 日本色播在线视频| 中文乱码字字幕精品一区二区三区 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 蜜桃亚洲精品一区二区三区| 免费观看av网站的网址| 国产黄片视频在线免费观看| 免费观看精品视频网站| 少妇被粗大猛烈的视频| 在线 av 中文字幕| 国产精品一区二区三区四区免费观看| 夫妻午夜视频| 尤物成人国产欧美一区二区三区| 亚洲欧美精品专区久久| 精品久久久精品久久久| 高清视频免费观看一区二区 | 亚洲人成网站在线观看播放| 一级毛片 在线播放| 国产中年淑女户外野战色| 国产毛片a区久久久久| 亚洲精品色激情综合| 非洲黑人性xxxx精品又粗又长| 久久久久国产网址| 夫妻午夜视频| 观看免费一级毛片| 自拍偷自拍亚洲精品老妇| 美女内射精品一级片tv| 中国美白少妇内射xxxbb| 久久久成人免费电影| 国产人妻一区二区三区在| 国产一区二区三区综合在线观看 | 26uuu在线亚洲综合色| 大香蕉97超碰在线| 色综合色国产| 精品久久国产蜜桃| 十八禁国产超污无遮挡网站| av国产免费在线观看| 好男人在线观看高清免费视频| av在线亚洲专区| 日韩大片免费观看网站| 啦啦啦啦在线视频资源| 久久久久久九九精品二区国产| 免费观看无遮挡的男女| 国产伦理片在线播放av一区| 天堂√8在线中文| 一级av片app| 免费观看在线日韩| 亚洲精品久久午夜乱码| 波多野结衣巨乳人妻| 在线观看一区二区三区| 日韩 亚洲 欧美在线| 白带黄色成豆腐渣| 视频中文字幕在线观看| 秋霞伦理黄片| 国产精品国产三级国产av玫瑰| 国产69精品久久久久777片| 免费观看av网站的网址| 嫩草影院入口| 国产亚洲精品av在线| 青春草视频在线免费观看| 午夜福利高清视频| 欧美高清性xxxxhd video| 中文字幕av成人在线电影| 国产成人a∨麻豆精品| 亚洲人成网站在线观看播放| 七月丁香在线播放| 18禁动态无遮挡网站| 我的老师免费观看完整版| 亚洲精品久久午夜乱码| 日韩在线高清观看一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲美女搞黄在线观看| 日日啪夜夜撸| 国产黄色视频一区二区在线观看| 国产精品久久久久久久电影| 亚洲高清免费不卡视频| 久久亚洲国产成人精品v| 久久久久久久久久黄片| 欧美极品一区二区三区四区| 日韩欧美精品v在线| 男人和女人高潮做爰伦理| 少妇熟女aⅴ在线视频| 男女边吃奶边做爰视频| 国产高清三级在线| 97人妻精品一区二区三区麻豆| 五月天丁香电影| 黄色日韩在线| 最近中文字幕高清免费大全6| 国产欧美日韩精品一区二区| 国产成人免费观看mmmm| 欧美极品一区二区三区四区| 深夜a级毛片| 国产精品久久久久久精品电影小说 | 国产精品蜜桃在线观看| 免费高清在线观看视频在线观看| 久久久色成人| 国产av码专区亚洲av| 日韩成人av中文字幕在线观看| 成人欧美大片| 国产精品人妻久久久影院| 精品99又大又爽又粗少妇毛片| 日日啪夜夜撸| 亚洲精品第二区| 全区人妻精品视频| 在线观看免费高清a一片| 国产又色又爽无遮挡免| 国产成人福利小说| 久久久久久久亚洲中文字幕| www.色视频.com| 又爽又黄无遮挡网站| 纵有疾风起免费观看全集完整版 | 亚洲在久久综合| 色综合站精品国产| 国产一区有黄有色的免费视频 | 国产精品国产三级国产专区5o| 肉色欧美久久久久久久蜜桃 | 亚洲久久久久久中文字幕| 亚洲av国产av综合av卡| 国产精品熟女久久久久浪| 人人妻人人澡人人爽人人夜夜 | 日韩一区二区视频免费看| 亚洲国产精品专区欧美| 久久精品国产鲁丝片午夜精品| 美女大奶头视频| 国产麻豆成人av免费视频| 可以在线观看毛片的网站| 色综合亚洲欧美另类图片| 又爽又黄无遮挡网站| 亚洲精品视频女| 高清欧美精品videossex| 亚洲经典国产精华液单| 亚洲av一区综合| 国产午夜精品一二区理论片| 成人特级av手机在线观看| 男女那种视频在线观看| 亚洲av.av天堂| 欧美日韩亚洲高清精品| 国产av在哪里看| 成人性生交大片免费视频hd| 亚洲国产高清在线一区二区三| 性插视频无遮挡在线免费观看| 亚洲,欧美,日韩| 亚洲欧洲日产国产| 免费人成在线观看视频色| 国产精品一及| 国产在视频线在精品| eeuss影院久久| 干丝袜人妻中文字幕| 亚洲精品自拍成人| 狂野欧美白嫩少妇大欣赏| 国产一级毛片七仙女欲春2| 性色avwww在线观看| 一个人看的www免费观看视频| 搡老妇女老女人老熟妇| 国产一区二区三区综合在线观看 | 久久久久久国产a免费观看| 又黄又爽又刺激的免费视频.| 精品国产三级普通话版| 高清毛片免费看| av免费观看日本| 免费看av在线观看网站| 亚洲精品久久久久久婷婷小说| 熟女人妻精品中文字幕| 亚洲精品乱码久久久久久按摩| 久久久色成人| 最近中文字幕2019免费版| 69av精品久久久久久| 久久久久久久久久成人| 69人妻影院| 91狼人影院| a级毛片免费高清观看在线播放| 日韩电影二区| 日韩av免费高清视频| 日本一本二区三区精品| 97超视频在线观看视频| 国产乱人偷精品视频| 国产午夜精品久久久久久一区二区三区| 熟女人妻精品中文字幕| 神马国产精品三级电影在线观看| 麻豆精品久久久久久蜜桃| 婷婷色麻豆天堂久久| 精品人妻熟女av久视频| 国产精品1区2区在线观看.| 色播亚洲综合网| www.av在线官网国产| 一个人免费在线观看电影| 赤兔流量卡办理| 久久久久久久午夜电影| 十八禁国产超污无遮挡网站| 又粗又硬又长又爽又黄的视频| 亚洲在线观看片| 国产一区亚洲一区在线观看| 久久精品熟女亚洲av麻豆精品 | 伊人久久国产一区二区| 国产成人91sexporn| 亚洲av在线观看美女高潮| 国产午夜精品久久久久久一区二区三区| 建设人人有责人人尽责人人享有的 | 三级毛片av免费| 久久久久久九九精品二区国产| 亚洲欧美中文字幕日韩二区| 91午夜精品亚洲一区二区三区| 简卡轻食公司| 蜜臀久久99精品久久宅男| 在现免费观看毛片| 成人毛片60女人毛片免费| 日本一二三区视频观看| 1000部很黄的大片| 97精品久久久久久久久久精品| 久久热精品热| 大香蕉97超碰在线| 全区人妻精品视频| 国产一区二区在线观看日韩| 国产视频内射| 亚洲国产精品sss在线观看| 欧美3d第一页| 日本午夜av视频| 久久草成人影院| 国产 一区精品| 女人被狂操c到高潮| 久久国产乱子免费精品| 免费黄网站久久成人精品| 亚洲精品久久久久久婷婷小说| 国产乱来视频区| 国产午夜精品久久久久久一区二区三区| 成人鲁丝片一二三区免费| 日本色播在线视频| 国产在线一区二区三区精| 高清日韩中文字幕在线| 国产探花在线观看一区二区| 搡老乐熟女国产| 久久久久免费精品人妻一区二区| 国产av国产精品国产| 国产成人精品福利久久| 18禁裸乳无遮挡免费网站照片| 国产av码专区亚洲av| 亚洲最大成人手机在线| 男女啪啪激烈高潮av片| 成年免费大片在线观看| 看非洲黑人一级黄片| 亚洲图色成人| 亚洲自偷自拍三级| 午夜免费男女啪啪视频观看| 亚洲精品久久午夜乱码| 三级男女做爰猛烈吃奶摸视频| 最后的刺客免费高清国语| 夜夜看夜夜爽夜夜摸| 成年人午夜在线观看视频 | 亚洲高清免费不卡视频| 亚洲欧美一区二区三区黑人 | 男女边摸边吃奶| 日本黄色片子视频| 五月伊人婷婷丁香| 精品亚洲乱码少妇综合久久| 国产欧美另类精品又又久久亚洲欧美| 国产在线男女| 久久综合国产亚洲精品| 亚洲成色77777| 91久久精品国产一区二区成人| 80岁老熟妇乱子伦牲交| 国产亚洲5aaaaa淫片| 国产亚洲5aaaaa淫片| 久久综合国产亚洲精品| 国产色爽女视频免费观看| 精品国产三级普通话版| 午夜激情欧美在线| 亚洲成人中文字幕在线播放| videos熟女内射| 五月玫瑰六月丁香| 国产精品不卡视频一区二区| 两个人视频免费观看高清| 老司机影院毛片| 在现免费观看毛片| 自拍偷自拍亚洲精品老妇| 18禁动态无遮挡网站| 日日撸夜夜添| 欧美变态另类bdsm刘玥| 国产视频内射| 国产一区二区三区av在线| 久久久久久久国产电影| 国产精品人妻久久久影院| 久久这里有精品视频免费| 美女xxoo啪啪120秒动态图| 内射极品少妇av片p| 免费观看性生交大片5| 亚洲成人av在线免费| 亚洲av免费高清在线观看| 国产乱人视频| 欧美极品一区二区三区四区| 亚洲自偷自拍三级| 免费黄色在线免费观看| 日韩欧美精品v在线| 蜜桃亚洲精品一区二区三区| 国产黄频视频在线观看| videossex国产| 国产人妻一区二区三区在| 人人妻人人澡人人爽人人夜夜 | 男插女下体视频免费在线播放| 网址你懂的国产日韩在线| 观看美女的网站| a级一级毛片免费在线观看| 日日撸夜夜添| 国产乱人偷精品视频| 日本wwww免费看| 亚洲国产成人一精品久久久| av免费观看日本| 亚洲成人一二三区av| 自拍偷自拍亚洲精品老妇| 最近的中文字幕免费完整| 成人毛片60女人毛片免费| 一本一本综合久久| 一级毛片黄色毛片免费观看视频| 国产精品久久视频播放| 成人亚洲精品一区在线观看 | 特级一级黄色大片| 狂野欧美激情性xxxx在线观看| 国产欧美日韩精品一区二区| 欧美精品国产亚洲| 丝袜喷水一区| 一边亲一边摸免费视频| 内射极品少妇av片p| 一本久久精品| 午夜日本视频在线| av线在线观看网站| 丰满乱子伦码专区| 少妇高潮的动态图| 在线 av 中文字幕| 精华霜和精华液先用哪个| 日韩欧美国产在线观看| 激情 狠狠 欧美| .国产精品久久| 国产成人aa在线观看| 亚洲精品日韩av片在线观看| 国产 一区 欧美 日韩| 成人特级av手机在线观看| 女人久久www免费人成看片| 国产片特级美女逼逼视频| 亚洲国产精品sss在线观看| 最近最新中文字幕大全电影3| 亚洲最大成人av| 免费av不卡在线播放| 欧美xxxx性猛交bbbb| 日韩精品有码人妻一区| 联通29元200g的流量卡| 国产av国产精品国产| 国产精品一区二区性色av| 美女高潮的动态| 国产 亚洲一区二区三区 | 一级毛片aaaaaa免费看小| 国产黄片视频在线免费观看| 欧美xxⅹ黑人| 国产v大片淫在线免费观看| 亚洲一级一片aⅴ在线观看| 中国美白少妇内射xxxbb| 精品久久久精品久久久| 国产黄频视频在线观看| 国产高清三级在线| 亚洲精品久久午夜乱码| 人妻夜夜爽99麻豆av| 99久国产av精品| 又粗又硬又长又爽又黄的视频| 国产成人freesex在线| 禁无遮挡网站| 日本色播在线视频| 免费人成在线观看视频色| 欧美zozozo另类| 1000部很黄的大片| 国产91av在线免费观看| 久久久精品94久久精品| 久久久久久久久久人人人人人人| 国产精品国产三级国产专区5o| 亚洲av.av天堂| 91狼人影院| 国内精品宾馆在线| 亚洲国产精品sss在线观看| 黄色欧美视频在线观看| 一夜夜www| 国产女主播在线喷水免费视频网站 | 一区二区三区高清视频在线| 亚洲久久久久久中文字幕| 一级毛片电影观看| 成人亚洲欧美一区二区av| 99热6这里只有精品| 亚洲伊人久久精品综合| 日韩成人伦理影院| 日本欧美国产在线视频| 看非洲黑人一级黄片| 亚洲av中文字字幕乱码综合| 欧美日韩国产mv在线观看视频 | 国产永久视频网站| 久久精品国产鲁丝片午夜精品| 舔av片在线| 熟女电影av网| 91精品国产九色| 色哟哟·www| av天堂中文字幕网| 麻豆av噜噜一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人亚洲精品一区在线观看 | 亚洲精品久久午夜乱码| 国产高潮美女av| 国产久久久一区二区三区| 丰满乱子伦码专区| 国产熟女欧美一区二区| 极品教师在线视频| 色综合站精品国产| 一级毛片黄色毛片免费观看视频| 蜜桃亚洲精品一区二区三区| 一级a做视频免费观看| h日本视频在线播放| 99热这里只有是精品在线观看| 五月玫瑰六月丁香| 成年版毛片免费区| 国产视频内射| 亚洲熟妇中文字幕五十中出| 国产成人精品福利久久| 麻豆乱淫一区二区| 亚洲av免费在线观看| 色综合色国产| a级毛色黄片| 777米奇影视久久| 三级经典国产精品| 美女高潮的动态| 国产黄a三级三级三级人| av天堂中文字幕网| 国产精品伦人一区二区| 国产黄a三级三级三级人| 在线a可以看的网站| 啦啦啦啦在线视频资源| 精品久久久久久久久av| 人体艺术视频欧美日本| 一级av片app| 丰满乱子伦码专区| 成年女人看的毛片在线观看| 精品久久久久久成人av| 深爱激情五月婷婷| av网站免费在线观看视频 | 精品久久久久久久久亚洲| 亚洲天堂国产精品一区在线| 国产成人福利小说| 一个人看视频在线观看www免费| 欧美xxxx黑人xx丫x性爽| 边亲边吃奶的免费视频| 成人鲁丝片一二三区免费| 美女国产视频在线观看| 国产欧美日韩精品一区二区|