含有并行机器组的串行短生产线的性能分析

第２２卷第ｌＯ期　２０１２年１０月　计算机技术与发展　ＣＯＭＰＵＴＥＲ　ＴＥＣＨＮＯＬＯＧＹ　ＡＮＤ　ＤＥＶＥＬＯＰＭＥＮＴ　Ｖ０Ｉ．２２　Ｎｏ．１０　Ｏｃｔ．　２０１２　含有并行机器组的串行短生产线的性能分析　穆华灵，汪峥　（东南大学自动化学院，江苏南京２１００９６）　摘要：文中主要研究含有并行机器组的串行短生产线的生产性能，此生产线含有一个并行机器组及其上下游缓冲区。　由概率守恒原理推导出短生产线的概率平衡方程，求得到短生产线的稳态概率，由此计算串行短生产线的性能指标：平均　生产率和平均在制品数量。在Ｍａｔｌａｂ平台下分析系统平均生产率与部分系统参数之间的关系。同时，建立短生产线的　Ｐｅｔｒｉ网模型，采用事件调度法对短生产线进行仿真，将数值结果与仿真结果进行对比分析，说明文中以机器为中心的分析　方法的有效性。　关键词：并行机器组；短生产线；稳态概率；生产性能；事件调度法　中图分类号：ＴＰ３９　文献标识码：Ａ　文章编号：１６７３—６２９Ｘ（２０１２）ｌ０—０１９４一ｏ４　Ａｎａｌｙｓｉｓ　ｏｆ　Ｓｈｏｒｔ－ｓｅｒｉａｌ　Ｐｒｏｄｕｃｔｉｏｎ　Ｌｉｎｅ　ｗｉｔｈ　Ｍｕｌｔｉｐｌｅ　Ｐａｒａｌｌｅｌ－ｍａｃｈｉｎｅ　Ｗｏｒｋｓｔａｔｉｏｎｓ　ＭＵ　Ｈｕａ－ｌｉｎｇ，ＷＡＮＧ　Ｚｈｅｎｇ　（Ｓｃｈｏｏｌ　ｏｆ　Ａｕｔｏｍａｔｉｏｎ，Ｓｏｕｔｈｅａｓｔ　Ｕｎｉｖｅｒｓｉｔｙ，Ｎａｎｊｉｎｇ　２１００９６，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｉｔ　ａｎａｌｙｚｅｓ　ｔｈｅ　ｐｅｒｆｏｒｍａｎｃｅ　ｏｆ　ａ　ｓｈｏｒｔ　ｓｅｒｉａｌ　ｐｒｏｄｕｃｔｉｏｎ　ｌｉｎｅ　ｗｉｔｈ　ａ　ｐａｒａｌｌｅｌ——ｍａｃｈｉｎｅ　ｗｏｒｋｓｔａｔｉｏｎ　ａｎｄ　ｉｔｓ　ｕｐｓｔｒｅａｍ　ａｎｄ　ｄｏｗｎ・－　ｓｔｒｅａｍ．Ｃｏｎｓｔｒｃｕｔｅｄ　ｉｔｓ　ｐｒｏｂａｂｉｌｉｔｙ　ｂａｌａｎｃｅ　ｅｑｕａｔｉｏｎ　ａｎｄ　ｅｖａｌｕａｔｅｄ　ｉｔｓ　ｐｅｒｆｏｒｍａｎｃｅｓ　ｓｕｃｈ　ａｓ　ｐｒｏｄｕｃｔｉｏｎ　ｒａｔｅ　ａｎｄ　ｗｏｒｋ—ｉｎ—ｐｒｏｃｅｓｓ　ｂｙ　ａ　ｎｕ—　ｍｅｒｉｃ　ｍｅｔｈｏｄ．Ａｎａｌｙｚｅ　ｔｈｅ　ｒｅｌａｔｉｏｎ　ｂｅｔｗｅｅｎ　ａｖｅｒａｇｅ　ｐｒｏｄｕｃｔｉｏｎ　ｒａｔｅ　ａｎｄ　ｐａｒｔｉｃｌｅ　ｓｙｓｔｅｍ　ｐａｒａｍｅｔｅｒｓ　ｕｎｄｅｒ　Ｍａｔｌａｂ　ｐｌａｔｆｏｒｍ．Ｍｅａｎｗｈｉｌｅ，ｅｓ－　ｔａｂｌｉｓｈ　ｔｈｅ　Ｐｅｔｒｉ　ｍｏｄｅｌ　ｆｏｒ　ｈｉｔｓ　ｓｈｏｒｔ　ｓｅｒｉａｌ　ｐｒｏｄｕｃｔｉｏｎ　ｌｉｎｅ　ａｎｄ　ｔｈｅ　ｓｈｏｒｔ　ｌｉｎｅ　ｉｓ　ｓｉｍｕｌａｔｅｄ　ｂｙ　ｅｖｅｎｔ　ｓｃｈｅｄｕｌｉｎｇ．Ｃｏｍｐａｒｉｎｇ　ｔｈｅ　ｎｕｍｅｒｉｃａｌ　ｒｅ—　ｓｕｌｔ　ａｎｄ　ｓｉｍｕｌａｔｉｏｎ　ｒｅｓｕｌｔ　ｓｈｏｗｓ　ｔｈａｔ　ｔｈｅ　ｐｒｏｐｏｓｅｄ　ｍｅｔｈｏｄ　ｉｓ　ｆｅａｓｉｂｌｅ．　Ｋｅｙ　ｗｏｒｄｓ：ｐａｒａｌｌｅｌ—ｍａｃｈｉｎｅ　ｗｏｒｋｓｔａｔｉｏｎ；ｓｈｏｒｔ　ｌｉｎｅ；ｓｔｅａｄｙ－ｓｔａｔｅ　ｐｒｏｂａｂｉｌｉｔｙ；ｐｒｏｄｕｃｔｉｏｎ　ｐｅｒｆｏｒｍａｎｃｅ；ｅｖｅｎｔ　ｓｃｈｅｄｕｌｉｎｇ　Ｕ　引　吾　相同。分解算法的一般步骤为：先将一个长度为Ⅳ的　长生产线分解成Ｎ一１个两机器的生产线；然后推导　出每个短生产线的方程，这些方程中含有短生产线的　在生产制造系统中，生产能力以及平均在制品数　量的分析对于生产线的设计、管理都起着至关重要的　作用　。对生产性能的分析可分为对可靠机器和不可　靠机器、有限缓冲区和无限缓冲区的分析。常见的生　产线类型有串行生产线、并行生产线、可重人生产线、　参数，如故障率、修复率、生产率等，并计算这些参数。　总体来说，这些近似的分析方法都以两机器单缓冲区　的分析为基础。但在实际的生产制造系统中，一个工　装配生产线、环形生产线　。　目前，已经解决了对于只含有两个机器一个缓冲　区的串行生产线的性能的精确分析　。对于多机器　串行生产线的分析目前大多采用分解／聚合方法来近　作站中通常包含多个机器（文中称其为并行机器组），　在以往的分析中，将每一个并行机器组等效成为一个　机器　，这样等效之后的分析方法就与传统的串行生　产线的分析方法相同。然而这种简化的分析方法，在　似分析，聚合算法主要是基于Ｍｅｅｒｋｏｖ的聚合方法　，　分解算法主要是基于Ｇｅｒｓｈｗｉｎ分解方法　。聚合方　法的主要思想是将一个两机器的生产线用一个等效机　来代替，这个等效机的生产能力和两机器的生产能力　收稿日期：２０１２－０２—１９：修回日期：２０１２－０５—２３　基金项目：国家自然科学基金资助项目（６０９７４０９６）　作者简介：穆华灵（１９９０一），女，硕士研究生，研究方向为计算机集成　制造；汪峥，副教授，研究方向为离散事件动态系统、计算机集成　全状态的状态转移图中近似忽略了正常工作状态之间　的状态转移，这样对于缓冲区容量比小的系统，带来的　近似误差比较大。　为了提高分析的精度，对于含有并行机器组的串　行生产线性能的分析，同样采取分解方法，并且考虑到　并行机器组对生产性能的影响，将把一个含有并行机　器组的串行生产线分解为含有两个缓冲区一个机器组　的短生产线，如图１所示。文中对于这样的生产单元　即两个缓冲区一个机器组的短生产线的性能进行精确　制造、制造系统分析与控制、产品开发过程分析与规划。　第１０期　穆华灵等：含有并行机器组的串行短生产线的性能分析　・１９５・　分析，并建立２个缓冲区１个工作站的Ｐｅｔｒｉ网模型，　进行仿真。　一（５）当一个工件从日　进入机器组被加工，同时有　个工件加工好后进入　，即（ｎ　，ｎ　，　）一＞（ｎ　一　１，ｎＤ＋１，　）ｏ　１短串行生产线的分析　根据概率平衡方程建立方程组并由此计算分解后　的短生产线中各个机器的参数，从而计算平均生产率　以及平均在制品数量ＩＶＩＰ。　为分析问题简便，先给出以下定义：　但这些事件的发生会受系统本来状态的限制，比　如当ｎ　＝０时，事件（２）不会发生，所以每个事件的发　生都有一定的约束条件，对应于上面五个事件的约束　条件分别为：　（１）ｎ　＝　，上游缓冲区中工件达到缓冲区最大　值时，工件无法再进入缓冲区；　（１）工作站中机器的数量为固定的　，每个机器都　是不可靠的，工件的平均加工时间ｌ　，机器单位时间　的故障率Ｐ，单位时间的维修率ｒ，工件平均到达时间　１／ｈ．，需求间隔时间１／ｄ均服从负指数分布；　（２）在短生产线中，只加工一种类型的工件；　（３）机器只在加工工件时才会发生故障，机器相　互独立的；　（４）上下游缓冲区大小分别为Ｊ７、ｒ　Ｎ　（５）系统的状态：（ｎ　ｎ。，　），其中ｎ　，＝０，１，…，　Ｎ　，ｎ。＝０，１，…，Ｎ。，　＝０，１，…，　，ｎ　表示上游缓冲区　中工件的数目，ｎ　下游缓冲区中工件的数目，　表示能　够工作的机器数目；　（６）缓冲区在存储和传递工件的过程中是完全可　靠的，不会出现故障。　由于有多台机器，因此每个机器组的生产率与等　待加工的零件数及机器数有关，即：　＝／ｘｍｉｎ｛／７，　，　｝。　图１为含有并行机器组的短串行生产线。　Ｉ　｝ｊＵ　Ｌ　ｌ　图１含有并行机器组的短串行生产线　１．１稳态概率平衡方程　定义Ｐ（　，／７，。，　）为系统状态（ｎ　，ｎ　，　）的稳态　概率。　●以下事件的发生会使系统离开（ｎ　ｎ　，ｘ）这个　状态：　（１）当有一个工件到达口　，即（ｎ　ｎ。，　）一＞（ｎ　＋１，ｎＤ，　）；　（２）当有一个工件离开Ｂ　，即（ｎ　ｎ。，　）一＞　（ｎｆ，，ｒ／，。一１，　）；　（３）当有一个正在工作的机器发　故障，即（ｎ　，　ｎＤ，　）一＞（ｎ“，ｎＤ，ｘ一１）；　（４）当有一个故障的机器恢复正常，即（ｎ　，ｎ。，　）——＞（ｎ“，ｎＤ，　＋１）；　（２）ｎ　：０，下游缓冲区中工件为０时，无法从中取　得工件；　（３）　＝０或ｎ　＝０或　’。＝Ｎ。，能够正常工作的机　器数量为０或机器处于被饥饿／阻塞时，机器不可能发　生故障；　（４）　＝　，所有机器都可正常工作时，没有机器可　被修复；　（５）ｎ　，＝０或ｎ。＝Ｎｏ，上游缓冲区为空或者下游缓　冲区为满时，机器组处于饥饿或阻塞状态，不能正常工　作。　●以下事件的发生会使系统进入（ｎ　，ｎ　，　）这个　状态（同时考虑到约束条件）：　（１）（ｎ　＋１，ｎ　一１，　）—　＞（ｎ　，ｎ　，　），系统中有　一个工件从上游缓冲区被加工的同时有一个加工好的　工件放到下游缓冲区，且当一　＝Ｎ　或ｎ。＝０时，不存在　状态（ｎｆ，＋１，ｎ，Ｊ一１，　）；　（２）（ｎ“一１，ｎＤ，　）一＞（ｎ　ｎ。，　），系统中有一　个工件到达日　，且当ｎ　，＝０时，不存在状态（ｎ　一１，ｎ　，　）；　（３）（ｎ　，ｎ。＋１，Ｘ）一＞（ｎ　ｎ。，　），系统中有一　个工件离开该系统，且当ｎ。＝Ｎｏ时。不存在状态（ｎ　，　ｎＤ＋１，　）；　（４）（ｎ“，ｎＤ，　一１）一＞（ｎｕ，ｎＤ，　），有一个机器　从故障状态下恢复正常工作，当　＝０时不存在状态　（ｎⅣ，ｎＤ，　一１）；　（５）（ｎｕ，ｎＤ，　＋１）一＞（ｎⅣ，ｎＤ，　），有一个机器　发生故障，当　＝Ｘ时，不存在状态（ｎ　ｎ　，　＋１），当ｎ　＝０或１７，。＝Ｎｏ时，机器被饥饿／阻塞，机器不能正常工　作，也就不会发生故障。　根据概率平衡方程的构造原理　：离开状态（ｎ　ｎ　，　）的概率等于进入状态（ｎ　ｎ。，　）的概率，所以得　到以下稳态概率平衡方程：　Ｐ（ｎｕ，ｎＤ，　）［Ａ，（ｎⅣ¨－　，　）＋ｄｌ（凡　＝０）＋ｒＩ（　：　）＋ｐｌ（　＝０　ＬＪ　＝０　ｕ凡Ｄ＝ＮＤ）＋ｂ￣ｍｉｎ（　，　ｕ），（凡　＝０　ｔ＿Ｊ　ｎＤ＝ＮＤ）］＝Ｐ（ｎｕ＋１，ｎＤ一１，　）／ｘｍｉｎ（　，凡　＋１）　・，（凡　＝ＮＵ　ｕ凡Ｄ＝０）＋Ｐ（ｎ　一１，ｎＤ，　）Ａ，（ｎｕ＝０）＋　Ｐ（凡　，　Ｄ＋１，　）ｄｌ（ｎＤ＝Ｎｏ）＋Ｐ（凡　，ｎＤ，　一Ｉ）ｒｌ（　：　・１９６・　计算机技术与发展　第２２卷　０）＋Ｐ（Ｉｔｕ，１１，ｏ，　＋１）ｐｌ（　＝Ｘ　ｕ　ｎｕ＝０　ｕ　ｎＤ＝ＮＤ）　（１）　模型包括某些变迁为立即变迁的情形¨　。　文中采用随机Ｐｅｔｒｉ网建模，具有一个并行机器组　和上下游缓冲区的短生产线Ｐｅｔｒｉ网模型如图２所示。　其中，Ｉｔｕ＝０，１，…，Ｎｕ，／９，Ｄ＝０，１，…，ＮＤ，　＝０，　ｌ，…，　，，（ｃｏｎｄｉｔｉｏｎｓ）＝Ｌ０ｃｏｎｄｉ‘ｉ。ｎｓ满足　，其中各符号的定义如表１所示。　由文中第二部分的假设可知，工作站中机器是否　可正常工作受上游缓冲　与下游缓冲区曰。的制约，　这一假设条件可以通过将抑制弧连接到机器的ｔ　来表　示。　表示库所，ｔ　表示变迁。　ｌ，其他概率Ｐ（ｎ　，ｎ　，　）满足如下的归一化条件：　∑∑∑Ｐ（ｎ　，　）＝１　（２）　由公式（１）和（２）组成的方程组可求出稳态概率　Ｐ（ｎ　，ｎ　，　），即可得到短串行生产线的平均生产率以　及平均在制品数量ＷＩＰ。　１．２短生产线平均生产率及平均在制品数量　尸　平均生产率分析：　从上游缓冲区的角度看，短生产线的稳态生产率　为：　刑　＝Ａ∑∑∑Ｐ（ｎ　ｎ　，　）　（３）　从机器的角度得到的短生产线稳态生产率为：　删　＝　∑∑∑ｍｉｎ（ｎ　）Ｐ（ｎ　，　，　）　（４）　从下游缓冲区的角度分析短生产线的稳态生产率　为：　删　＝∑∑∑ｍｉｎ（ｎ。，ｄ）Ｐ（Ｉｔ　，ｎ。，　）　（５）　当生产线处于平衡状态时，要满足ＴＨ　＝ＴＨ　＝　ＴＨｏ。　平均在制品数量ＷＩＰ分析：　上游缓冲区的平均在制品数量：　ＷＩＰ　＝∑∑∑ｒｇＵＰ（ｎ　，ｎ。，　）　（６）　下游缓冲区的平均在制品数量：　ＷＩＰ。＝∑∑∑１７＂０Ｐ（ｎ　，　）　（７）　短生产线的平均在制品数量：　Ｎ“　Ｎ¨ｘ　ＷＩＰ＝∑∑∑（ｎ＝０　』　＝０　＝０　　＋ｎ　）Ｐ（ｎ　，ｎ　，　）　（８）　２仿真实验　Ｐｅｔｒｉ网作为一种图形化的数学建模方法，其直观　性和严谨性可以描述系统和过程活动间复杂的逻辑关　系，所以成为离散事件动态系统有效的分析和建模工　具　。随机Ｐｅｔｒｉ网（Ｓｔｏｃｈａｓｔｉｃ　Ｐｅｔｒｉ　Ｎｅｔｓ）由Ｍｏｌｌｙ首　先提出的，他研究的是有界随机Ｐｅｔｒｉ网，随机变量服　从负指数分布，同时证明了Ｐｅｔｒｉ网状态可达图同构于　一个马氏链，借助马氏理论便可得到系统的相关性能　指标；Ｃｈｏｉ等人推广了Ｍｏｌｌｙ的工作，采用了一种广义　随机Ｐｅｔｒｉ网（Ｇｅｎｅｒａｌｉｚｅｄ　Ｓｔｏｃｈａｓｔｉｃ　Ｐｅｔｒｉ　Ｎｅｔｓ）模型，此　从仿真模型中可以看出，单工件单品种类型短生　产线的离散事件仿真中包含六种事件类型，分别是：上　游工件到达事件、机器开始加工事件、机器完成加工事　件、机器故障事件、机器修复完成事件、工件离开下游　缓冲区事件。　文中采用事件调度法对系统进行仿真，即用事件　的观点来分析真实系统，通过定义事件及每个事件发　生引起系统状态的变化，按时间顺序确定并执行每个　事件发生时有关的逻辑关系－一　。所有事件均放在事　件表中。事件调度法算法主要包括：初始化操作、事件　表操作、推进仿真时钟操作　。事件的选择与处理不　断地进行，直到仿真终止的条件产生为止。下一事件　发生的时间总是保存在事件表中，当从事件表中取出　最先发生时间的事件，此时将仿真钟推进到该事件发　生时间，转向该事件的处理子程序执行。　ｎ　］　ＮＤ　－Ｖ　Ｃ　ｔｌ　ｐ２　１３６　ｔ　上游缓冲区　Ｆ游缓冲区　图２　短生产线的Ｐｅｔｒｉ网模型　事件中的对象实体有两种，分别是：机器和工件。　每个事件的要素有：事件的发生时间、事件的主体以及　事件类型。　图２中各个符号的意义如下：　Ｐ．：上游缓冲区的输入库所（假设一直提供原料）；　Ｐ：：上游缓冲区容量为Ⅳ　；ｐ　：机器组工作；　：机器组中有，个机器空闲；　Ｐ　：下游缓冲区容量为Ｎ。；　Ｐ，：下游缓冲区的输出库所（假设容量为无限大）；　ｔ　：新工件到达生产系统；ｔ　：机器故障修复；　ｔ　：机器开始加工；ｔ　：机器发生故障；　ｔ　：工件加工完毕；ｔ　：工件离开系统。　第ｌＯ期　穆华灵等：含有并行机器组的串行短生产线的性能分析　・１９７・　３数值结果与仿真结果　３．１　Ｍａｔｌａｂ计算出的数值结果　（１）Ｎｕ＝４，ＮＤ＝３，Ｘ＝２，　＝１，Ｐ＝０．０１，ｒ＝０．１，　在制品数量的仿真结果记为　Ｐ　ｓ（ｉ），数值结果记为　Ｐｒ＿Ｔ（ｉ），平均在制品数量的平均误差记为ｅ＿ＷＩＰ。　（１）Ｎｕ＝４，ⅣＤ＝３，Ｘ＝２，ｕ：１，Ｐ＝０．Ｏｌ，ｒ＝０．１，　ｄ＝１，Ａ∈｛１，２，…，ｌ０｝，步长是１。机器端的生产率　ＴＨ　，上游缓冲区的生产率ＴＨ　下游机器端的生产率　ｄ＝１，Ａ∈｛１，２，…，ｌ０｝，数值分析结果与仿真实验结　果对比见表１。平均生产率与平均在制品数量的平均　误差分别为：　见图３，横坐标代表Ａ，纵坐标分别代表ＴＨ　，ＴＨ　，　ＴＨｎ　ｏ　至　Ｚ　１．－　擗　糌　｛Ｌ　｛Ｌ　ｇ　窖　罂　裴　髂　一＜　晕　勰　工件到达率ｖ　工件到达宰Ｖ　图３　Ａ∈｛１，２，…，ｌ０｝，步长为１－，对应的　ＴＨ　，ＴＨ”ＴＨ　Ｄ　（２）Ｎｏ：３，Ｘ＝２，　：１，Ｐ＝０．０１，ｒ：０．１，ｄ＝１，Ａ　＝１０，Ｎ　Ｅ｛１，２，…，ｌ０｝，步长为１，机器端的生产率　ＴＨ　，上游缓冲区的生产率ＴＨ，，，下游机器端的生产率　见图４，横坐标代表　，纵坐标分别代表ＴＨ　，　ｊ　王　ｏ　卜　Ｚ　卜　褂　糌　盔　ｇ　ｇ　释　释　一＜　羽　舞　簿　０　５　ｌｏ　ｏ　５　１０　上游缓冲区Ｎｕ　上游缓冲区Ｎｕ　图４　Ｎｕ∈｛１，２，…，ｌ０｝，步长为１对应的　ＴＨ　，ＴＨ￡，，ＴＨＤ　由图３和图４可知，删　＝ＴＨ　＝ＴＨ。。与分析的　结果一致。且系统的平均生产率是关于参数Ａ，Ⅳ，，的　单调递增有限函数。　３．２数值结果与仿真结果的对比　系统平均生产率的仿真结果记为ＴＨ＿Ｓ（ｉ），数值　结果记为ＴＨ＿Ｔ（ｉ），平均误差记为ｅ＿ＴＨ；系统的平均　１０　∑（ＴＨ＿Ｓ（　）一ＴＨ＿Ｔ（　））／ＴＨ＿Ｔ（ｉ）　ｅ删＝　——一　———一×　一１００％＝１．０４％　１０　∑（ＷＩＰ＿Ｓ（　）一一　　Ｐ－ｅＷＩＰ：兰。＿——一—————一　（ｆ））／ＷＩＰ＿Ｔ（　）　　—Ｘ】００％＝一２．０５％　表ｌ　到达率变化系统的数值值与仿真结果　（２）Ｎｄ＝３，Ｘ＝２，Ｈ＝１，Ｐ＝０．０１，ｒ＝０．１，ｄ＝１，Ａ　＝ｌ０，Ｎｕ＝｛１，２，…，ｌ０｝。理论分析结果与仿真实验　结果对比见表２。平均生产率与平均在制品数量的平　均误差分别为：　ｌ０　∑（ｒＨ＿Ｓ（　）一ｒＨ　ｒ（ｉ））／ＴＨ　Ｔ（　）　册＝旦——一　———一×　一１（）０％＝０．（）０６％　１０　ｅ唧＝旦∑（ｒ—ａ—ｍ＿ｓ（一　）一　尸＿，—　—，（—　））—／ｍ／一Ｐ＿ｒ（　）　　一×１００％＝０．３３％　４结束语　文中主要对含有一并行机器组及上下游缓冲区的　短串行生产线的性能进行数值分析与仿真，构建这种　短串行生产线的稳态概率平衡方程，及其Ｐｅｔｒｉ网模　型，根据Ｐｅｔｒｉ网模型进行离散事件仿真。由表１和表　２数据计算得到仿真结果的误差在允许误差范围之　内，由此验证了文中分析方法的可行性。　（下转第２０１页）　第ｌ０期　王瑛：ＧＩＳ在土地评价中的应用研究　国土资源遥感，２００８（３）：１０・・ｌ４．　・２Ｏｌ・　况进行监督施工，或者先行开挖保护，或者不允许开　发，总之依据情况必须事先征求文物行政部门的意见。　２）对大块地下文物埋藏的地区，国家政府可以组　织考古文物部门主动考察勘探，整批出土保护。　［４］　Ｚｈａｎｇ　Ｊ。Ｇｏｏｄｃｈｉｌｄ　Ｍ　Ｆ．Ｕｎｃｅｒｔａｉｎｔｙ　ｉｎ　Ｇｅｏｇｒａｐｈｉｃ　Ｉｆｕｏｒｍａ－　ｔｉｏｎ［Ｍ］．Ｎｅｗ　Ｙｏｒｋ：Ｔａｙｌｏｒ＆＿Ｆｒａｎｃｉｓ，２００２．　［５］涂超．ＧＩＳ在历史文化资源保护中的应用研究［Ｊ］．计算　平．基：于遥感和ＧＩＳ的青岛市土地利　机技术与发展，２００６，１６（７）：１６５—１６７．　３）苏州城市规划部门作古城规划的时候需要考　虑地下文物保护区的范围，有地下古文物需要保护的　地块给予保护性规划，比如古城区南端以及北端的靠　近护城河附近区域。　［６］　李瑞俊，孙希华，李用动态变化分析［Ｊ］．山东师范大学学报（自然科学版），　２００４，１９（４）：５０－５２．　［７］石常蕴，周慧珍．ＧＩＳ技术在土地质量评价中的应用一以苏　州市水田为例［Ｊ］．土壤学报，２００１，３８（３）：２４８—２５４．　［８］　闫保银，孙在宏．地理信息系统技术在土地评价中的应用一　３结束语　ＧＩＳ技术运用于各类土地评价研究是一种趋势，　结合城镇土地定级估价信息系统［Ｊ］．农机化研究，２００５　（２）：１９２－１９４．　先进的空间信息处理、统计、分析、归纳、决策技术给土　地评价以更具科学化、技术化，全面化的发展，为各项　［９］李永香，赵俊三，李洪玉，等．ＧＩＳ空间分析在土地定级中　的应用［Ｊ］．矿山测量，２００６（２）：１０－１２．　［１Ｏ］王红扬．８０年代以来苏州地区城镇土地演化特征与机制　［Ｊ］．地理科学，１９９９，１９（２）：１２８－１３４．　土地决策管理做出重要的参考性指导性服务¨　’”］。　本研究通过ＧＩＳ技术，建立时空模型，进行土地利用评　价，在土地评价的研究方式上作了部分研究应用。　参考文献：　［１］杨国清．土地评价的哲学思考［Ｊ］．广东工业大学学报（社　会科学版），２ｏｏ５，５（３）：９２—９５．　［２］Ｐａｎ　Ｊ　Ｊ，Ｚｈａｎｇ　Ｔ　Ｌ，Ｚｈａｏ　Ｑ　Ｇ．Ｄｙｎａｍｉｃｓ　ｏｆ　ｓｏｉｌ　ｅｒｏｓｉｏｎ　ｉｎ　Ｘｉｎｇｇｕｏ　ｃｏｕｎｔｙ，Ｃｈｉｎａ，ｄｅｔｅｒｍｉｎｅｄ　ｕｓｉｎｇ　ｒｅｍｏｔｅ　ｓｅｎｓｉｎｇ　ａｎｄ　［１１］乔伟峰，孙在宏．ＧＩＳ辅助下的城市土地集约利用潜力评　价方法研究一以江苏省苏州市为例［Ｊ］．国土资源科技管　理，２００４（１）：３４－３７．　［１２］黄小芳．ＧＩＳ在城市土地利用生态适宜性评价中的应用一　以上海市浦东新区为例［Ｊ］．科学技术与工程，２０１ｌ，ｌ１　（３１）：７８４１－７８４６．　［１３］Ｙａｎｇ　Ｈ　Ｄ，Ｈｕ　Ｙ　Ｍ，Ｄｅｎｇ　Ｆ　Ｑ，ｅｔ　ａ１．Ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ｉｍｍｕｎｅ　ＧＩＳ［Ｊ］．Ｐｅｄｏｓｐｈｅｒｅ，２００５，１５（３）：３５６—３６２．　ａｌｇｏｒｉｔｈｍ　ｔｏ　ｅｖａｌｕａｔｉｏｎ　ｏｆ　ｓｏｉｌ　ｒｅｓｏｕｒｃｅ　ｑｕａｌｉ￣［Ｊ］．Ｐｅｄｏ．　ｓｐｈｅｒｅ，２００５，１５（５）：６６０—６６８．　［３］陈华，孙丹峰．基于ＧＩＳ技术的土地评价研究进展［Ｊ］．　（上接第１９７页）　表２　Ｎｕ变化系统的数值结果与仿真结果　ｒｅｖｉｅｗ　ｏｆ　ｍｏｄｅｌｓ　ａｎｄ　ａｎａｌｙｔｉｃａｌ　ｒｅｓｕｌｔｓ［Ｊ］．Ｑｕｅｕｅｉｎｇ　Ｓｙｓ—　ｔｅｒｎｓ。１９９２，１２（１）：３－３４．　［４］“Ｊ，Ｂｌｕｍｅｎｆｅｌｄ　Ｄ　Ｅ，Ａｌｄｅｎ　Ｊ　Ｍ．Ｃｏｍｐａｉｒｓｏｎｓ　ｏｆ　ｔｗｏ—ｍａ．　ｃｈｉｎｅ　ｌｉｎｅ　ｍｏｄｅｌｓ　ｉｎ　ｔｈｒｏｕｇｈｐｕｔ　ａｎａｌｙｓｉｓ［Ｊ］．Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｊｏｕｒｎａｌ　ｏｆ　Ｐｒｏｄｕｃｔｉｏｎ　Ｒｅｓｅａｒｃｈ，２００６，４４（７）：１３７５—１３９８．　［５］　Ｍｅｅｒｋｏｖ　Ｓ　Ｍ，Ｚｈａｎｇ　Ｌ．Ｔｒａｎｓｉｅｎｔ　ｂｅｈａｖｉｏｒ　ｏｆ　ｅｒｓｉａｌ　ｐｒｏｄｕｃｔｉｏｎ　ｌｉｎｅｓ［Ｊ］．１ｉｅ　Ｔｒａｎｓａｃｔｉｏｎｓ，２００８，４（３）：２９７—３１２．　［６］　Ｇｅｒｓｈｗｉｎ　Ｓ　Ｂ．Ｍａｎｕｆａｃｔｕｉｒｎｇ　ｓｙ￣ｅｍｓ　ｅｎｇｉｎｅｅｒｉｎｇ［Ｍ］．Ｅｎ－　ｇｅｌｗｏｏｄ　Ｃｌｉｆｓ，ＮＪ：ＰＴＲ　Ｐｒｅｎｔｉｃｅ　Ｈａｌｌ，１９９４：５９—１７２．　［７］Ｂｕｒｍａｎ　Ｍ　Ｈ．Ｎｅｗ　ｒｅｓｕｌｔｓ　ｉｎ　ｆｌｏｗ　ｌｉｎｅ　ａｎａｌｙｓｉｓ［Ｄ］．Ｍａｓｓａ．　ｅｌｌｕｓｅｔｔｓ：ＭＩＴ，１９９５．　［８］梁之舜，邓集贤．概率论及数理统计［Ｍ］．北京：高等教育　出版社。１９８８．　参考文献：　［１］郑力，陈恳，张伯鹏．制造系统［Ｍ］．北京：清华大学出　版社，２００１：７７—１２０．　［９］　江志斌．Ｐｅ　网及其在制造系统建模与控制中的应用　［Ｍ］．ｊｔ京：机械工业出版社，２００４：９７一ｌ１７．　［１Ｏ］孟宪刚，严洪森．基于自学习模糊Ｐｅｔｒｉ网的知识化制造系　统采购预测［Ｊ］．控制与决策，２００９，２４（３）：３７１—３７６．　［Ｉ１］郑大钟，赵千Ｊｉｌ．离散事件动态系统［Ｍ］．北京：清华大学　出版社，２００１．　［２］Ｌｉ　Ｊｉｎｇｓｈａｎ，Ｂｌｕｍｅｎｆｅｌｄ　Ｅ，Ａｌｄｅｎ　Ｍ．ｅｔ　ａ１．Ｔｈｒｏｕｇｈｐｕｔ　ａｎａｌｙ・　ｓｉｓ　ｏｆ　ｐｒｏｄｕｃｔｉｏｎ　ｓｙｓｔｅｍｓ：ｒｅｃｅｎｔ　ａｄｖａｎｃｅｓ　ａｎｄ　ｆｕｔｕｒｅ　ｔｏｐｉｃｓ　［Ｊ］．Ｉｎｔｅｎａｔｒｉｏｎａｌ　Ｊｏｕｒｎａｌ　ｏｆ　Ｐｒｏｄｕｃｔｉｏｎ　Ｒｅｓｅａｒｃｈ，２００９．４７　（１４）：３８２３－３８５１．　［１２］刘勇，王德才，冯正超．离散事件系统仿真建模与仿真策　［３］Ｄａｌｌｅｒｙ　Ｙ，Ｇｅｒｓｈｗｉｎ　Ｓ　Ｂ．Ｍａｎｕｆａｃｔｕｒｉｎｇ　ｆｌｏｗ　ｌｉｎｅ　ｓｙ￣ｅｍｓ：ａ　略［Ｊ］．西南师范大学学报，２００５，３０（６）：１０１９—１０２５．