基于查算表的中尺度大气光化学模拟分段线性回归算法的制作方法

本发明属于大气环境污染物检测技术领域，尤其涉及中尺度大气光化学污染的改进型模拟预测模型。

背景技术：

空气质量模式，是在对污染物排入大气环境后传输、扩散、转化和清除等一系列物理和化学过程的认识基础上，利用气象、环境、物理、化学等学科的研究方法和计算机技术，实现模拟和预报不同空间尺度上空气污染物浓度分布状况及变化趋势的方法，在空气质量预报、大气污染控制、环境规划与管理、城市建设及公共卫生等方面均有重要的实际应用价值，具有广阔的发展前景。

数值模拟在当前对流层光化学研究的中得到广泛应用，它对于了解O3的生消机理以及制定O3控制对策具有重要作用。目前，为模拟区域对流O3及其前体物NOX、VOC的分布和演化，大气化学模式大多将化学反应模块(见一的(五))耦合入传输模式进行计算。CALGRID也是这样，包含的化学机制有SAPRC90和CBM4。在求解(1)式中(P-L)GAS项时，CALGRID默认方法是利用QSSA或Hybrid方法求解刚性的浓度变化方程组。应用CALGRID模式进行模拟计算时，大量的时间将花费在每个格点的光化学反应的计算上。因此，在当前的计算条件下，CALGRID模式这种直接耦合的计算方法对于污染物变化趋势的模拟还没有达到污染物预测的时效要求。于是，如何缩短化学模式的计算时间，以便于光化学污染事件的快速预测和响应，成为实际工作中有待解决的问题。

技术实现要素：

针对现有技术存在的问题，本发明提供了一种基于查算表的中尺度大气光化学模拟预测分段线性回归算法，能对CALGRID模型中(P-L)GAS化学变换项的处理兼具有效性和快速性。

为解决上述技术问题，本发明采用了以下技术方案：一种基于查算表的中尺度大气光化学模拟分段线性回归算法，包括以下步骤：

步骤1：以CALGRID模式为基础，考虑大气化学反应、大气输送与扩散、沉降、底面积高架排放源的影响，并对化学物种原始浓度中物理量的平均量分解为脉动量，得到化学物种浓度变化方程如式(1)，

式中，C是化学物种平均浓度，V是三维风矢量的平均量，K是湍流扩散系数，E是污染源排放，是由于沉降引起的物种浓度变化，PCHEM是化学产生率，LCHEM是化学损失率；式(1)中二阶湍流扩散项通过湍流扩散系数K理论闭合转化得到；式(1)中右侧多项式依次分别为平流项、扩散项、源项、沉降项、化学变化项；

步骤2：以三维空气质量模式WRF-CHEM中的气相化学模块对CALGRID模式中化学变化项(PCHEM-LCHEM)建立零维化学机理模式，并考虑初始O3浓度、NOX浓度、VOC/NOX、太阳天顶角、高度、云水混合比、温度七个因子对O3浓度的影响；

步骤3：对步骤2中七个因子进行分档，见表1：

表1建立O3、NOX、VOC查算表时考虑的因子及其分档

步骤4：对步骤3中的表1中各种不同情况下的O3、NOX和VOC的浓度进行计算，最终得到查算表。

进一步的，将步骤2中所述初始O3浓度、NOX浓度、VOC/NOX、太阳天顶角、高度、云水混合比、温度七个因子的取值进行分段取值，各因子的取值分段内的取值与对流层臭氧浓度建立线性拟合，以各影响因子为自变量，臭氧浓度y为因变量建立多元线性回归模型，其中：

将初始臭氧浓度设为x1，由于其与最终臭氧浓度高度线性相关，因此不设置特定的区间；将NOX浓度设为x2，取2—10、10—50、50—80、80—120四个区间；用VOC浓度代表VOC/NOX，设为x3，取VOC/NOX值为0.5—2、2—10两个区间；将太阳天顶角(°)设为x4，取0—30、30—90两个区间；将高度设为x5，取0—4000、4000—17500两个区间；将云水混合比，设为x6，取0.000025—0.0002、0.0002—0.0005、0.0005—0.00125三个区间；将温度设为x7，取268—293、293—313两个区间，对以上七个主要因子的各取值区间进行排列组合，得到七元线性回归模型；

最后将查算表和分段线性回归模型引用语区域空气质量模式中，模拟区域光化学污染物的变化和分布。

进一步的，所述七个主要因子进行排列组合，得到的七元线性回归模型如下：(1)y＝-3.240208E-01+9.964009E-01*x1-2.065380E-02*x2+5.990107E-02*x3-5.931150E-04*x4+5.583519E-05*x5-585.583800*x6+1.901724E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(2)y＝-2.723271E-01+9.960434E-01*x1-3.568400E-02*x2+1.028046E-01*x3-1.463320E-03*x4+7.182861E-05*x5-974.030000*x6+1.929903E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(3)y＝-3.272779E-02+9.964944E-01*x1-2.671353E-02*x2+2.294145E-02*x3-7.859087E-04*x4+5.055262E-05*x5+73.771790*x6+6.419718E-04*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(4)y＝1.309971E-01+9.956155E-01*x1-3.668855E-02*x2+4.392214E-02*x3-7.471496E-04*x4+6.450880E-05*x5-59.405400*x6+3.592733E-04*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(5)y＝5.821121E-01+9.963047E-01*x1-2.874934E-02*x2+9.012221E-03*x3-1.427647E-04*x4+4.905372E-05*x5-59.390620*x6-1.278407E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(6)y＝1.959712E-01+9.954416E-01*x1-3.973377E-02*x2+2.169921E-02*x3-4.534774E-04*x4+6.300749E-05*x5-72.632220*x6+2.979054E-04*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(7)y＝-1.824551+1.000130*x1-1.821033E-02*x2+1.606262E-01*x3+1.350742E-03*x4+6.068224E-05*x5-1563.824000*x6+5.499422E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(8)y＝-1.850930+1.000620*x1-2.396976E-02*x2+1.859919E-01*x3+5.090907E-04*x4+6.428771E-05*x5-1959.192000*x6+5.397780E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(9)y＝-6.568277E-01+9.988475E-01*x1-1.576972E-02*x2+1.314013E-01*x3+3.311812E-03*x4+4.718762E-05*x5-493.532100*x6+1.668974E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(10)y＝-9.148871E-01+9.988598E-01*x1-1.815725E-02*x2+1.318464E-01*x3+2.217415E-03*x4+5.257754E-05*x5-589.966300*x6+2.543386E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(11)y＝-3.072163E-01+9.988124E-01*x1-1.546677E-02*x2+6.877820E-02*x3+1.716454E-03*x4+3.392976E-05*x5-110.801900*x6+6.493933E-04*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(12)y＝-3.726518E-01+9.986563E-01*x1-1.816652E-02*x2+6.200048E-02*x3+1.061970E-03*x4+3.902526E-05*x5-121.240200*x6+8.591420E-04*x7(x2∈(2,10),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(13)y＝3.081077E-01+9.962174E-01*x1-2.381927E-02*x2+2.289869E-02*x3-2.000344E-03*x4+5.104330E-05*x5-228.536400*x6-8.294025E-05*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(14)y＝3.124634E-02+9.956563E-01*x1-3.774594E-02*x2+3.971992E-02*x3-3.319472E-03*x4+6.447939E-05*x5-401.404300*x6+1.272312E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(15)y＝2.790606E-01+9.964704E-01*x1-3.118108E-02*x2+7.124420E-03*x3-1.110823E-03*x4+5.087587E-05*x5+124.614000*x6-3.404965E-04*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(16)y＝2.116826E-01+9.955676E-01*x1-4.177957E-02*x2+1.276735E-02*x3-1.994362E-03*x4+6.592003E-05*x5-8.116969*x6+3.596415E-04*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(17)y＝7.329265E-01+9.963008E-01*x1-3.148335E-02*x2-5.452790E-03*x3-5.851859E-04*x4+5.065716E-05*x5-27.441320*x6-1.791743E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(18)y＝3.548505E-01+9.953990E-01*x1-4.457491E-02*x2-3.247681E-03*x3-1.413143E-03*x4+6.860931E-05*x5-42.275850*x6-5.781889E-05*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(19)y＝-7.467231E-01+9.985067E-01*x1+1.758507E-04*x2+7.397251E-02*x3-4.808137E-03*x4+2.150459E-05*x5-448.531900*x6+3.220839E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(20)y＝-9.694228E-01+9.988482E-01*x1-8.947784E-03*x2+6.914812E-02*x3-7.142411E-03*x4+2.723278E-05*x5-569.478800*x6+4.124503E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(21)y＝-4.678397E-01+9.985849E-01*x1-5.550351E-03*x2+4.727684E-02*x3-2.991841E-03*x4+1.881227E-05*x5-25.732930*x6+1.870373E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(22)y＝-6.108454E-01+9.986129E-01*x1-1.183719E-02*x2+4.064662E-02*x3-4.458526E-03*x4+2.416343E-05*x5-123.106600*x6+2.509716E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(23)y＝3.601069E-02+9.985658E-01*x1-7.485946E-03*x2+2.058266E-02*x3-2.205337E-03*x4+1.678560E-05*x5-46.871980*x6+1.430738E-04*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(24)y＝-2.478765E-01+9.985423E-01*x1-1.413486E-02*x2+1.163396E-02*x3-3.303614E-03*x4+2.205447E-05*x5-56.847990*x6+1.230728E-03*x7(x2∈(2,10),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(25)y＝-4.368167+9.922656E-01*x1+6.788243E-02*x2+3.744743E-02*x3-3.311060E-03*x4+1.379768E-04*x5-1634.702000*x6+1.619045E-02*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(26)y＝-2.471550E-02+9.892216E-01*x1+1.110628E-01*x2+1.041566E-01*x3+2.299108E-03*x4+1.454969E-04*x5-2180.613000*x6-2.214575E-04*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(27)y＝-2.345732+9.937021E-01*x1+2.603346E-02*x2+7.284725E-03*x3-2.369582E-03*x4+1.032368E-04*x5-322.408100*x6+9.295531E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(28)y＝-7.304223E-01+9.918218E-01*x1+4.502220E-02*x2+4.798276E-02*x3+1.123742E-03*x4+1.061126E-04*x5-573.156600*x6+2.935073E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(29)y＝-7.697839E-01+9.945830E-01*x1+1.698890E-03*x2+5.369201E-03*x3-1.579162E-03*x4+7.510532E-05*x5-109.837900*x6+3.670128E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(30)y＝3.451672E-01+9.923933E-01*x1+1.700697E-02*x2+2.844298E-02*x3+2.227075E-04*x4+8.482892E-05*x5-266.547800*x6-1.435603E-04*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(31)y＝-2.887084+9.975139E-01*x1+8.601198E-02*x2+1.315123E-01*x3+6.905591E-03*x4+1.185616E-04*x5-3223.169000*x6+6.233539E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(32)y＝-3.641623+9.959201E-01*x1+1.396759E-01*x2+1.710466E-01*x3+2.253887E-03*x4+9.390853E-05*x5-3501.076000*x6+8.641128E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(33)y＝-2.185009+9.974869E-01*x1+5.556955E-02*x2+5.737754E-02*x3+1.145451E-03*x4+9.912318E-05*x5-973.314000*x6+5.725153E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(34)y＝-3.443886+9.965073E-01*x1+7.965098E-02*x2+8.619615E-02*x3+2.980210E-03*x4+9.343903E-05*x5-1138.395000*x6+9.333369E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(35)y＝-7.970152E-01+9.976070E-01*x1+2.907395E-02*x2+4.144833E-02*x3+3.184299E-04*x4+7.061394E-05*x5-401.295600*x6+1.806758E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(36)y＝-1.118565+9.959723E-01*x1+5.649757E-02*x2+6.056751E-02*x3+2.781465E-04*x4+6.650959E-05*x5-476.712600*x6+2.790875E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(37)y＝-1.664885+9.926918E-01*x1+4.470778E-02*x2+1.192070E-02*x3-8.350832E-03*x4+1.114390E-04*x5-374.144100*x6+7.782010E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(38)y＝2.384221+9.889320E-01*x1+9.486073E-02*x2+4.406160E-02*x3-1.341317E-02*x4+1.392843E-04*x5-556.936100*x6-5.802624E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(39)y＝-9.304208E-01+9.949607E-01*x1+1.771931E-03*x2+5.474444E-03*x3-4.258139E-03*x4+7.577182E-05*x5+25.288380*x6+4.452014E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(40)y＝4.312990E-01+9.932446E-01*x1+1.509342E-02*x2+2.704296E-02*x3-7.816507E-03*x4+9.210577E-05*x5-19.379630*x6-2.417344E-05*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(41)y＝1.800170E-01+9.952010E-01*x1-1.269084E-02*x2+2.098616E-03*x3-2.475425E-03*x4+6.712468E-05*x5-60.639230*x6+5.135286E-04*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(42)y＝1.633596+9.947525E-01*x1-5.643301E-03*x2+3.988956E-02*x3-5.384440E-03*x4+4.026228E-07*x5-142.120200*x6-4.129409E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(43)y＝-1.211134+9.952862E-01*x1+1.040021E-01*x2+7.374608E-02*x3-1.656395E-02*x4+3.953829E-05*x5-743.448900*x6+5.159724E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(44)y＝1.159495+9.933154E-01*x1+1.414483E-01*x2+1.059156E-01*x3-2.061150E-02*x4+3.011642E-05*x5-762.916600*x6-2.612516E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(45)y＝-1.761401+9.978502E-01*x1+4.146809E-02*x2+3.555228E-02*x3-1.001503E-02*x4+3.728897E-05*x5+252.215300*x6+6.186068E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(46)y＝-9.576707E-01+9.971114E-01*x1+5.710833E-02*x2+5.901777E-02*x3-1.372914E-02*x4+3.614406E-05*x5+175.402100*x6+3.423160E-03*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(47)y＝1.003763E-01+9.976429E-01*x1+2.734632E-02*x2+2.377641E-02*x3-7.450061E-03*x4+2.935615E-05*x5-236.897200*x6+5.483131E-04*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(48)y＝3.493710E-01+9.965521E-01*x1+4.329569E-02*x2+3.648369E-02*x3-1.015974E-02*x4+3.073803E-05*x5-304.258400*x6+7.001265E-05*x7(x2∈(2,10),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(49)y＝-4.068297+9.793764E-01*x1+1.151110E-02*x2+7.595756E-01*x3+1.235805E-03*x4+2.484733E-04*x5-3853.707000*x6+1.560803E-02*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(50)y＝-9.385123E-01+9.807823E-01*x1+1.911071E-02*x2+1.441576*x3+2.773124E-03*x4+2.258778E-04*x5-5662.412000*x6+1.532608E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(51)y＝-3.152284+9.802212E-01*x1-8.397728E-03*x2+4.473833E-01*x3+2.582128E-03*x4+2.291173E-04*x5-711.046600*x6+1.341236E-02*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(52)y＝-9.400721E-01+9.784378E-01*x1+1.437047E-03*x2+9.855671E-01*x3+4.323706E-03*x4+2.428531E-04*x5-266.570100*x6+2.377880E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(53)y＝-9.009330E-01+9.801690E-01*x1-1.615795E-02*x2+2.702809E-01*x3+3.475548E-04*x4+2.225798E-04*x5-314.630900*x6+6.281761E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(54)y＝-1.904919+9.758961E-01*x1-6.886203E-03*x2+7.410426E-01*x3+2.983172E-03*x4+2.597334E-04*x5-195.206100*x6+7.180118E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(55)y＝-2.325166+9.932447E-01*x1+2.016241E-02*x2+1.213590*x3+5.513672E-03*x4+9.500864E-05*x5-5423.164000*x6+4.525600E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(56)y＝-7.933130E-01+9.949995E-01*x1+2.141758E-02*x2+1.574441*x3-5.372633E-03*x4+8.326916E-05*x5-6831.569000*x6-9.892448E-04*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(57)y＝-1.603997E-02+9.904087E-01*x1+8.389562E-03*x2+7.700873E-01*x3+1.881665E-02*x4+1.022214E-04*x5-2295.905000*x6-1.701525E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(58)y＝-1.606499E-01+9.914483E-01*x1+1.362824E-02*x2+1.140721*x3+3.043819E-03*x4+7.996368E-05*x5-1199.129000*x6-2.514759E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(59)y＝-8.960595E-01+9.910848E-01*x1+2.019272E-03*x2+4.667926E-01*x3+1.015497E-02*x4+8.518482E-05*x5-342.326900*x6+1.424380E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(60)y＝-3.055983E-01+9.909601E-01*x1+9.919812E-03*x2+8.338131E-01*x3+9.918531E-04*x4+6.013386E-05*x5-98.426550*x6-1.720040E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(61)y＝-1.983608+9.786157E-01*x1-1.020248E-03*x2+4.558848E-01*x3-2.140059E-02*x4+2.885153E-04*x5-1433.667000*x6+1.218415E-02*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(62)y＝-1.145338E-01+9.769415E-01*x1+7.079334E-03*x2+9.532245E-01*x3-2.540609E-02*x4+3.299398E-04*x5-2970.766000*x6+4.133362E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(63)y＝-9.222234E-01+9.802819E-01*x1-1.549873E-02*x2+3.102884E-01*x3-1.376449E-02*x4+2.451195E-04*x5-362.136300*x6+7.860989E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(64)y＝-1.448364+9.753239E-01*x1-6.751831E-03*x2+5.875958E-01*x3-2.248026E-02*x4+3.408139E-04*x5-1385.056000*x6+1.055895E-02*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(65)y＝1.162606+9.803939E-01*x1-2.227140E-02*x2+1.426101E-01*x3-8.088112E-03*x4+2.462770E-04*x5-80.102230*x6-8.188606E-05*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(66)y＝1.252874E-01+9.746158E-01*x1-1.825038E-02*x2+2.976048E-01*x3-1.596521E-02*x4+3.530837E-04*x5-18.288860*x6+4.183286E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(67)y＝-1.247715+9.917399E-01*x1+1.903525E-02*x2+5.325450E-01*x3-2.494083E-02*x4+7.179968E-05*x5-2323.908000*x6+6.390039E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(68)y＝-2.152172+9.926344E-01*x1+2.250471E-02*x2+9.150219E-01*x3-2.604114E-02*x4+5.518390E-05*x5-3385.534000*x6+8.231821E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(69)y＝-4.612176E-01+9.913595E-01*x1+8.873804E-03*x2+3.426107E-01*x3-1.768040E-02*x4+6.938956E-05*x5-1153.404000*x6+3.982377E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(70)y＝-2.077077+9.911726E-01*x1+1.436401E-02*x2+6.348427E-01*x3-2.161581E-02*x4+6.287773E-05*x5-1144.878000*x6+8.541320E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(71)y＝-3.199580E-01+9.917509E-01*x1+1.850327E-03*x2+1.795252E-01*x3-1.332499E-02*x4+7.191153E-05*x5-132.073500*x6+2.276957E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(72)y＝-9.026651E-01+9.909828E-01*x1+7.062571E-03*x2+3.947388E-01*x3-1.556330E-02*x4+6.711970E-05*x5+169.633400*x6+2.783605E-03*x7(x2∈(10,50),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(73)y＝-67.277380+9.526719E-01*x1+2.150776E-01*x2+5.926964E-01*x3-7.657846E-02*x4+2.157669E-04*x5-23710.240000*x6+2.425835E-01*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(74)y＝-8.144266E-02+9.177592E-01*x1+4.016503E-01*x2+1.280748*x3-2.653963E-02*x4+3.457272E-04*x5-25139.850000*x6-1.373116E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(75)y＝-37.575570+9.568909E-01*x1+1.022116E-01*x2+1.755540E-01*x3-2.705634E-02*x4+4.951081E-04*x5-5718.874000*x6+1.428326E-01*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(76)y＝-13.179510+9.432965E-01*x1+2.066247E-01*x2+5.014628E-01*x3-3.672522E-02*x4+5.000925E-04*x5-9456.723000*x6+4.944349E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(77)y＝-20.995860+9.654289E-01*x1+5.134772E-02*x2+1.254074E-01*x3-2.736752E-02*x4+3.203315E-04*x5-1267.697000*x6+8.155224E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(78)y＝-3.095740+9.533170E-01*x1+1.150660E-01*x2+3.162320E-01*x3-3.560353E-02*x4+2.957570E-04*x5-2319.914000*x6+1.660067E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(79)y＝-37.598520+9.648107E-01*x1+2.345120E-01*x2+8.868453E-01*x3-3.590833E-02*x4+3.450520E-04*x5-24431.410000*x6+1.204598E-01*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(80)y＝7.471283E-01+9.347926E-01*x1+4.131311E-01*x2+1.759082*x3-2.945104E-02*x4+2.495652E-05*x5-26058.020000*x6-2.419468E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(81)y＝-23.804860+9.662620E-01*x1+1.530079E-01*x2+4.063675E-01*x3-2.134277E-03*x4+3.757846E-04*x5-10554.820000*x6+8.349140E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(82)y＝-10.703670+9.454387E-01*x1+2.962156E-01*x2+9.829471E-01*x3-3.761961E-03*x4+1.709403E-04*x5-10859.940000*x6+2.730107E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(83)y＝-21.028100+9.729833E-01*x1+9.942949E-02*x2+2.700489E-01*x3-2.090836E-02*x4+2.957499E-04*x5-2387.810000*x6+7.314803E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(84)y＝-2.027840+9.527982E-01*x1+2.088033E-01*x2+6.149850E-01*x3+7.802768E-03*x4+2.190025E-04*x5-2798.272000*x6-4.952510E-04*x7(x2∈(10,50),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(85)y＝-26.144640+9.474962E-01*x1+1.633582E-01*x2+2.916913E-01*x3-1.130467E-01*x4+4.863362E-04*x5-3504.274000*x6+1.094477E-01*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(86)y＝18.844900+9.044617E-01*x1+3.445308E-01*x2+5.393649E-01*x3-1.841617E-01*x4+1.106130E-03*x5-1268.477000*x6-4.462981E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(87)y＝-14.782700+9.684497E-01*x1+6.124426E-02*x2+1.431382E-01*x3-6.062901E-02*x4+3.115574E-04*x5-1485.088000*x6+6.435586E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(88)y＝4.516720+9.546168E-01*x1+1.319053E-01*x2+3.093738E-01*x3-1.032822E-01*x4+4.815080E-04*x5-2369.073000*x6-6.053004E-04*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90)v,x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(89)y＝-3.817863+9.718370E-01*x1+2.321394E-02*x2+9.504198E-02*x3-3.749003E-02*x4+2.754935E-04*x5-669.811400*x6+2.298088E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(90)y＝9.645473+9.586214E-01*x1+6.172916E-02*x2+1.673288E-01*x3-6.508395E-02*x4+4.466397E-04*x5-1082.474000*x6-1.899072E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(91)y＝-18.997940+9.593664E-01*x1+2.240199E-01*x2+4.594425E-01*x3-1.440167E-01*x4+2.284849E-04*x5-4680.951000*x6+7.764503E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(92)y＝15.481980+9.294597E-01*x1+4.155883E-01*x2+9.893557E-01*x3-2.085239E-01*x4+1.107711E-04*x5-5259.125000*x6-4.212879E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(93)y＝-17.594770+9.753751E-01*x1+1.256366E-01*x2+2.183914E-01*x3-9.094041E-02*x4+2.023214E-04*x5-1989.589000*x6+7.082541E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(94)y＝4.659071+9.568983E-01*x1+2.535890E-01*x2+5.442517E-01*x3-1.450103E-01*x4+1.721430E-04*x5-1866.580000*x6-7.468818E-03*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(95)y＝-7.603465+9.781001E-01*x1+8.452740E-02*x2+1.309597E-01*x3-7.049275E-02*x4+1.910222E-04*x5-1563.986000*x6+3.641737E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(96)y＝9.080085+9.613805E-01*x1+1.797981E-01*x2+3.144742E-01*x3-1.152370E-01*x4+2.081061E-04*x5-2463.687000*x6-1.733872E-02*x7(x2∈(10,50),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(97)y＝-41.304200+9.530057E-01*x1+4.776467E-02*x2+3.208994*x3-6.524942E-03*x4+2.731961E-04*x5-11693.470000*x6+1.444051E-01*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(98)y＝5.204725+9.395179E-01*x1+9.188972E-02*x2+5.835634*x3-9.410817E-03*x4+1.234824E-04*x5-10758.030000*x6-3.190811E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(99)y＝-29.224990+9.532951E-01*x1+1.595931E-02*x2+2.251197*x3+7.499541E-03*x4+3.379353E-04*x5-3433.718000*x6+1.062357E-01*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(100)y＝3.796019E-01+9.376501E-01*x1+6.784893E-02*x2+4.624281*x3+1.455173E-02*x4+1.731985E-04*x5-5112.771000*x6-9.759421E-03*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(101)y＝-20.527810+9.585254E-01*x1-5.665152E-03*x2+1.405822*x3-1.018896E-02*x4+2.883298E-04*x5-306.282300*x6+7.734944E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(102)y＝-1.328678+9.410257E-01*x1+4.607726E-02*x2+3.365378*x3-3.279839E-03*x4+1.617774E-04*x5-786.168200*x6-1.319220E-03*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(103)y＝-15.007470+9.661753E-01*x1+6.987780E-02*x2+3.628845*x3+1.280424E-02*x4+1.612308E-04*x5-8891.931000*x6+3.945967E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(104)y＝-4.072728+9.663899E-01*x1+7.518981E-02*x2+4.807811*x3+2.443779E-03*x4+8.927848E-05*x5-9903.953000*x6-1.009444E-03*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(105)y＝-7.749023+9.658133E-01*x1+5.160597E-02*x2+2.700959*x3+4.388245E-02*x4+2.076135E-04*x5-6759.989000*x6+1.888338E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(106)y＝-6.372048+9.668870E-01*x1+4.975396E-02*x2+3.906536*x3+2.338038E-02*x4+1.219001E-04*x5-5385.154000*x6+1.227501E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(107)y＝-11.734570+9.734503E-01*x1+2.217565E-02*x2+1.827746*x3+1.576758E-02*x4+1.928362E-04*x5-895.860300*x6+3.519030E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(108)y＝-5.621375+9.710658E-01*x1+3.366676E-02*x2+3.083864*x3+1.170058E-02*x4+1.066690E-04*x5-909.787000*x6+1.035425E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125)x7∈(293,313))

(109)y＝-22.936570+9.579544E-01*x1+5.724278E-02*x2+2.871880*x3-9.870274E-02*x4+3.801531E-05*x5-6902.101000*x6+8.881824E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(110)y＝9.807717+9.443876E-01*x1+8.880144E-02*x2+4.567663*x3-1.013640E-01*x4+1.585254E-04*x5-7090.578000*x6-3.302963E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(111)y＝-13.004060+9.588403E-01*x1+3.139330E-02*x2+2.101687*x3-6.846342E-02*x4+8.158744E-05*x5-2087.987000*x6+5.324405E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(112)y＝8.251905E-01+9.461744E-01*x1+5.958668E-02*x2+3.308819*x3-9.395129E-02*x4+2.312203E-04*x5-5585.741000*x6+5.799404E-03*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(113)y＝-7.229942+9.612972E-01*x1-1.303386E-02*x2+1.105971*x3-3.490436E-02*x4+2.487584E-04*x5+151.219100*x6+3.404372E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(114)y＝2.940816+9.488205E-01*x1+2.135809E-02*x2+1.913224*x3-6.864620E-02*x4+3.160611E-04*x5-108.112300*x6-2.163213E-03*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(115)y＝-11.936830+9.776380E-01*x1+6.193003E-02*x2+2.304178*x3-8.485042E-02*x4+1.382874E-04*x5-6347.165000*x6+4.329282E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(116)y＝-3.759118+9.768391E-01*x1+6.314319E-02*x2+3.390205*x3-7.563021E-02*x4+7.818541E-05*x5-7006.673000*x6+1.072720E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(117)y＝-8.380909+9.774119E-01*x1+4.517767E-02*x2+1.826621*x3-6.218132E-02*x4+1.452567E-04*x5-3694.888000*x6+3.168184E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(118)y＝-5.121690+9.758981E-01*x1+4.465432E-02*x2+2.733763*x3-6.214803E-02*x4+1.134381E-04*x5-4856.670000*x6+1.930363E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))(119)y＝-8.438271+9.809998E-01*x1+1.334219E-02*x2+1.199683*x3-4.444391E-02*x4+1.722316E-04*x5-474.692800*x6+3.173314E-02*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(120)y＝-2.193199+9.771178E-01*x1+2.837287E-02*x2+1.985922*x3-4.883602E-02*x4+1.484692E-04*x5+-36.9098*x6+5.128405E-03*x7(x2∈(50,80),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(121)y＝-260.765900+9.711961E-01*x1+2.150301E-01*x2+1.687296*x3-2.830020E-01*x4+5.204327E-05*x5-81260.560000*x6+9.487989E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(122)y＝60.973270+9.514270E-01*x1+4.959467E-01*x2+3.576864*x3-3.251543E-01*x4-2.070621E-03*x5-81029.090000*x6-2.347771E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(123)y＝-139.579400+9.323776E-01*x1+1.212930E-01*x2+3.945718E-01*x3-1.036026E-01*x4+1.486933E-03*x5-13164.720000*x6+5.192564E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(124)y＝7.690414+9.308496E-01*x1+2.452191E-01*x2+1.228685*x3-1.879364E-01*x4+3.050077E-04*x5-20069.220000*x6-1.170580E-02*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(125)y＝-70.658500+9.505330E-01*x1+3.101810E-02*x2+3.995057E-01*x3-8.973462E-02*x4+7.469725E-04*x5-4941.216000*x6+2.776060E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(126)y＝-1.217375E-01+9.402534E-01*x1+9.291220E-02*x2+7.223848E-01*x3-1.410149E-01*x4+2.824067E-04*x5-7602.890000*x6+3.011519E-02*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(127)y＝-196.651900+9.627561E-01*x1+3.195279E-01*x2+2.105033*x3-1.969038E-01*x4+7.079219E-04*x5-55868.230000*x6+6.688890E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(128)y＝-18.965620+9.443374E-01*x1+4.645689E-01*x2+4.124652*x3-1.348671E-01*x4+2.493041E-04*x5-45530.960000*x6-1.443615E-03*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(129)y＝-153.037700+9.626841E-01*x1+2.135167E-01*x2+1.160856*x3-1.355926E-01*x4+8.474919E-04*x5-15493.460000*x6+5.277170E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(130)y＝-36.322600+9.531107E-01*x1+3.297354E-01*x2+2.409677*x3-1.151161E-01*x4+6.805047E-04*x5-19335.200000*x6+8.560818E-02*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(131)y＝-119.181200+9.703280E-01*x1+9.997182E-02*x2+8.234596E-01*x3-1.439857E-01*x4+7.214918E-04*x5-4789.802000*x6+4.230369E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(132)y＝-37.210080+9.579915E-01*x1+1.549139E-01*x2+1.627510*x3-1.085771E-01*x4+8.137726E-04*x5-9912.122000*x6+1.225420E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(133)y＝-107.435400+9.590725E-01*x1+5.004289E-02*x2+9.231682E-01*x3-3.163382E-01*x4+5.380478E-04*x5-25310.420000*x6+4.530263E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(134)y＝91.958580+9.387231E-01*x1+2.922581E-01*x2+1.463459*x3-4.862408E-01*x4+1.608543E-04*x5-30287.100000*x6-2.481259E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(135)y＝-57.047200+9.596692E-01*x1+4.970643E-02*x2+2.923747E-01*x3-1.808815E-01*x4+6.931752E-04*x5-3302.231000*x6+2.401877E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(136)y＝35.977620+9.491119E-01*x1+1.629910E-01*x2+4.414748E-01*x3-2.940354E-01*x4+4.907094E-04*x5-5422.053000*x6-7.877889E-02*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(137)y＝-20.477960+9.621268E-01*x1-2.040702E-02*x2+2.692859E-01*x3-1.007270E-01*x4+5.317654E-04*x5-1864.671000*x6+1.052656E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(138)y＝28.030350+9.450977E-01*x1+4.358010E-02*x2+2.066310E-01*x3-1.745209E-01*x4+6.787182E-04*x5-2388.090000*x6-5.422955E-02*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(139)y＝-119.219000+9.749209E-01*x1+2.392741E-01*x2+9.175242E-01*x3-4.105793E-01*x4+5.578739E-04*x5-16466.030000*x6+4.604685E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(140)y＝66.254810+9.604233E-01*x1+4.113871E-01*x2+2.156885*x3-5.221239E-01*x4+4.180325E-04*x5-11748.850000*x6-2.023202E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(141)y＝-102.137200+9.844571E-01*x1+1.594711E-01*x2+3.747348E-01*x3-2.865466E-01*x5+6.113834E-04*x5+2115.782000*x6+3.793078E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(142)y＝30.316650+9.728167E-01*x1+3.262946E-01*x2+9.050711E-01*x3-4.082671E-01*x4+7.042955E-04*x5+3238.043000*x6-9.871565E-02*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(143)y＝-60.855910+9.817285E-01*x1+6.019156E-02*x2+3.002397E-01*x3-1.831116E-01*x4+5.305468E-04*x5-1513.607000*x6+2.380477E-01*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(144)y＝8.981695+9.754149E-01*x1+1.382817E-01*x2+5.870919E-01*x3-2.953852E-01*x4+7.110565E-04*x5-4829.298000*x6+1.628661E-03*x7(x2∈(50,80),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(145)y＝-92.838850+9.227836E-01*x1+6.425527E-02*x2+6.028917*x3-1.806976E-02*x4+5.382141E-04*x5-20763.350000*x6+3.247094E-01*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(146)y＝2.388206+8.954878E-01*x1+1.495277E-01*x2+11.846850*x3-4.805484E-03*x4+4.449294E-04*x5-18165.670000*x6-4.770325E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(147)y＝-67.442600+9.242344E-01*x1+1.762759E-02*x2+4.265764*x3+1.070559E-03*x4+6.385496E-04*x5-6414.220000*x6+2.464900E-01*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(148)y＝-5.830338+8.929380E-01*x1+1.181507E-01*x2+9.242343*x3+3.508816E-02*x4+4.624882E-04*x5-8761.871000*x6-7.505241E-03*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(149)y＝-46.162940+9.322760E-01*x1-1.123576E-03*x2+2.728999*x3-2.740139E-02*x4+4.889521E-04*x5-1404.580000*x6+1.742756E-01*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(150)y＝-8.258425+8.992592E-01*x1+8.998610E-02*x2+6.470732*x3-9.643477E-04*x4+3.061728E-04*x5-1600.246000*x6+8.861413E-03*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(151)y＝-36.584650+9.407218E-01*x1+9.209546E-02*x2+6.446417*x3+2.790491E-02*x4+2.181304E-04*x5-13796.740000*x6+1.058400E-01*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(152)y＝-4.526431+9.349814E-01*x1+1.488102E-01*x2+9.967759*x3-2.090372E-03*x4-9.052463E-05*x5-17042.830000*x6-2.122961E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(153)y＝-26.345320+9.450253E-01*x1+5.580748E-02*x2+4.841273*x3+7.355987E-02*x4+3.202475E-04*x5-11482.310000*x6+8.215021E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(154)y＝-7.838434+9.394884E-01*x1+1.101672E-01*x2+8.102083*x3+2.875557E-02*x4+2.306845E-05*x5-8195.557000*x6-4.924939E-04*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(155)y＝-30.220320+9.553940E-01*x1+3.044287E-02*x2+3.432111*x3+2.255534E-02*x4+2.899454E-04*x5-2304.941000*x6+9.926283E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(156)y＝-6.240188+9.453992E-01*x1+9.693360E-02*x2+6.418663*x3+1.322339E-02*x4+2.785376E-05*x5-621.651600*x6-8.352687E-03*x7(x2∈(80,120),x3∈(0.5,2),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(157)y＝-56.166280+9.313858E-01*x1+6.434479E-02*x2+5.197959*x3-1.747978E-01*x4+2.568768E-04*x5-10905.080000*x6+2.140545E-01*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(158)y＝15.262020+9.052225E-01*x1+1.988708E-01*x2+9.159864*x3-2.287593E-01*x4+3.195819E-04*x5-12574.970000*x6-7.208276E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(159)y＝-36.66480+9.352257E-01*x1+2.278437E-02*x2+3.681256*x3-1.154942E-01*x4+3.139853E-04*x5-3029.993000*x6+1.463580E-01*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(160)y＝4.473351E-01+9.125183E-01*x1+1.508374E-01*x2+6.265718*x3-2.092544E-01*x4+3.331003E-04*x5-6671.407000*x6-6.449669E-03*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(161)y＝-19.893200+9.382225E-01*x1-1.516939E-02*x2+2.079074*x3-6.232395E-02*x4+4.315766E-04*x5-72.262220*x6+8.734944E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(162)y＝3.092788+9.176164E-01*x1+7.248460E-02*x2+3.429343*x3-1.411566E-01*x4+4.488971E-04*x5-1.948532*x6-6.627731E-03*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(163)y＝-31.419360+9.633741E-01*x1+8.711424E-02*x2+4.264078*x3-1.618496E-01*x4+1.968414E-04*x5-10502.140000*x6+1.126269E-01*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(164)y＝1.087515+9.543275E-01*x1+1.587685E-01*x2+7.399195*x3-1.713219E-01*x4-1.797065E-05*x5-12978.470000*x6-2.181080E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(165)y＝-26.774160+9.658952E-01*x1+4.654263E-02*x2+3.384149*x3-1.125798E-01*x4+2.274321E-04*x5-5968.216000*x6+1.015292E-01*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(166)y＝-1.290382+9.552841E-01*x1+1.231650E-01*x2+5.813636*x3-1.421108E-01*x4+9.779599E-05*x5-7434.211000*x6-6.114227E-03*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(167)y＝-22.733760+9.691562E-01*x1+1.827365E-02*x2+2.364750*x3-8.333100E-02*x4+2.696480E-04*x5-1251.194000*x6+8.630653E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(168)y＝4.615942+9.563699E-01*x1+9.893165E-02*x2+4.106521*x3-1.048091E-01*x4+1.759027E-04*x5+772.164400*x6-3.478598E-02*x7(x2∈(80,120),x3∈(0.5,2),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(169)y＝-396.933600+9.526166E-01*x1+2.525258E-01*x2+2.359423*x3-4.925722E-01*x4+2.039524E-05*x5-134947.800000*x6+1.453570*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(170)y＝108.318000+8.828685E-01*x1+7.554229E-01*x2+4.076378*x3-4.258380E-01*x4-2.183178E-03*x5-133877.500000*x6-4.469640E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(171)y＝-199.800700+8.927680E-01*x1+1.548989E-01*x2+4.112165E-01*x3-1.653029E-01*x4+2.343168E-03*x5-22553.510000*x6+7.463322E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(172)y＝15.899720+8.824918E-01*x1+3.443186E-01*x2+1.083500*x3-2.492890E-01*x4+1.090580E-03*x5-31172.740000*x6-4.351455E-02*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(173)y＝-96.727840+9.163131E-01*x1+8.519273E-02*x2+3.823752E-01*x3-1.368241E-01*x4+1.171222E-03*x5-9083.279000*x6+3.788665E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(174)y＝5.162071+8.935758E-01*x1+1.869899E-01*x2+4.364303E-01*x3-1.916598E-01*x4+8.639818E-04*x5-11472.850000*x6+1.389372E-02*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(175)y＝-376.946100+9.637045E-01*x1+1.896301E-01*x2+3.915342*x3-5.681339E-01*x4+5.002005E-04*x5-95770.010000*x6+1.375055*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(176)y＝45.467280+8.759269E-01*x1+6.952113E-01*x2+6.200656*x3-3.094572E-01*x4-3.327774E-05*x5-76868.910000*x6-2.665758E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(177)y＝-289.691300+9.540589E-01*x1+1.421749E-01*x2+2.124764*x3-3.772774E-01*x4+1.048473E-03*x5-24551.100000*x6+1.051429*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(178)y＝-20.584040+9.101976E-01*x1+4.357477E-01*x2+3.221443*x3-2.537760E-01*x4+1.138324E-03*x5-33570.140000*x6+1.762219E-02*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(179)y＝-210.721700+9.535397E-01*x1+9.784803E-02*x2+1.291639*x3-3.158618E-01*x4+1.057248E-03*x5-11632.280000*x6+7.765433E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(180)y＝-27.843060+9.174219E-01*x1+2.474206E-01*x2+1.786803*x3-2.424458E-01*x4+1.535123E-03*x5-19111.410000*x6+9.709703E-02*x7(x2∈(80,120),x3∈(2,10),x4∈(0,30),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))

(181)y＝-150.850200+9.438012E-01*x1+4.992467E-03*x2+1.467345*x3-4.841207E-01*x4+5.618173E-04*x5-51561.840000*x6+6.676195E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(268,293))

(182)y＝129.848500+9.095764E-01*x1+3.875103E-01*x2+1.547508*x3-8.350937E-01*x4+3.741471E-04*x5-60389.210000*x6-3.354945E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.000025,0.0002),x7∈(293,313))

(183)y＝-76.339510+9.403742E-01*x1+1.456072E-02*x2+4.535122E-01*x3-2.650149E-01*x4+1.139739E-03*x5-7988.424000*x6+3.431920E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(268,293))

(184)y＝46.027450+9.214661E-01*x1+2.053119E-01*x2+2.631423E-01*x3-4.862009E-01*x4+1.073730E-03*x5-8441.940000*x6-8.808275E-02*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0002,0.0005),x7∈(293,313))

(185)y＝-25.274350+9.422040E-01*x1-3.142700E-02*x2+3.649581E-01*x3-1.508127E-01*x4+8.275011E-04*x5-3300.267000*x6+1.459524E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(268,293))

(186)y＝35.286630+9.171444E-01*x1+3.988732E-02*x2+1.064684E-01*x3-2.726866E-01*x4+1.209301E-03*x5-3790.051000*x6-5.075740E-02*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(0,4000),x6∈(0.0005,0.00125),x7∈(293,313))

(187)y＝-224.357800+9.639255E-01*x1+1.029272E-01*x2+1.989228*x3-6.078703E-01*x4+6.183745E-04*x5-36700.150000*x6+9.069020E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(268,293))

(188)y＝155.940400+9.173321E-01*x1+5.615374E-01*x2+3.024570*x3-9.138510E-01*x4+6.172678E-04*x5-28181.020000*x6-4.881493E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.000025,0.0002),x7∈(293,313))

(189)y＝-178.120300+9.768384E-01*x1+9.404033E-02*x2+8.618064E-01*x3-4.338468E-01*x4+8.617358E-04*x5-165.373000*x6+6.873459E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(268,293))

(190)y＝64.612730+9.551004E-01*x1+3.925678E-01*x2+1.107687*x3-6.915261E-01*x4+1.199616E-03*x5+5390.041000*x6-2.096809E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0002,0.0005),x7∈(293,313))

(191)y＝-104.615300+9.739983E-01*x1-1.266501E-03*x2+6.242684E-01*x3-2.733380E-01*x4+8.120720E-04*x5-3839.086000*x6+4.300113E-01*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(268,293))

(192)y＝16.755670+9.678883E-01*x1+7.977404E-02*x2+8.320093E-01*x3-4.597599E-01*x4+1.180964E-03*x5-9162.614000*x6+1.993333E-02*x7(x2∈(80,120),x3∈(2,10),x4∈(30,90),x5∈(4000,17500),x6∈(0.0005,0.00125),x7∈(293,313))。

附图说明

图1为本发明所述基于查算表的中尺度大气光化学模拟分段线性回归算法取代复杂耦合化学机制算法的原理框图；

图2(a)为现有技术直接耦合法计算地面O3平均浓度的分布图；

图2(b)为本发明基于查算表的快速预测法计算的地面O3平均浓度的分布图；

图3(a)为现有技术直接耦合法计算的污染过程中各天地面O3平均浓度的分布图；

图3(b)为本发明基于查算表的快速预测法计算的污染过程中各天地面O3平均浓度的分布图；

图4为现有技术直接耦合法与本发明的基于查算表的快速预测法计算的香港地区地面O3浓度逐时变化图。

具体实施方式

下面结合附图并以具体实施例，进一步阐明本发明。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

本研究主要对(1)式中(P-L)GAS项的处理方法进行改进。

首先，利用具有国际先进水平的三维空气质量模式WRF－CHEM中的气相化学模块建立箱模式，并应用其探讨了一些可能影响光化学过程的物理和化学因子进行分析，找出对O3浓度影响大的因子。利用光化学反应的机理分析，以及具有国际先进水平的零维化学机理模式(仅考虑气相化学机制，没有考虑排放、输送、沉降的作用)的化学敏感性分析，发现影响对流层臭氧及其前体物浓度的主要因子有：初始O3浓度、NOX浓度、VOC/NOX、太阳天顶角、高度、云水混合比、温度。

其次，对这些因子进行分档计算，建立各种情况下O3、NOX、VOC浓度查算表。然后，经过分析可发现，各影响因子在不同的数值区间内对臭氧浓度的影响呈现一定的线性化特征，因此可建立各影响因子在不同数值区间内的臭氧浓度分段线性回归方程。

在充分调研我国O3、NOX浓度水平，VOC/NOX的实际状况以及各种真实的气象条件的基础上，对上述重要的影响因子进行分档，如表1所示。分档不仅力求包括各种真实的情形，而且尽量做到反映因子之间相互关系的非线性性。表1共计1406250种不同的情况。应用零维化学机理模式，对不同组合情况下O3、NOX和VOC的浓度进行计算，最终得到查算表1。

表1建立O3、NOX、VOC查算表时考虑的因子及其分档

在对影响对流层臭氧浓度的主要因子进行分析的时候，可以发现，虽然各因子与臭氧浓度之间的关系大多是非线性的，但在特定的数值区间内，却大都可以将二者的关系近似看作线性的，因此可以先找出这些存在这样线性关系的区间，再将各影响因子的取值进行分段组合，以各影响因子为自变量，臭氧浓度为因变量建立多元线性回归模型。

将初始臭氧浓度(ppb)设为x1，由于其与最终臭氧浓度高度线性相关，因此不设置特定的区间；将NOX浓度(ppb)设为x2，取2—10、10—50、50—80、80—120四个区间；用VOC浓度(ppb)代表VOC/NOX，设为x3，取VOC/NOX值为0.5—2、2—10两个区间；将太阳天顶角(°)设为x4，取0—30、30—90两个区间；将高度(m)设为x5，取0—4000、4000—17500两个区间；将云水混合比(kg/kg)设为x6，取0.000025—0.0002、0.0002—0.0005、0.0005—0.00125三个区间；将温度(K)设为x7，取268—293、293—313两个区间。将以上七个主要因子的各取值区间进行排列组合，一共可以得到192个七元线性回归模型。

最后，将查算表和分段线性回归算法应用于区域模式中，模拟区域光化学污染物的变化和分布。如图1所示，完成对CALGRID模型(1)式中(P-L)GAS项的处理方法的改进，使其兼具有效性和快速性。

具体案例：

在三维区域空气质量模式CALGRID中应用基于查算表数据库的分段线性回归算法模拟上述研究个例，并验证其计算结果是否与直接耦合法的一致。模拟区域为18.6°N～53.5°N，90.6°E～151.8°E，大致覆盖中国地区；水平方向分66×56网格，格距为75km。模拟时间为2000年3月25日～31日。气象场由WRF提供，顶气压100hPa；垂直方向31层；采用Mellor-Yamada-Janjic TKE边界层方案、时变边界条件和Grell-Devenyi ensemble scheme积云对流参数化方案。

图2是分别利用直接耦合法和新方法计算的3月27日～31日地面O3平均浓度的分布。由图可见，直接耦合法得到的浓度大值区在四川和华北地区；在华南，珠江三角洲地区也有相对较大的浓度分布(30～40ppb)。新方法得到的臭氧浓度大小和分布同直接耦合法基本一致，特别是大值区分布非常相似。可见，新的快速预测方法对于区域光化学污染的有较好的模拟能力。但是，在个别地区(如甘肃、西藏的部分地区)，两种方法还是存在一定差异，这可能是因为基于查算表数据库的分段线性回归算法的基础是RADM2机制，而CALGRID直接耦合法应用的是SAPRC机制，两种机制的细微差别所致。

图3是直接耦合法和新方法计算的污染过程中各天(28～31日)地面O3平均浓度的分布O3平均浓度的分布。可见，直接耦合法和新方法都较好的模拟了这次污染事件的特征，即28日、29日受上风方(珠三角)污染物影响，污染较为严重(40～90ppb)，30日、31日受海洋上清洁空气的影响，O3浓度较低(小于20ppb)。此外，两种方法模拟的各天臭氧浓度值大小和区域分布都非常相似。

图4给出了两种方法计算的香港地区地面臭氧浓度的逐时变化。可见，两种方法的模拟结果，无论是值的大小，还是曲线的日变化形态，峰值、谷值的出现时间都非常的一致。虽然查表法模拟值略微偏高，比直接耦合法高9ppb，但是，总体来讲两种方法的结果还是符合的较好，相关系数为0.902。

综上所述，本研究建立的快速预测方法同直接耦合法一样，较好的反映了光化学反应的非线性性，对光化学污染事件有较强的模拟能力，可以应用于区域空气质量模式中对区域污染问题进行模拟。另外，由于本方法可以大大节省时间，比直接耦合法的计算速度高3～5倍，因此有利于光化学污染事件的快速预测以及污染控制的快速决策。