Matlab编程天然气压缩因子计算模型

1程序目的利用AGA8-92DC模型计算天然气的压缩因子，该程序主要应用于在输气 和配气正常进行的压力P和温度T范围内的管输气的压缩因子计算2数学模型：AGA8-92DC模型2.1模型介绍此模型是已知气体详细的摩尔分数组成和相关压力、温度來计算气体压缩 因子。输入变量包括绝对压力、热力学温度和摩尔组成。摩尔组成是以摩尔分数表示下列组分:C02、际比、CO、CH,曲、CsHs、i-C】Hio、 n-CiHio i-CsH、 n-CsHi2 n-CbHni、n-CsHw o2.2模型适用条件绝对压力：0MPaVPV12MPa热力学温度：263KWTW338K高位发热量：30MJm3WHsW45 MJm3相对密度：天然气中各组分的摩尔分数应在以下范围内：CH4： N2： 0Wxn?W0.20CO2： 0Wxco2W0.20C2H6： OWxc2H6WO.1OC3H8： 0Wxc3H8W0.035C4H10： OWxc4HioW0.015C5H12： 0Wxc5H12W0.005C6H14： OWXC6H14WO.OOIC7H16： 0Wxc7H16W0.0005C8H18和更高碳数炷类：C8H18： 0Wxc8H18W0.0005H2： OWxh/WO.IOCO： 0WxcoW0.03如果已知体积分数组成，则应将其换算成摩尔分数组成。所有摩尔分数大于0.00005的组分都不可忽略。2.3模型描述已知条件绝对压力P、热力学温度T、组分数N；各组分的摩尔分数治，i = lN；查附表1、2、3得到的以下数据：58种物质的状态方程参数a】，bn, cn, kn, un, gn, qn, sn, wn ；14种识别组分的特征参数Mi，Ei，Ki，Gi，Qi，Fi，W；14种识别组分的二元交互作用参数场,U”，埼,G；j。待求量压缩因子Z计算步骤a）第二维利系数B的计算：ISN N3E =工工兀 xR（KK，n=li=l j=l1丄b*j = 卩+1 - 乳产（QQj +1 - q严 呼+1 - fn）:（筮+i-sn）5（wj+i-wn）w 二元参数和Gj,由以下两式计算：E厂 E；（EE,q=q；（q+Gj/2b）计算系数G，Il = 1358C； = J（G+1 - 弘严（Q，+ Q - q严（F+1- 严IT叮用以下方程求解混合方程，计算混合物参数U, G, QoN5N-l N5U5 = QE + 2工 X （U；-1）（EE护i=li=l j =1+1NN-l NG =工 g + 2工工家（G*- l）（Gx + GJi=l1=1 j=i+l1=1c) 计算混合物体积参数K；N5N-l N5K5 = 肩 J2 + 2工 初曲-1)(KK护1=11=1 j=i-Hd) 计算对比密度摩尔密度为：pm = P/(ZRT)式中，P为绝对压力，Mpa： R为摩尔气体常数；T为热力学温度，Ko对比密度p r同摩尔密度P M相关：R = K九e) 利用AGA8-92DC方程，对压缩因子进行迭代计算IS5SZ = 1 + E久-久工 C： + 工 C：(bx - cnknexp(-cn P)n=13n=13迭代过程：给出ZO的初始值为1,先计算出Pm，将Pm K和己知量带入AGA8-92DC方程方程，得到新的Z值，当(Z-Z0)的绝对值小于0.000001时, 停止迭代，得到Z值。3程序代码fiinction Z = YSYZ( T,p“)%计算天然气给定组分的压缩因子% x 为天然气组分，按照 CO2 N2 H2 CO CH4 C2H6 C3H8 i-C4H10 n-C4H10 i-C5H12 11-C5H12%n-C6H14 11-C7H16 11-C8H18 的顺序输入%T为温度，单位为K%P为压力，单位为兆帕N=14;R=8.314;%状态参数值b=l,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,4,4,4,4, 4,4,4,5,5,5,5,5,60,7,7,8,8,8,9,9;k=0,0,0,0,0,0,0,0,0,0,0,0,3,2,2,2,4,4,0,0,2,2,2,4,4,4,4,0,1,1,2,2,3,3,4,4,4,0,02,2, 2,4,4,02,2,4,4,0,2,0,2,1,2,2,2,2;c=0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,0,0,1,0,1,1,1,1,1,1,1,1,1,0,0,1,1,1g=0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,1,0,0,1,1,1,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,04,0,0,1,0,1,0,0;b0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,1,0,0,0,0,1,0,0,0,0,0,0,0 ,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;q=0,0,0,0,0,0,l ,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,1,0,0,0,0,0,1;s=0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 ,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;w=0,0,0,0,0,0,0,0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;a=0.153832600,1.341953000,-2.998583000,-0.048312280,0.375796500,-1.589575000,.0053588470,0.886594630,0710237040,J.471722000,1.321850350,- 0.786659250,.2.291290e-9,0.157672400,-0.436386400,-0.044081590,-0.003433888,0.032059050,.0.024873550,0.073322790, 0.001600573,0.642470600,0416260100,0.066899570,.0.279179500,-0.696605100,-0.002860589,-0.008098836,3.150547000,0.007224479,.-0.705752900,0.534979200,-0.079314910,-1.418465000,-5.99905e 17,0.105840200,.0.034317290, 0.007022847,0.024955870,0.042968180,0.746545300,0.291961300,.7.294616000,-9.936757000,-0.005399808,-0.243256700,0.049870160,0.033733797,.1.874951000,0.002168144,-0.658716400,0.000205518,0.009776195,-0.020487080,.0.015573220,0.006862415,0001226752,0.002850908;u=0.0,0.5,1.0,3.5,0.5,4.5,0.5,7.5,9.5,6.0,12.0,12.5,6.0,2.0,3.0,2.0,2.0,11.0,0.5,.0.5,00406021023.022.0,10，05,7.0,10,6040，1090,1302108.0, 0.5,0.02.0,7.0.022.0,23.0,1.03.0,3.0,8.0,23.0,1.5,5.0,0.5,4O7O3O0O1O0.0;%特征参数值M=44.0100,28.0135,2.0159,2&0100,16.0430,30.0700,44.0970,58.1230,58.1230 ,72.1500,72.1500,86.1770,100.2024,114.2310;E=241.960600,99.737780,26.957940,105.534800,151.318300,244.166700,298. 118300,324.068900,337.638900,365.599900,370.682300,402.636293,427.722630,45 0.325022;G=0.189065,0.027815,0.034369,0.038953,0.0,0.079300,0.141239,0.256692,0.281835,0.332267,0.366911,0.289731,0.337542,0.383381;Q=0.690,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0;K=0.4557489,0.4479153,0.3514916,0.4533894,0.4619255,0.5279209,0.5837490,0.6406937,0.6341423,0.6738577,0.6798307,0.7175118,0.7525189,0.7849550;F=zeros(N);S=zeros(N);W=zeros(N);F =1；%交互作用参数值Ex=1.0,1.022740,1.281790,1.5,0.960644,0.925053,0.960237,0.906849,0.89736 2,0.726255,0.859764,0.855134,0.831229,0.808310;.1.022740,1.0,1.086320,1.005710,0.971640,0.970120,0.945939,0.946914,0.973384,0. 959340,0.945520,1.0,1.0,1.0;.1.281790,1.086320,1.0,1.1,1.170520,1.164460,1.034787,1.3,1.3,1.0,1.0,1.0,1.0,1.0;.1.5,1.005710,1.1,1.0,0.990126,1.0,1.0,1.0,1.0049,1.0,1.0,1.0,1.0,1.0;.0.960644,0.971640,1.170520,0.990126,1.0,1.0,0.994635,1.019530,0.989844,1 .00235,0.999268,1.107274,0.88080,0.880973;.0.925053,0.970120,1.164460,1.0,1.0,1.0,1.022560,1.0,1.013060,1.0,1.005320,1.0,1.0,1.0;.0.960237,0.945939,1.034787,1.0,0.994635,1.022560,1.0,1.0,1.0049,1.0,1.0,1.0,1.0,1.0;.0.906849,0.946914,1.3,1.0,1.019530,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.897362,0.973384,1.3,1.0,0.989844,1.01306,1.0049,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.726255,0.959340,1.0,1.0,1.00235,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.859764,0.945520,1.0,1.0,0.999268,1.00532,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.855134,1.0,1.0,1.0,1.107274,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.831229,1.0,1.0,1.0,0.88088,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.808310,1.0,1.0,1.0,0.880973,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;Gx=1.0,0.982746,1.0,1.0,0.807653,0.370296,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.982746,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.95731,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.807653,1.0,1.95731,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.370296,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;Ux=1.0,0.835058,1.0,0.9,0.963827,0.969870,1.0,1.0,1.0,1.0,1.0,1.066638,1.077 634,1.088178;.0.835058,1.0,0.408838,1.0,0.886106,0.816431,0.915502,1.0,0.993556,1.0,1.0,1.0,1.0,1.0;.1.0,0.408838,1.0,1.0,1.156390,1.616660,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0；.0.9,1.0,1.0,1.0,1.156390,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.963827,0.886106,1.156390,1.0,1.0,1.0,0.990877,1.0,0.992291,1.0,1.00367,1.302576,1.191904,1.205769;.0.969870,0.816431,1.616660,1.0,1.0,1.0,1.065173,1.25,1.25,1.25,1.25,1.0,1.0,1.0;.1.0,0.915502,1.0,1.0,0.990877,1.065173,1.0,1.0,1.0,1.0,1.0,1.0,1.0,10;.1.0,1.0,1.0,1.0,1.0,1.25,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.1.0,0.993556,1.0,1.0,0.992291,1.25,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.1.0,1.0,1.0,1.0,1.0,1.25,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.1.0,1.0,1.0,1.0,1.00367,1.25,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.066638,1.0,1.0,1.0,1.302576,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.1.077634,1.0,1.0,1.0,1.191904,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.1.088178,1.0,1.0,1.0,1.205769,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;Kx=1 0,0.982361,1.0,1.0,0.995933,1.00851,1.0,1.0,1.0,1.0,1.0,0.910183,0.8953 62,0.881152;0.982361,1.0,1.03227,1.0,1.00363,1.00796,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.03227,1.0,1.0,1.02326,1.02034,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.995933,1.00363,1.02326,1.0,1.0,1.0,1.007619,1.0,0.997596,1.0,1.002529,0.982962,0.983565,0.982707;1.00851,1.00796,1.02034,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.007619,0.986893,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,0.997596,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.002529,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.910183,1.0,1.0,1.0,0.982962,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.895362,1.0,1.0,1.0,0.983565,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.881152,1.0,1.0,1.0,0.982707,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;Z0=l;M0=0;B=0;%参数Eij、Gij以及第二维利系数的计算siiml=0;siim2=0;for n=l:18siim=0;ZJCS=(TA-u(n);for i=l:14for j=l:14Eij(ij)=Ex(ij)*(E(i)*E(j)A0.5);Gij(ij)=Gx(ij)*(G(i)+G(j)/2;Bij(ij)=(Gij(ij)+l-g(n)Ag(n)*(Q(i)*QG)+l - q(n)Aq(n)*(F(i)*FG)A0.5)+l-f(n)Af(n)*(S(i)*S(j)+l-s(n)As(n)*(W(i)：1：W(j)+l- w(n)八 w(i】)；siim=siim4-x(i)*x(j)!,:Bij(ij):l:(K(i)*K(j)/x1.5);endendE=B+a(n)*ZJCS 木 sum;end%K值的计算F0=0;G0=0;Q0=0;U0=0;for i=l:NF0=F0+(x(i)八2)网)；Q0=Q0+x(i)*Q(i);sum 1=sum 1 +x(i)水 G(i);sum2=sum2+x(i) * (E(i)/X2.5);endfor i=l:(N-l)for j=2:NGO=GO+x(i)*xG)*(Gx(ij)-l)*(G(i)+G(j);U0=U0+(Ux(ij)A5)-l)*(E(i)*E(j)A2.5);endendG0=siiml+2*G0;U0=(sum2 八 2)+2 水 U0F0.2;siiml=O;for i=l:Nsum 1=suml +x(i) * (K(i)A2.5);endsiim2=0;for i=l:(N-l)for j=2:N sum2=sum2+x(i) *x(j) * (Kx(i j)八5) 1)水(K(i) *K(j )八2.5);endendK0=(suml 八 2)+2 木 siim2 )八0.2);%对比密度的计算pr=pM： (K0 八 3)/(Z0 *R+T);SUM1=O;SUM2=0;%计算Cii的值，共有46个值for n=13:18Cn=a (1】)冷(GO+1 欽11)八欽11)尸(QO 八 2)+Q0 q(n)/q(n)水(FO+1 f(n)Af(n) * (U0Au(n)* (TA-u(n);SUM1=SUM1+Cn;endfor i=13:58Cn=a(i)*(GO+l-g(i)Ag(i)5,:(QO*QO+QO q(i)Fq(i):l:(F0+l - f(i)Af(i)*(UOAu(i)*(TA-u(i);SUM2=SUM24-Cn* (b(i) - c(i) *k(i) (prAk(i) (pr Ab(i) * exp(-c(i) * (prAk(i); end%Z的计算Z=l+B* (pr/(K0A3)-pr+SUM 1+SUM2;%迭代过程wliile abs(Z-Z0)=0.000001Z0=(Z+Z0)/2;pr=p*(K0A3)/(Z0+R*T);SUM2=0;for n=13:58Cn=a(i)*(G0+l-g(i)Ag(i)*(Q0*Q0+Q0-q(i)Aq(i)*(F0+l-f(i)Af(i)* (UO Au(i)* (TA-u(i);SUM2=SUM2+Cn,:(b(i)-c(i)!,:k(i)*(prAk(i)*(prAb(i)*exp(-c(i)*(prAk(i);endZ=l+B* (pi7(K0A3)-pr+SUM 1+SUM2;endZend4运行结果对气样如下表的天然气计算在300K、O.IMPa下的压缩因子表4.1气样组分气体组成摩尔分数CO20.006N20.003H20CO0CH40.965C2H60.018C3H80.0045i-C4H100.00111-C4H100.001i-C5H120.000511-C5H120.0003C6H140.0007C7H160C8H180 P=0. 1; 1=300; x= 0. 006 0.003 0.00 0.00 0.965 0.018 0.0045 0.001 0.001 0.0005 0.0003 0.0007 0 0： z=YSYZ(T, P,x)0. 9995z =0.9995附录附表1状态参数方程表0.153832600,1.341953000, -2.998583000, -0.048312280, 0.375796500, 1.589575000, -0.053588470, 0.886594630, -0.710237040, -1.471722000, 1.321850350, 0786659250, 2.291290e-9, 0.157672400, -0.436386400,- 0.044081590, -0.003433888, 0.032059050, 0.024873550, 0.073322790, 0.001600573, 0.642470600, -0.416260100, -0.066899570, 0.279179500,- 0.696605100, -0.002860589, -0.008098836, 3.150547000, 0.007224479,- 0.705752900, 0.534979200, -0.079314910, -1.418465000, -5.99905e-17, 0.105840200,0.034317290,0007022847,0.024955870, 0.042968180, 0.746545300, -0.291961300, 7.294616000, -9.936757000, -0.005399808,- 0.243256700, 0.049870160, 0.033733797,1.874951000, 0.002168144,- 0.658716400, 0.000205518, 0.009776195, -0.020487080,0.015573220, 0.006862415, -0.001226752, 0.002850908bn1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,12,22,22,22,22,3,3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,5,5,5,5,5,6,6,7,7,8,8,8,9Cn1,1,0,1,1,1,1,0,1,0,1,1,1,1,14kn0,0,0,0,0,0,0,0,0,0,0,0,3,22,2,4,4,0,02,22,4,4,4,4,0,1,12,2,3,3,4,4,4,0,0,22,2,4,4,0,22,4,4,02,02,1222,2Un0.0,0.5,1.0,3.5 ,0.5,4.5,0.5,7.5,9.5,6.0,12.0,12.5, 6.0,2.0,3.0,2.0,2.0,11.0, 0.5,0.5,0.0,4.0,6.0,21.0,23.0,22.0, 1.0, 0.5,7.0,-1.0,6.0,4.0,1.0.0,13.0,21.0,8.0,-0.5,0.02.0,7.0.022.0,23.0,1.0夕.0,3.0,8.0,23.0,1.5,5.0,0.5,4.0,7.0, 3.0,0.0, 1.0, 0.0Sn0,0,0,0，l,l，0,0,0,0,0O0,0,0,0,0O0,0,0,0,0Ol,0,0Ol,0,011,l，0,0,0,0,0,0,0,0,0,0,0,0，0,0,0,0 J ,0,0 J，0 J ,0,0Qn0O0,0,0,0，l,0,0,0,0O0,0,010O0,0,0,0,0,0,010,l，0,0,0,0,0,0,0,0，l,0,0,0,0,l,0,0,0,0 J ,0 J ,0,0 J ,0,0,0,0,0,1缶0O0,0,0,0,0,0,0,0,0Ol,0,0,0,0O0,0,0,0,0,0,0,0，lO0,l，0,0,0,0，l,0,0,0,0,0,0,0,0,0,0,0，0,0,0,0,0,0,0,0,0,0,0,0Sil0,0，0,0,0,0,011,0,0O0,0,0,0,0,0,0,0,0,0,0,0,0,0,0O0,0,0,0,0,0,0,0,0 0000,0，0,0,0,0，0,0,0,0,0,0,0,0,0,0,0,0Wn0,0，0,0,0,0,0,0,0,l，l,l，0,0,0,0,0,0,0,0,0,0,0,0,0,0,0O0,0,0,0,0,0,0,0,0 0000,0，0,0,0,0，0,0,0,0,0,0,0,0,0,0,0,0附表2特征参数表M44.01002&0135,2.0159,2&0100，160430,300700,440970,5 8.1230,5&1230,72.1500,72.1500,86.1770,100.2024,114.2310E241.960600,99.737780,26.957940,105.534800,151.318300,244.166700,298.118300,324.068900,337.638900,365.599900,370.682300,402.636293,427.722630,450.325022G0.189065,0.027815,0.034369,0.038953,0.0,0.079300,0.141239,0.256692,0.281835,0.332267,0.366911,0.289731,0.337542,0.383381Q0690,00,0O00,00,0O00,00,0O00,0O00,00,00K0.4557489,0.4479153,0.3514916,0.4533894,0.4619255,0.5279209,0.5837490,0.6406937,0.6341423,0.6738577,0.6798307,0.7175118,0.7525189,0.7849550FS00000000000000W00000000000000附表3二元交互作用参数值表Eij*1.0,1.022740,1.281790,1.5,0.960644,0.925053,0.960237,0.906849,0.897362,0.726255,0.859764,0.855134,0.831229,0.808310;1.022740,1.0,1.086320,1.005710,0.971640,0.970120,0.945939,0.946914,0.973384,0.959340,0.945520,1.0,1.0,1.0;1.281790,1.086320,1.0,1.1,1.170520,1.164460,1.034787,1.3,1.3,1.0,1.0,1.0,1.0 ,10;1.5,1.005710,1.1,1.0,0.990126,1.0,1.0,1.0,1.0049,1.0,1.0,1.0,1.0,1.0;0.960644,0.971640,1.170520,0.990126,1.0,1.0,0.994635,1.019530,0.989844,1.00235,0.999268,1.107274,0.88080,0.880973；0.925053,0.970120,1.164460,1.0,1.0,1.0,1.022560,1.0,1.013060,1.0,1.005320,1.0,1.0,1.0;0.960237,0.945939,1.034787,1.0,0.994635,1.022560,1.0,1.0,1.0049,1.0,1.0,1.0 ,1.0,1.0;0.906849,0.946914,13,1.0,1.019530,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.897362,0.973384,1.3,1.0,0.989844,1.01306,1.0049,1.0,1.0,1.0,1.0,1.0,1.0,1.0 0.726255,0.959340,1.0,1.0,1.00235,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.859764,0.945520,1.0,1.0,0.999268,1.00532,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.855134,1.0,1.0,1.0,1.107274,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.831229,1.0,1.0,1.0,0.88088,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.808310,1.0,1.0,1.0,0.880973,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0Uij1.0,0.835058,1.0,0.9,0.963827,0.969870,1.0,1.0,1.0,1.0,1.0,1.066638,1.077634 ,1.088178;.0.835058,1.0,0.408838,1.0,0.886106,0.816431,0.915502,1.0,0.993556,1.0,1.0,1.0,1.0,1.0;.1.0,0.408838,1.0,1.0,1.156390,1.616660,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.9,1.0,1.0,1.0,1.156390,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;.0.963827,0.886106,1.156390,1.0,1.0,1.0,0.990877,1.0,0.992291,1.0,1.00367,1.302576,1.191904,1.205769;0.969870,0.816431,1.616660,1.0,1.0,1.0,1.065173,1.25,1.25,1.25,1.25,1.0,1.0,1.0;1.0,0.915502,1.0,1.0,0.990877,1.065173,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.25,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,0.993556,1.0,1.0,0.992291,1.25,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.25,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.00367,1.25,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.066638,1.0,1.0,1.0,1.302576,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.077634,1.0,1.0,1.0,1.191904,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.088178,1.0,1.0,1.0,1.205769,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0Kij1.0,0.982361,1.0,1.0,0.995933,1.00851,1.0,1.0,1.0,1.0,1.0,0.910183,0.895362, 0.881152;0.982361,1.0,1.03227,1.0,1.00363,1.00796,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.03227,1.0,1.0,1.02326,1.02034,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.995933,1.00363,1.02326,1.0,1.0,1.0,1.007619,1.0,0.997596,1.0,1.002529,0.9 82962,0.983565,0.982707;1.00851,1.00796,1.02034,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.007619,0.986893,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,0.997596,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;1.0,1.0,1.0,1.0,1.002529,1.0,1.0,10,1.0,1.0,1.0,1.0,1.0,1.0；0.910183,1.0,1.0,1.0,0.982962,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.895362,1.0,1.0,1.0,0.983565,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.881152,1.0,1.0,1.0,0.982707,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0Gij水1.0,0.982746,1.0,1.0,0.807653,0.370296,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;0.982746,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0;O to to To to ro to To ro to to To to To tO to to To to ro to To ro to to To to To tO to to To to ro to To ro to to To to To tO to to To to ro to To ro to to To to To tO to to To to ro to To ro to to To to To tO to to To to ro to To ro to to To to To tO to to To to ro to To ro to to To to To tO to to To to ro to To ro to to To to To t0 To ro To To ro To To ro To ro To ro T96ZOZ,roSToToToroToToToToToToTlg6ToT$9O2OO to to To to ro to To ro to to To to To to ro roto to ro to roto rifz.?6TorororoT