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

    Simultaneous optimization of even flow and land and timber value in forest planning:a continuous approach

    2021-02-28 09:13:22JosGonzlezGonzlezMiguelzquezndezandUlisesDiguezAranda
    Forest Ecosystems 2021年4期

    José M.González-González,Miguel E.Vázquez-Méndez and Ulises Diéguez-Aranda*

    Abstract

    Keywords: Even-aged forest management,Forest-level optimization,Continuous optimization,Gradient-type algorithms

    Background

    A forest can be considered a set of stands, which are contiguous communities of trees sufficiently uniform in composition, age-class distribution and site quality,to distinguish each community from its adjacent ones(Helms 1998). In forest planning, decision-making processes often use optimisation approaches for developing optimal harvest schedules that will best meet the objectives of the landowners or land managers (Kaya et al.2016).The main issue consists in determining when(and how) to harvest each stand during the planning horizon(nextPyears),to achieve the proposed objectives.To formulate this problem, it is necessary to have a simulator that allows predicting the development of forest stands efficiently and accurately.In the last decades,many papers have been devoted to this topic (see, for instance García 1994; Castedo-Dorado et al. 2007; álvarez-González et al. 2010), and for any well-know species can be assumed that output functions are available,given,for example,the timber volumeVj(y)(volume/area) of the standjalong the time. To clearly show the novelty and interest of this work,we begin by analyzing different formulations of the forest planning problem in the optimization framework.For simplicity,we assume that maximizing timber volume is the main objective.Without any additional hypothesis,denotingAjthe area of standj, the problem is to determine the harvest ratehrj(y)(area/time)of each stand over time(Heaps 1984;2015),such that

    Within the framework of the binary linear programming(BLP), the exact solution of (4) can be obtained by the branch and bound or cutting plane methods,as long as the problem does not exceed a certain relatively small dimension (Bettinger et al. 1999). For large problems, a spatial application of the RC method can be also useful(Pukkala et al.2009).Additionally,problem(4)also accepts a combinatorial formulation in the frame of non-linear integer programming(NLIP),by considering the harvest instants y =(y1,...,yns)as the decision variables.In this case,(4)is equivalent to find vector y such that

    This formulation is more appropriate if a heuristic method is used to solve the problem (which generates a solution within a reasonable amount of time, although of uncertain quality(Davis et al.2001,p.741;Bettinger et al.2009,p.172).

    If hypothesis H1 is not considered,then H2 is equivalent to assume thathrj(y)=Ajδ(y-yj),and without any additional assumption,problem(1)changes into the following non-linear problem(NLP)Now,hypothesis H1 can be understanding as a discretization of the continuous variable of this NLP, and if it is assumed (what may be necessary if the dimension of the problem –number of stands–, is high and the available methods do not work in a reasonable time), problem (6) becomes straightly the previous NLIP problem(5). Equivalences and transformations between all these formulations are summarized in Fig.1.

    Fig.1 Equivalences and transformations between different formulations of the forest planning basic problem(harvest scheduling problem)

    Problems(3)-(6)can be improved by considering more sophisticated goals,such as land and timber value(LTV).They can also be completed with additional objectives and/or constraints related to even flow of products (EF),size of contiguous area harvested(spatial constraints),etc.There are many papers dealing with models based on formulations (3)-(5), and Table 1 includes some examples classified by the formulation type,and the objectives and constraints considered. On the contrary, due to the inheritance of the original LP formulation, and to the apparent absence of suitable methods for numerical resolution,models based on the continuous formulation given in problem (6) have been much less explored (see, for instance, Roise 1990). The aim of this work is to help fill this gap in the forestry literature, by presenting a novel continuous bi-objective model in forest planning.

    In recent papers,Arias-Rodil et al.(2017)and González-González et al. (2020) presented rigorous mathematical analysis of two continuous metrics for measuring,respectively,land expectation value(LEV)and EF.In this work,without losing their good properties for working with gradient-type methods, these metrics are appropriately modified to deal with the scheduling harvest problem of a forest, assuming future rotations based on economically optimal management prescriptions at stand level. These modifications include a novel way to estimate the planning horizon in an automatic manner, and the possibility that a stand is harvested several times.The combination of these new metrics leads to a biobjective model,which can be efficiently solved by using gradient-type techniques, does not require to set management prescriptions in advance,and avoids the division of the planning horizon into periods, therefore resulting an useful tool for the decision-making processes in forest planning.

    The paper is organized as follows.In Methods section,the new metrics are detailed in a suitable mathematical framework, and the forest planning problem is formulated by means of the new continuous bi-objective model,which is studied from a cooperative point of view,proposing a numerical method to obtain its Pareto-optimal frontier.The efficiency and accuracy of the numerical method is shown in Results and discussion section,where results in someEucalyptus globulusLabill. forests are presented and discussed. Finally, some interesting conclusions are summarized in Conclusions section.

    Methods

    Mathematical framework

    To formulate the forest planning problem, we need a growth simulator which allows prediction of stands development. To that end, we used the dynamic systemsbased framework (frequently referred to as “state-space”approach),first used in forestry by García(1994).Particularly,our approach is based on:

    (i) Each standjis characterized by three state variables:mean height of dominant trees in the stand(Hj,inmeters),number of stems per hectare(Nj),and stand basal area(Gj,in m2/ha).

    Table 1 Literature review:some references dealing with optimization models for forest planning.They are grouped by the formulation type and classified according to the objectives(columns 3-5)and constraints(columns 6-8)considered

    (ii) Outputs can be obtained from the predicted values of these state variables.For example,the timber volume(in m3/ha)is given by a known species-specific functionvj(H,N,G).

    Economic metric:the land and timber value(LTV)

    The economic value of the forest is the sum of the value of each standj,which depends on:

    wherer∈(0,1]is the interest rate,andRj(ˉtj)andCj(ˉtj)represent,respectively,functions of present values of revenues and costs.

    The land and timber value (LTV) can be defined as the present value of the sum of the timber value at the first clearcutting and the LEV just after clearcutting. For each standj, it depends not only on the time of the first clearcutting (yj), but also on the future rotation age (ˉtj).We only seek foryjand, consequently, ˉtjhas to be set in advance.We assume future rotations based on the optimal management prescription at stand level,in such a way thatˉtjis the solution of the following optimization problem

    The regularity properties ofJLTVare given by the properties of functionshj,nj,gjandvj,and we can assume that it is smooth enough(Arias-Rodil et al.2017).

    Even flow(EF)metric

    Even flow (EF) and sustained yield are two of the oldest objectives of forest management. The first has been related to the concept of maintaining a stable timber supply (nowadays provisioning ecosystem services) through appropriate forest planning and management techniques.There are numerous examples of how to address the even flow issue in harvest scheduling. Some studies use constraints that ensure volume levels in each period within some upper and lower limits (O’Hara 1989; Murray and Church 1995) or that allow a certain variation in volume between periods(McDill et al.2002;Constantino et al. 2008), while others include even flow as an objective in which volume variations between periods should be minimized (Kao and Brodie 1979; Brumelle et al. 1998;Ducheyne et al.2004).

    In this study we recover the concept of harvest rate(now volume/time) and seek for a constant rate during a period ofTyears, named tentative planning horizon for even-flow. It represents the period for which, clearcutting all stands at least once time, a stable timber supply is desired. The value ofTis a technical decision, which depends on the species,site qualities,stand ages...of the forest.

    We consider that interventions may not begin at the present time if, for example, the forest is too young. We denote bylj >0, in years, the minimum harvest age for each standj, in such a way thattj≥ljshould be verified. From these bounds, the instanta≥0 (years) when harvests can begin in the forest is given:it will be zero if any stand can be harvested now, otherwise it will be the minimum time guaranteeingtj≥lj.Specifically:

    Taking into account that a management prescription can be greater thana+T, and assuming that all stands have to be harvested at least once during the planning horizon,its upper limit must be automatically computed by

    whereH(y)is the Heaviside function(Oldham et al.2009,Ch.9).Therefore,from the EF point of view,the objective is to achieve the next goal harvest rate function(see Fig.2)

    Fig.2 Example of harvest rate(top)and volume harvested(bottom).The grey solid lines represent real instant harvests(top,Eq.(10))and real volumes harvested(bottom,Eq.(12)),while the black dashed lines represent goal harvest rate functions(top,Eq.(9))and goal volume harvested functions(bottom,Eq.(11))

    The Dirac delta is not a function in the traditional sense,so the comparison between expressions(9)and(10)is not appropriate. Instead of comparing harvest rates we propose to compare the volumes harvested (goal and real),which are obtained by integrating those expressions (9)and(10):

    The EF metric must measure the distance between both functions for the planning horizon[a,b(y)].By using the associated distance to theL2norm(Adams 1975),for any planning strategy y ∈Yad,it is given by

    where the minus sign is included precisely so that a higher value ofJEF(y)corresponds to a greater similarity of the functionV(y,.)to the goal functionVg(y,.),which is what is meant by a higher even flow.

    Following (González-González et al. 2020), it can be proved that functionJEFhas continuous derivatives in almost all points,and an explicit expression for its gradient can be obtained.

    Bi-objective model

    In this work, it is assumed that the main goal of forest planning is to determine the management prescriptions to be applied to each stand (i.e., to determine y ∈Yad) to get the best profitability(i.e.,to maximizeJLTV)with the highest possible even flow(i.e.,to maximizeJEF).Hence,the following Bi-objective Optimization Problem(BOP)is considered:

    This problem must be studied from a cooperative point of view, looking for non-dominated admissible vectors,also known as efficient solutions or Pareto optima. An admissible vector y ∈Yadis said dominated if there exists a better admissible vector, in the sense that it improves at least one objective, without making the other worse.If y ∈Yadis non-dominated (Pareto optimum), the corresponding objective vector J(y)∈R2is also known as Pareto-optimal, and the set of Pareto-optimal objective vectors is known as Pareto frontier (see Fig. 3). A more detailed definition can be seen,for instance,in Miettinen(1998,p.11).

    As commented above,the final aim is to determine the management prescription to be applied to each stand.Therefore, a satisfactory solution must be chosen within the set of Pareto optima.A satisfactory solution is the best compromise for the preferences of the decision maker.These preferences may be articulated prior to the analysis(allowing to reformulate the BOP into a single objective problem through weights or/and treating some objectives as constraints),after graphing the Pareto frontier,or interactively in a progressive way (requiring to compute only some Pareto optima) (Miettinen 1998, pp. 61-65). In the next section we propose a numerical method for obtaining the set of Pareto optima and graphing the Pareto frontier.It allows to articulate preferences a posteriori,resulting in a very interesting tool for the decision-making process in forest planning.

    Numerical resolution

    Regarding BOP (13), to take advantage of the regularity properties ofJLTVandJEF, we propose to obtain the Pareto frontier by the weighting method (Miettinen 1998, pp. 78-84), solving each single-optimization problem by a gradient-type method, combined with a random multi-start if local minima are detected. It is also advisable to normalize the objective functions so that their objective values are of approximately the same magnitude. For the results presented in the next section, both objectives were normalized using the ideal objective vector (that obtained by maximizing each of the objective functions individually subject to the constraints of the problem) and an estimation of the nadir objective vector (given by the value of each objective function at the point where the other objective function reaches its maximum) (Miettinen 1998, pp. 15-19),weights were taking equally spaced, and the L-BFGS-B algorithm was used again without any multi-start procedure. To check this approach, BOP (13) was also solved by the evolutionary algorithm NSGA-II(Deb et al.2002), widely used in multi-objective optimization, and already implemented in the Python library Inspyred 1.0(Tonda 2020).

    Fig.3 Example of an admissible set Yad ?R3(corresponding with ns =3)and an image set J(Yad)?R2,where the Pareto frontier is highlighted

    Results and discussion

    We analyzed the usefulness of our approach in a real forest ofEucalyptus globulusLabill. located in the municipality of Xove (Galicia, NW Spain), considering a tentative planning horizon for even flow ofT= 13.5 years. The forest comprisesns= 51 stands whose area and inventory data can be seen in Table 2. From these data, the state variables until clearcutting are computed from the specific transition functions given by (García-Villabrille 2015).Specifically:

    After clearcutting,we assume no changes on site quality (site index) and a constant plantation density of 1333 trees/ha for all stands.Under this hypotheses,by using an appropriate function for initializing the stand basal area(García-Villabrille 2015),we predicted the data at one year after clearcutting(see Table 2).These data were used for computing the state variables (and outcomes) after the first clearcutting. To further test our approach in a forest with a different structure,we also used these data for simulating a hypothetical forest(hereinafter the simulated young forest)with 51 stands of the same species,area and site index,but all at one-year age.

    For the real forest of the case study, Fig. 4 shows the Pareto fronts obtained with the weighting method and with the NSGA-II, for a population size ofP=200 individuals andG=500 generations. Additionally, this figure also shows the comparison between real and goal volumes corresponding to all (ten) Pareto optima obtained with the weighting method. The corresponding results for the simulated young forest can be seen in Fig. 5. In view of these results, it can be highlighted that:

    · Both methods,although using different techniques,provide the same Pareto fronts,which warrants the appropriateness of our approach.

    · The continuous approach works well for forest with very different stand structures,as suggested from the results of the real forest(which has many mature stands of different ages)and the simulated young forest(in which all stands are of the same age).

    · Figures 4 and 5 are very useful for the decision-making process:

    – In the top graphs,the normalized-JLTVisJLTV/JLTV(y*),i.e.,the LTV divided by its best value.This way of proceeding allows an intuitive comparison of the different optimal management strategies:for the real forest,there was a reduction of about 8%in LTV between the best solutions from the economic and even flow perspectives(points 1 and 10 in the Pareto front,respectively),while for the simulated young forest it was of only 5%.

    Table 2 Inventory data and hypothetical predicted data at one year after clearcutting for the forest of the case study.Aj:stand area in hectares,t0j:stand age in years,H0j:dominant height in meters,N0j:number of trees per hectare,G0j:stand basal area in m2/ha(Continued)

    – The even flow metric has not such an easy interpretation in terms of numerical values,but the grey area of the bottom graphs of the figures clearly indicates the achievement of the even flow objective as well as the moment in which the clearcuts are done.Additionally,these graphs also show the total harvest volume of each Pareto optimum,which is reduced as the even flow objective is improved.

    Finally, to analyze our continuous approach in larger forests, we simulated two new hypothetical forests: one with 204 stands and other with 1632 stands. They were obtained by replicating the stands of the real forest 4 and 32 times,respectively,and randomly generating their area between the minimum and maximum real-area values.Then,we solved the problem on the simulated forest with 204 stands and observed that the weighting method captured the complete Pareto front in a very reasonable time (about half a minute), but the NSGA-II needed a greater number of generations than the previous experiments with the 51-stand forests and, consequently, a higher computation time.This fact was confirmed by the results obtained on the simulated forest with 1632 stands,where the weighting method kept providing the complete Pareto front in an acceptable time(about seven minutes),but the NSGA-II withG=2000 generations only captured a small arc(see Fig.6),with a computation time of about eleven hours. Table 3 summarizes the computation time of all performed experiments.

    The use of non-linear continuous models to formulate the forest management problems has been scarce.In the framework of optimal control theory, Heaps(1984; 2015) studied the general problem of determining the harvest rate of each stand over time, but did not address with examples how the continuous forestry age class model he proposed could be solved numerically. Tahvonen (2004) discussed the use of discrete or continuous variables in describing time and age class structure in models for optimization of forest harvesting.Roise (1990) proposed and solved a continuous model,but without performing a mathematical analysis of the objective and constraint functions,applying optimization techniques that may not be the most adequate to solve the problem.On the contrary,the continuous model proposed in this paper uses functions with appropriate regularity properties, that guarantee a successful numerical resolution using gradient-type methods.

    In this study we have dealt with a species whose main use is pulp production, therefore the land and timber value has been calculated on the basis of the stumpage price times the volume at the harvest instant. Nevertheless,for other species which require merchantable volume estimates,a more sophisticated method must be used.For example, in a study of optimization at stand level, Arias-Rodil et al. (2015) evaluated two alternative methods of estimating the merchantable volume of a stand at a given age: (i) a disaggregation system, which involves aggregating the merchantable volume previously predicted by diameter classes,and(ii)a stand volume ration function.The latter considers a top diameter limit for classifying the volume by timber assortments,while the former also allows specification of long length requirements.The use of alternative (i) may cause discontinuities in theJLTVfunction, which makes a gradient-based search challenging.However,alternative(ii)does not present any problem and,although will provide slightly higher optimal values ofJLTV,may even be more accurate(Arias-Rodil et al.2016).

    Fig.4 Results for the real forest:normalized Pareto front obtained with the weighting method and the NSGA-II(top),and comparison between real and goal volumes for the points obtained with the weighting method(bottom)

    Fig.5 Results for the simulated young forest:normalized Pareto front obtained with the weighting method and the NSGA-II(top),and comparison between real and goal volumes for the points obtained with the weighting method(bottom)

    Fig.6 Normalized Pareto fronts obtained with the weighting method(·)and the NSGA-II,for two simulated forest comprising 204 stands(top)and 1632 stands(bottom).The Pareto fronts obtained with NSGA-II correspond with a population size of P=200 individuals and a number of generations of G=500(dashed line)and G=2000(solid lines)

    Our proposal provides a graphical tool to analyze the relationship(trade-off)between the two objectives,helping the landowners to choose a suitable weighting after seeing the effect of even-flow requirement on economic profitability. This methodology can include constraints and be extended to other objectives, in pairwise comparison or considering three or more objectives. In the latter case,despite a posteriori techniques combined with special visualization methods have been used in forest management problems (Borges et al. 2014), interactive techniques can also be useful (Diaz-Balteiro 1998). In any case, the features of the problem and the capabilities and type of the decision maker have to be considered before selecting the most appropriate solution method(Miettinen 1998,pp.227-231).

    Table 3 Computation times(seconds)for obtaining the Pareto fronts shown in Figures 4,5 and 6 on an Intel?CoreTM i7 Processor 4 GHz iMac.Comparison,on forest with different structure or number of stands,between the weighting method and the NSGA-II with a population size of P=200 individuals and different number of generations(G)

    Conclusions

    In this paper we present a novel continuous bi-objective model (economic vs. even flow) in forest planning.According to the results obtained,it is worth highlighting the following aspects:

    · The planning horizon is not set in advance and only a tentative value is required.It is involved in the optimization process and obtained as an output of the model.

    · The LTV and EF functions have smooth enough properties,which allows the use of gradient-type methods for its optimization,behaving in a more efficient,effective and robust way than many other gradient-free methods commonly applied in forestry.

    · The continuous formulation does not require to set management prescriptions in advance and avoids the division of the planning horizon into periods,providing solutions which may be better than those obtained with traditional combinatorial formulations.

    · This model allows forest-level problems to be treated with continuous decision variables,in a way similar to that used generally with stand-level problems.As envisaged in previous works,the performance of the continuous approach combined with gradient-type methods improved as the forest size(the number of decision variables)increased.

    · With this approach,it is not necessary to articulate the preferences of the decision maker a priori.Furthermore,the graphical outputs of our numerical method provide the Pareto front and other useful trade-off information in an intuitive way,which makes the proposed continuous bi-objective model a very interesting tool for the decision-making process in forest planning.

    · This model was successfully used in forests comprising single-species,even-aged stands which are intended to be managed with the rotation forest management system.The applicability of the model in forests with more complex structures and dynamics(forests consisting of mixed stands,which have gradual advance regeneration,and are treated with partial cuttings)is the topic of an ongoing research.

    Abbreviations

    BLP:Binary linear programming;BOP:Bi-objective optimization problem;EF:Even flow;LEV:Land expectation value;LP:Linear programming;LTV:Land and timber value;NLIP:Non-linear integer programming;NLP:Non-linear problem;NW:North West

    Acknowledgements

    Not applicable.

    Authors’contributions

    All authors have contributed similarly in research planning,analysis,and writing of the manuscript.

    Authors’information

    JMG-G is a PhD Student focused in forest management planning with optimization methods.MEV-M is Associate Profesor and his research interests are optimal control of partial differential equations and numerical optimization with applications in engineering and environment.UD-A is Associate Profesor interested on quantitative methods for forest management(growth modelling,forest inventory and forest management planning).Currently,the three are working together on forest management planning within the framework of continuous optimization.

    Funding

    Not applicable.

    Availability of data and materials

    All data generated or analyzed during this study are included in this published article.

    Declarations

    Ethics approval and consent to participate

    Not applicable.

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Author details

    1 Departamento de Enxe?aría Agroforestal,Universidade de Santiago de Compostela.Escola Politécnica Superior de Enxe?aría,R/Benigno Ledo,27002 Lugo,Spain.2Departamento de Matemática Aplicada,Instituto de Matemáticas,Universidade de Santiago de Compostela.Escola Politécnica Superior de Enxe?aría,R/Benigno Ledo,27002 Lugo,Spain.

    Received:9 April 2021 Accepted:1 June 2021

    日日干狠狠操夜夜爽| 国产欧美日韩一区二区三| 我的老师免费观看完整版| 757午夜福利合集在线观看| 波多野结衣高清无吗| 久久久久免费精品人妻一区二区| 少妇熟女aⅴ在线视频| 婷婷丁香在线五月| 欧美成人性av电影在线观看| 老司机福利观看| 精品日产1卡2卡| 亚洲av成人一区二区三| 国内精品久久久久精免费| 午夜激情福利司机影院| 一本综合久久免费| 成人手机av| 成人高潮视频无遮挡免费网站| 一二三四在线观看免费中文在| 999精品在线视频| 此物有八面人人有两片| 婷婷六月久久综合丁香| 露出奶头的视频| 久久精品91无色码中文字幕| 真人一进一出gif抽搐免费| 黄片大片在线免费观看| 午夜精品在线福利| 亚洲av日韩精品久久久久久密| 国产精品99久久99久久久不卡| 欧洲精品卡2卡3卡4卡5卡区| 国产v大片淫在线免费观看| 国产精品av久久久久免费| 最近最新中文字幕大全免费视频| 久久中文字幕人妻熟女| 亚洲精品中文字幕一二三四区| 国产精品 国内视频| 日韩高清综合在线| 久久香蕉激情| 久久精品国产99精品国产亚洲性色| 精品欧美一区二区三区在线| 人妻久久中文字幕网| 免费看十八禁软件| 日日干狠狠操夜夜爽| АⅤ资源中文在线天堂| 久久精品91无色码中文字幕| 两性夫妻黄色片| 这个男人来自地球电影免费观看| 三级毛片av免费| 久久天躁狠狠躁夜夜2o2o| 99热6这里只有精品| 后天国语完整版免费观看| 国产熟女午夜一区二区三区| 又粗又爽又猛毛片免费看| 在线观看免费视频日本深夜| 午夜福利欧美成人| 国产91精品成人一区二区三区| 国产91精品成人一区二区三区| 国产私拍福利视频在线观看| 国产精品九九99| 18禁裸乳无遮挡免费网站照片| 亚洲色图 男人天堂 中文字幕| 国产成人aa在线观看| 亚洲国产精品久久男人天堂| 成人18禁高潮啪啪吃奶动态图| 亚洲男人的天堂狠狠| 啦啦啦韩国在线观看视频| 久久 成人 亚洲| 亚洲欧美精品综合一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 深夜精品福利| 日韩大码丰满熟妇| 成人亚洲精品av一区二区| 正在播放国产对白刺激| 久久久国产成人精品二区| 亚洲中文日韩欧美视频| 精品福利观看| 此物有八面人人有两片| 久久中文看片网| 两个人免费观看高清视频| 99国产极品粉嫩在线观看| 国产一区二区激情短视频| 成人av一区二区三区在线看| 黄色成人免费大全| 无限看片的www在线观看| 老司机在亚洲福利影院| 悠悠久久av| 国产高清激情床上av| 亚洲成人国产一区在线观看| 黄色女人牲交| 老汉色av国产亚洲站长工具| 免费高清视频大片| 国内毛片毛片毛片毛片毛片| 日韩大码丰满熟妇| 午夜精品久久久久久毛片777| 欧美日韩一级在线毛片| 欧美另类亚洲清纯唯美| 国产熟女xx| 人妻夜夜爽99麻豆av| 香蕉av资源在线| 琪琪午夜伦伦电影理论片6080| 久久精品91无色码中文字幕| 桃色一区二区三区在线观看| 久久久久精品国产欧美久久久| 日日爽夜夜爽网站| 熟女电影av网| 男女午夜视频在线观看| 波多野结衣高清作品| 国产一区二区在线av高清观看| 黄片小视频在线播放| 黄色女人牲交| 国内毛片毛片毛片毛片毛片| 国产精品自产拍在线观看55亚洲| 午夜激情福利司机影院| 中文字幕熟女人妻在线| 中国美女看黄片| 两人在一起打扑克的视频| 国产精品av视频在线免费观看| 他把我摸到了高潮在线观看| 国产av不卡久久| 午夜福利在线观看吧| 在线观看美女被高潮喷水网站 | 午夜福利高清视频| 麻豆成人av在线观看| 久久国产精品影院| 一级黄色大片毛片| 蜜桃久久精品国产亚洲av| 亚洲欧洲精品一区二区精品久久久| 成人高潮视频无遮挡免费网站| 国产亚洲精品久久久久久毛片| 99国产极品粉嫩在线观看| 精品电影一区二区在线| 国产av在哪里看| 2021天堂中文幕一二区在线观| 婷婷六月久久综合丁香| 精品国产美女av久久久久小说| 国产高清激情床上av| 18禁国产床啪视频网站| 欧美性猛交黑人性爽| 久久久久久久午夜电影| 欧美精品亚洲一区二区| 熟女电影av网| 久久久久免费精品人妻一区二区| 亚洲五月婷婷丁香| 青草久久国产| 熟女少妇亚洲综合色aaa.| 日日爽夜夜爽网站| 亚洲欧美精品综合一区二区三区| 欧美三级亚洲精品| 亚洲全国av大片| 国产精品一及| 黄色视频,在线免费观看| 欧美成狂野欧美在线观看| 中文字幕精品亚洲无线码一区| 又爽又黄无遮挡网站| 午夜久久久久精精品| 真人做人爱边吃奶动态| 熟女电影av网| 国产真人三级小视频在线观看| 亚洲av电影在线进入| 后天国语完整版免费观看| 黑人欧美特级aaaaaa片| 两性午夜刺激爽爽歪歪视频在线观看 | 男人的好看免费观看在线视频 | 国产乱人伦免费视频| 国产高清视频在线播放一区| 99国产极品粉嫩在线观看| 亚洲,欧美精品.| 在线观看午夜福利视频| 亚洲国产中文字幕在线视频| 看片在线看免费视频| 精品一区二区三区视频在线观看免费| 久久人人精品亚洲av| АⅤ资源中文在线天堂| 国产精品自产拍在线观看55亚洲| 可以在线观看的亚洲视频| 国内久久婷婷六月综合欲色啪| 久久久久久久久中文| 曰老女人黄片| 亚洲av成人不卡在线观看播放网| 国产蜜桃级精品一区二区三区| 国产伦人伦偷精品视频| 欧美丝袜亚洲另类 | 国产精品九九99| 国产午夜精品久久久久久| 亚洲精品久久国产高清桃花| 三级毛片av免费| 91大片在线观看| 麻豆成人av在线观看| 日韩欧美国产一区二区入口| 长腿黑丝高跟| 丝袜美腿诱惑在线| 男插女下体视频免费在线播放| 日韩欧美在线二视频| 熟妇人妻久久中文字幕3abv| ponron亚洲| 国产高清视频在线观看网站| 亚洲黑人精品在线| 老司机福利观看| 国产午夜精品论理片| 久久婷婷成人综合色麻豆| 久久人人精品亚洲av| 亚洲精品粉嫩美女一区| 香蕉国产在线看| 夜夜夜夜夜久久久久| 国产av一区在线观看免费| 一区二区三区国产精品乱码| 黑人巨大精品欧美一区二区mp4| 日韩成人在线观看一区二区三区| 美女黄网站色视频| 日韩av在线大香蕉| 国产单亲对白刺激| 久久久久免费精品人妻一区二区| 国产亚洲精品第一综合不卡| 久久久水蜜桃国产精品网| 女警被强在线播放| 国产激情欧美一区二区| 美女 人体艺术 gogo| 国产欧美日韩一区二区精品| 欧洲精品卡2卡3卡4卡5卡区| 午夜日韩欧美国产| 日韩精品青青久久久久久| 欧美黄色片欧美黄色片| 欧美色视频一区免费| 国产久久久一区二区三区| 国产精品98久久久久久宅男小说| 国产成年人精品一区二区| 中亚洲国语对白在线视频| 最新在线观看一区二区三区| 成人一区二区视频在线观看| 999久久久国产精品视频| 亚洲av成人av| 久久精品91无色码中文字幕| 性欧美人与动物交配| 9191精品国产免费久久| 91字幕亚洲| 琪琪午夜伦伦电影理论片6080| 日韩高清综合在线| 国产又黄又爽又无遮挡在线| 亚洲一区中文字幕在线| 午夜a级毛片| 久久精品国产亚洲av高清一级| 亚洲成av人片免费观看| 亚洲一区高清亚洲精品| 亚洲午夜理论影院| 欧美日韩黄片免| 黄色视频不卡| 国产精品久久视频播放| 美女大奶头视频| 88av欧美| 午夜福利视频1000在线观看| 午夜成年电影在线免费观看| 亚洲国产中文字幕在线视频| av福利片在线| 亚洲在线自拍视频| 国产午夜福利久久久久久| 波多野结衣高清作品| 国模一区二区三区四区视频 | 成人av在线播放网站| 久久久久久久久免费视频了| 精品福利观看| 欧美黑人精品巨大| 18美女黄网站色大片免费观看| 日本在线视频免费播放| 亚洲欧美精品综合一区二区三区| 99国产极品粉嫩在线观看| 亚洲中文av在线| 成人手机av| 亚洲第一欧美日韩一区二区三区| 中文字幕高清在线视频| 伊人久久大香线蕉亚洲五| 搡老妇女老女人老熟妇| 亚洲一区高清亚洲精品| 99国产精品一区二区三区| 国产成人系列免费观看| 日日干狠狠操夜夜爽| 成人三级做爰电影| 国产精品亚洲一级av第二区| 在线永久观看黄色视频| 男插女下体视频免费在线播放| 久久久久精品国产欧美久久久| 老熟妇乱子伦视频在线观看| 久久精品国产亚洲av香蕉五月| 成人午夜高清在线视频| 亚洲av电影不卡..在线观看| 成人国产综合亚洲| 久久精品国产亚洲av香蕉五月| 中亚洲国语对白在线视频| 一本精品99久久精品77| 精品不卡国产一区二区三区| 日本精品一区二区三区蜜桃| 99在线视频只有这里精品首页| 又黄又粗又硬又大视频| 欧美性猛交黑人性爽| 淫妇啪啪啪对白视频| 国产又色又爽无遮挡免费看| 91麻豆av在线| 在线视频色国产色| av片东京热男人的天堂| 欧美又色又爽又黄视频| 日本一区二区免费在线视频| 国内少妇人妻偷人精品xxx网站 | 国产av一区在线观看免费| 搡老岳熟女国产| 狂野欧美白嫩少妇大欣赏| 久久精品91蜜桃| 亚洲精华国产精华精| 不卡av一区二区三区| 精品久久久久久久毛片微露脸| 久久久久久久久免费视频了| 久久久久久久久久黄片| 日本一本二区三区精品| 他把我摸到了高潮在线观看| 在线播放国产精品三级| 中文字幕人成人乱码亚洲影| 99久久综合精品五月天人人| 97人妻精品一区二区三区麻豆| 精品一区二区三区四区五区乱码| 欧美+亚洲+日韩+国产| 午夜福利在线观看吧| 精品国内亚洲2022精品成人| www.自偷自拍.com| 亚洲精品国产精品久久久不卡| 久久久精品国产亚洲av高清涩受| 国产精品一区二区三区四区免费观看 | 亚洲国产精品999在线| 露出奶头的视频| 亚洲av电影不卡..在线观看| 三级国产精品欧美在线观看 | a级毛片在线看网站| 看免费av毛片| 99re在线观看精品视频| 日本一二三区视频观看| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美一区二区三区在线观看| 国产真实乱freesex| 国内精品一区二区在线观看| 男女床上黄色一级片免费看| 女人被狂操c到高潮| 亚洲熟女毛片儿| 久久中文字幕一级| 最新美女视频免费是黄的| 欧美+亚洲+日韩+国产| 老司机午夜福利在线观看视频| 日韩有码中文字幕| 亚洲av中文字字幕乱码综合| 午夜福利高清视频| 午夜亚洲福利在线播放| 日韩欧美在线二视频| 亚洲成人久久性| 91成年电影在线观看| 蜜桃久久精品国产亚洲av| 午夜激情av网站| АⅤ资源中文在线天堂| 亚洲欧美精品综合久久99| 国产成人av激情在线播放| 国产精品一及| 国产在线精品亚洲第一网站| 美女高潮喷水抽搐中文字幕| 免费观看精品视频网站| 国产精品日韩av在线免费观看| 国产激情欧美一区二区| 欧美黑人欧美精品刺激| 国产99久久九九免费精品| 国产男靠女视频免费网站| 18美女黄网站色大片免费观看| 十八禁网站免费在线| 久久性视频一级片| 国产精品电影一区二区三区| 美女免费视频网站| 十八禁人妻一区二区| 日日摸夜夜添夜夜添小说| 一a级毛片在线观看| 成人av在线播放网站| 亚洲专区字幕在线| 他把我摸到了高潮在线观看| 很黄的视频免费| 每晚都被弄得嗷嗷叫到高潮| cao死你这个sao货| 人人妻人人看人人澡| 国产精品日韩av在线免费观看| 免费一级毛片在线播放高清视频| 亚洲精品美女久久av网站| 精品一区二区三区av网在线观看| 国产午夜精品久久久久久| 黄色毛片三级朝国网站| 国产一区二区三区在线臀色熟女| 无限看片的www在线观看| 国产黄色小视频在线观看| 日本熟妇午夜| 亚洲国产精品成人综合色| 日韩精品中文字幕看吧| 免费看十八禁软件| 亚洲精品中文字幕在线视频| 国产麻豆成人av免费视频| 久久精品人妻少妇| 免费在线观看亚洲国产| 国产在线观看jvid| 成人18禁在线播放| 五月玫瑰六月丁香| 看片在线看免费视频| 一个人免费在线观看的高清视频| 国产伦一二天堂av在线观看| 久久精品国产亚洲av高清一级| 色老头精品视频在线观看| 一个人观看的视频www高清免费观看 | 免费看日本二区| 国产真人三级小视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 日韩 欧美 亚洲 中文字幕| 中出人妻视频一区二区| 两个人免费观看高清视频| 国内久久婷婷六月综合欲色啪| 18禁国产床啪视频网站| 亚洲 欧美一区二区三区| 国产一区二区激情短视频| 日韩欧美在线乱码| 精品免费久久久久久久清纯| 国产精品99久久99久久久不卡| 国模一区二区三区四区视频 | 中文字幕人妻丝袜一区二区| 波多野结衣巨乳人妻| 99国产综合亚洲精品| 99国产精品一区二区三区| 免费人成视频x8x8入口观看| 亚洲国产高清在线一区二区三| 一本一本综合久久| 久久精品人妻少妇| 亚洲第一欧美日韩一区二区三区| 国产蜜桃级精品一区二区三区| 国产精品久久久久久精品电影| 高清在线国产一区| 99在线人妻在线中文字幕| 十八禁人妻一区二区| 变态另类成人亚洲欧美熟女| 久久午夜亚洲精品久久| 久久欧美精品欧美久久欧美| 精品午夜福利视频在线观看一区| 国产av一区二区精品久久| 男女那种视频在线观看| 九九热线精品视视频播放| 久久精品亚洲精品国产色婷小说| 免费人成视频x8x8入口观看| 精品国产超薄肉色丝袜足j| 欧美日韩福利视频一区二区| 久久久国产欧美日韩av| 国产激情久久老熟女| 欧美黑人精品巨大| 人妻丰满熟妇av一区二区三区| 国产亚洲欧美在线一区二区| 欧美 亚洲 国产 日韩一| cao死你这个sao货| 妹子高潮喷水视频| 亚洲成人久久爱视频| 男女午夜视频在线观看| 午夜成年电影在线免费观看| 亚洲国产欧洲综合997久久,| 在线观看舔阴道视频| 岛国在线免费视频观看| 一边摸一边做爽爽视频免费| 在线观看日韩欧美| 日韩成人在线观看一区二区三区| 国产片内射在线| 亚洲熟妇中文字幕五十中出| 女人爽到高潮嗷嗷叫在线视频| 欧美zozozo另类| 国产精品久久久久久久电影 | 俺也久久电影网| 12—13女人毛片做爰片一| av片东京热男人的天堂| 黄色女人牲交| 国产av麻豆久久久久久久| 精品国产乱码久久久久久男人| 老司机在亚洲福利影院| 久久久国产成人免费| 午夜福利在线在线| 午夜成年电影在线免费观看| 亚洲精品美女久久av网站| 亚洲av成人一区二区三| 国产成人精品无人区| 麻豆成人av在线观看| 欧美日韩一级在线毛片| 97人妻精品一区二区三区麻豆| 丰满人妻熟妇乱又伦精品不卡| 天堂√8在线中文| 一本综合久久免费| 又大又爽又粗| 成人国产一区最新在线观看| 亚洲avbb在线观看| 岛国视频午夜一区免费看| 亚洲av片天天在线观看| 日本 欧美在线| 12—13女人毛片做爰片一| 欧美3d第一页| 1024香蕉在线观看| 国产亚洲精品综合一区在线观看 | 黄色 视频免费看| 日本撒尿小便嘘嘘汇集6| 国内揄拍国产精品人妻在线| 免费观看精品视频网站| 国产黄色小视频在线观看| 熟女电影av网| 99久久精品国产亚洲精品| 国产高清激情床上av| 欧美最黄视频在线播放免费| 搡老妇女老女人老熟妇| 久久久久久免费高清国产稀缺| 99久久精品国产亚洲精品| 两个人免费观看高清视频| 亚洲精品av麻豆狂野| 黄频高清免费视频| 亚洲成人免费电影在线观看| xxxwww97欧美| 精品国产美女av久久久久小说| 真人做人爱边吃奶动态| 男女视频在线观看网站免费 | 狂野欧美激情性xxxx| 黄片小视频在线播放| 亚洲国产欧美网| 色播亚洲综合网| 日韩欧美 国产精品| 精品欧美一区二区三区在线| 深夜精品福利| 好看av亚洲va欧美ⅴa在| 丝袜人妻中文字幕| 日韩欧美国产一区二区入口| √禁漫天堂资源中文www| 女生性感内裤真人,穿戴方法视频| 亚洲专区字幕在线| 男人的好看免费观看在线视频 | 高清在线国产一区| 久久久久久大精品| 777久久人妻少妇嫩草av网站| 国产亚洲av嫩草精品影院| 身体一侧抽搐| 无人区码免费观看不卡| 午夜免费成人在线视频| 少妇的丰满在线观看| 国产在线观看jvid| 国产69精品久久久久777片 | 欧美黑人精品巨大| 国产片内射在线| 美女黄网站色视频| 中文字幕久久专区| 九九热线精品视视频播放| 黑人操中国人逼视频| 性色av乱码一区二区三区2| 国内精品一区二区在线观看| 丝袜人妻中文字幕| 国产99白浆流出| 欧美人与性动交α欧美精品济南到| 日韩免费av在线播放| 国产日本99.免费观看| 美女大奶头视频| 校园春色视频在线观看| 久久亚洲精品不卡| 久久国产精品人妻蜜桃| 亚洲国产欧美人成| 国产精品1区2区在线观看.| 18禁黄网站禁片午夜丰满| 国产免费av片在线观看野外av| 亚洲av片天天在线观看| 日韩高清综合在线| 色在线成人网| 精品国产乱子伦一区二区三区| 亚洲av熟女| 久久久久久九九精品二区国产 | 久久久水蜜桃国产精品网| 曰老女人黄片| 美女 人体艺术 gogo| 香蕉av资源在线| 国产伦在线观看视频一区| 久久九九热精品免费| 国产真人三级小视频在线观看| 日日爽夜夜爽网站| 中亚洲国语对白在线视频| av中文乱码字幕在线| 男女那种视频在线观看| 久久精品国产亚洲av高清一级| 国产爱豆传媒在线观看 | 国产精品日韩av在线免费观看| 亚洲欧美精品综合一区二区三区| 久久久久性生活片| 美女高潮喷水抽搐中文字幕| 欧美日韩瑟瑟在线播放| 黄色片一级片一级黄色片| 日韩欧美一区二区三区在线观看| 国产熟女xx| tocl精华| 亚洲国产看品久久| 亚洲国产欧美网| 给我免费播放毛片高清在线观看| 国产精品久久久av美女十八| 日日夜夜操网爽| 一边摸一边抽搐一进一小说| 日韩欧美精品v在线| 啦啦啦观看免费观看视频高清| 一个人观看的视频www高清免费观看 | 欧美成人午夜精品| 欧美3d第一页| 19禁男女啪啪无遮挡网站| 一级黄色大片毛片| 久久精品国产99精品国产亚洲性色| 黄色视频不卡| 国产亚洲av高清不卡| 中出人妻视频一区二区| 欧美一区二区国产精品久久精品 | 成人午夜高清在线视频| 国产精品综合久久久久久久免费| 国产精品 欧美亚洲| 香蕉丝袜av| 久久久久久九九精品二区国产 | 精品久久久久久久毛片微露脸| 久久久久久免费高清国产稀缺| 国产精品久久电影中文字幕| 丝袜美腿诱惑在线| ponron亚洲|