本發(fā)明屬于多能源耦合的系統(tǒng)建模領(lǐng)域，特別涉及一種計(jì)及管網(wǎng)儲(chǔ)熱效益的集中供熱系統(tǒng)建模方法。
背景技術(shù)：
：供熱網(wǎng)絡(luò)的管網(wǎng)具有良好的儲(chǔ)熱特性，即對(duì)供熱管網(wǎng)的熱量改變具有很強(qiáng)的延時(shí)性，這在宏觀時(shí)間上體現(xiàn)為熱量?jī)?chǔ)存特性。在我國(guó)北方集中供熱地區(qū)，冬季供暖的主要方式為大型熱電聯(lián)產(chǎn)集中供熱。目前，我國(guó)北方的熱電聯(lián)產(chǎn)機(jī)組主要采用“以熱定電”的模型運(yùn)行，即熱電聯(lián)產(chǎn)機(jī)組的發(fā)電出力取決于供熱負(fù)荷，這制約了熱-電的出力，體現(xiàn)出了熱-電強(qiáng)制約關(guān)系。在“以熱定電”運(yùn)行模式下，熱電聯(lián)產(chǎn)機(jī)組受到最小供電出力的限制來(lái)保證供熱量滿足熱負(fù)荷需求。在冬季供熱期，為了保證滿足供熱需求，大量的熱電聯(lián)產(chǎn)機(jī)組運(yùn)行于最小技術(shù)出力狀態(tài)。受“以熱定電”模式的限制，熱電聯(lián)產(chǎn)機(jī)組的調(diào)峰能力十分有限，導(dǎo)致系統(tǒng)向下旋轉(zhuǎn)備用容量不足，難以為消納風(fēng)電提供向下調(diào)節(jié)空間。在傳統(tǒng)的電力系統(tǒng)經(jīng)濟(jì)調(diào)度、機(jī)組組合問(wèn)題中，集中供熱系統(tǒng)并沒(méi)有被建?？紤]。供熱管網(wǎng)的溫度變化具有高延時(shí)性，即熱量在管網(wǎng)中被儲(chǔ)存。應(yīng)用供熱管網(wǎng)的儲(chǔ)熱效益，可以松弛“以熱定電”模式下的熱-電強(qiáng)制約關(guān)系。作為眾多北方城市的基礎(chǔ)公共設(shè)施，區(qū)域集中供熱網(wǎng)絡(luò)能夠?yàn)闊犭娐?lián)產(chǎn)機(jī)組的運(yùn)行提供額外的靈活性。實(shí)際規(guī)模的區(qū)域集中供熱網(wǎng)由成千上萬(wàn)絕熱管道組成，具有巨大的儲(chǔ)熱能力。為了對(duì)管網(wǎng)的儲(chǔ)能特性合理利用，對(duì)集中供熱系統(tǒng)供熱管網(wǎng)的詳細(xì)模型的研究十分必要。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的目的是是為填補(bǔ)已有技術(shù)的空白之處，提出一種計(jì)及管網(wǎng)儲(chǔ)熱效益的集中供熱系統(tǒng)建模方法。本發(fā)明方法很好地反映出熱網(wǎng)管網(wǎng)的儲(chǔ)熱特性。所建立的模型適宜應(yīng)用在熱網(wǎng)與電網(wǎng)的聯(lián)合調(diào)度決策中。本發(fā)明提出的一種計(jì)及管網(wǎng)儲(chǔ)熱效益的集中供熱系統(tǒng)建模方法，其特征在于，集中供熱系統(tǒng)模型包括熱電聯(lián)產(chǎn)機(jī)組模型，循環(huán)水泵模型，換熱站模型和供熱網(wǎng)絡(luò)模型四個(gè)子模型，每個(gè)子模型由各自的約束條件構(gòu)成；該方法具體包括以下步驟：1)構(gòu)建熱電聯(lián)產(chǎn)機(jī)組模型熱電聯(lián)產(chǎn)機(jī)組模型包含以下約束：1.1)供熱與供電關(guān)系約束熱電聯(lián)產(chǎn)機(jī)組的供電出力和供熱出力的約束采用多邊形區(qū)域極點(diǎn)的凸組合來(lái)描述，如式(1)所示：其中，式(1)中，pi,t為第i臺(tái)機(jī)組在第t個(gè)調(diào)度時(shí)段的供電出力，hi,t為第i臺(tái)熱電聯(lián)產(chǎn)機(jī)組在第t個(gè)調(diào)度時(shí)段的供熱出力，為第i臺(tái)熱電聯(lián)產(chǎn)機(jī)組運(yùn)行可行域近似多邊形的第k個(gè)極點(diǎn)，為第i臺(tái)熱電聯(lián)產(chǎn)機(jī)組在第t個(gè)調(diào)度時(shí)段的運(yùn)行點(diǎn)對(duì)應(yīng)第k個(gè)極點(diǎn)的凸組合系數(shù)，NKi為第i臺(tái)熱電聯(lián)產(chǎn)機(jī)組的運(yùn)行可行域近似多邊形的極點(diǎn)個(gè)數(shù)，為熱電聯(lián)產(chǎn)機(jī)組下標(biāo)集合，為調(diào)度時(shí)段下標(biāo)集合；1.2)循環(huán)水供熱約束熱電聯(lián)產(chǎn)機(jī)組的供熱出力用于加熱供熱網(wǎng)絡(luò)中的循環(huán)水流，如式(3)所示：其中，c為水的比熱容，為在第t個(gè)調(diào)度時(shí)段內(nèi)流經(jīng)第j個(gè)熱力站的循環(huán)水質(zhì)量，即循環(huán)水流量；分別為供熱網(wǎng)絡(luò)中第n個(gè)節(jié)點(diǎn)在第t個(gè)調(diào)度時(shí)段的供水、回水溫度，為熱力站下標(biāo)集合，為供熱網(wǎng)絡(luò)中與第j個(gè)熱力站連接的節(jié)點(diǎn)下標(biāo)；1.3)供熱網(wǎng)絡(luò)節(jié)點(diǎn)溫度約束供熱網(wǎng)絡(luò)節(jié)點(diǎn)溫度控制在合理范圍內(nèi)以保證供熱質(zhì)量并防止循環(huán)水汽化，如式(4)所示：其中，和分別為供熱網(wǎng)絡(luò)第n個(gè)節(jié)點(diǎn)的溫度下限和上限；2)構(gòu)建循環(huán)水泵模型循環(huán)水泵模型包含以下約束：2.1)循環(huán)水泵消耗的電功率約束循環(huán)水泵消耗的電功率正比于供熱網(wǎng)絡(luò)節(jié)點(diǎn)的供回水壓強(qiáng)差以及循環(huán)水流量，如式(5)所示：其中，為循環(huán)水泵消耗的電功率，分別表示供熱網(wǎng)絡(luò)第n個(gè)節(jié)點(diǎn)在第t個(gè)調(diào)度時(shí)段的供水、回水壓強(qiáng)，為循環(huán)水泵的工作效率，ρ為水密度；2.2)循環(huán)水電功率的上下限約束，如式(6)所示：其中，與分別代表循環(huán)水泵電功率的上限和下限；3)構(gòu)建換熱站模型換熱站模型包含以下約束：3.1)供回水溫度與換熱量的關(guān)系約束換熱站的供回水溫度與換熱量的關(guān)系如式(7)所示：其中，為第l個(gè)換熱站在第t個(gè)調(diào)度時(shí)段的循環(huán)水流量，為第l個(gè)換熱站在第t個(gè)調(diào)度時(shí)段的用熱負(fù)荷功率，為換熱站下標(biāo)集合，為供熱網(wǎng)絡(luò)中與第l個(gè)換熱站連接的節(jié)點(diǎn)下標(biāo)；3.2)節(jié)點(diǎn)供回水壓強(qiáng)約束換熱站處的節(jié)點(diǎn)供回水壓強(qiáng)差高于一定水平以維持循環(huán)水流，如式(8)所示：其中，為第l個(gè)換熱站的最小供回水壓強(qiáng)差；3.3)回水溫度約束換熱站的回水溫度保持在一定的溫度范圍之內(nèi)，如式(9)所示：其中，與分別代表?yè)Q熱站回水溫度的上限和下限；4)構(gòu)建供熱網(wǎng)絡(luò)模型供熱網(wǎng)絡(luò)中：表示以第i個(gè)供熱網(wǎng)絡(luò)節(jié)點(diǎn)為終點(diǎn)的供熱管道下標(biāo)集合，表示以第i個(gè)供熱網(wǎng)絡(luò)節(jié)點(diǎn)為起點(diǎn)的供熱管道下標(biāo)集合，分別表示第b條供水管道在第t個(gè)調(diào)度時(shí)段的首端、末端溫度，分別表示第b條回水管道在第t個(gè)調(diào)度時(shí)段的首端、末端溫度，分別表示第b條供水、回水管道在第t個(gè)調(diào)度時(shí)段的水流量；供熱網(wǎng)絡(luò)模型包含以下約束：4.1)流量連續(xù)性約束進(jìn)入同一節(jié)點(diǎn)的水流量之和為零，如式(10)和(11)所示：其中，為供熱網(wǎng)絡(luò)節(jié)點(diǎn)下標(biāo)集合；4.2)溫度混合約束來(lái)自不同管道的水流在同一網(wǎng)絡(luò)節(jié)點(diǎn)混合后的溫度滿足以下方程如式(12)和(13)所示：4.3)網(wǎng)絡(luò)節(jié)點(diǎn)溫度約束從網(wǎng)絡(luò)節(jié)點(diǎn)流出的水流溫度等于該網(wǎng)絡(luò)節(jié)點(diǎn)的溫度，如式(14)和(15)所示：4.4)流量限制約束循環(huán)水流量限制在一定的范圍以內(nèi)以防止管道振動(dòng)，如式(16)和(17)所示：其中，為第b條供熱管道的水流流速上限，為供熱管道下標(biāo)集合；4.4)壓強(qiáng)損失約束由于水流與管道內(nèi)壁摩擦引起的沿管壓強(qiáng)損失與流速的平方成正比，如式(18)所示：其中，μb為第b條供熱管道的壓強(qiáng)損失系數(shù)，分別為第b條供熱管道的首端、末端節(jié)點(diǎn)下標(biāo)；4.5)水溫變化延時(shí)約束與沿管熱損耗約束此約束分為兩個(gè)步驟：4.5.1)利用過(guò)去時(shí)段的管道入口溫度估計(jì)忽略沿管熱損耗的管道出口溫度，如式(19)和式(20)所示：其中，和分別為第b條供熱管道在第t個(gè)調(diào)度時(shí)段忽略沿管熱損耗的管道出口溫度；變量Kb,t,k的取值由循環(huán)水流速?zèng)Q定，如式(21)所示：Kb,t,k=(msb,tpipe·Δt-Sb,t+ρAbLb)/(msb,tpipe·Δt),k=t-φb,t(msb,tpipe·Δt)/(msb,tpipe·Δt),k=t-φb,t+1,...,t-γb,t-1(Rb,t-ρAbLb)/(msb,tpipe·Δt),k=t-γb,t0,otherwise---(21)]]>其中，Δt為相鄰調(diào)度時(shí)段的時(shí)間間隔，Ab為第b條供熱管道的截面積，Lb為第b條供熱管道的長(zhǎng)度；式(19)～式(21)中，整數(shù)變量φb,t和γb,t分別表示與水溫變化延時(shí)相關(guān)的調(diào)度時(shí)段數(shù)，如式(22)和式(23)所示：式(21)中的Rb,t和Sb,t表達(dá)式分別如式(24)和式(25)所示：Rb,t=Σk=0γb,t(msb,t-kpipe·Δt)---(24)]]>Sb,t=Σk=0φb,t-1(msb,t-kpipe·Δt),ifφb,t≥γb,t+1Rb,t,otherwise---(25)]]>4.5.2)對(duì)管道出口溫度進(jìn)行熱損耗修正，如式(26)和式(27)所示：τb,tPS,out=τtam+(τb,t′PS,out-τtam)exp[-λbΔtAbρc(γb,t+12+Sb,t-Rb,tmsb,t-γb,tpipeΔt)]---(26)]]>τb,tPR,out=τtam+(τb,t′PR,out-τtam)exp[-λbΔtAbρc(γb,t+12+Sb,t-Rb,tmsb,t-γb,tpipeΔt)]---(27)]]>其中，為第t個(gè)調(diào)度時(shí)段的外界溫度，λb為第b條供熱管道的沿管熱損耗系數(shù)；5)集中供熱系統(tǒng)模型構(gòu)建完畢.本發(fā)明的特點(diǎn)及有益效果：本發(fā)明所提出的計(jì)及管網(wǎng)儲(chǔ)熱效益的集中供熱系統(tǒng)建模方法，針對(duì)區(qū)域熱網(wǎng)特性，依次對(duì)熱電聯(lián)產(chǎn)機(jī)組、循環(huán)水泵、換熱站以及供熱網(wǎng)絡(luò)分別進(jìn)行建模。在所建立的模型中，考慮了熱網(wǎng)管網(wǎng)的時(shí)間擬動(dòng)態(tài)特性，針對(duì)集中供熱系統(tǒng)的儲(chǔ)熱特性描述準(zhǔn)確，適宜應(yīng)用在電力系統(tǒng)的運(yùn)行調(diào)度中，以充分利用集中供熱系統(tǒng)引入的調(diào)度靈活性。附圖說(shuō)明圖1為本發(fā)明方法的流程框圖。圖2(a)為背壓式熱電聯(lián)產(chǎn)機(jī)組出力運(yùn)行可行域示意圖。圖2(b)為抽凝式熱電聯(lián)產(chǎn)機(jī)組出力運(yùn)行可行域示意圖。具體實(shí)施例本發(fā)明提出的一種計(jì)及管網(wǎng)儲(chǔ)熱效益的集中供熱系統(tǒng)建模方法，下面結(jié)合附圖和具體實(shí)施例對(duì)本發(fā)明進(jìn)一步詳細(xì)說(shuō)明如下。本發(fā)明提出的一種計(jì)及管網(wǎng)儲(chǔ)熱效益的集中供熱系統(tǒng)建模方法，流程框圖如圖1所示，集中供熱系統(tǒng)模型包括熱電聯(lián)產(chǎn)機(jī)組模型，循環(huán)水泵模型，換熱站模型和供熱網(wǎng)絡(luò)模型四個(gè)子模型，每個(gè)子模型由各自的約束條件構(gòu)成；該方法具體包括以下步驟：1)構(gòu)建熱電聯(lián)產(chǎn)機(jī)組模型熱電聯(lián)產(chǎn)機(jī)組模型包含以下約束：1.1)供熱與供電關(guān)系約束熱電聯(lián)產(chǎn)機(jī)組分為背壓式和抽凝式兩種類型，兩種機(jī)組的出力運(yùn)行可行域，即供熱與供電出力關(guān)系，如圖2所示，其中圖2(a)為背壓式熱電聯(lián)產(chǎn)機(jī)組，圖2(b)為抽凝式熱電聯(lián)產(chǎn)機(jī)組。兩個(gè)圖的橫坐標(biāo)均為機(jī)組的供熱出力，縱坐標(biāo)均為機(jī)組的供電出力。其中圖2(a)的背壓式熱電聯(lián)產(chǎn)機(jī)組供電出力與供熱出力為線性關(guān)系，圖2(b)的抽凝式機(jī)組供電出力與供熱出力關(guān)系為凸四邊形可行域。熱電聯(lián)產(chǎn)機(jī)組的供電出力和供熱出力的約束采用多邊形區(qū)域極點(diǎn)的凸組合來(lái)描述，如式(1)所示：其中，式(1)中，pi,t為第i臺(tái)機(jī)組在第t個(gè)調(diào)度時(shí)段的供電出力，hi,t為第i臺(tái)熱電聯(lián)產(chǎn)機(jī)組在第t個(gè)調(diào)度時(shí)段的供熱出力，為第i臺(tái)熱電聯(lián)產(chǎn)機(jī)組運(yùn)行可行域近似多邊形的第k個(gè)極點(diǎn)，為第i臺(tái)熱電聯(lián)產(chǎn)機(jī)組在第t個(gè)調(diào)度時(shí)段的運(yùn)行點(diǎn)對(duì)應(yīng)第k個(gè)極點(diǎn)的凸組合系數(shù)，NKi為第i臺(tái)熱電聯(lián)產(chǎn)機(jī)組的運(yùn)行可行域近似多邊形的極點(diǎn)個(gè)數(shù)，為熱電聯(lián)產(chǎn)機(jī)組下標(biāo)集合，為調(diào)度時(shí)段下標(biāo)集合；1.2)循環(huán)水供熱約束熱電聯(lián)產(chǎn)機(jī)組的供熱出力用于加熱供熱網(wǎng)絡(luò)中的循環(huán)水流如式(3)所示：其中，c為水的比熱容，其具體值根據(jù)實(shí)際水質(zhì)實(shí)驗(yàn)測(cè)出，為在第t個(gè)調(diào)度時(shí)段內(nèi)流經(jīng)第j個(gè)熱力站的循環(huán)水質(zhì)量(為方便起見(jiàn)，下文稱之為循環(huán)水流量)，分別為供熱網(wǎng)絡(luò)中第n個(gè)節(jié)點(diǎn)在第t個(gè)調(diào)度時(shí)段的供水、回水溫度，為熱力站下標(biāo)集合，為供熱網(wǎng)絡(luò)中與第j個(gè)熱力站連接的節(jié)點(diǎn)下標(biāo)；1.3)供熱網(wǎng)絡(luò)節(jié)點(diǎn)溫度約束供熱網(wǎng)絡(luò)節(jié)點(diǎn)溫度控制在合理范圍內(nèi)以保證供熱質(zhì)量并防止循環(huán)水汽化，如式(4)所示：其中，和分別為供熱網(wǎng)絡(luò)第n個(gè)節(jié)點(diǎn)的溫度下限和上限；2)構(gòu)建循環(huán)水泵模型循環(huán)水泵模型包含以下約束：2.1)循環(huán)水泵消耗的電功率約束循環(huán)水泵消耗的電功率正比于供熱網(wǎng)絡(luò)節(jié)點(diǎn)的供回水壓強(qiáng)差以及循環(huán)水流量，如式(5)所示：其中，為循環(huán)水泵消耗的電功率，分別表示供熱網(wǎng)絡(luò)第n個(gè)節(jié)點(diǎn)在第t個(gè)調(diào)度時(shí)段的供水、回水壓強(qiáng)，為循環(huán)水泵的工作效率，ρ為水密度；2.2)循環(huán)水電功率的上下限約束，如式(6)所示：其中，與分別代表循環(huán)水泵電功率的上限和下限；3)構(gòu)建換熱站模型換熱站模型包含以下約束：3.1)供回水溫度與換熱量的關(guān)系約束在輸熱系統(tǒng)中，換熱站可以被視為一個(gè)用熱負(fù)荷；換熱站的供回水溫度與換熱量的關(guān)系如式(7)所示：其中，為第l個(gè)換熱站在第t個(gè)調(diào)度時(shí)段的循環(huán)水流量，為第l個(gè)換熱站在第t個(gè)調(diào)度時(shí)段的用熱負(fù)荷功率，為換熱站下標(biāo)集合，為供熱網(wǎng)絡(luò)中與第l個(gè)換熱站連接的節(jié)點(diǎn)下標(biāo)；3.2)節(jié)點(diǎn)供回水壓強(qiáng)約束換熱站處的節(jié)點(diǎn)供回水壓強(qiáng)差高于一定水平以維持循環(huán)水流，如式(8)所示：其中，為第l個(gè)換熱站的最小供回水壓強(qiáng)差；3.3)回水溫度約束換熱站的回水溫度保持在一定的溫度范圍之內(nèi)，如式(9)所示：其中，與分別代表?yè)Q熱站回水溫度的上限和下限；4)構(gòu)建供熱網(wǎng)絡(luò)模型定義供熱網(wǎng)絡(luò)中的如下變量：表示以第i個(gè)供熱網(wǎng)絡(luò)節(jié)點(diǎn)為終點(diǎn)的供熱管道下標(biāo)集合，表示以第i個(gè)供熱網(wǎng)絡(luò)節(jié)點(diǎn)為起點(diǎn)的供熱管道下標(biāo)集合，分別表示第b條供水管道在第t個(gè)調(diào)度時(shí)段的首端、末端溫度，分別表示第b條回水管道在第t個(gè)調(diào)度時(shí)段的首端、末端溫度，分別表示第b條供水、回水管道在第t個(gè)調(diào)度時(shí)段的水流量。值得注意的是，供熱網(wǎng)絡(luò)中第i個(gè)節(jié)點(diǎn)在第t個(gè)調(diào)度時(shí)段的供水、回水溫度指的是流進(jìn)該節(jié)點(diǎn)的水流相互混合后的穩(wěn)態(tài)溫度，供水和回水管道在第t個(gè)調(diào)度時(shí)段的首末端溫度指的是在管道內(nèi)的相應(yīng)位置上混合前的水流溫度。供熱網(wǎng)絡(luò)模型包含以下約束：4.1)流量連續(xù)性約束根據(jù)質(zhì)量守恒定律，由于水是不可壓縮流體，進(jìn)入同一節(jié)點(diǎn)的水流量之和為零，如式(10)和(11)所示：其中，為供熱網(wǎng)絡(luò)節(jié)點(diǎn)下標(biāo)集合；4.2)溫度混合約束根據(jù)能量守恒定律，來(lái)自不同管道的水流在同一網(wǎng)絡(luò)節(jié)點(diǎn)混合后的溫度滿足以下方程如式(12)和(13)所示：4.3)網(wǎng)絡(luò)節(jié)點(diǎn)溫度約束從網(wǎng)絡(luò)節(jié)點(diǎn)流出的水流溫度等于該網(wǎng)絡(luò)節(jié)點(diǎn)的溫度，如式(14)和(15)所示：4.4)流量限制約束循環(huán)水流量限制在一定的范圍以內(nèi)以防止管道振動(dòng)，如式(16)和(17)所示：其中，為第b條供熱管道的水流流速上限，為供熱管道下標(biāo)集合；4.4)壓強(qiáng)損失約束根據(jù)達(dá)西-魏斯巴赫公式，由于水流與管道內(nèi)壁摩擦引起的沿管壓強(qiáng)損失與流速的平方成正比，如式(18)所示：其中，μb為第b條供熱管道的壓強(qiáng)損失系數(shù)，分別為第b條供熱管道的首端、末端節(jié)點(diǎn)下標(biāo)；4.5)水溫變化延時(shí)約束與沿管熱損耗約束采用節(jié)點(diǎn)法進(jìn)行描述，此約束分為兩個(gè)步驟：4.5.1)利用過(guò)去時(shí)段的管道入口溫度估計(jì)忽略沿管熱損耗的管道出口溫度，如式(19)和式(20)所示：其中，和分別為第b條供熱管道在第t個(gè)調(diào)度時(shí)段忽略沿管熱損耗的管道出口溫度；變量Kb,t,k的取值由循環(huán)水流速?zèng)Q定，如式(21)所示：Kb,t,k=(msb,tpipe·Δt-Sb,t+ρAbLb)/(msb,tpipe·Δt),k=t-φb,t(msb,tpipe·Δt)/(msb,tpipe·Δt),k=t-φb,t+1,...,t-γb,t-1(Rb,t-ρAbLb)/(msb,tpipe·Δt),k=t-γb,t0,otherwise---(21)]]>其中，Δt為相鄰調(diào)度時(shí)段的時(shí)間間隔，Ab為第b條供熱管道的截面積，Lb為第b條供熱管道的長(zhǎng)度；式(19)～式(21)中，整數(shù)變量φb,t和γb,t分別表示與水溫變化延時(shí)相關(guān)的調(diào)度時(shí)段數(shù)，如式(22)和式(23)所示：式(21)中的Rb,t和Sb,t表達(dá)式分別如式(24)和式(25)所示：Rb,t=Σk=0γb,t(msb,t-kpipe·Δt)---(24)]]>Sb,t=Σk=0φb,t-1(msb,t-kpipe·Δt),ifφb,t≥γb,t+1Rb,t,otherwise---(25)]]>4.5.2)對(duì)管道出口溫度進(jìn)行熱損耗修正，如式(26)和式(27)所示：τb,tPS,out=τtam+(τb,t′PS,out-τtam)exp[-λbΔtAbρc(γb,t+12+Sb,t-Rb,tmsb,t-γb,tpipeΔt)]---(26)]]>τb,tPR,out=τtam+(τb,t′PR,out-τtam)exp[-λbΔtAbρc(γb,t+12+Sb,t-Rb,tmsb,t-γb,tpipeΔt)]---(27)]]>其中，為第t個(gè)調(diào)度時(shí)段的外界溫度，λb為第b條供熱管道的沿管熱損耗系數(shù)。5)集中供熱系統(tǒng)模型構(gòu)建完畢。當(dāng)前第1頁(yè)1 2 3