Vol.23 朱红青等：黄金分割法求解火源特性×间临界点 ·503· 表1 Fibonacci数 点H,,其位置在区间长度的0.618处： Table 1 Fibonacci table taba H=0.618×H,kJ/min (4) n0123456789101112 然后应用火灾模拟软件进行模拟计算，若关键 F.1233581321345589144233 分支风向反转，模拟区间变为[H。,H,];若关键 式(1)中的缩短率数列为： 分支风向不变，模拟区间为[H,H]. F- (3)针对子区间[H,H(或[H,,H]),确定 第2个试点H2,其位置在子区间长度的0.618处： F1 H=0.618×H,kJ/min (4) 然后应用火灾模拟软件进行模拟计算，若关键 F 分支风向反转，模拟区间变为[H,H]:若关键 图11个搜家试点的位置 分支风向不变，模拟区间为[H2,H]. Fig.1 The location of two searching points (4)按上述步骤反复进行，进一步缩少特性 会兵 F 区间，在n步后，其区间为[H,H。],区间长度 为Ln,ALn=H。-Hn1,若： 将上.述数列分为奇数项F-F和偶数项 Ln≤6 (6) Fw/Fw,因为这2个数列都收敛于同一极限.即： 则H,或Hn为火源特性区间的临界点.式(S)中 F业=μ im5=入imF (2) F4 6为预设的精度，其值一般定为10、 (S)根据临界点H,或Hn,结合定性分析法 且 =4=50.618 2 (3) 确定相应的控风措施，形成控风方案 所以，以不变的区间缩短率0.618，代替Fb- onacci搜索法每次不同的缩短率，即得到了黄 3结论 金分割法（或为0.618法）.对于某一矿井或高层 建筑物，若某一位置发生火灾，首先根据可燃物 通过理论分析和以往的现场应用实践证 的种类、数量、分布范围，估计火源燃烧强度的 明，与用一端逼近搜索法搜索火源特性区间临 大小，其火源强度范围的大小可结合相应的火 界点相比，黄金分割法的最大优点是每进行一 灾试验来推算.对于某一具体位置火源，假设根 次搜索，其区间范围减少0.618倍或0.382倍，大 据推算，估计火源燃烧释放的最大热量为H, 大提高了搜索效率，相应地减少了全风网动态 kJ/min,释放的最小热量为H,kJ/min,则火源特 模拟的次数，对于减少作量和火灾事故救灾 性区间为[H,H],搜索区间也为[H。,H],则利 很有意义. 用黄金分割法搜索临界点的步骤如下， 参考文献 (1)按H.和以、利用火灾非稳态模拟软件模 1钱颂迪.运筹学.北京：清华大学出版社，北京：1999.149 拟，若某些关键分支风流流向同向，则在搜索区 2 Zhou Xinquan,Zhu Hongging.The Application of the 间[H。,H]内不存在特性区间临界点：若反向则 Qualitative and Quantitative Synthetical Analysis Method 在搜索区间[H.,H,]内存在一特性区间临界点. to Decision Making of Mine Fire Fighting.in:Mining Sci- (2)在[H.,H]区间内按照黄金分割法取试 ence and Technology '99.1999.83 Searching the Critical Point of Fire Characteric Zone by Golden Partition Method ZHU Hongqing.ZHOU Xinquan.WU Yangyiu China University of Mining and Technology-Beijing.100083.Chian ABSTRACT Because of fire intensity changing,the critical point exists within fire characteric zone when using the qualitative and quantitative synthetical analysis method and fire characteric zone method to deter- mine the airflow control scheme.The golden partition method is applied to search the critical point,enhancing the searching efficiency. KEY WORDS golden partition method;critical point;fire characteric zone.匕 朱 红 青等 黄金 分割 法 求 解 火源特 性 区 间 临 界点 衰 数 地 皿 式 中的缩短 率数列 为 凡 凡 一 九 一 七一仁二付二州 ” ‘’ 圈 个搜索试 点的位 柱 一 ， 凡 一 凡 ’ ， ” 凡 ’ 将 王 述 数 列 分 为 奇 数 项 凡 一 ， 凡 和 偶 数 项 凡阴小 ，，因为这 个数列都收敛 于 同一极 限 ， 即 点 ， ， 其位置 在 区 间长度 的 处 从 然后 应 用 火 灾 模拟 软件进行 模拟 计算 ， 若关键 分支风 向反 转 ， 模拟 区 间 变 为〔 ， 」若关键 分支 风 向不 变 ， 模拟 区 间 为〔 ， ， 从 针 对 子 区 间 〔 ， 从」或〔 ， ， 〕 ， 确定 第 个 试点从 ， 其位 置在 子 区 间 长 度 的 处 ， 然 后 应 用 火灾 模拟 软件进 行模拟计算 ， 若关键 分支风 向反 转 ， 模拟区 间 变 为 ， 从 若关键 分支风 向不变 ， 模拟 区 间 为〔从 ， 从 〕 按 匕述 步骤 反 复进 行 ， 进 一 步缩 少 特性 区 间 ， 在 步 后 ， 其 区 间 为 〔从 ， 从 ， ， 区 间 长度 为 。 ， 几 ， 从 一 从 一 ， 若 ， ‘ 占 则从 或从 一 为火 源 特性 区 间 的临 界点 式 中 咨为预 设 的精度 ， 其值一 般定 为 一 ’ 根 据临界点从 或从 一 ，， 结合定性分析法 确定 相 应 的控 风措施 ， 形 成控 风 方案 、 ， 、 少产、夕 才 ‘，︼、 ‘龟 ， ‘、了 且 又二 产 以 卿合 一 ， 粤 、 所 以 ， 以不 变 的区 间缩短 率 ， 代 替 搜索法 每次不 同的缩 短 率 ， 即得 到 了 黄 金 分割法 或 为 法 对 于 某一 矿 井或 高层 建筑物 ， 若某一 位置发生 火 灾 ， 首先根据 可燃 物 的 种类 、 数量 、 分布范 围 ， 估计 火源燃烧强 度 的 大小 ， 其火源强 度范围的大小可结合相应 的火 灾试验来推算 对于某一具体位 置火 源 ， 假设根 据推算 ， 估计火源燃烧 释放 的最大 热量 为 ， 曰 ， 释放 的最 小热量 为 从 ， ， 则火 源 特 性 区 间 为 ， 从〕 ， 搜索区 间也 为〔 ， 从〕 ， 则利 用 黄 金 分割 法搜 索临 界 点 的 步 骤 如 下 按 和 从 利用 火 灾 非 稳 态模 拟 软件模 拟 ， 若某些 关键分支风流 流 向 同向 ， 则在搜索区 间 〔 ， 从 内不 存在 特性 区 间 临 界点 若反 向则 在搜索区 间 〔 ， 从 内存 在 一 特性 区 间临 界点 在 ， 从飞区 间 内按 照 黄金 分割法 取试 结论 通 过 理 论 分 析 和 以 往 的 现 场 应 用 实 践 证 明 ， 与用 一 端逼 近搜索法搜索火源 特性 区 间临 界点相 比 ， 黄金 分割法 的最 大 优点 是 每进 行 一 次搜 索 其 区 间 范 围减 少 倍或 倍 ， 大 大提 高了搜索效率 ， 相 应 地 减 少 了 全风 网动 态 模拟 的 次 数 ， 对 于 减 少 ’ 作量 和 火 灾 事故 救灾 很 有 意 义 参 考 文 献 钱颂迪 运筹学 ， 北 京 清华大学 出版社 ， 北京 ， 介 ’ ” ， 体 舒 ” 一 ， ，