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

    基于模糊聚類篩選和交通可達(dá)性的心血管醫(yī)院選址布局研究

    2016-10-14 12:56:42韓振軍焦建彬
    管理工程學(xué)報(bào) 2016年3期
    關(guān)鍵詞:合格心血管聚類

    趙 韡,韓振軍,焦建彬

    ?

    基于模糊聚類篩選和交通可達(dá)性的心血管醫(yī)院選址布局研究

    趙 韡,韓振軍,焦建彬

    (中國科學(xué)院大學(xué),北京,100049)

    心血管病(CVD)正在成為當(dāng)今影響我國居民健康的主要威脅,如何通過構(gòu)建重點(diǎn)心血管病防治一體化網(wǎng)絡(luò),從而有效地對其進(jìn)行預(yù)防和控制成為當(dāng)務(wù)之急。該問題可歸結(jié)為選址問題,以前的研究成果要么定性的從多屬性中評估候選地址并對其進(jìn)行排序,要么定量建立目標(biāo)函數(shù)求解最優(yōu)結(jié)果。然而地址的選擇問題,尤其對于醫(yī)院這類不僅要考慮醫(yī)院自身各項(xiàng)指標(biāo),同時(shí)也要考慮其他影響(如群眾方便性等)的復(fù)雜問題,求解方式應(yīng)該是定性和定量綜合評估和優(yōu)化的過程。所以,本文在理論上首次把定性和定量(交通可達(dá))模型融合在一起用于解決醫(yī)院選址的問題。首先使用模糊聚類的方法,定性的根據(jù)醫(yī)院的各項(xiàng)指標(biāo)篩選出能夠滿足目標(biāo)的合格醫(yī)院;其次,在從空間分布(即交通可達(dá)性)角度建立目標(biāo)函數(shù)優(yōu)化模型,定量的從所有合格的醫(yī)院當(dāng)中確定最終的入圍醫(yī)院。通過兩步篩選既保證了醫(yī)療資源的優(yōu)質(zhì)性,又能從空間上兼顧醫(yī)院地理上的布局,同時(shí)還引入了權(quán)重參數(shù),方便決策者對結(jié)果進(jìn)行修正。文章最后結(jié)合我國的實(shí)際情況進(jìn)行案例分析,為心血管病醫(yī)院的選址布局提供決策依據(jù)和方法參考。

    心血管?。会t(yī)院選址;模糊聚類

    0引言

    目前,心血管疾病已經(jīng)成為全球面臨的問題,其“發(fā)病率高、致殘率高、死亡率高、復(fù)發(fā)率高,并發(fā)癥多”—“四高一多”的特點(diǎn)[1],使之成為全世界面臨的重大公共衛(wèi)生乃至社會問題,且在發(fā)展中國家亦日益突出。在我國,心血管疾病死亡率長期居于首位,高于腫瘤及其它疾病[2]。近幾年來由中國醫(yī)學(xué)科學(xué)院阜外心血管醫(yī)院和首都醫(yī)科大學(xué)附屬北京安貞醫(yī)院-北京市心肺血管疾病研究生與美國哥倫比亞大學(xué)合作完成的中國冠心病預(yù)測模型是目前針對我國國內(nèi)心血管發(fā)病預(yù)測模型中功能最全面的一個(gè)模型。報(bào)告中的模型指出,中國國內(nèi)由于人口老齡化、高血壓、高血脂、吸煙等諸多因素,到2030年中國男性的心血管發(fā)病率會增加13%,女性增加14%[2]。

    由于心血管疾病的嚴(yán)重性、普遍性和高病死率,衛(wèi)生保健體系亟待采取各種措施降低心血管危險(xiǎn)的發(fā)生和加大相關(guān)醫(yī)療資源的投入[4]。現(xiàn)代心血管醫(yī)院的選址布局,對于構(gòu)建完整合理的心血管病防治一體化網(wǎng)絡(luò)、優(yōu)化地區(qū)醫(yī)療資源的配置至關(guān)重要。心血管醫(yī)院的防治一體化網(wǎng)絡(luò)布局建設(shè)指的是把上級大中型醫(yī)院與社區(qū)衛(wèi)生服務(wù)中心醫(yī)院(簡稱社區(qū)醫(yī)院)二者結(jié)合起來的二級網(wǎng)絡(luò)模式,主要用于防范和治療心血管疾病。該網(wǎng)絡(luò)主要以防治基本原則為依據(jù),“發(fā)展優(yōu)勢, 改善劣勢、抓住機(jī)遇、避免威脅”為指導(dǎo)原則,全面構(gòu)建心血管病防治戰(zhàn)略體系,打通防治各個(gè)環(huán)節(jié),建立防治網(wǎng)絡(luò),充分發(fā)揮各級醫(yī)療機(jī)構(gòu)在防治體系中的作用,提高心血管病防治醫(yī)療服務(wù)能力。具體來說,應(yīng)該由頂層根據(jù)各區(qū)域的發(fā)病率、人口分布特點(diǎn)、醫(yī)院接診能力及交通便利程度等因素,遴選出區(qū)域重點(diǎn)心中心醫(yī)院,劃分社區(qū)醫(yī)療機(jī)構(gòu)歸屬,明確各級醫(yī)療機(jī)構(gòu)防治職責(zé)與任務(wù),形成國家主抓重點(diǎn)醫(yī)院、重點(diǎn)醫(yī)院指導(dǎo)社區(qū)、社區(qū)負(fù)責(zé)居民的層級管理網(wǎng)絡(luò)。大中型醫(yī)院是網(wǎng)絡(luò)的中心,可以從眾多現(xiàn)有醫(yī)院當(dāng)中挑選出合適的幾家[5]。社區(qū)衛(wèi)生服務(wù)中心是對接單位,向上可以對接中心醫(yī)院的資源,向下可以直接面向居民用戶。網(wǎng)絡(luò)建設(shè)的成本主要是中心醫(yī)院的資源和人員配置管理的費(fèi)用。因而,采用科學(xué)有效的方法對大中型醫(yī)院進(jìn)行選址布局,可以降低建設(shè)成本,高效公平的為居民提供醫(yī)療資源。同時(shí)對發(fā)揮重點(diǎn)心血管醫(yī)院的作用、保障患者生命健康以及促進(jìn)醫(yī)療衛(wèi)生系統(tǒng)的平衡快速發(fā)展也有重要意義。

    然而,構(gòu)建心血管病防治一體化網(wǎng)絡(luò),對心血管醫(yī)院進(jìn)行選址布局是一個(gè)復(fù)雜的問題,既要考慮各個(gè)醫(yī)院的客觀指標(biāo),也需要考慮群眾出行的方便性,結(jié)合定性和定量方法進(jìn)行綜合的分析。對于決策者來說,他們既希望客觀快速的篩選醫(yī)院,又常常需要人為的調(diào)節(jié)結(jié)果。于是找到一個(gè)融合定性定量分析,同時(shí)又能動態(tài)可調(diào)節(jié)的模型成為當(dāng)務(wù)之急。

    1 相關(guān)研究

    對于選址布局問題而言,傳統(tǒng)的做法大致可以分成2類,即定性和定量。定性主要通過對地址的主觀或客觀打分,最后獲得最優(yōu)的方案,如以下幾篇文章。文獻(xiàn)[6]使用了要素權(quán)重排序的方法,采用了多種不同要素評估的測量,一些要素是使用主觀評分,另一些是主觀排序,最后通過這些評判結(jié)果對所有地址進(jìn)行打分排序挑選出排在前面的地址。文獻(xiàn)[7]中所有的要素根據(jù)主觀和客觀特性進(jìn)行分組,關(guān)鍵的要素權(quán)重更高,最后通過對這些地址的得分選出候選方案。文獻(xiàn)[8]使用對不同屬性評估其效用值的方法,并給出核電選址的例子。文獻(xiàn)[9]使用對目標(biāo)分配權(quán)重的方法,不同的目標(biāo)根據(jù)模糊集合理論進(jìn)行定義,把一個(gè)線性問題轉(zhuǎn)換成了非線性問題。文獻(xiàn)[10]使用模糊數(shù)學(xué)的相關(guān)理論對不同的屬性進(jìn)行評估,該模型可以評估一些非精確描述的屬性。文獻(xiàn)[11]使用語言描述的方案來評估主觀屬性,并給出了一個(gè)簡單的基于模糊聚類的方法。文獻(xiàn)[12]使用三角模糊隸屬度來評估不同權(quán)威決策者的分析結(jié)果,最后根據(jù)模糊適性指數(shù)排序不同的地址。文獻(xiàn)[13]使用了三重標(biāo)準(zhǔn)評估方案,即最大最小、最小最大、最小和,用于評估凸型區(qū)域的選擇。文獻(xiàn)[14]使用模糊多目標(biāo)的方法建模,并引入遺傳算法進(jìn)行求解,這篇文章首次將模糊集合和遺傳算法結(jié)合在一起求解選擇問題。文獻(xiàn)[15]根據(jù)語言描述來排序不同的候選地址,使用了模糊集合理論。文獻(xiàn)[16]使用了三個(gè)偏好模型來評估不同的候選地址。另外,文獻(xiàn)[17][18][19][20][21]在決策的目標(biāo)函數(shù)中加入了不確定性因素進(jìn)行了研究。文獻(xiàn)[22] [23] [24] [25]基于博弈論對選擇問題進(jìn)行決策。

    定量分析主要采用數(shù)值優(yōu)化的方案,根據(jù)目標(biāo)函數(shù)求解最優(yōu)值,而目標(biāo)函數(shù)通常為成本最小函數(shù)。文獻(xiàn)[26]的整型規(guī)劃是對以醫(yī)療資源的可接近性與資源在地區(qū)內(nèi)的平均分布性作為約束的雙目標(biāo)函數(shù)進(jìn)行求解。文獻(xiàn)[27]中P-中值問題是指選定p個(gè)設(shè)施的位置,使最壞的情況達(dá)到情況最優(yōu)。文獻(xiàn)[28]中選址的目標(biāo)函數(shù)主要基于商品價(jià)格和到市場的距離,這是首次提出這種基于花費(fèi)函數(shù)的選址決策理論。文獻(xiàn)[29]地址的選址考慮了運(yùn)輸、人力和工業(yè)分布等成本,由于這個(gè)模型并沒有考慮地價(jià)等因素,因而相對比較簡單。文獻(xiàn)[30]選址模型考慮了利潤而并非只有成本花費(fèi),最終目標(biāo)是讓利潤最大化。文獻(xiàn)[31]使用優(yōu)化模型來評估廠站的選址,考慮的因素包括交通和資源等,不過文章的前提假設(shè)過于簡單。文獻(xiàn)[32]選址使用了電勢點(diǎn)的概念。

    以上的選址方案理論研究成果要么定性的從多屬性中評估候選地址并對其進(jìn)行排序,要么定量建立目標(biāo)函數(shù)求解最優(yōu)結(jié)果。其實(shí)地址的選擇問題,尤其對于醫(yī)院這類不僅要考慮醫(yī)院自身各項(xiàng)指標(biāo),同時(shí)也要考慮其他影響(如群眾方便性等)的復(fù)雜問題,求解方式應(yīng)該是定性和定量綜合評估和優(yōu)化的過程。所以,本文在理論上首次把定性和定量(交通可達(dá))模型融合在一起用于解決醫(yī)院選址的問題。具體來說,本文采取兩步?jīng)Q策的步驟。首先使用模糊聚類的方法,定性的根據(jù)醫(yī)院的各項(xiàng)指標(biāo)篩選出能夠滿足目標(biāo)的合格醫(yī)院;其次,在從空間分布(即交通可達(dá)性)角度建立目標(biāo)函數(shù)優(yōu)化模型,定量的從所有合格的醫(yī)院當(dāng)中確定最終的入圍醫(yī)院。通過這兩個(gè)過程,我們一方面可以保證優(yōu)秀的醫(yī)院資源得到合理利用,另一方面又考慮了交通便利,方便了各社區(qū)醫(yī)院到中心醫(yī)院的便利性。同時(shí)還引入了權(quán)重參數(shù),方便決策者對結(jié)果進(jìn)行修正。在目前我國醫(yī)患關(guān)系問題突出的背景下,使用本文的兩步方法決策,可以充分體現(xiàn)對人民群眾的重視,也展現(xiàn)了以人為本的決策思路。文章最后用真實(shí)數(shù)據(jù)求解模型并得出了結(jié)果。

    2 問題描述和建模思路

    2.1 問題描述

    圖1 心血管醫(yī)院選址問題空間描述

    問題的描述如圖1所示,該城市有候選的n家醫(yī)院(十字),問題的目標(biāo)就是要從這n家醫(yī)院當(dāng)中選出b家來構(gòu)建心血管防治一體化網(wǎng)絡(luò)的中心醫(yī)院,同時(shí)該網(wǎng)絡(luò)布局使得社區(qū)醫(yī)院(居民)交通上最為便利,使各社區(qū)醫(yī)院(圓圈)可以方便到達(dá)網(wǎng)絡(luò)里某一家中心醫(yī)院即可。于是,本文選址的方法主要分成兩步,第一步根據(jù)醫(yī)院的各項(xiàng)指標(biāo)用模糊C均值聚類方法把醫(yī)院分成兩類,一類是合格,另一類是不合適,然后從合格的醫(yī)院里篩選取得分最高的h個(gè)醫(yī)院作為第二步的輸入;第二步通過交通可達(dá)模型確定醫(yī)院最終選址布局(從h家合格醫(yī)院里選出b家,b

    2.2 建模思路

    醫(yī)院地址的選擇問題需要把定性和定量模型結(jié)合在一起。首先,我們假設(shè)候選的醫(yī)院是m維空間的對象,每一維空間代表某一屬性或指標(biāo)。首先我們采用定性的方法,根據(jù)醫(yī)院本身的屬性或指標(biāo)進(jìn)行聚類,找出合格的醫(yī)院集合。接著,對于空間分布這個(gè)比較中性的屬性,可以采用定量的方法,建立交通可達(dá)目標(biāo)函數(shù),再從這些合格醫(yī)院當(dāng)中找出最優(yōu)的醫(yī)院布局。

    2.3 定性篩選原則

    考慮哪些因素來篩選醫(yī)院非常重要,一般來說,選址的因素可以分成兩類,即主觀因素和客觀因素。對于主觀因素,我們主要用定性的方式來描述,如差、一般、非常好等。而對于客觀因素則可以使用定量的方式來描述,如年門診人數(shù)等。醫(yī)院的篩選應(yīng)該綜合考慮主觀因素和客觀因素。由于本文是針對心血管疾病的醫(yī)院進(jìn)行選擇,因而客觀上需要考慮心血管專科的排名和心血管床位數(shù),這是最能代表該??颇芰Φ膬蓚€(gè)指標(biāo)。另外,也需要考慮醫(yī)院的綜合影響,所以把日門診量和綜合排名也考慮進(jìn)來。最后,患者的滿意度也是醫(yī)院評估的重要考量,這是一個(gè)主觀的因素。綜上分析,我們主要考量以下5個(gè)因素

    1)心血管??婆琶?/p>

    2)心血管床位數(shù)

    3)患者滿意度

    4)日門診量

    5)醫(yī)院綜合排名

    考慮到不同因素的值區(qū)間不一致,我們最終通過歸一化處理,把這5個(gè)因素的值映射到數(shù)值0-100之間。

    另外,考慮到各要素對醫(yī)院的評估貢獻(xiàn)是不一樣的,我們引入要素分配權(quán)重,區(qū)間為[0,1],值越大表示越重要。不同的權(quán)重使得各要素歸一化處理的過程中可以得到不同的權(quán)值。這樣我們分配好各權(quán)重值后,便可以在后面的距離值計(jì)算模型中對各要素直接相加和對比。歸一化公式如下

    max為該指標(biāo)樣本數(shù)據(jù)的最大值,min為該樣本數(shù)據(jù)的最小值。就是樣本數(shù)據(jù),為歸一化后的值。為權(quán)重取值[0,1]。

    2.4 定性篩選邏輯

    通常對于數(shù)據(jù)的聚類依托于對象屬性的內(nèi)在關(guān)系,不同的內(nèi)在關(guān)系可能會把數(shù)據(jù)分成1到n類。本文聚類篩選的目的主要是挑選出合格的醫(yī)院,因而可以人為定義成兩類,即“合格”與“不合格”。為了可以明確的區(qū)分哪一類是“合格”,哪一類是“不合格”,我們又人為的引入兩個(gè)虛擬醫(yī)院,即“最好醫(yī)院”和“最差醫(yī)院”。這兩個(gè)虛擬醫(yī)院的屬性值就是從所有的候選醫(yī)院對應(yīng)的值中選擇最好的值或最壞的值。經(jīng)過模糊聚類后,包含“最好醫(yī)院”醫(yī)院的組就是合格的組。

    2.5 定量優(yōu)化機(jī)制

    定量優(yōu)化主要考慮的是醫(yī)院的空間分布,對決策者來說就是交通可達(dá)性。交通可達(dá)性的不同對救治的及時(shí)性和人們的出行有很大的影響。而本文交通可達(dá)性的主要衡量方法就是計(jì)算社區(qū)醫(yī)院到達(dá)某中心醫(yī)院所花費(fèi)的時(shí)間,時(shí)間越短我們認(rèn)為交通可達(dá)性越好。衡量某家中心醫(yī)院的交通可達(dá)性,我們可以計(jì)算不同社區(qū)醫(yī)院到這家醫(yī)院所花費(fèi)的時(shí)間平均值。于是,建立一個(gè)目標(biāo)函數(shù),使得該方案下所得到的時(shí)間花費(fèi)最小來選擇出最后方案。

    3 模型建立和求解

    3.1模糊C均值聚類定性篩選模型

    根據(jù)上一章節(jié)的分析,我們首先目標(biāo)是把候選醫(yī)院首先分成合格醫(yī)院與不合格醫(yī)院。具體方法就是加入兩個(gè)虛擬的醫(yī)院,一個(gè)是理想的合格醫(yī)院(BestHospital),一個(gè)是最差的不合格醫(yī)院(WorstHospital)。對于BestHospital,我們設(shè)置其所有的考慮因素得分是候選醫(yī)院里最好的值,相反,對于WorstHospital,我們設(shè)置所有的得分都是最差的值。

    根據(jù)模糊C均值模型[33],目標(biāo)最小化以下目標(biāo)函數(shù)

    (3)

    而每個(gè)集群的中心v用以下公式獲得

    n表示待分類的對象數(shù)目,m表示對象的維度(考慮因素?cái)?shù)目),里面U表示聚類矩陣,是在聚類k的隸屬度,是第k個(gè)對象,表示分類i聚類中心,c是聚類的數(shù)目,z是加權(quán)指數(shù)。

    然后,根據(jù)合格醫(yī)院集合,計(jì)算每個(gè)醫(yī)院與BestHospital的距離值,并選取最高的h個(gè)醫(yī)院作為第二步的輸入。具體的距離值計(jì)算模型如下

    3.2交通可達(dá)定量優(yōu)化模型

    根據(jù)上面模型的結(jié)果,我們選取了h家得分最高的醫(yī)院,但并沒有考慮他們的空間布局。所以接下來,我們通過交通可達(dá)模型從這h家中選出b家作為最終的方案??紤]到選出的這h家合格醫(yī)院也有不同,決策者有時(shí)需要人為調(diào)節(jié)結(jié)果或人們有時(shí)也不管距離多遠(yuǎn)都會去某家中心醫(yī)院看相關(guān)疾病,于是我們引入權(quán)重值來表示醫(yī)院j的權(quán)重以區(qū)分不同醫(yī)院的重要性。的本質(zhì)就是把重要醫(yī)院的距離(可達(dá)時(shí)間)縮短,這樣可以統(tǒng)一歸結(jié)到可達(dá)時(shí)間這一個(gè)目標(biāo)上來。

    首先,我們以最短可達(dá)時(shí)間作為交通可達(dá)模型的評價(jià)指標(biāo),于是各社區(qū)醫(yī)院到中心醫(yī)院的平均可達(dá)時(shí)間的模型為:

    p為統(tǒng)計(jì)的社區(qū)醫(yī)院數(shù)目

    根據(jù)公式(6),我們就可以評估某方案X(b個(gè)醫(yī)院集合)的交通可達(dá)指標(biāo)值,具體如下

    Max(8)

    3.3模型求解

    3.3.1模糊C均值聚類模型求解

    我們假設(shè)候選的醫(yī)院為n家,評價(jià)的要素為m個(gè),第一步的目的是從這n家醫(yī)院選出h(h<=n)家合格的中心醫(yī)院。于是對每個(gè)醫(yī)院每個(gè)因素評分,可以形成一個(gè)n*m的矩陣。第k個(gè)醫(yī)院的第j因素的評分用表示,0

    1) 初始化矩陣U

    3) 利用如下公式(9)更新聚類矩陣

    5) 根據(jù)公式(5)的結(jié)果計(jì)算合格的醫(yī)院評分,然后根據(jù)評分選出前h家醫(yī)院作為結(jié)果

    3.3.2交通可達(dá)模型求解

    求解通過O-D矩陣分析獲得,我們借助ArcGIS的網(wǎng)絡(luò)分析模塊。具體做法就是通過公式(8),從可行的候選方案集合X中選擇一個(gè)方案,并把這方案中的合格醫(yī)院作為目的點(diǎn)設(shè)置在地圖上;另外設(shè)置典型的p個(gè)社區(qū)醫(yī)院作為起點(diǎn)。同時(shí)設(shè)置每條道路的通行時(shí)間(公共交通和其他方式分別設(shè)置)。經(jīng)過計(jì)算,我們可以得到p個(gè)社區(qū)醫(yī)院到h個(gè)合格醫(yī)院的通行時(shí)間,用矩陣p*h來表示。根據(jù)該矩陣,我們就可通過公式(6)計(jì)算交通可達(dá)指標(biāo)值。最后求出最優(yōu)的方案。

    4 結(jié)果分析與討論

    為驗(yàn)證模型和求解算法的有效性,本節(jié)根據(jù)前述分析,以北京的10所醫(yī)院作為候選集為例(n=10),來計(jì)算最終方案。本文數(shù)據(jù)綜合醫(yī)院排名和??婆虐鎭碜訦C3i中國數(shù)字醫(yī)療網(wǎng),其他數(shù)據(jù)來自各醫(yī)院的官方網(wǎng)站。實(shí)例當(dāng)中,候選的醫(yī)院有10家,我們的目標(biāo)是從這10家醫(yī)院當(dāng)中選出3家作為最終的方案。b=3表示問題的目標(biāo)是選出3家,由專家規(guī)劃決定。

    第一步,根據(jù)不同因素的值篩選出5家合格的醫(yī)院,h=5表示總的集合的一半,因?yàn)槲覀冏罱K把醫(yī)院分組2組,best和worst,所以就取了中間值5。第二步再從這5家合格醫(yī)院里選出交通可達(dá)的最優(yōu)3家(b=3)。為了驗(yàn)證我們第二步模型中權(quán)重的有效性,我們先假設(shè)所有醫(yī)院的權(quán)重值都為1,計(jì)算出對應(yīng)結(jié)果。然后,根據(jù)實(shí)際情況認(rèn)為某一家醫(yī)院是必須入選的,我們把其權(quán)重改為10,然后再計(jì)算對應(yīng)結(jié)果,并對比兩次結(jié)果的不同。其他參數(shù)m=5,z=2,,因素表見表1(經(jīng)過公式1歸一化處理,權(quán)重設(shè)置都為1)。實(shí)質(zhì)是醫(yī)院的平均可達(dá)時(shí)間的模型中的時(shí)間倍數(shù),等于1是表示時(shí)間不變,等于10時(shí)表示10倍的時(shí)間。這里設(shè)置10表示把該醫(yī)院在平均可達(dá)時(shí)間修改為原來10倍,主要是為了顯示出該系數(shù)改變后的效果,從地圖上,設(shè)置到10倍距離肯定可以排除該醫(yī)院。m=5表示5個(gè)評分要素。z=2是根據(jù)[35]中的試驗(yàn)結(jié)果得到,建議取值2。是算法結(jié)束的精度。

    表1 10所醫(yī)院候選集及因素評分

    首先加入BestHostpital和WorstHopstal,見表2。

    表2 BestHostpital和worstHopstal的因素評分

    通過公式2,計(jì)算出聚類值,見表3。

    根據(jù)聚類2排除不合格醫(yī)院(聚類2中值>0.5,該值表示相似度,區(qū)間是[0,1],我們選擇中間值0.5)1,3,8,10 醫(yī)院。從排除結(jié)果來看,所有的沒有心血管專科的醫(yī)院(心血管專科排名=100,床位=0),即1,3,10這幾家醫(yī)院都被排除,而且他們的聚類值都非常高(大于0.9),符合常理,也進(jìn)一步表明模型的正確性。

    接著,根據(jù)聚類1選擇合格的6(聚類1中值>0.5, 聚類2中值>0.5,該值表示相似度,區(qū)間是[0,1],我們選擇中間值0.5)家醫(yī)院,然后根據(jù)公式5計(jì)算合格醫(yī)院的打分值如下表4所示。

    表3 聚類評分結(jié)果

    表4 合格醫(yī)院打分結(jié)果

    選擇得分最小的5家醫(yī)院2, 醫(yī)院4, 醫(yī)院5, 醫(yī)院6, 醫(yī)院9進(jìn)入下一步。醫(yī)院7被排除也符合預(yù)期,從專業(yè)排名上看,7的排名在這里面是最靠后的,同時(shí),他的床位數(shù)也是最少的。

    全北京(5環(huán)內(nèi))的社區(qū)醫(yī)院大概有600家,而社區(qū)醫(yī)院的分布一定程度上代表了患者的需求分布。為了考慮需求分布及簡化模型計(jì)算,我們把時(shí)間精度控制在15分鐘內(nèi),按照市區(qū)時(shí)速35公里計(jì)算,15分鐘大概可以走9公里。因而我們把5環(huán)內(nèi)區(qū)域(大概直徑27公里)劃分為3*3=9個(gè)區(qū)域,每個(gè)區(qū)域直徑大概9公里。我們選取15家(平均大概50平方公里內(nèi)一家)作為典型的社區(qū)醫(yī)院作為計(jì)算依據(jù),按照5環(huán)內(nèi)大概有600家社區(qū)醫(yī)院計(jì)算(社區(qū)醫(yī)院的分布通過百度地圖搜索社區(qū)醫(yī)院得到),也就是每40家社區(qū)醫(yī)院里面選取一家代表。這樣我們可以根據(jù)9個(gè)區(qū)域包含的社區(qū)醫(yī)院數(shù)目,決定每個(gè)區(qū)域有幾家典型社區(qū)醫(yī)院。

    通過GIS,我們根據(jù)上面的說明,計(jì)算出典型的15(p=15)個(gè)社區(qū)醫(yī)院作為起點(diǎn),具體分布見圖2,同時(shí)計(jì)算出到這5家中心醫(yī)院的平均可達(dá)時(shí)間,=0.5(表示乘坐公共交通的比例為50%),t=45(表示到中心醫(yī)院時(shí)間如果超過45分鐘就不計(jì)入模型計(jì)算),如表5所示。

    圖2 五家合格醫(yī)院及典型社區(qū)醫(yī)院分布地圖

    表5 典型社區(qū)醫(yī)院到中心醫(yī)院平均可達(dá)時(shí)間

    接著,利用公式(7)計(jì)算可能組合的交通可達(dá)值,如表6所示。

    可見組合方案9是最優(yōu)方案(得分86.78322最低),最后選擇醫(yī)院2,醫(yī)院9和醫(yī)院5,分布結(jié)果見圖3。從圖上看,比較集中的幾家醫(yī)院(6,2,4)我們只選取了一家,這樣在整體上,2,9,5的分布性最好,可以更好的方便群眾出行。

    表6 交通可達(dá)得分

    圖3 不考慮權(quán)重指數(shù)的最終醫(yī)院選取方案

    但是,實(shí)際情況是醫(yī)院6是必須入選的醫(yī)院,因?yàn)槠渑琶谝?,是全國最好的醫(yī)院,大家通常不考慮距離也會去,決策者也把它列入必選項(xiàng)。于是,我們改變模型中的權(quán)重系數(shù)。繼續(xù)計(jì)算,結(jié)果如下:

    方案序號醫(yī)院組合得分 12,6,465.928024 22,6,961.402046 32,6,562.402046 46,4,962.781174 56,4,563.781174 66,9,559.255196 72,4,990.3092 82,4,591.3092 92,9,586.78322 104,9,588.16234

    5 結(jié)論

    隨著人們生活水平的提高,心血管疾病問題越發(fā)突出,因而構(gòu)建心血管防治一體化網(wǎng)絡(luò),篩選出重點(diǎn)心血管醫(yī)院成為亟待解決的問題。本文通過運(yùn)用模糊聚類篩選合格醫(yī)院,然后通過交通可達(dá)模型計(jì)算出最優(yōu)的空間布局方案。從計(jì)算結(jié)果來看,該模型給出的選擇結(jié)果與我們的預(yù)期基本一致,能給實(shí)際的選址工作提供一種全新的方法。主要優(yōu)點(diǎn)如下

    1)模型采用模糊聚類的方法,可以很方便的擴(kuò)展候選對象的屬性空間,因而對于需要考慮大規(guī)模維度的對象選擇問題,使用該方法可以快速收斂得到有效值。

    2)模型采用的分步計(jì)算的方法,因而可以很方便的對每一步的結(jié)果進(jìn)行修正和驗(yàn)證。比如第一步計(jì)算的結(jié)果沒有包括想要的醫(yī)院,我們也可以人為加入想要的醫(yī)院直接進(jìn)入第二步優(yōu)化計(jì)算。

    3)模型引入了權(quán)重指數(shù),使得可以更方便的人為修正結(jié)果成為可能。

    雖然該模型取得了一定的效果,當(dāng)將來還是可以從以下幾個(gè)方面進(jìn)行擴(kuò)展或更深度的研究

    1)文中心血管醫(yī)院評價(jià)因素的選擇是人工挑選的,在一定程度上存在主觀性,因而將來可以對候選醫(yī)院因素的選擇也建立模型,使模糊聚類的方法更加有效。

    2)文中第二步優(yōu)化問題的計(jì)算是窮盡所有結(jié)果得到,對于大規(guī)模的計(jì)算有一定影響,可以在將來的研究當(dāng)中引入遺傳算法或其他算法。

    總之,我們的模型對于實(shí)際工作有借鑒意義,能為國家和醫(yī)療機(jī)構(gòu)的決策提供更精細(xì)化的參考工具。

    [1] 楊捷. 社區(qū)居民心血管疾病相關(guān)危險(xiǎn)因素調(diào)查研究[J]. 長春醫(yī)學(xué), 2007, 5(2): 11

    [2] 衛(wèi)生部. 中國衛(wèi)生統(tǒng)計(jì)年鑒[M]. 中國協(xié)和醫(yī)科大學(xué), 2003.

    [3] Moran A, Gu D, Zhao D, et al. Future Cardiovascular Disease in China Markov Model and Risk Factor Scenario Projections From the Coronary Heart Disease Policy Model–China[J]. Circulation: Cardiovascular Quality and Outcomes, 2010, 3(3): 243-252.

    [4] 胡盛壽, 孔靈芝. 中國心血管病報(bào)告 2007[M]. 北京:衛(wèi)生部心血管病防治研究中心, 2008: 2-34

    [5] 王雪飛. 心血管病防治一體化進(jìn)程漸漸提速[J]. 中國保健營養(yǎng), 2005, 第2期:10-11.

    [6] Keeney R L, Raiffa H. Decisions with multiple objectives: preferences and value trade-offs [M]. Cambridge university press, 1993.

    [7] Brown P A, Gibson D F. A quantified model for facility site selection-application to a multiplant location problem [J]. AIIE transactions, 1972, 4(1): 1-10.

    [8] Kirkwood C W. A case history of nuclear power plant site selection [J]. Journal of the operational research society, 1982: 353-363.

    [9] Yager R R. Fuzzy decision making including unequal objectives [J]. Fuzzy sets and systems, 1978, 1(2): 87-95.

    [10] Narasimhan R. A FUZZY SUBSET CHARACTERIZATION OF A SITE‐SELECTION PROBLEM[J]. Decision Sciences, 1979, 10(4): 618-628.

    [11] Raoot A D, Rakshit A. A ‘fuzzy ‘a(chǎn)pproach to facilities lay-out planning [J]. The International Journal of Production Research, 1991, 29(4): 835-857.

    [12] Liang G S, Wang M J J. A fuzzy multi-criteria decision-making method for facility site selection [J]. The International Journal of Production Research, 1991, 29(11): 2313-2330.

    [13] Bhattacharya U, Rao J R, Tiwari R N. Fuzzy multi-criteria facility location problem [J]. Fuzzy Sets and Systems, 1992, 51(3): 277-287..

    [14] Tzeng G H, Chen Y W. The optimal location of airport fire stations: a fuzzy multi‐objective programming and revised genetic algorithm approach[J]. Transportation Planning and Technology, 1999, 23(1): 37-55.

    [15] Chu T C. Facility location selection using fuzzy TOPSIS under group decisions [J]. International journal of uncertainty, fuzziness and knowledge-based systems, 2002, 10(06): 687-701.

    [16] CINAR D. FACILITY LOCATION SELECTION USING A FUZZY OUTRANKING METHOD[C]//Applied Artificial Intelligence: Proceedings of the 7th International FLINS Conference, Genova, Italy, 29-31 August 2006. World Scientific, 2006: 359.

    [17] 陳青豐, 魯建廈, 劉敏. 非對稱 AHP 方法在物流中心選址中的應(yīng)用[J]. 工業(yè)工程, 2005, 8(1): 75-78.

    [18] 劉曉峰, 陳通, 柳錦銘, 等. 基于 ANP 的物流中心選址[J]. 工業(yè)工程, 2007, 10(5): 136-140.

    [19] LIU Xiaofeng, CHEN Tong, LIU Jinming, et al .The ANP based evaluation model and application of the site selection of logistics center [J]. Industrial Engineering, 2007, 10 (5): 136-140.

    [20] 談英姿, 沈炯, 呂震中. 免疫優(yōu)化算法及其前景展望[J]. 信息與控制, 2002, 31(5): 385-390.

    [21] TAN Yingzi, SHEN Jiong, LU Zhenzhong. Study of immune based optimization algorithm and its future development [J]. Information and Control, 2002, 31(5): 385-390.

    [22] 周華珍, 王花蘭, 盧柏蓉, 等. 基于靜態(tài)博弈論的物流配送中心選址優(yōu)化研究[J]. 山東交通學(xué)院學(xué)報(bào), 2014, 22(3): 31-34.

    [23] 李東, 晏湘濤, 匡興華. 考慮設(shè)施失效的軍事物流配送中心選址模型[J]. 計(jì)算機(jī)工程與應(yīng)用, 2010, 46(11): 3-6.

    [24] 周洪超, 李海鋒. 基于博弈論的電動汽車充電站選址優(yōu)化模型研究[J]. 科技和產(chǎn)業(yè), 2011, 11(2): 51-54.

    [25] 李如琦, 蘇浩益. 基于博弈論的核電廠選址優(yōu)化模型[J]. 電網(wǎng)技術(shù), 2011, 35(2): 216-220.

    [26] P. Mitropoulos, I. Mitropoulos, I. Giannikos, A. Sissouras, A Bi-objective model for the location planning of hospitals and health centers [J]. Health Care and Management Science, 2006(9): P171-179

    [27] Brandeau M, Sanfort F, Pierskallas W P. Operations research and health care: A handbook of methods and applications [M]. United States: Springer Publication, 2004.

    [28] von Thünen J H, Braeuer W. Der isolierte staat in beziehung auf landwirtschaft und national?konomie[M]. Fischer, 1960.

    [29] Giersch H. Economic union between nations and the location of industries [J]. The Review of Economic Studies, 1949, 17(2): 87-97.

    [30] Martin P, Rogers C A. Industrial location and public infrastructure [J]. Journal of International Economics, 1995, 39(3): 335-351.

    [31] Koopmans T C, Beckmann M. Assignment problems and the location of economic activities [J]. Econometrica: journal of the Econometric Society, 1957: 53-76.

    [32] Muther R, Webster D B. Plant Layout and Materials Handling [J]. Materials Handling Handbook, 2nd Edition, 1985: 19-77.

    [33] 高新波, 裴繼紅, 謝維信. 模糊 c-均值聚類算法中加權(quán)指數(shù) m 的研究[J]. 2000.

    [34] 毛罕平, 張艷誠, 胡波. 基于模糊 C 均值聚類的作物病害葉片圖像分割方法研究[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2008, 24(9): 136-140.

    [35] Pal N R, Bezdek J C. On cluster validity for the fuzzy c-means model [J]. Fuzzy Systems, IEEE Transactions on, 1995, 3(3): 370-379.

    Fuzzy C Mean and Commute TimeBased Optimizing Locations of Cardiovascular Hospitals

    ZHAO Wei ,HAN Zhen-jun,JIAO Jian-bin

    ( University of Chinese Academy of Science (UCAS),Beijing 100049,China)

    Cardiovascular diseases have become a global problem, especially in developing countries. The long-term cardiovascular mortality is always one of top problems in China. In recent years, the Chinese Academy of Medical Sciences Fu Wai Hospital and other hospitals design a comprehensive heart disease prediction model. It is important to pay more attention to the problem because of the severity of cardiovascular diseases. To solve this problem, we must build an integrated model to prevent cardiovascular diseases. The integrated network consists of the centered hospitals and community health centers. Specifically, the centered hospital is chosen from a few hospitals in a city dealing with cardiovascular disseases . Community health center is the docking unit. How to use scientific and effective method to choose a right centered hospital is critical to the integrated network.The first part is the background and introduction. As for how to choose centered hospitals, the traditional approach can be divided into two types, namely qualitative or quantitative. Few articles integrate these two approaches. Therefore, this is the first article in the current literature trying to integrate qualitative and quantitative models in order to solve the problem for the chosen hospitals. Specifically, a two-step decision model is built in this article. First, we use fuzzy clustering methods, and qualitatively choose qualified hospitals to address our research questions. As a result, the utilization of hospital resources can be achieved. Finally, we solve the model with real data. The second part describes our modeling framework. First, we assume that candidate hospitals are located is m-dimensional space, with each dimension of the space representing one hospital. We then use qualitative methods to identify qualified hospitals. Next, quantitative method is used to find the best hospitals from the first step.The third part provides detail information about the model. For qualitative model based on fuzzy C-means clustering, our first goal is to choose qualified candidates. The method is to add two virtual hospitals: Best Hospital and Worst Hospital. For Best Hospital, we set all of its score from the best candidates. On the contrary, for Worst Hospital, we set all of its scores as the worst values. For transportation quantitative optimization model, we select final hospitals from the first step based on space distribution. Sometimes we need manual adjustments so that we can include weights into the hospital model. The fourth part is the result analysis and discussion. To verify the validity of the model, this section introduces 10 hospitals in Beijing as the candidate set. The ranking data is from HC3i Chinese digital medical network. Other data are collected from the official website of each hospital. Our goal is to choose 3 hospitals as the final solution from these 10 hospitals. In the first step, we identify 5 qualified hospitals. Afterwards, we select the final 3 hospitals from based on the analysis in the second step. After checking the result with the real map, we believe this model is as good as we expected.The fifth part is the conclusion. The main advantages are summarized as follows. First, the model uses fuzzy clustering which can be easily extended to analyzing similar problems in large-scale dimensions. Secondly, the method uses a two-step model so that it is convenient to adjust and verify each step’s result. For example, the results of the first step in the calculation do not include hospitals we want. Instead, we artificially add one hospital directly to the hospital list of the second step. Third, the model introduces weights parameters, making it easier for us to adjust the model.

    cardiovascular hospital; optimizing location; fuzzy c mean

    中文編輯:杜 ??;英文編輯:Charlie C. Chen

    O223

    A

    1004-6062(2016)03-0202-07

    10.13587/j.cnki.jieem.2016.03.025

    2015-10-09

    2016-01-04

    國家科技支撐計(jì)劃課題(2011BAI11B02)

    趙韡(1974—),男,北京人;中國科學(xué)院大學(xué)在讀博士,高級工程師,研究方向:生物醫(yī)學(xué)工程,醫(yī)學(xué)信息學(xué)。

    猜你喜歡
    合格心血管聚類
    COVID-19心血管并發(fā)癥的研究進(jìn)展
    基于DBSACN聚類算法的XML文檔聚類
    電子測試(2017年15期)2017-12-18 07:19:27
    我是合格的小會計(jì)
    誰不合格?
    做合格黨員
    大江南北(2016年8期)2016-02-27 08:22:46
    基于改進(jìn)的遺傳算法的模糊聚類算法
    lncRNA與心血管疾病
    胱抑素C與心血管疾病的相關(guān)性
    一種層次初始的聚類個(gè)數(shù)自適應(yīng)的聚類方法研究
    自適應(yīng)確定K-means算法的聚類數(shù):以遙感圖像聚類為例
    国产精品综合久久久久久久免费| 国产精品嫩草影院av在线观看 | 欧美又色又爽又黄视频| www.色视频.com| 岛国在线免费视频观看| 亚洲av免费高清在线观看| 欧美中文综合在线视频| 99久久99久久久精品蜜桃| 我的老师免费观看完整版| www日本黄色视频网| 国产欧美日韩精品亚洲av| 亚洲精品久久国产高清桃花| 欧美国产日韩亚洲一区| 波野结衣二区三区在线 | 亚洲国产色片| 成人一区二区视频在线观看| 欧美三级亚洲精品| 亚洲精品美女久久久久99蜜臀| 成人午夜高清在线视频| 亚洲国产中文字幕在线视频| h日本视频在线播放| 免费看a级黄色片| 亚洲国产欧美人成| 午夜激情欧美在线| 午夜视频国产福利| 免费在线观看影片大全网站| 亚洲国产欧美人成| 免费在线观看亚洲国产| 中文资源天堂在线| 国产三级中文精品| 欧美性猛交╳xxx乱大交人| 国产精品美女特级片免费视频播放器| 丁香六月欧美| 97碰自拍视频| 露出奶头的视频| 精品熟女少妇八av免费久了| 成人精品一区二区免费| 12—13女人毛片做爰片一| 97超视频在线观看视频| 亚洲精华国产精华精| 国产69精品久久久久777片| 免费一级毛片在线播放高清视频| 亚洲 国产 在线| 亚洲18禁久久av| 日韩欧美精品免费久久 | 一个人看的www免费观看视频| 丁香欧美五月| 欧美激情在线99| 精品日产1卡2卡| 搞女人的毛片| 国产不卡一卡二| 亚洲av第一区精品v没综合| 亚洲成a人片在线一区二区| 午夜免费男女啪啪视频观看 | 成年女人毛片免费观看观看9| 韩国av一区二区三区四区| 亚洲av成人不卡在线观看播放网| 欧美性猛交╳xxx乱大交人| 老司机午夜十八禁免费视频| 无限看片的www在线观看| 国产真实伦视频高清在线观看 | 长腿黑丝高跟| 久久午夜亚洲精品久久| 天堂网av新在线| 国产午夜福利久久久久久| 国产精品,欧美在线| 国模一区二区三区四区视频| 啦啦啦观看免费观看视频高清| 成人特级黄色片久久久久久久| 午夜免费观看网址| 精品一区二区三区视频在线 | 日韩欧美免费精品| 精品久久久久久成人av| 网址你懂的国产日韩在线| av欧美777| 夜夜爽天天搞| 夜夜躁狠狠躁天天躁| 身体一侧抽搐| 国内毛片毛片毛片毛片毛片| 精品99又大又爽又粗少妇毛片 | 免费av观看视频| 国产精品久久久久久亚洲av鲁大| 极品教师在线免费播放| 亚洲国产中文字幕在线视频| 欧美区成人在线视频| 香蕉av资源在线| 午夜老司机福利剧场| 制服人妻中文乱码| 91久久精品电影网| 国产精品女同一区二区软件 | 欧美日韩福利视频一区二区| 久久精品亚洲精品国产色婷小说| 亚洲欧美日韩无卡精品| 亚洲av熟女| 亚洲人成网站在线播| 又黄又爽又免费观看的视频| 高清日韩中文字幕在线| 熟女人妻精品中文字幕| 色老头精品视频在线观看| eeuss影院久久| 99国产精品一区二区蜜桃av| 天堂网av新在线| 日本免费一区二区三区高清不卡| 免费大片18禁| www.熟女人妻精品国产| 亚洲成人免费电影在线观看| 欧美日韩中文字幕国产精品一区二区三区| a在线观看视频网站| 国产色爽女视频免费观看| 精品国产亚洲在线| 俄罗斯特黄特色一大片| 99久久成人亚洲精品观看| 国产视频一区二区在线看| 97超视频在线观看视频| 国内精品久久久久久久电影| 久久精品国产综合久久久| 亚洲乱码一区二区免费版| 国产黄色小视频在线观看| 特大巨黑吊av在线直播| 国产亚洲精品久久久久久毛片| 丰满的人妻完整版| 国产成人aa在线观看| 波多野结衣巨乳人妻| 亚洲乱码一区二区免费版| 女人十人毛片免费观看3o分钟| 熟女少妇亚洲综合色aaa.| 婷婷精品国产亚洲av| 91麻豆精品激情在线观看国产| 免费观看的影片在线观看| 午夜精品一区二区三区免费看| 一进一出抽搐gif免费好疼| 欧美性猛交黑人性爽| 国产真实乱freesex| 欧美高清成人免费视频www| 色在线成人网| 在线观看午夜福利视频| 波多野结衣巨乳人妻| 亚洲成人精品中文字幕电影| 桃红色精品国产亚洲av| 好看av亚洲va欧美ⅴa在| 国产一区在线观看成人免费| 精品欧美国产一区二区三| 小蜜桃在线观看免费完整版高清| 国产真人三级小视频在线观看| 久久精品91无色码中文字幕| 黑人欧美特级aaaaaa片| 亚洲美女黄片视频| 久久久久九九精品影院| 看免费av毛片| 一a级毛片在线观看| 亚洲精品久久国产高清桃花| 亚洲国产欧美网| 精华霜和精华液先用哪个| 亚洲精华国产精华精| 午夜精品一区二区三区免费看| 久久欧美精品欧美久久欧美| 欧美激情在线99| 99热这里只有是精品50| 亚洲专区中文字幕在线| 国产精品久久视频播放| 午夜免费男女啪啪视频观看 | 久久国产精品人妻蜜桃| 国产午夜精品久久久久久一区二区三区 | 51国产日韩欧美| 欧美日韩国产亚洲二区| 国产免费av片在线观看野外av| 久9热在线精品视频| 亚洲成av人片在线播放无| 久久国产乱子伦精品免费另类| 亚洲欧美一区二区三区黑人| 国产精品电影一区二区三区| 国产不卡一卡二| 午夜精品在线福利| 91九色精品人成在线观看| 搡女人真爽免费视频火全软件 | 婷婷精品国产亚洲av| 一二三四社区在线视频社区8| 色视频www国产| 波野结衣二区三区在线 | 一本精品99久久精品77| 国产男靠女视频免费网站| 日日干狠狠操夜夜爽| 久久久久久久亚洲中文字幕 | 亚洲精品国产精品久久久不卡| 又爽又黄无遮挡网站| 国产亚洲精品久久久久久毛片| 久久久久久大精品| 蜜桃久久精品国产亚洲av| 不卡一级毛片| 制服丝袜大香蕉在线| 国内久久婷婷六月综合欲色啪| 亚洲中文字幕一区二区三区有码在线看| 天天添夜夜摸| 看免费av毛片| 久久久成人免费电影| 在线观看日韩欧美| 九九在线视频观看精品| 一个人免费在线观看的高清视频| 观看免费一级毛片| 日本五十路高清| 一个人观看的视频www高清免费观看| 亚洲av成人精品一区久久| 亚洲精品在线观看二区| 美女cb高潮喷水在线观看| 99久久99久久久精品蜜桃| 久久99热这里只有精品18| 可以在线观看的亚洲视频| 天天一区二区日本电影三级| 亚洲无线在线观看| 日韩欧美国产在线观看| 国产午夜精品久久久久久一区二区三区 | 日本三级黄在线观看| 每晚都被弄得嗷嗷叫到高潮| 欧美黑人巨大hd| 看片在线看免费视频| eeuss影院久久| 在线观看日韩欧美| 九九在线视频观看精品| 91久久精品国产一区二区成人 | 久久99热这里只有精品18| 三级毛片av免费| 老熟妇乱子伦视频在线观看| 最近视频中文字幕2019在线8| 欧美国产日韩亚洲一区| 最好的美女福利视频网| 99精品久久久久人妻精品| 成人国产综合亚洲| 最近视频中文字幕2019在线8| 中文字幕人妻熟人妻熟丝袜美 | 亚洲精品乱码久久久v下载方式 | 真实男女啪啪啪动态图| 国产av在哪里看| 国产精华一区二区三区| 少妇高潮的动态图| 黄色视频,在线免费观看| 观看免费一级毛片| 国产精品av视频在线免费观看| 非洲黑人性xxxx精品又粗又长| 国产精品98久久久久久宅男小说| 欧美一区二区亚洲| 美女被艹到高潮喷水动态| 国产不卡一卡二| 久99久视频精品免费| 激情在线观看视频在线高清| 亚洲性夜色夜夜综合| 高清在线国产一区| 3wmmmm亚洲av在线观看| 可以在线观看的亚洲视频| 可以在线观看毛片的网站| 成人三级黄色视频| 欧美日韩一级在线毛片| 国产亚洲精品久久久久久毛片| 观看美女的网站| 欧美最黄视频在线播放免费| 国产精品久久久人人做人人爽| 夜夜夜夜夜久久久久| 亚洲 欧美 日韩 在线 免费| 18美女黄网站色大片免费观看| 9191精品国产免费久久| 男女床上黄色一级片免费看| 久久伊人香网站| 午夜免费成人在线视频| 久久午夜亚洲精品久久| 久久香蕉精品热| 在线观看美女被高潮喷水网站 | 麻豆一二三区av精品| av欧美777| 国产视频内射| 国产精品影院久久| 精品欧美国产一区二区三| 国产乱人伦免费视频| av在线天堂中文字幕| 淫妇啪啪啪对白视频| 18禁黄网站禁片免费观看直播| 久久久国产成人免费| 老司机在亚洲福利影院| av视频在线观看入口| 看黄色毛片网站| 国产精品,欧美在线| 色综合欧美亚洲国产小说| 日韩人妻高清精品专区| 日韩亚洲欧美综合| 国产精品乱码一区二三区的特点| 天天添夜夜摸| 人人妻人人看人人澡| 黄片小视频在线播放| 欧美成人免费av一区二区三区| ponron亚洲| 男插女下体视频免费在线播放| 伊人久久精品亚洲午夜| 制服丝袜大香蕉在线| 一区二区三区国产精品乱码| 成年女人永久免费观看视频| 精品国内亚洲2022精品成人| 国产欧美日韩一区二区三| 国产单亲对白刺激| 婷婷精品国产亚洲av在线| 十八禁人妻一区二区| 久久精品国产亚洲av香蕉五月| 精品无人区乱码1区二区| 久久久久久久亚洲中文字幕 | 淫妇啪啪啪对白视频| 神马国产精品三级电影在线观看| 成人午夜高清在线视频| 在线播放国产精品三级| 一区二区三区高清视频在线| 乱人视频在线观看| 人妻丰满熟妇av一区二区三区| 日本 av在线| 免费看日本二区| 久久精品人妻少妇| 亚洲人成电影免费在线| 悠悠久久av| 一本久久中文字幕| 国语自产精品视频在线第100页| 午夜亚洲福利在线播放| 国产伦在线观看视频一区| 欧美性猛交黑人性爽| 欧美最黄视频在线播放免费| 19禁男女啪啪无遮挡网站| 国产aⅴ精品一区二区三区波| 3wmmmm亚洲av在线观看| 国产v大片淫在线免费观看| 婷婷精品国产亚洲av在线| 在线观看免费午夜福利视频| 一进一出抽搐gif免费好疼| 国产亚洲精品久久久com| 婷婷精品国产亚洲av| 亚洲av熟女| 国产探花极品一区二区| 午夜亚洲福利在线播放| 老司机深夜福利视频在线观看| 午夜福利免费观看在线| 中文字幕av在线有码专区| 国产高清视频在线观看网站| 一级毛片女人18水好多| 日本精品一区二区三区蜜桃| 真实男女啪啪啪动态图| 男女床上黄色一级片免费看| 欧美一级a爱片免费观看看| 高清日韩中文字幕在线| 男女那种视频在线观看| 精品乱码久久久久久99久播| 亚洲精品乱码久久久v下载方式 | 看片在线看免费视频| 黄片大片在线免费观看| 亚洲一区二区三区不卡视频| 亚洲性夜色夜夜综合| 人人妻人人看人人澡| 久久香蕉国产精品| 真人做人爱边吃奶动态| 床上黄色一级片| 国产探花在线观看一区二区| 变态另类丝袜制服| 美女被艹到高潮喷水动态| 天堂√8在线中文| 国产在线精品亚洲第一网站| 午夜免费观看网址| 国产精品一区二区免费欧美| 亚洲熟妇中文字幕五十中出| 国产激情欧美一区二区| 久久久久久大精品| 国产午夜精品论理片| 18禁在线播放成人免费| 午夜视频国产福利| 天天添夜夜摸| 特级一级黄色大片| 丝袜美腿在线中文| 老熟妇乱子伦视频在线观看| 国产精品久久视频播放| 日韩中文字幕欧美一区二区| 国产亚洲精品久久久com| 动漫黄色视频在线观看| 一本一本综合久久| 亚洲专区中文字幕在线| 免费在线观看成人毛片| 色在线成人网| 亚洲人成网站在线播| 亚洲片人在线观看| 国产精品 欧美亚洲| 久久这里只有精品中国| 高清在线国产一区| 岛国视频午夜一区免费看| 精品一区二区三区av网在线观看| 长腿黑丝高跟| 午夜日韩欧美国产| 一个人免费在线观看电影| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美一区二区三区黑人| 黄色日韩在线| www.www免费av| 又紧又爽又黄一区二区| 高清毛片免费观看视频网站| 每晚都被弄得嗷嗷叫到高潮| 国产美女午夜福利| 亚洲自拍偷在线| 亚洲精品456在线播放app | 我的老师免费观看完整版| 别揉我奶头~嗯~啊~动态视频| 国产精品综合久久久久久久免费| 老司机午夜十八禁免费视频| 一级a爱片免费观看的视频| 亚洲第一欧美日韩一区二区三区| 日本 av在线| 精品久久久久久久久久久久久| 国产精品久久久人人做人人爽| 亚洲人成网站在线播| 久久久成人免费电影| 国产av一区在线观看免费| 国产91精品成人一区二区三区| 观看免费一级毛片| 久久精品影院6| av专区在线播放| 国产精品爽爽va在线观看网站| 人人妻人人看人人澡| 18禁在线播放成人免费| 国产精品一区二区三区四区免费观看 | www.熟女人妻精品国产| 天堂动漫精品| 久久国产精品影院| 久久久久国内视频| 国产高清视频在线播放一区| 国产亚洲欧美98| 国产真人三级小视频在线观看| 国产精品 欧美亚洲| 国产私拍福利视频在线观看| 国产精品久久久久久精品电影| 欧美极品一区二区三区四区| www.色视频.com| 日本一本二区三区精品| 有码 亚洲区| 欧美高清成人免费视频www| 制服丝袜大香蕉在线| 99久久九九国产精品国产免费| 久久天躁狠狠躁夜夜2o2o| 三级男女做爰猛烈吃奶摸视频| 亚洲成av人片免费观看| 国产精品影院久久| 99久国产av精品| 男女下面进入的视频免费午夜| 99久久成人亚洲精品观看| 在线免费观看的www视频| 最新美女视频免费是黄的| 亚洲18禁久久av| 中文字幕av在线有码专区| 免费观看人在逋| 国语自产精品视频在线第100页| 国产精华一区二区三区| 亚洲av成人不卡在线观看播放网| 狂野欧美激情性xxxx| 一a级毛片在线观看| 欧美大码av| 国产欧美日韩精品亚洲av| 99久久99久久久精品蜜桃| 91在线精品国自产拍蜜月 | 麻豆国产av国片精品| 好男人电影高清在线观看| 国产单亲对白刺激| 国产精品一及| 国产av麻豆久久久久久久| 国产一区二区三区在线臀色熟女| 亚洲最大成人中文| 午夜福利视频1000在线观看| 日韩欧美国产一区二区入口| 一卡2卡三卡四卡精品乱码亚洲| 免费电影在线观看免费观看| av专区在线播放| 亚洲av一区综合| 亚洲av五月六月丁香网| 亚洲欧美日韩高清专用| 国产野战对白在线观看| 国产av一区在线观看免费| 亚洲内射少妇av| 搡老熟女国产l中国老女人| 久久亚洲真实| 两个人的视频大全免费| www国产在线视频色| 亚洲美女视频黄频| 国产色爽女视频免费观看| 90打野战视频偷拍视频| 岛国在线观看网站| 久久性视频一级片| 欧美日本亚洲视频在线播放| 国产在视频线在精品| 床上黄色一级片| 内射极品少妇av片p| 午夜老司机福利剧场| 精品人妻偷拍中文字幕| 欧美极品一区二区三区四区| 亚洲精品美女久久久久99蜜臀| 国产熟女xx| 色精品久久人妻99蜜桃| 亚洲专区中文字幕在线| 熟女少妇亚洲综合色aaa.| 国产精品一区二区三区四区久久| 国产亚洲精品av在线| 亚洲国产精品成人综合色| 成人特级黄色片久久久久久久| x7x7x7水蜜桃| 国产高清激情床上av| 天堂网av新在线| 欧美一区二区亚洲| 中文亚洲av片在线观看爽| 国产美女午夜福利| av欧美777| 脱女人内裤的视频| 两个人视频免费观看高清| 国产极品精品免费视频能看的| 中文字幕精品亚洲无线码一区| 精华霜和精华液先用哪个| 国产精品一区二区免费欧美| 超碰av人人做人人爽久久 | 69人妻影院| 精品国产亚洲在线| 国产97色在线日韩免费| 丰满人妻一区二区三区视频av | 天堂√8在线中文| 伊人久久精品亚洲午夜| tocl精华| 在线视频色国产色| 天堂网av新在线| 午夜亚洲福利在线播放| 国产一级毛片七仙女欲春2| 国产69精品久久久久777片| 国产精品一及| 香蕉久久夜色| a在线观看视频网站| 人人妻,人人澡人人爽秒播| 亚洲精品乱码久久久v下载方式 | 手机成人av网站| 亚洲国产色片| 最新中文字幕久久久久| a级一级毛片免费在线观看| 精品免费久久久久久久清纯| av在线蜜桃| 色综合欧美亚洲国产小说| АⅤ资源中文在线天堂| 757午夜福利合集在线观看| 麻豆国产av国片精品| 国产真人三级小视频在线观看| 国产探花在线观看一区二区| 99国产精品一区二区蜜桃av| 日韩免费av在线播放| 色av中文字幕| 精品乱码久久久久久99久播| 欧美+日韩+精品| 成年人黄色毛片网站| 看黄色毛片网站| 亚洲av免费在线观看| 亚洲人成网站高清观看| 色吧在线观看| 成人18禁在线播放| 亚洲美女视频黄频| 最新美女视频免费是黄的| 天堂√8在线中文| 麻豆一二三区av精品| av在线天堂中文字幕| 内射极品少妇av片p| 欧美在线一区亚洲| 老司机在亚洲福利影院| www.www免费av| 国产免费一级a男人的天堂| 黄色片一级片一级黄色片| 亚洲真实伦在线观看| 最新在线观看一区二区三区| 亚洲精品成人久久久久久| 亚洲精品乱码久久久v下载方式 | 成人亚洲精品av一区二区| 国产蜜桃级精品一区二区三区| 亚洲人成网站高清观看| 天堂√8在线中文| 欧美日本视频| 两个人视频免费观看高清| 久久国产精品影院| 久久午夜亚洲精品久久| 又粗又爽又猛毛片免费看| 亚洲精品456在线播放app | 免费av毛片视频| 亚洲成av人片在线播放无| 男人的好看免费观看在线视频| 美女大奶头视频| 啦啦啦观看免费观看视频高清| 91久久精品国产一区二区成人 | 日韩亚洲欧美综合| av欧美777| 亚洲成av人片在线播放无| 日本一二三区视频观看| 少妇人妻一区二区三区视频| 国产欧美日韩精品一区二区| 成人无遮挡网站| 国产真实伦视频高清在线观看 | 亚洲天堂国产精品一区在线| 极品教师在线免费播放| 亚洲一区二区三区色噜噜| 内地一区二区视频在线| 淫秽高清视频在线观看| 欧美三级亚洲精品| 欧美日韩亚洲国产一区二区在线观看| 女生性感内裤真人,穿戴方法视频| 又紧又爽又黄一区二区| 又爽又黄无遮挡网站| 国产一区二区激情短视频| 日本成人三级电影网站| 蜜桃亚洲精品一区二区三区| 国产精品三级大全| 国产亚洲精品久久久久久毛片| 蜜桃亚洲精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 国产视频一区二区在线看| 欧美一级毛片孕妇| 淫秽高清视频在线观看| 亚洲精品在线美女| 成熟少妇高潮喷水视频| 久久久精品欧美日韩精品| 一区二区三区激情视频|