油井維護(hù)的設(shè)置調(diào)度與原油運輸?shù)木W(wǎng)絡(luò)設(shè)計（數(shù)學(xué)建模）范文22

1、 油井維護(hù)的設(shè)置調(diào)度與原油運輸?shù)木W(wǎng)絡(luò)設(shè)計摘要 本文對某油田設(shè)置維護(hù)班組和原油運輸?shù)南嚓P(guān)問題，應(yīng)用Floyd算法、貪心法、01整數(shù)規(guī)劃等算法，針對不同情況分別建立了相應(yīng)數(shù)學(xué)模型，給出了合理的解決方案。首先根據(jù)某油田A區(qū)的交通道路和現(xiàn)有的20個維護(hù)班組的設(shè)置情況示意圖及相關(guān)的數(shù)據(jù)信息，整理并運用Floyd算法在matlab軟件下求出20個維護(hù)班組到各采油點的最短距離，再運用貪心法思想得出分配方案。然后由于實際中一個班組最多維修一個采油點，一個采油點也僅需一個班組維修，所以我們采用0-1整數(shù)規(guī)劃，利用lingo強大的求解器功能求得分配給各漏油點的維護(hù)班組。為了設(shè)計一個合理的公路運輸網(wǎng)絡(luò)，使得修路和運

2、輸本錢最小并將各采油點的原油快速運輸?shù)絆點坐標(biāo)原點，經(jīng)分析，為使總費用降低，我們應(yīng)該確立一個修建直接通向O運輸?shù)缆返挠吞?，然后以這個油田為中心，向四周慢慢擴散，我們采用窮舉法從92個油田中選取一個適宜的油田修建直接通向O點的道路從而使費用最優(yōu)?？紤]到管道運輸比擬道路運輸較為經(jīng)濟，故假設(shè)全部使用管道運輸。而采取管道運輸只存在修建費用，故總費用與總路程成正比，于是把問題化簡為求最短路程。故我們采用圖論中的最小生成樹算法，計算最短路程。最后對文中所建模型進(jìn)行了深刻探討，并加以評價，文中所列某些方法雖不能得到最優(yōu)解，但是是經(jīng)過多方面考慮顧全大局所的最終結(jié)果。關(guān)鍵詞：Floyd算法，貪心法，0-1整數(shù)規(guī)

3、劃，窮舉法，最小生成樹1問題重述 試就某油田設(shè)置維護(hù)班組和原油運輸?shù)南嚓P(guān)情況建立數(shù)學(xué)模型分析研究下面問題： 1附件1中的附圖給出了某油田A區(qū)的交通道路和現(xiàn)有的20個維護(hù)班組的設(shè)置情況示意圖，相關(guān)的數(shù)據(jù)信息見附件2。維護(hù)班組每天的工作是巡視其管轄范圍內(nèi)的采油井，請為20個維護(hù)班組合理分配管轄范圍，使各維護(hù)班組的工作量盡量均衡，且在其所管轄油井發(fā)生事故時能盡快到達(dá)。 2假設(shè)13個采油點同時發(fā)生漏油事故，需要調(diào)度20個維護(hù)班組實現(xiàn)快速維修，恢復(fù)生產(chǎn)。實際中一個班組最多維修一個采油點，一個采油點也僅需一個班組維修，為了盡量減少漏油損失，請給出20個維護(hù)班組的合理調(diào)度方案。注：假設(shè)采油點的原油開采速度與

4、可采儲量成正比 3為了將各采油點的原油快速運輸?shù)絆點坐標(biāo)原點，方案將局部現(xiàn)有道路修建為高等級公路，作為原油的運輸專線。假設(shè)公路修建及采油期內(nèi)的養(yǎng)護(hù)本錢合計為1百萬/公里，公路運輸本錢為0.1百萬/萬噸公里，請設(shè)計一個合理的公路運輸網(wǎng)絡(luò)，使得修路和運輸本錢最小。再假設(shè)經(jīng)過勘查，各采油點的可采儲量增長了10倍，請重新設(shè)計一個合理的公路運輸網(wǎng)絡(luò)。注：假設(shè)各采油點均可修建直線到達(dá)O點的公路 4管道輸油具有建設(shè)本錢高、輸油本錢低的特點。假設(shè)輸油管道的修建及采油期內(nèi)的養(yǎng)護(hù)本錢合計為1千萬/公里，輸油本錢忽略不計，且管道只能沿著現(xiàn)有道路修建。請比照公路運輸分析修建輸油管道的經(jīng)濟性，設(shè)計一個合理的原油運輸方案

5、。再假設(shè)經(jīng)過勘查，各采油點的可采儲量增長了10倍，請重新設(shè)計一個合理的原油運輸方案。注：假設(shè)各采油點均可修建直線到達(dá)O點的輸油管道2、模型假設(shè)1、假設(shè)運輸過程是理想化的，不考慮任何意外情況。2、假設(shè)兩節(jié)點的路線是直線。符號說明為1時表示第i個維護(hù)點維護(hù)第j個漏油點，為0時反之表示第i個維護(hù)站到第j個漏油點的距離所有的維護(hù)站到分配的漏油點的距離和4、模型分析、建立與求解 1為了使各維護(hù)班組的工作量盡量均衡，且在其所管轄油井發(fā)生事故時能盡快到達(dá)，根據(jù)假設(shè)，我們選擇了經(jīng)典的Floyd算法附錄1，運用matlab軟件實現(xiàn)得出各維護(hù)站依次到各采油點的最短距離表見附件1，代碼見附錄2，然后利用貪心法思想用

6、matlab編程代碼見附錄3實現(xiàn)依次求得20個維護(hù)站到92個采油井的分配方案如下表 按上述分配方式所得各維護(hù)站的工作量如下：維護(hù)站序號工作量149.3358226284584253.055197326951163122.74657190046571449.88174974576336546.03045655454906654.73875750698796756.89234232643393866.35530147660691983.7670213043798810150.569137541325211147.0021719149844412192.2526163255014613246.384

7、8783616538214178.723956589628951577.200061641378541651.3233173060420317240.46693111852281886.706493981951671964.093585281005322056.914000469498745由上表所得結(jié)果運用matlab畫出如下所示圖畫圖程序代碼見附錄3：為了區(qū)分20個維護(hù)站所管轄的范圍，我們分別不同的顏色和形狀的組合來表示，同一顏色和形狀代表同一管轄范圍。綜上可得結(jié)果如下：維護(hù)站1: 1,68,69,73,74,75,78維護(hù)站2: 2,43,70,71，72維護(hù)站3: 3,44,55,65

8、,91維護(hù)站4: 4,57,60,62,63維護(hù)站5: 5,49,50,53,56維護(hù)站6: 6,47,58,59維護(hù)站7： 7,30,32,51維護(hù)站8: 8,33,46,54維護(hù)站9： 9，34,35，45,52維護(hù)站10: 10,26,61維護(hù)站11: 11,27,48.維護(hù)站12： 12，25,39維護(hù)站13: 13,22,23,24，64維護(hù)站14: 14,21，66維護(hù)站15: 15,28,31維護(hù)站16: 16,36,37,38維護(hù)站17: 17,29,40,41,42維護(hù)站18: 18，81,82,83,90,92維護(hù)站19: 19，67,76，77,79,80維護(hù)站20: 2

9、0,84,85,86,87,88,89 2由于一個班組最多維修一個采油點，一個采油點也僅需一個班組維修，為了盡量減少漏油損失，實現(xiàn)快速維修盡快恢復(fù)生產(chǎn)，所以維護(hù)站到達(dá)漏油點的距離應(yīng)該越短越好，故我們延續(xù)Floyd算法思想得出20個維護(hù)站分別到13個漏油點的最小距離代碼見附錄4，結(jié)合0-1整數(shù)規(guī)劃模型，利用lingo軟件計算出最優(yōu)的分配方案，使得目標(biāo)函數(shù)到達(dá)最小。設(shè) 其中 并假設(shè)s為所有維護(hù)站到分配的漏油點的距離和 維護(hù)站對應(yīng)維護(hù)漏油點如下表：維護(hù)站序號對應(yīng)漏油點序號2114135961274810105117121136142158163為了設(shè)計一個合理的公路運輸網(wǎng)絡(luò)，使得修路和運輸本錢最小并

10、將各采油點的原油快速運輸?shù)絆點坐標(biāo)原點，我們通過計算得出如下結(jié)論：1.公路的修路費用占據(jù)了主要的費用，故為使總費用降低，應(yīng)該讓絕大多數(shù)油田的運輸?shù)缆范急M可能多的沿用其他運輸?shù)缆?，這樣就可以減少修路費用2.通過程序的測試，修建直接通向O點的運輸?shù)缆窞?條時修路費用最低。其原因經(jīng)過分析為：O點距離其他采油點過遠(yuǎn)，修建多條直接通向O點的運輸?shù)缆窌蟠笤黾涌傎M用。因此，根據(jù)這兩條結(jié)論，為使總費用降低，我們應(yīng)該確立一個修建直接通向O運輸?shù)缆返挠吞铮缓笠赃@個油田為中心，向四周慢慢擴散，這樣就可以使得絕大多數(shù)油田的運輸?shù)缆范急M可能多的沿用其他運輸?shù)缆罚瑥亩菇Y(jié)果接近最優(yōu)值。我們采用窮舉法從92個油田中選取

11、一個適宜的油田修建直接通向O點的道路從而使費用最優(yōu)，并通過程序(見附錄3)計算得出應(yīng)采用編號為11的油田。并作出運輸?shù)缆沸藿ňW(wǎng)路如下列圖：圖中各線段兩端點對應(yīng)的采油點標(biāo)號依次為：0，1111,2211,2611,25)(22,13)(22,21)(26,27)(26,10)(25,12)(25,24)(13,23)(21,14)(10,34)(14,16)(34,9)(34,31)(34,33)(16,38)(16,36)(9,35)(9,8)(31,32)(31,15)(38,19)(38,41)(36,37)(35,45)(8,47)(8,46)(32,7)(15,28)(39,40)(3

12、9,4)(41,17)(41,92)(45,3)(47,48)(47,6)(47,5)(46,55)(7,30)(28,29)(40,2)(4,63)(4,57)(4,62)(17,42)(17,81)(92,87)(92,91)(3,65)(3,44)(48,61)(6,59)(5,49)(5,50)(55,54)(2,43)(2,70)(63,64)(57,58)(57,60,)(57,56)(62,85)(81,82)(81,18)(87,88)(87,86)(91,90)(65,66)(44,67)(59,51)(49,53)(43,72)(70,69)(64,76)(56,52)(8

13、5,20)(85,84)(82,83)(18,73)(18,80)(88,89)(67,68)(72,71)(69,1)(76,77)(76,75)(73,74)(80,79)(1,78)(77,19)按上述路徑計算其所需的最小總費用為176659.1512。當(dāng)各采油點的可采儲量增長了10倍后，依然用上述方法，但我們發(fā)現(xiàn)雖然可采儲量增長十倍，但最終路徑與未增長前相同，原因為道路運輸修路費用過高，故所得最終結(jié)果差異不大，為179376.895。4考慮到管道運輸比擬道路運輸較為經(jīng)濟，故假設(shè)全部使用管道運輸。而采取管道運輸只存在修建費用，故總費用與總路程成正比，于是把問題化簡求最短路程。我們將140

14、條道路和各采油點到O點的直接道路看做一個連通圖，并用最小生成樹算法求最短路徑，而所求路徑就是使整體路程最小的路徑，從而得到使總體費用最低的修建方案。經(jīng)程序計算，該運輸網(wǎng)路如下：5、模型評價模型一的求解首先運用了Floyd算法，容易理解，代碼編寫簡單，求得20個維護(hù)點到各采油井的最短距離問題，然后再用貪心法合理分配出20個維護(hù)班組的管轄范圍，為使各維護(hù)班組的工作量盡量均衡，可能最后結(jié)果不是最優(yōu)解。3模型二的求解運用了0-1整數(shù)規(guī)劃，得出20個維護(hù)班組的合理調(diào)度方案。模型較為簡單，不適宜數(shù)據(jù)過多的情況。3模型三的求解是根據(jù)此題情況分析，其中局部引用窮舉法，得到一個近似最優(yōu)解，雖有些復(fù)雜但易于理解。

15、4模型四在求解過程中雖然通過假設(shè)簡化了問題，但是也造成結(jié)果與實際問題存在偏差，故所求結(jié)果較為片面。6、參考文獻(xiàn)1王紅梅，?算法設(shè)計與分析?，清華大學(xué)出版社2楊桂元，黃己立，?數(shù)學(xué)建模?，中國科學(xué)技術(shù)大學(xué)出版社，3嚴(yán)蔚敏，?數(shù)據(jù)結(jié)構(gòu)?，清華大學(xué)出版社4劉浩，?MATLAB一本通?,電子工業(yè)出版社附錄附錄1Floyd算法：functionD,R=floyd(a)n=size(a,1);D=a;for i=1:n for j=1:n R(i,j)=j; endendfor k=1:n for i=1:n for j=1:n if D(i,k)+D(k,j)D(i,j) D(i,j)=D(i,k)+D

16、(k,j); R(i,j)=R(i,k); end end endendend分配函數(shù)：function C,x=distribution(B)%B矩陣的行代表維護(hù)班組，列代表油田，矩陣中元素代表維護(hù)班組到油田的最小距離，C矩陣中行代表維護(hù)班組，矩陣中元素所在行代表哪些油田歸屬于哪些維護(hù)班組h=size(B,1);l=size(B,2);C=zeros(h,l);%C矩陣初始化y=ones(h);x=zeros(1,h);%h代表維護(hù)班組,x代表各個班組的工作量for ii=1:l a,b=min(x); c,d=min(B(b,:); x(b)=x(b)+c; B(:,d)=inf; C(b

17、,y(b)=d; y(b)=y(b)+1;end各采油點間直接距離 C=zeros(92,92);%C矩陣初始化for ii=1:140; C(B(ii,1),B(ii,2)=sqrt(A(B(ii,1),2)-A(B(ii,2),2)2+(A(B(ii,1),3)-A(B(ii,2),3)2);%A矩陣?yán)锎鎯α瞬捎忘c的標(biāo)號，橫縱坐標(biāo)和可開采量，B矩陣?yán)镌卮砹藗€采油點之間的道路 C(B(ii,2),B(ii,1)=sqrt(A(B(ii,1),2)-A(B(ii,2),2)2+(A(B(ii,1),3)-A(B(ii,2),3)2);endC(C+eye(92)=0)=inf%C矩陣為各采

18、油點之間存在的直線距離畫圖for ii=1:20 plot(A(ii,2),A(ii,3),o); hold on;endL=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20;for ii=1:20 text(A(ii,2)+0.2,A(ii,3)+0.2,Lii);endj=1;while F(1,j)=0 plot(A(F(1,j),2),A(F(1,j),3),*r); hold on; j=j+1;endj=1;while F(2,j)=0 plot(A(F(2,j),2),A(F(2,j),3),*g); hold on; j=j+

19、1;endj=1;while F(3,j)=0 plot(A(F(3,j),2),A(F(3,j),3),*b); hold on; j=j+1;endj=1;while F(4,j)=0 plot(A(F(4,j),2),A(F(4,j),3),*c); hold on; j=j+1;endj=1;while F(5,j)=0 plot(A(F(5,j),2),A(F(5,j),3),*m); hold on; j=j+1;endj=1;while F(6,j)=0 plot(A(F(6,j),2),A(F(6,j),3),*y); hold on; j=j+1;endj=1;while F

20、(7,j)=0 plot(A(F(7,j),2),A(F(7,j),3),*k); hold on; j=j+1;endj=1;while F(8,j)=0 plot(A(F(8,j),2),A(F(8,j),3),sr); hold on; j=j+1;endj=1;while F(9,j)=0 plot(A(F(9,j),2),A(F(9,j),3),sg); hold on; j=j+1;endj=1;while F(10,j)=0 plot(A(F(10,j),2),A(F(10,j),3),sb); hold on; j=j+1;endj=1;while F(11,j)=0 plot

21、(A(F(11,j),2),A(F(11,j),3),sc); hold on; j=j+1;endj=1;while F(12,j)=0 plot(A(F(12,j),2),A(F(12,j),3),sm); hold on; j=j+1;endj=1;while F(13,j)=0 plot(A(F(13,j),2),A(F(13,j),3),sy); hold on; j=j+1;endj=1;while F(14,j)=0 plot(A(F(14,j),2),A(F(14,j),3),sk); hold on; j=j+1;endj=1;while F(15,j)=0 plot(A(F

22、(15,j),2),A(F(15,j),3),r); hold on; j=j+1;endj=1;while F(16,j)=0 plot(A(F(16,j),2),A(F(16,j),3),g); hold on; j=j+1;endj=1;while F(17,j)=0 plot(A(F(17,j),2),A(F(17,j),3),b); hold on; j=j+1;endj=1;while F(18,j)=0 plot(A(F(18,j),2),A(F(18,j),3),c); hold on; j=j+1;endj=1;while F(19,j)=0 plot(A(F(19,j),2

23、),A(F(19,j),3),m); hold on; j=j+1;endj=1;while F(20,j)=0 plot(A(F(20,j),2),A(F(20,j),3),y); hold on; j=j+1;end附錄220個維護(hù)站到13個漏油點的最短距離G=zeros(20,13);%G為20個維護(hù)站到13個漏油點的最短距離for ii=1:20 for j=1:13 G(ii,j)=D(ii,y(j); endendLingo過程代碼Min=x1_1*222.3615+x1_2*160.2847+x1_3*92.86812+x1_4*192.9344+x1_5*210.9621+x1

24、_6*225.0175+x1_7*228.932+x1_8*190.0116+x1_9*195.1581+x1_10*120.8344+x1_11*58.80935+x1_12*118.5011+x1_13*48.85217+x2_1*204.6392+x2_2*141.2972+x2_3*73.88063+x2_4*173.9469+x2_5*191.9747+x2_6*206.03+x2_7*211.2097+x2_8*172.2893+x2_9*177.4358+x2_10*103.1121+x2_11*39.82186+x2_12*103.0954+x2_13*60.35068+x3_

25、1*183.5227+x3_2*127.6723+x3_3*60.25566+x3_4*160.3219+x3_5*178.3497+x3_6*192.4051+x3_7*190.0932+x3_8*151.1728+x3_9*156.3192+x3_10*81.9956+x3_11*60.9384+x3_12*81.97884+x3_13*43.93385+x4_1*219.9738+x4_2*150.0851+x4_3*82.66853+x4_4*182.7348+x4_5*200.7626+x4_6*214.8179+x4_7*226.5443+x4_8*162.2691+x4_9*15

26、5.3534+x4_10*81.02976+x4_11*48.60976+x4_12*73.95869+x4_13*3.5+x5_1*176.2819+x5_2*129.6963+x5_3*62.27967+x5_4*162.3459+x5_5*177.4952+x5_6*191.5506+x5_7*182.8524+x5_8*113.0687+x5_9*106.1529+x5_10*31.82933+x5_11*94.21119+x5_12*24.75826+x5_13*52.55075+x6_1*176.5878+x6_2*130.0021+x6_3*62.58552+x6_4*162.6

27、518+x6_5*177.8011+x6_6*191.8565+x6_7*183.1583+x6_8*113.3745+x6_9*106.4588+x6_10*32.13518+x6_11*94.51704+x6_12*25.06411+x6_13*53.37332+x7_1*149.1494+x7_2*109.0122+x7_3*41.59559+x7_4*14.6619+x7_5*150.3627+x7_6*164.4181+x7_7*155.7199+x7_8*85.70218+x7_9*80.15457+x7_10*5.830952+x7_11*73.52711+x7_12*12.90

28、202+x7_13*79.91722+x8_1*140.9251+x8_2*94.33943+x8_3*26.92282+x8_4*126.9891+x8_5*142.1384+x8_6*156.1938+x8_7*147.4956+x8_8*102.2803+x8_9*104.9318+x8_10*30.6082+x8_11*58.85434+x8_12*30.99467+x8_13*86.77283+x9_1*130.1071+x9_2*82.74202+x9_3*15.3254+x9_4*115.3917+x9_5*131.3205+x9_6*145.3759+x9_7*136.6776

29、+x9_8*97.75722+x9_9*107.2441+x9_10*34.92304+x9_11*47.25692+x9_12*41.9941+x9_13*93.36668+x10_1*75.86585+x10_2*127.7566+x10_3*69.5667+x10_4*95.10693+x10_5*77.07918+x10_6*91.13456+x10_7*82.43635+x10_8*141.9486+x10_9*151.4355+x10_10*79.11446+x10_11*101.4982+x10_12*86.18553+x10_13*147.608+x11_1*37.91353+

30、x11_2*83.37298+x11_3*113.9503+x11_4*50.72332+x11_5*32.69557+x11_6*46.75095+x11_7*38.05274+x11_8*186.3323+x11_9*195.8191+x11_10*123.4981+x11_11*145.8818+x11_12*130.5691+x11_13*191.9916+x12_1*0+x12_2*119.5028+x12_3*145.4326+x12_4*86.85316+x12_5*68.82541+x12_6*64.77002+x12_7*35.9163+x12_8*217.8145+x12_

31、9*227.3013+x12_10*154.9803+x12_11*177.3641+x12_12*162.0514+x12_13*223.4738+x13_1*59.77002+x13_2*59.7328+x13_3*127.1494+x13_4*27.08314+x13_5*9.055385+x13_6*5+x13_7*23.85372+x13_8*228.0832+x13_9*237.57+x13_10*165.249+x13_11*161.2082+x13_12*172.3201+x13_13*213.3179+x14_1*119.5028+x14_2*0+x14_3*67.41662

32、+x14_4*32.64966+x14_5*50.67741+x14_6*64.7328+x14_7*83.5865+x14_8*180.4992+x14_9*189.1668+x14_10*144.8432+x14_11*101.4754+x14_12*121.9142+x14_13*153.5851+x15_1*170.2961+x15_2*132.9808+x15_3*65.52421+x15_4*165.6305+x15_5*171.5094+x15_6*185.5648+x15_7*176.8666+x15_8*47.51842+x15_9*57.00525+x15_10*44.01

33、472+x15_11*97.49573+x15_12*51.08579+x15_13*118.101+x16_1*145.4325+x16_2*67.41662+x16_3*0+x16_4*100.0663+x16_5*118.094+x16_6*132.1494+x16_7*151.0031+x16_8*113.0826+x16_9*121.7502+x16_10*47.42655+x16_11*34.05877+x16_12*54.49761+x16_13*86.16853+x17_1*218.9211+x17_2*149.0324+x17_3*81.6158+x17_4*181.6821

34、+x17_5*199.7098+x17_6*213.7652+x17_7*225.4916+x17_8*186.5712+x17_9*195.2387+x17_10*120.9151+x17_11*47.55703+x17_12*127.9862+x17_13*78.20525+x18_1*242.4718+x18_2*185.1448+x18_3*117.7282+x18_4*217.7945+x18_5*235.8223+x18_6*249.8776+x18_7*249.0423+x18_8*210.1219+x18_9*215.2683+x18_10*140.9447+x18_11*83

35、.66946+x18_12*136.9926+x18_13*67.34362+x19_1*225.46523+x19_2*169.6148+x19_3*102.1982+x19_4*202.2645+x19_5*220.2923+x19_6*234.3476+x19_7*232.0358+x19_8*193.1153+x19_9*198.2618+x19_10*123.9382+x19_11*76.39281+x19_12*119.9861+x19_13*50.33709+x20_1*269.458+x20_2*212.1311+x20_3*144.7145+x20_4*244.7808+x2

36、0_5*262.8085+x20_6*276.8639+x20_7*276.0286+x20_8*230.1082+x20_9*223.1925+x20_10*148.8689+x20_11*110.6557+x20_12*141.7978+x20_13*64.4888;x1_1+x1_2+x1_3+x1_4+x1_5+x1_6+x1_7+x1_8+x1_9+x1_10+x1_11+x1_12+x1_13=1;x2_1+x2_2+x2_3+x2_4+x2_5+x2_6+x2_7+x2_8+x2_9+x2_10+x2_11+x2_12+x2_13=1;x3_1+x3_2+x3_3+x3_4+x3

37、_5+x3_6+x3_7+x3_8+x3_9+x3_10+x3_11+x3_12+x3_13=1;x4_1+x4_2+x4_3+x4_5+x4_5+x4_6+x4_7+x4_8+x4_9+x4_10+x4_11+x4_12+x4_13=1;x5_1+x5_2+x5_3+x5_4+x5_5+x5_6+x5_7+x5_8+x5_9+x5_10+x5_11+x5_12+x5_13=1;x6_1+x6_2+x6_3+x6_4+x6_5+x6_6+x6_7+x6_8+x6_9+x6_10+x6_11+x6_12+x6_13=1;x7_1+x7_2+x7_3+x7_4+x7_5+x7_6+x7_7+x7_

38、8+x7_9+x7_10+x7_11+x7_12+x7_13=1;x8_1+x8_2+x8_3+x8_4+x8_5+x8_6+x8_7+x8_8+x8_9+x8_10+x8_11+x8_12+x8_13=1;x9_1+x9_2+x9_3+x9_4+x9_5+x9_6+x9_7+x9_8+x9_9+x9_10+x9_11+x9_12+x9_13=1;x10_1+x10_2+x10_3+x10_4+x10_5+x10_6+x10_7+x10_8+x10_9+x10_10+x10_11+x10_12+x10_13=1;x11_1+x11_2+x11_3+x11_4+x11_5+x11_6+x11_7+x11_8+x11_9+x11_10+x11_11+x11_12+x11_13=1;x12_1+x12_2+x12_3+x12_4+x12_5+x12_6+x12_7+x12_8+x12_9+x12_10+x12_11+x12_12+x12_13=1;x13_1+x13_2+x13_3+x13_4+x13_5+x13_6+x13_7+x13_8+x13_9+x13_10+x13_11+x13_12+x13_13=1;x14_1+x14_2+x14_3+x14_4+x14_5+x14_6+x14_7+x14_8+x14_9+x14_10