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

    基于Excel VBA的漁具選擇性分析SELECT模型實(shí)現(xiàn)

    2014-05-17 01:34:48金宇鋒
    實(shí)驗(yàn)室研究與探索 2014年3期
    關(guān)鍵詞:拖網(wǎng)漁獲漁具

    金宇鋒, 張 健

    (上海海洋大學(xué) 大洋漁業(yè)資源可持續(xù)開(kāi)發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室;

    國(guó)家遠(yuǎn)洋工程技術(shù)研究中心;農(nóng)業(yè)部大洋漁業(yè)資源環(huán)境科學(xué)觀測(cè)實(shí)驗(yàn)站,上海 201306)

    0 引言

    漁具的選擇性是指漁具在作業(yè)過(guò)程中只能捕獲特定種類(lèi)或者具有特定生物特征捕撈對(duì)象的特性,它是漁具設(shè)計(jì)和漁業(yè)管理的重要依據(jù)。在拖網(wǎng)漁具選擇性研究中,早期對(duì)于比較試驗(yàn)法實(shí)驗(yàn)結(jié)果的處理方法因缺乏嚴(yán)謹(jǐn)統(tǒng)計(jì)模型而被質(zhì)疑[1-2],而SELECT模型的建立和提出解決了這一難題,并被廣泛運(yùn)用于各類(lèi)漁具的選擇性研究中[3-4]。SELECT模型作為廣義線性模型,在具體的漁具選擇性分析中有多種實(shí)現(xiàn)手段,除了使用SAS、Splus等統(tǒng)計(jì)分析軟件的線性函數(shù)[5-6]、或者使用Fortran或者R語(yǔ)言編寫(xiě)的處理程序包以外[7-8],還有使用專(zhuān)門(mén)編寫(xiě)的軟件,例如Constat的拖網(wǎng)漁具選擇性分析軟件CC2000[9-12]以及Herrmann等開(kāi)發(fā)的分析軟件SELNET[13-14]。此外,基于表格處理軟件Excel的“規(guī)劃求解”(SOLVER Add-in)功能實(shí)現(xiàn)模型解析的方法自運(yùn)用于拖網(wǎng)選擇性分析以來(lái)也得到了廣泛的運(yùn)用[15-18],在科研和教學(xué)實(shí)踐中發(fā)現(xiàn),要充分利用這一方法實(shí)現(xiàn)拖網(wǎng)等漁具的選擇性分析不僅需要對(duì)SELECT模型的充分認(rèn)識(shí),還需要具備一定的Excel函數(shù)知識(shí)。本文根據(jù)對(duì)SELECT模型的解析,結(jié)合Excel內(nèi)置的VBA(Visual Basic for Application)程序語(yǔ)言,開(kāi)展模型在拖網(wǎng)等過(guò)濾性漁具選擇性分析中的自動(dòng)實(shí)現(xiàn)的研究和探討,為今后開(kāi)展?jié)O具選擇性研究或教學(xué)提供參考。

    1 模型簡(jiǎn)介

    1.1 SELECT 模型原理

    SELECT 模型[3,19]從捕撈對(duì)象的分布特點(diǎn)出發(fā),認(rèn)為在漁具選擇性研究中,特定漁具i(研究對(duì)象)所捕獲的捕撈對(duì)象的數(shù)量Ci服從一個(gè)多項(xiàng):

    Ci~Multi(∑C;Φ1,…,Φi,…,ΦI)

    式中:Φi為單元概率,表示一尾被捕獲的個(gè)體被i漁具捕獲的的期望概率,即Pi為i漁具的相對(duì)作業(yè)強(qiáng)度;Sij為i漁具的選擇率;θ為模型參數(shù),即包含選擇性參數(shù),也包括相對(duì)作業(yè)強(qiáng)度。

    1.2 SELECT模型在拖網(wǎng)選擇性分析中的解析

    在拖網(wǎng)選擇性試驗(yàn)中,若使用比較試驗(yàn)法,對(duì)比對(duì)照網(wǎng)和試驗(yàn)網(wǎng)的漁獲進(jìn)行分析,SELECT模型轉(zhuǎn)化為二項(xiàng)分布,即

    式中:C1、C2分別為試驗(yàn)網(wǎng)和對(duì)照網(wǎng)捕獲的漁獲數(shù)量;單元概率為試驗(yàn)網(wǎng)的相對(duì)作業(yè)強(qiáng)度。

    若使用套網(wǎng)法開(kāi)展選擇性試驗(yàn),則C1、C2分別為網(wǎng)囊和套網(wǎng)中的漁獲數(shù)量,且Φ=S。

    1.3 SELECT模型在拖網(wǎng)選擇性分析中實(shí)現(xiàn)的方法

    在拖網(wǎng)選擇性分析中,通常使用極大似然估計(jì)法估算模型參數(shù),似然函數(shù)為

    對(duì)于選擇性曲線,除了Logistic曲線以外,考慮到特定情況下曲線的非對(duì)稱(chēng)性,還可選擇Richards曲線(通過(guò)假設(shè)檢驗(yàn)可驗(yàn)證是否可以簡(jiǎn)化為L(zhǎng)ogistic曲線)。在對(duì)比試驗(yàn)法下,既可將相對(duì)作業(yè)強(qiáng)度假設(shè)為P=1/2(稱(chēng)為等分模型)、也可將相對(duì)作業(yè)強(qiáng)度作為參數(shù)進(jìn)行估算以此比較捕撈對(duì)象對(duì)實(shí)驗(yàn)漁具和對(duì)照漁具的行為差異(稱(chēng)為估計(jì)分隔模型)[19],當(dāng)實(shí)驗(yàn)網(wǎng)和對(duì)照網(wǎng)捕撈努力量(例如作業(yè)時(shí)間或者漁具數(shù)量)不同時(shí),也可直接給定P(定義為固定分隔模型)[20]。那么,在拖網(wǎng)選擇性分析中,SELECT模型的拖網(wǎng)選擇性分析中的實(shí)現(xiàn)可歸納為8種模型,表1中分別用1~8來(lái)表示,模型簡(jiǎn)化可歸納為7類(lèi),分別為表中的11-17。

    表1 拖網(wǎng)選擇性分析中SELECT模型實(shí)現(xiàn)方法和簡(jiǎn)化方法

    2 拖網(wǎng)選擇性分析的Excel VBA自動(dòng)實(shí)現(xiàn)

    2.1 自動(dòng)實(shí)現(xiàn)流程圖

    在不同的拖網(wǎng)選擇性分析模型和模型簡(jiǎn)化,建立如圖1所示的主程序流程圖。

    2.2 交互界面

    2.2.1 實(shí)驗(yàn)參數(shù)交互界面

    圖1 主程序流程圖

    圖2 實(shí)驗(yàn)參數(shù)錄入和曲線選擇交互界面

    實(shí)驗(yàn)參數(shù)錄入和曲線選擇交互界面如圖2所示。除可設(shè)置網(wǎng)次數(shù)量、選擇試驗(yàn)方法、選擇性曲線以外,為了后續(xù)可以開(kāi)展網(wǎng)次間差異的分析,同時(shí)設(shè)置了試驗(yàn)參數(shù)(既包括定性,也包含定量參數(shù))。通過(guò)不同的設(shè)置方法和參數(shù)選擇,向主程序傳遞參數(shù)。

    2.2.2 選擇性分析參數(shù)選擇交互界面

    選擇性分析中的模型簡(jiǎn)化、初始值設(shè)定、標(biāo)準(zhǔn)差計(jì)算選擇交互界面如圖3所示。其中,嵌套模型(Nested Model)功能是通過(guò)假設(shè)檢驗(yàn)以對(duì)模型進(jìn)行簡(jiǎn)化,包括對(duì)選擇性曲線和對(duì)估算分隔模型的簡(jiǎn)化判斷(具體簡(jiǎn)化的方法見(jiàn)表1)。

    2.3 選擇性分析程序和調(diào)用

    圖3 選擇性分析參數(shù)選擇交互界面

    從功能上程序可分為2類(lèi),其中選擇率、單元概率和似然值計(jì)算、模型選擇和曲線簡(jiǎn)化以及使用“規(guī)劃求解”功能求解模型參數(shù)為主程序,其他部分包括參數(shù)標(biāo)準(zhǔn)差計(jì)算、圖表顯示等為輔助程序。模型實(shí)現(xiàn)的完 整 程 序 可 從http://hyxy.shou.edu.cn/keylab/admin/UploadFiles/2013581059115.rar處下載并查看。

    程序采用了RibbonX代碼的方式,在功能區(qū)制作了調(diào)用控件以實(shí)現(xiàn)程序調(diào)用。而主程序則采用加載宏的方式調(diào)用,當(dāng)然直接程序打開(kāi)亦可使用。

    3 實(shí)例分析

    分別選用文獻(xiàn)Fryer(1991)中套網(wǎng)法拖網(wǎng)選擇性試驗(yàn)數(shù)據(jù)[21](僅以第1種漁具的第2網(wǎng)次數(shù)據(jù))和文獻(xiàn)Millar(1992)中對(duì)比試驗(yàn)拖網(wǎng)選擇性試驗(yàn)數(shù)據(jù)[3]為例,使用上述代碼(結(jié)合相應(yīng)的輔助代碼),選擇不同分析模型并選擇模型簡(jiǎn)化,進(jìn)行數(shù)據(jù)擬合,所得結(jié)果如表2所示。

    對(duì)比原文獻(xiàn)中的實(shí)驗(yàn)數(shù)據(jù)擬合結(jié)果發(fā)現(xiàn),在相同的模型下,擬合結(jié)果是完全一致的。不僅如此,在程序中,我們內(nèi)置了不同模型的比較和模型簡(jiǎn)化的方法,根據(jù)優(yōu)劣性指標(biāo)選擇合適的模型,例如在表2中,根據(jù)套網(wǎng)法的擬合和簡(jiǎn)化結(jié)果,可使用 Logistic曲線代替Richards曲線以?xún)?yōu)化模型(P>0.05);而在對(duì)比試驗(yàn)法中,使用等分模型的假設(shè)不能被接受(P<0.01),但可接受Logistic代替Richards曲線的假設(shè)(P>0.05)。

    4 討論

    目前SELECT模型在拖網(wǎng)漁具選擇性分析可通過(guò)多種途徑實(shí)現(xiàn),但是無(wú)論是使用大型統(tǒng)計(jì)軟件的相關(guān)模塊、還是使用專(zhuān)門(mén)的選擇性分析軟件,除了需要具備一定的軟件(或語(yǔ)言)基礎(chǔ)外,商業(yè)軟件本身價(jià)格不菲。通過(guò)使用VBA實(shí)現(xiàn)SELECT模型在漁具選擇性分析雖然受到Excel平臺(tái)和VBA語(yǔ)言自身?xiàng)l件限制,例如無(wú)法定義單元格的數(shù)據(jù)類(lèi)型、在浮點(diǎn)計(jì)算、矩陣計(jì)算等方面的先天不足、代碼執(zhí)行效率低和計(jì)算時(shí)間較長(zhǎng)等,但是在漁具選擇性分析中,因?yàn)橛?jì)算量不大,這些缺點(diǎn)并不突出。此外,使用VBA實(shí)現(xiàn)的方法與其它幾種實(shí)現(xiàn)方法相比,仍存在一些優(yōu)勢(shì):

    表2 不同試驗(yàn)方法下模型擬合和簡(jiǎn)化結(jié)果

    (1)繼承Excel平臺(tái)和VBA語(yǔ)言的優(yōu)點(diǎn),例如無(wú)需運(yùn)行庫(kù),無(wú)需配置環(huán)境參數(shù),僅依靠Excel平臺(tái)而忽略操作系統(tǒng)和數(shù)據(jù)格式;可調(diào)用Excel內(nèi)置的函數(shù);Excel自帶的打印、文件處理、格式化或者文本編輯功能,且VBA支持調(diào)用Windows底層API;VBA與VB(VisalBasci)代碼格式一致,可調(diào)用一切基于VB的函數(shù)等;Excel的加載項(xiàng)“規(guī)劃求解”功能強(qiáng)大,除可設(shè)置迭代算法以外,還可設(shè)置精度等條件;此外,可對(duì)Excel自帶的圖表格式進(jìn)行修改以滿(mǎn)足不同的需求;Excel自帶的宏錄制功能,不僅方便處理一些不熟悉的實(shí)現(xiàn)步驟,也便于對(duì)程序代碼的理解;這些都為今后的繼續(xù)開(kāi)展?jié)O具選擇性分析的VBA自動(dòng)實(shí)現(xiàn)提供了便利。

    (2)Excel軟件普及度較高,通用性強(qiáng),而且其中“規(guī)劃求解”功能也屬于免費(fèi)加載項(xiàng);此外,多數(shù)情況下前期的漁獲數(shù)據(jù)處理會(huì)使用Excel軟件,因此前后期數(shù)據(jù)銜接較好,便于對(duì)數(shù)據(jù)的統(tǒng)計(jì)分析和分類(lèi)保存,數(shù)據(jù)格式比較統(tǒng)一;

    (3)因?yàn)樾枰褂谩耙?guī)劃求解”功能求模型參數(shù),因此在SELECT模型的實(shí)現(xiàn)過(guò)程中,使用了單元格公式編輯(.formula或者.formulaR1C1語(yǔ)法)的方法,而不是直接計(jì)算單元格的值,這一方法不僅有助于了解模型實(shí)現(xiàn)的步驟進(jìn)而掌握模型,而且有助于利用Excel自帶的查錯(cuò)功能迅速定位代碼中的邏輯錯(cuò)誤或者計(jì)算過(guò)程中的數(shù)值錯(cuò)誤,一方面可方便對(duì)模型的修改或者完善補(bǔ)充;另一方面也有助于在教學(xué)中的理解和應(yīng)用。

    表1基本涵蓋了拖網(wǎng)等過(guò)濾性漁具選擇性分析中對(duì)套網(wǎng)法和對(duì)比試驗(yàn)法的試驗(yàn)數(shù)據(jù)分析的基本模型,但在此次SELECT模型的實(shí)現(xiàn)過(guò)程中,仍對(duì)部分模型和試驗(yàn)條件進(jìn)行了限制。例如,雖然過(guò)濾性漁具的選擇性曲線模型有多種,但是從目前多數(shù)研究來(lái)看,基本使用Logistic和Richards曲線,因此在此次VBA實(shí)現(xiàn)中僅考慮了這2個(gè)曲線模型;考慮到分析結(jié)果與后續(xù)使用混合影響模型進(jìn)行網(wǎng)次間差異分析的相互銜接,限定了定性的試驗(yàn)漁具種類(lèi)不超過(guò)3種,定量的試驗(yàn)參數(shù)不超過(guò)4種(實(shí)際中這些限定條件能滿(mǎn)足絕大部分分析要求)。除了對(duì)模型和試驗(yàn)條件限制以外也應(yīng)該看到,模型的實(shí)現(xiàn)的主要流程也可被應(yīng)用于其它漁具、其它試驗(yàn)方法或者其它選擇性模型,張健等就利用了類(lèi)似的實(shí)現(xiàn)手段對(duì)分隔裝置的分隔性能進(jìn)行了分析[22],因此利用Excel完成實(shí)現(xiàn)模型這一思路具有較為廣泛的擴(kuò)展性。

    在模型的擬合中也應(yīng)注意,“規(guī)劃求解”加載項(xiàng)功能雖然強(qiáng)大,但其對(duì)初始值的選取比較敏感,當(dāng)始值選取不合理的時(shí)候,迭代過(guò)程會(huì)出現(xiàn)不收斂的情況,進(jìn)而無(wú)法得出合理的試驗(yàn)結(jié)果,特別是在使用Richards曲線作為選擇性曲線(存在指數(shù)計(jì)算)的時(shí)候。因此在進(jìn)行自動(dòng)化求解選擇性參數(shù)之前,須根據(jù)經(jīng)驗(yàn)或者制作漁獲分布圖對(duì)漁獲數(shù)據(jù)進(jìn)行簡(jiǎn)單分析判斷出選擇性參數(shù)的大致范圍,并可在選擇性分析參數(shù)選擇交互界面(見(jiàn)圖3)中預(yù)先設(shè)定。

    過(guò)濾性漁具作業(yè)效率較高,因此會(huì)出現(xiàn)對(duì)捕獲的漁獲進(jìn)行取樣的情況。此次模型實(shí)現(xiàn)中,根據(jù)有關(guān)的理論[23]也將這一因素考慮到模型之中,因此模型的適用性得到了進(jìn)一步的加強(qiáng)。

    [1]CADIGAN N G,MILLAR R B.Reliability of selection curves obtained from trouser trawl or alternate hual experiments[J].Canadian Journal of Fisheries and Aquatic Sciences,1992,49:1625-1632.

    [2]梁振林,葛長(zhǎng)字,劉英光.國(guó)外漁具選擇性研究進(jìn)展[J].青島海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,31(6):835-841.

    LIANG Zhen-lin,GE Chang-zi,LIU Ying-guang.Review of study on grid separators for beam trawl the gear selectivityabroad[J].Journal of ocean university of china(Natural Science),2001,31(6):835-841.

    [3]MILLAR R B.Estimation the size-selectivity of fishing gear by conditioning on the total catch[J].Journal of the American Statistical Association,1992,87:962-968.

    [4]MILLAR R B,F(xiàn)RYER R J.Estimating the size-selection curves of towed gears,traps,nets and hooks[J].Reviews in Fish Biology and Fisheries,1999,9(1):89-116.

    [5]MILLAR R B.Analysis of trawl selectivity studies(addendum):implementation in SAS[J].Fisheries Research,1993,17(3-4):373-377.

    [6]MILLAR R B,Holst R.Estimation of gillnet and hook selectivity using log-linear models[J].ICES Journal of Marine Science,1997,54:471-477.

    [7]MILLAR R B,CADIGAN N G.A FORTRAN program for fitting selectivity curves to trouser trawl data[M].St.John’s:Dept.of Fisheries and Oceans, Newfoundland Region, Science Branch,1990.

    [8]MILLAR R B.Selectivity software and documentation[EB/OL].The University of Auckland,New Zealand,2010 [2013-02-03];http://www.stat.auckland.ac.nz/~ millar/selectware/code.html.

    [9]AYDIN C,TOSUNOGLU Z.Selectivity of square and hexagonal mesh codendsforthe deep waterroseshrimp,Parapenaeus longirostris(Lucas,1846)(Decapoda,Penaeidae)in the Aegean Sea[J].Crustaceana,2009,82(1):89-98.

    [10]AYDIN C,TSUNGLU Z,TOKAC A.Sorting grid trials to improve size selectivity of red mullet(Mullus barbatus)and annular sea bream(Diplodus annularis)in Turkish bottom trawl fishery[J].Journal of Applied Ichthyology,2008,24(3):306-310.

    [11]CONSTAT,CC2000Selectivityprogram[CP], ConStatDK,Hjorring,Denmark,2000.

    [12]FRANDSEN R P,MADSEN N,KRAG L A.Selectivity and escapement behaviour of five commercial fishery species in standard square-and diamond-mesh codends[J].ICES Journal of Marine Science:Journal du Conseil,2010,67(8):1721-1731.

    [13]HERRMANN B,MADSEN N,SISTIAGA M,et al.Modelling escapement during the fishing process as a dual sequence-Introducing SELNE.Poster session presented at ICES Working Group on Fishing Technology and Fish Behaviour(WGFTFB).Ancona,Italy:ICES;2009.

    [14]SISTIAGA M,HERRMANN B,GRIMALDO E,et al.Assessment of dual selection in grid based selectivity systems[J].Fisheries Research,2010,105(3):187-199.

    [15]TOKAI T.Maximum likelihood parameter estimates of a mesh selectivity logistic model through SOLVER on MS-Excel[J].Bulletin of The Japanese Society of Fisheries Oceanography,1997,61(3):288-298.

    [16]ZUUR G,F(xiàn)RYER R J,F(xiàn)ERRO R S T,et al.Modelling the size selectivities of a trawl codend and an associated square mesh panel[J].ICES Journal of Marine Science,2001,58(3):657-671.

    [17]WILEMAN D A,F(xiàn)ERRO R S T,F(xiàn)ONTEYNE R,et al.Manual of Methods of Measuring the Selectivity of Towed Fishing Gears[R].Copenhagen Denmark:ICES,1996.No.215.

    [18]孫滿(mǎn)昌,張 健,許柳雄,等.漁具漁法選擇性[M].北京:中國(guó)農(nóng)業(yè)出版社,2004.

    [19]MILLAR R B,WALSH S J.Analysis of trawl selectivity studies with an application to trouser trawl[J].Fisheries Research,1992,13(3):205-220.

    [20]XU X,MILLAR R B.Estimation of trap selectivity for male snow crab(Chionoecetes opilio)using the SELECT modeling approach with unequal sampling effort[J].Canadian Journal of Fisheries and Aquatic Sciences,1993,50(11):2485-2490.

    [21]FRYER R J.A model of between-haul variation in selectivity[J].ICES J Mar Sci,1991,48(1):281-290.

    [22]張 健,石建高,張 鵬,等.分隔式桁拖網(wǎng)分隔裝置的分隔效率模型研究[J].水產(chǎn)學(xué)報(bào),2009,33(4):697-703.

    ZHANG Jian, SHI Jian-gao, ZHANG Peng, et al.Model for separation effeciency assessment of separator devices rigged in multilevel beam trawls[J].Journal of Fsheries of China,2009,33(4):697-703.

    [23]MILLAR R B.Sampling from trawl gears used in size selectivity experiments[J].ICES Journal of Marine Science,1994,51(3):293-298.

    猜你喜歡
    拖網(wǎng)漁獲漁具
    拖網(wǎng)的人(外一首)
    拖網(wǎng)的人(外一首)
    數(shù)據(jù)
    開(kāi)漁
    自制漁具
    小讀者之友(2020年8期)2020-09-01 07:08:23
    首批國(guó)家級(jí)海洋捕撈漁獲物定點(diǎn)上岸漁港廣西唯一入選名額落戶(hù)南漁港
    東海北部桁桿拖蝦漁獲組成分析
    南海北部深海區(qū)燈光罩網(wǎng)漁獲物組成及漁獲率的時(shí)空分布
    河北省框架拖網(wǎng)和桁桿拖網(wǎng)捕撈現(xiàn)狀調(diào)查
    漁獲
    寶藏(2017年4期)2017-05-17 03:34:52
    成人亚洲精品一区在线观看| 亚洲五月天丁香| 又紧又爽又黄一区二区| 国产极品粉嫩免费观看在线| 欧美激情极品国产一区二区三区| 欧美不卡视频在线免费观看 | 搡老岳熟女国产| 久久久久国产精品人妻aⅴ院| 国产欧美日韩综合在线一区二区| 久久久久久久午夜电影 | 99热只有精品国产| 亚洲精品国产区一区二| 黑人欧美特级aaaaaa片| 久久亚洲真实| 亚洲国产毛片av蜜桃av| 日韩精品中文字幕看吧| 国产精品 欧美亚洲| 久久精品亚洲av国产电影网| 亚洲视频免费观看视频| 国产视频一区二区在线看| 午夜久久久在线观看| 男女下面插进去视频免费观看| 国产又爽黄色视频| 女性被躁到高潮视频| 午夜成年电影在线免费观看| 久久精品成人免费网站| 久久久水蜜桃国产精品网| 成人国产一区最新在线观看| 日日夜夜操网爽| 97超级碰碰碰精品色视频在线观看| 狠狠狠狠99中文字幕| 99re在线观看精品视频| 侵犯人妻中文字幕一二三四区| 91大片在线观看| 青草久久国产| 免费高清在线观看日韩| 国产真人三级小视频在线观看| 在线观看免费午夜福利视频| 成年人免费黄色播放视频| 久久久久亚洲av毛片大全| 一边摸一边抽搐一进一出视频| 人人澡人人妻人| 久99久视频精品免费| 亚洲性夜色夜夜综合| 一进一出好大好爽视频| 乱人伦中国视频| 日韩精品中文字幕看吧| 国产精品乱码一区二三区的特点 | 欧美黑人欧美精品刺激| 国产精品免费视频内射| 久久影院123| 亚洲人成电影免费在线| 俄罗斯特黄特色一大片| 又紧又爽又黄一区二区| 国产精品 欧美亚洲| 精品国产乱码久久久久久男人| av在线播放免费不卡| 在线观看免费视频日本深夜| 久久久久久久久中文| 成人影院久久| 日本vs欧美在线观看视频| 久久久久久大精品| 午夜免费观看网址| 午夜影院日韩av| 超碰97精品在线观看| 757午夜福利合集在线观看| 少妇的丰满在线观看| 国产精品 国内视频| 日韩 欧美 亚洲 中文字幕| 一边摸一边抽搐一进一出视频| 精品久久久久久久久久免费视频 | 99热国产这里只有精品6| 午夜福利在线免费观看网站| 男女做爰动态图高潮gif福利片 | 日韩大码丰满熟妇| 国产黄a三级三级三级人| 80岁老熟妇乱子伦牲交| 精品国产亚洲在线| 老司机深夜福利视频在线观看| 91精品国产国语对白视频| 高潮久久久久久久久久久不卡| 欧美乱妇无乱码| 国产伦人伦偷精品视频| videosex国产| 欧美黄色片欧美黄色片| 啦啦啦免费观看视频1| 欧美日韩黄片免| 女性被躁到高潮视频| 一区二区三区国产精品乱码| 国产高清国产精品国产三级| 十分钟在线观看高清视频www| 母亲3免费完整高清在线观看| 国产一区二区激情短视频| 亚洲中文av在线| 国产成年人精品一区二区 | 国产精品一区二区三区四区久久 | 午夜福利,免费看| 亚洲成人精品中文字幕电影 | 黄色 视频免费看| 色婷婷久久久亚洲欧美| 黑人操中国人逼视频| 亚洲第一欧美日韩一区二区三区| 69精品国产乱码久久久| 亚洲av第一区精品v没综合| 久久天堂一区二区三区四区| 少妇粗大呻吟视频| 超碰成人久久| 又黄又粗又硬又大视频| 亚洲 欧美一区二区三区| av网站在线播放免费| 女人精品久久久久毛片| 久久久国产成人免费| 日本免费a在线| 日本五十路高清| 欧美+亚洲+日韩+国产| 成人特级黄色片久久久久久久| www日本在线高清视频| 国产精品香港三级国产av潘金莲| 久久久久久久久免费视频了| 日本vs欧美在线观看视频| 真人一进一出gif抽搐免费| 视频在线观看一区二区三区| 亚洲伊人色综图| 亚洲少妇的诱惑av| 韩国精品一区二区三区| 免费在线观看影片大全网站| 咕卡用的链子| 久久精品国产99精品国产亚洲性色 | 午夜亚洲福利在线播放| 老汉色av国产亚洲站长工具| 亚洲av美国av| 国产av一区二区精品久久| 天天躁夜夜躁狠狠躁躁| 久久狼人影院| 亚洲成人免费av在线播放| 久久婷婷成人综合色麻豆| 夜夜看夜夜爽夜夜摸 | 亚洲av成人不卡在线观看播放网| 变态另类成人亚洲欧美熟女 | 日韩大尺度精品在线看网址 | av网站在线播放免费| 色婷婷av一区二区三区视频| 91老司机精品| 一级毛片女人18水好多| 欧美乱码精品一区二区三区| 国产精品国产高清国产av| 热re99久久精品国产66热6| 亚洲精品在线美女| 日本免费一区二区三区高清不卡 | 免费一级毛片在线播放高清视频 | 精品国产一区二区久久| 日本一区二区免费在线视频| 欧美日韩一级在线毛片| 午夜福利在线观看吧| 亚洲国产欧美一区二区综合| 免费人成视频x8x8入口观看| 亚洲精华国产精华精| 欧美日韩视频精品一区| 在线观看免费日韩欧美大片| 91av网站免费观看| 亚洲精品粉嫩美女一区| 在线永久观看黄色视频| 成人三级黄色视频| 丝袜在线中文字幕| 美国免费a级毛片| 久久午夜综合久久蜜桃| 久久久精品国产亚洲av高清涩受| 制服诱惑二区| 国产成人精品久久二区二区91| 女生性感内裤真人,穿戴方法视频| 精品一区二区三卡| 亚洲精品一区av在线观看| 999久久久精品免费观看国产| a级毛片黄视频| 欧美 亚洲 国产 日韩一| 日日爽夜夜爽网站| 亚洲精品国产一区二区精华液| 热99re8久久精品国产| 久久精品成人免费网站| 少妇裸体淫交视频免费看高清 | 免费搜索国产男女视频| 香蕉丝袜av| av在线天堂中文字幕 | 日韩欧美一区二区三区在线观看| av免费在线观看网站| 国产成+人综合+亚洲专区| 亚洲精品久久午夜乱码| 亚洲在线自拍视频| 97超级碰碰碰精品色视频在线观看| 色精品久久人妻99蜜桃| а√天堂www在线а√下载| 久久亚洲真实| 日韩欧美免费精品| 久久久久久久午夜电影 | aaaaa片日本免费| 一进一出抽搐动态| 国产精品国产高清国产av| 亚洲色图 男人天堂 中文字幕| 香蕉国产在线看| 女人精品久久久久毛片| 亚洲avbb在线观看| 高清在线国产一区| 免费久久久久久久精品成人欧美视频| 久久精品国产99精品国产亚洲性色 | 老鸭窝网址在线观看| 这个男人来自地球电影免费观看| 99香蕉大伊视频| 99久久精品国产亚洲精品| 国内久久婷婷六月综合欲色啪| 色尼玛亚洲综合影院| 久久青草综合色| 亚洲欧美日韩另类电影网站| 午夜免费鲁丝| 国产片内射在线| 久久久精品国产亚洲av高清涩受| 性少妇av在线| 欧美日韩亚洲国产一区二区在线观看| 国产成人啪精品午夜网站| 欧美成人免费av一区二区三区| 精品欧美一区二区三区在线| 在线播放国产精品三级| 制服人妻中文乱码| 欧美久久黑人一区二区| 两个人看的免费小视频| 亚洲av成人一区二区三| 九色亚洲精品在线播放| 国产精品免费视频内射| 不卡一级毛片| 亚洲成人久久性| 色在线成人网| 五月开心婷婷网| 老鸭窝网址在线观看| 这个男人来自地球电影免费观看| 亚洲一区二区三区欧美精品| 桃红色精品国产亚洲av| e午夜精品久久久久久久| 9色porny在线观看| 成人永久免费在线观看视频| 不卡av一区二区三区| 色播在线永久视频| 色播在线永久视频| 黄色成人免费大全| 欧美日韩黄片免| 国产不卡一卡二| 亚洲,欧美精品.| 18禁国产床啪视频网站| 国产不卡一卡二| 久久人人精品亚洲av| 亚洲男人天堂网一区| 久久中文看片网| 51午夜福利影视在线观看| 久久性视频一级片| 午夜福利在线观看吧| 97超级碰碰碰精品色视频在线观看| 国产综合懂色| 日韩国内少妇激情av| 国产男靠女视频免费网站| 免费电影在线观看免费观看| 中亚洲国语对白在线视频| 一级av片app| 久久人妻av系列| 一区二区三区激情视频| 天堂动漫精品| 精品欧美国产一区二区三| 免费观看精品视频网站| 精品一区二区三区人妻视频| 日韩欧美一区二区三区在线观看| 国产精品爽爽va在线观看网站| 久久性视频一级片| 女人被狂操c到高潮| 老女人水多毛片| 亚洲无线观看免费| 日本成人三级电影网站| 神马国产精品三级电影在线观看| 久久久成人免费电影| 欧美性猛交黑人性爽| 欧美zozozo另类| 97热精品久久久久久| 在线免费观看不下载黄p国产 | 亚洲av第一区精品v没综合| 欧美成狂野欧美在线观看| 欧美丝袜亚洲另类 | 国产亚洲av嫩草精品影院| 国产成人av教育| 亚洲人成网站在线播| 亚洲三级黄色毛片| 亚洲国产精品合色在线| 国产高清视频在线播放一区| 亚洲片人在线观看| 我要看日韩黄色一级片| 日韩欧美精品免费久久 | 亚洲人成网站在线播放欧美日韩| 亚洲精品一区av在线观看| 在线免费观看的www视频| 丰满人妻一区二区三区视频av| 女生性感内裤真人,穿戴方法视频| 亚洲aⅴ乱码一区二区在线播放| 麻豆国产av国片精品| 午夜福利成人在线免费观看| 看片在线看免费视频| 国产精品久久久久久久电影| 亚洲欧美日韩无卡精品| 听说在线观看完整版免费高清| 国产69精品久久久久777片| 人妻丰满熟妇av一区二区三区| 成人国产综合亚洲| 日韩欧美 国产精品| 又粗又爽又猛毛片免费看| 白带黄色成豆腐渣| 免费观看精品视频网站| 精品久久久久久成人av| 乱码一卡2卡4卡精品| 欧美高清成人免费视频www| 日本一本二区三区精品| 亚洲成av人片免费观看| 女人十人毛片免费观看3o分钟| 亚洲片人在线观看| 国产精品久久久久久精品电影| 能在线免费观看的黄片| 精品不卡国产一区二区三区| 久久久久久久精品吃奶| avwww免费| 国产亚洲av嫩草精品影院| 国产精品亚洲美女久久久| 国产精品一区二区三区四区久久| 亚洲久久久久久中文字幕| 亚洲精品亚洲一区二区| 国产亚洲av嫩草精品影院| 日韩精品中文字幕看吧| 超碰av人人做人人爽久久| 国内精品久久久久精免费| avwww免费| 成年免费大片在线观看| 日韩人妻高清精品专区| 久久午夜亚洲精品久久| 久久久久久久久久成人| 琪琪午夜伦伦电影理论片6080| 久久精品国产亚洲av天美| 亚洲乱码一区二区免费版| 真人一进一出gif抽搐免费| 白带黄色成豆腐渣| 亚洲av熟女| 在线国产一区二区在线| 床上黄色一级片| 久久人妻av系列| 色噜噜av男人的天堂激情| 久久久久性生活片| 国产一级毛片七仙女欲春2| 欧美黑人巨大hd| 91九色精品人成在线观看| 一二三四社区在线视频社区8| 国产久久久一区二区三区| 国产午夜精品久久久久久一区二区三区 | 69av精品久久久久久| 99精品在免费线老司机午夜| 国产主播在线观看一区二区| 极品教师在线视频| 久久人人精品亚洲av| 1024手机看黄色片| 夜夜看夜夜爽夜夜摸| 久久亚洲真实| 黄色一级大片看看| 日本a在线网址| 亚洲午夜理论影院| 国产亚洲欧美在线一区二区| 国产激情偷乱视频一区二区| 哪里可以看免费的av片| 淫妇啪啪啪对白视频| 精品一区二区免费观看| 日韩欧美在线二视频| 97人妻精品一区二区三区麻豆| 一进一出抽搐动态| 一进一出抽搐gif免费好疼| 99riav亚洲国产免费| 成年版毛片免费区| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产日韩欧美精品在线观看| 久久久久久大精品| 午夜视频国产福利| 女生性感内裤真人,穿戴方法视频| 男女那种视频在线观看| 美女黄网站色视频| 中文在线观看免费www的网站| 欧美bdsm另类| 人妻久久中文字幕网| 久久人人精品亚洲av| aaaaa片日本免费| 午夜福利高清视频| 国产人妻一区二区三区在| 精品一区二区三区视频在线观看免费| 十八禁人妻一区二区| 一个人看的www免费观看视频| 搡老妇女老女人老熟妇| 老司机深夜福利视频在线观看| 亚洲狠狠婷婷综合久久图片| 小蜜桃在线观看免费完整版高清| 久久久久久国产a免费观看| 日韩大尺度精品在线看网址| 亚洲av免费在线观看| 99久久成人亚洲精品观看| 中文字幕人成人乱码亚洲影| 欧美成人免费av一区二区三区| 午夜影院日韩av| 亚洲欧美激情综合另类| 亚洲真实伦在线观看| 午夜激情欧美在线| 波多野结衣高清无吗| 国产成人aa在线观看| 人妻夜夜爽99麻豆av| 99久久精品国产亚洲精品| 美女 人体艺术 gogo| 他把我摸到了高潮在线观看| 九色成人免费人妻av| 极品教师在线免费播放| a级毛片a级免费在线| 亚洲人成电影免费在线| 国产精品永久免费网站| 天堂影院成人在线观看| 免费看光身美女| 午夜福利高清视频| 88av欧美| 国产精品嫩草影院av在线观看 | 高清在线国产一区| 黄色配什么色好看| 日韩欧美精品v在线| 琪琪午夜伦伦电影理论片6080| 欧美丝袜亚洲另类 | 在线播放无遮挡| 亚洲精品成人久久久久久| 嫁个100分男人电影在线观看| 怎么达到女性高潮| 三级男女做爰猛烈吃奶摸视频| 免费高清视频大片| 免费在线观看成人毛片| 伦理电影大哥的女人| 日韩成人在线观看一区二区三区| 中亚洲国语对白在线视频| 尤物成人国产欧美一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品在线观看二区| 亚洲熟妇熟女久久| 国内毛片毛片毛片毛片毛片| 99riav亚洲国产免费| 免费黄网站久久成人精品 | 淫秽高清视频在线观看| 久久99热6这里只有精品| 欧美潮喷喷水| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 我要搜黄色片| 国产久久久一区二区三区| 国产精品亚洲一级av第二区| 身体一侧抽搐| АⅤ资源中文在线天堂| 精品熟女少妇八av免费久了| 少妇的逼水好多| 免费看a级黄色片| 一级毛片久久久久久久久女| 老司机福利观看| 小说图片视频综合网站| 亚洲av成人av| 又爽又黄无遮挡网站| 精品福利观看| 亚洲国产欧洲综合997久久,| 一级黄片播放器| 天天一区二区日本电影三级| 久久中文看片网| 国产精品嫩草影院av在线观看 | av在线蜜桃| 别揉我奶头~嗯~啊~动态视频| 久久欧美精品欧美久久欧美| 国产精品亚洲av一区麻豆| 亚洲精品成人久久久久久| 亚洲国产精品久久男人天堂| 在线免费观看的www视频| 中文字幕久久专区| 午夜福利免费观看在线| 一本综合久久免费| 99国产精品一区二区蜜桃av| 全区人妻精品视频| 亚洲欧美日韩高清专用| 99国产极品粉嫩在线观看| 国产精品久久久久久久久免 | 亚洲人成伊人成综合网2020| 日韩欧美精品v在线| 99热这里只有是精品在线观看 | 麻豆国产97在线/欧美| 99久久精品国产亚洲精品| 99国产精品一区二区三区| 欧美午夜高清在线| 天天躁日日操中文字幕| 高清在线国产一区| 国产高潮美女av| 听说在线观看完整版免费高清| 国产精品久久久久久久久免 | 亚洲片人在线观看| 国产美女午夜福利| 一个人看视频在线观看www免费| 一区福利在线观看| 一级作爱视频免费观看| 国产探花极品一区二区| 天堂动漫精品| 久久精品国产99精品国产亚洲性色| 变态另类成人亚洲欧美熟女| 国产精品三级大全| 欧美中文日本在线观看视频| 91av网一区二区| 中文字幕人妻熟人妻熟丝袜美| 亚洲熟妇中文字幕五十中出| 91狼人影院| 色播亚洲综合网| 熟女电影av网| 日韩 亚洲 欧美在线| 亚洲国产精品999在线| 日韩欧美在线二视频| 国产精品一区二区三区四区久久| 91在线观看av| 男女视频在线观看网站免费| 亚洲成人久久性| 很黄的视频免费| 婷婷色综合大香蕉| 真实男女啪啪啪动态图| 国产精品日韩av在线免费观看| 啦啦啦韩国在线观看视频| 国产精品一区二区三区四区久久| 赤兔流量卡办理| 最近中文字幕高清免费大全6 | 一进一出抽搐gif免费好疼| 国产69精品久久久久777片| 草草在线视频免费看| 天堂网av新在线| 国产精品一区二区三区四区久久| 中文在线观看免费www的网站| 久久婷婷人人爽人人干人人爱| 国产成人av教育| 久久久国产成人精品二区| 老熟妇仑乱视频hdxx| 十八禁人妻一区二区| 国产探花在线观看一区二区| 精华霜和精华液先用哪个| 欧美xxxx性猛交bbbb| 99久久精品国产亚洲精品| 亚洲人成网站在线播放欧美日韩| 天美传媒精品一区二区| 国产成+人综合+亚洲专区| 久久精品影院6| 国产色爽女视频免费观看| 免费电影在线观看免费观看| 色尼玛亚洲综合影院| 一a级毛片在线观看| 亚洲男人的天堂狠狠| 午夜福利欧美成人| 亚洲avbb在线观看| 少妇被粗大猛烈的视频| 欧美日韩亚洲国产一区二区在线观看| 3wmmmm亚洲av在线观看| 校园春色视频在线观看| 午夜a级毛片| 欧洲精品卡2卡3卡4卡5卡区| 五月伊人婷婷丁香| 99久国产av精品| 久久精品久久久久久噜噜老黄 | 亚洲七黄色美女视频| 麻豆久久精品国产亚洲av| 特大巨黑吊av在线直播| 少妇的逼好多水| 欧美精品啪啪一区二区三区| 成人av在线播放网站| 日本 av在线| 成人欧美大片| 亚洲aⅴ乱码一区二区在线播放| 成人毛片a级毛片在线播放| 亚洲人成网站在线播| 色综合站精品国产| 亚洲经典国产精华液单 | 午夜福利欧美成人| 国产伦精品一区二区三区四那| 欧美在线一区亚洲| 97碰自拍视频| 级片在线观看| 极品教师在线免费播放| 成年女人看的毛片在线观看| 特大巨黑吊av在线直播| 国产精品日韩av在线免费观看| 欧美最新免费一区二区三区 | 亚洲精品在线美女| 性欧美人与动物交配| 亚洲国产欧洲综合997久久,| 激情在线观看视频在线高清| 少妇人妻精品综合一区二区 | 国产中年淑女户外野战色| 九九久久精品国产亚洲av麻豆| 国产亚洲精品av在线| 青草久久国产| 国产一区二区在线观看日韩| 成人精品一区二区免费| 俺也久久电影网| 亚洲成人精品中文字幕电影| 精品99又大又爽又粗少妇毛片 | 长腿黑丝高跟| 国产熟女xx| 免费搜索国产男女视频| 动漫黄色视频在线观看| 在线免费观看不下载黄p国产 | 亚洲美女视频黄频| 最新在线观看一区二区三区| 午夜两性在线视频| 伦理电影大哥的女人| a级毛片a级免费在线| 十八禁国产超污无遮挡网站| 村上凉子中文字幕在线| 18禁黄网站禁片午夜丰满| 一区二区三区激情视频| 亚洲精品久久国产高清桃花| 内射极品少妇av片p|