A Dynamic Programming Approach to Determining Optimal Forest Wildfire Initial Attack Responses
A Dynamic Programming Approach to Determining Optimal Forest Wildfire Initial Attack Responses
复制标题
确定最佳森林野火初始攻击响应的动态规划方法
DOI:
--
复制
发表时间:
2003
期刊:
影响因子:
--
通讯作者:
M. Wiitala
中科院分区:
文献类型:
--
作者:
M. Wiitala
A mathematical optimization model, based on the operations research technique of deterministic dynamic programming, is offered as a method to search quickly through available options to find the economically efficient set of initial attack resources to suppress a wildfire. Considerations in selecting initial attack resources for an efficient initial response include cost of transportation and use, line construction productivity, response times, resource complementaries, size of fire upon discovery, rate of fire spread, fire damage, mop-up cost, and fire benefits. The dispatch optimization model has several applications, such as developing pre planned and real-time dispatches and evaluating the economic efficiency of new initial attack technologies. For many years a quick and strong initial response was deemed paramount to successfully suppressing a forest fire. To achieve this objective, the criterion of closest forces was used to select the suppression resources deemed necessary for the initial attack effort. Much of the early supporting research focused on automating the process of finding closest forces (Mees 1978). Little attention was given to the cost of the initial attack response and the economic trade-offs that could be made with costs and resources losses associated with fire size. Over the past two decades the cost of maintaining an initial attack organization of sufficient size to make quick and strong initial attack responses increased steadily and significantly. Facing ever tightening budgets, dispatchers became more interested in balancing the cost of the initial suppression response against the benefit of reducing fire size. In response, researchers began to explore formal methods to help dispatchers select the most efficient initial attack suppression response. Even before the concern with economic efficiency, Parks (1964) presented an analytic solution for finding the economically optimal suppression effort. When the efficiency issue became more pressing, Parlar and Vickson (1982) revisited Parks' model. They offered an alternative solution using optimization techniques from control theory. Both models focused on the most efficient size and timing of the general suppression effort. Neither model considered the importance to the dispatch decision of differences among individual suppression resources in response times, line building rates, and other significant fire fighting traits. Because these considerations were and remain important aspects of the dispatch decision, neither model became operational. In subsequent research developments, resource differences important to the dispatch decision were considered. A foray in this area was initiated by Wiitala (1986) in pioneering the use of dynamic programming to find cost-effective dispatches. A similar approach was taken by Kourtz (1989) for dispatching water bombers and for delivering crews by helicopter. Although Wiitala (1986) attempted to account for all initial attack suppression costs, the dynamic programming algorithms by Kourtz (1989) minimized only transportation cost. Kourtz (1989) did not formally consider resource loss or other types of suppression related costs nor the economics of the strength of the attack. An abbreviated version of th i s paper was p resen ted at the Symposium on Fire Economics, Planning, and Policy: Bottom Lines, April 5-9, 1999, San Diego, California. Operations Research Analyst, Pacific Southwest Research Station, Forest Service, U.S. Department of Agriculture, 1221 SW Yamhill St., Suite 200, Portland, Oregon 97205. e -mai l : mwi i ta la / r6pnw_ Portland@fs.fed.us USDA Forest Service Gen. Tech. Rep. PSW-GTR-173. 1999. 115 Session III Optimal Forest Wildfire Responses---Wiitala In determining an appropriate suppression response, the efficiency-minded dispatcher considers all costs and losses. A balance must be struck between the cost of the suppression effort and those costs and losses associated with fire size. However, finding the cost efficient dispatch is a task made difficult by the availability of numerous initial attack resources dispersed over many locations that exhibit a wide range of characteristics. For each available suppression resource the dispatcher must know its location, transportation cost, line building cost, line construction rate, and response time. Also important are line building interactions between resources and any constraints affecting the intensity and duration of a resource's performance. Even with a modest 30 different resources from which to choose a dispatch, the number of combinations that can feasibly contain a fire can easily range into the millions. With so may choices, identifying the most efficient combination to dispatch to a fire is nearly impossible without assistance of computers and operations research techniques. This paper presents an operations research technique that can help the efficiency-minded dispatcher quickly identify an appropriate initial attack response. This objective is accomplished in two steps. The first step sets forth a general mathematical formulation of the wildfire dispatch optimization problem. The second step reformulates the general mathematical optimization problem to permit using the technique of deterministic dynamic programming to achieve an economically efficient dispatch. Model Formulation With many resources available to respond to a fire, the number of possible dispatches may be very large. Every dispatch will have its own fireline building trajectory, containment time, and cost profile. This profile will depend on the types, arrival times, and production rates of dispatched resources. Assuming resources will be sent immediately and will build fireline until the fire is contained, the cost of containment, J, can be formally stated as: J = ∑ x i (C0i + C1i (t − bi )) + C2(a(t)) + C4(t) [1]