Optimal Control: Theory and Application to Science, Engineering, and Social SciencesView this Special Issue
Optimal Control of Renewable Resources Based on the Effective Utilization Rate
The effective utilization rate of exploited renewable resources affects the final total revenue and the further exploitation of renewable resources. Considering the effective utilization rate, we propose an optimal control model for the exploitation of the renewable resources in this study. Firstly, we can prove that the novel model is nonsingular compared with the singular basic model. Secondly, we solve the novel model and obtain the optimal solution by Bang-Bang theory. Furthermore, we can determine the optimal total resources and the maximal total revenue. Finally, a numerical example is provided to verify the obtained theoretical results.
Renewable resources (such as fisheries resources) are considered to be “inexhaustible” at all times, but excessive exploitation will actually exhaust them. Thus, it is useful to study the reasonable exploitation of renewable resources and their effective utilization to obtain the maximum revenue.
Various optimum exploitation schemes for renewable resources have been studied in recent decades. In a pioneering study, Gordon [1, 2] proposed the Gordon-Schaefer bioeconomic model by introducing the concept of economic efficiency and cost management. Although this model included a large number of unrealistic assumptions, it exhibited a certain degree of concordance with the histories of empirical fisheries . Subsequently, further economic control models of renewable resources [4–10] were proposed based on the work of Gordon. Clark and Munro  extended the static version of the fisheries economics model to nonautonomous and nonlinear cases. Indeed, the models proposed by Gordon and other researchers were all static until the 1980s when Clark  established a dynamic bioeconomic model based on the Gordon-Schaefer biological model. In addition to the classical model of Clark, another study  provided a nonlinear optimal control bioeconomic model that used the variation in the fishing effort rate as the control. In the 1990s, Defeo and Seijo  developed a bioeconomic model using a yield-mortality model. Later, based on the Beverton-Holt age structure model , Beverton and Holt  constructed a dynamic bioeconomic model that considered interactions among populations, which was not based on the Gordon-Schaefer bioeconomic model.
In recent years, various bioeconomic models of fisheries resources have also been proposed. Smith  established a Bayesian bioeconomic dynamic model by introducing the Bayesian statistical method. Domíguez-Torreiro and Surís-Regueiro  introduced game theory into a bioeconomic model of fishery resources and proposed a management strategy that addressed fishery issues. Das et al.  proposed a new bioeconomic model that combined a predator-prey ecological model with marine environmental factors. Other studies [20–22] introduced the fundamentals of control parameterization, which is a popular numerical technique for solving optimal control problems. At the same time, a switched autonomous system was proposed to formulate a fed-batch culture process where the switching instants between the feed and batch processes were used as control variables, which is similar to the method proposed in the present study. The Food and Agriculture Organization also produced a series of bioeconomic models [23, 24] to provide a theoretical basis for policy to facilitate the sustainable use of fisheries resources.
However, all of these bioeconomic models and optimal exploitation schemes did not consider the effective utilization rate, especially for renewable resources. Indeed, if the exploited renewable resources are not utilized fully, they will be wasted, but they may also pollute the environment. Given the impact on the expected revenue of the effective utilization rate for renewable resources, we introduce the concept of effective utilization into a renewable resources development model and we propose an optimal control model to ensure that it approximates the actual situation.
In Section 2, we introduce the basic bioeconomic model and propose an optimal control bioeconomic model based on the effective utilization rate. In Section 3, the singularity of this novel model is analyzed based on the maximal principle. Using the Bang-Bang theory, we obtain the optimal exploitation scheme and the optimal total renewable resources subject to the maximum total revenue. In Section 4, a numerical example is given that verifies the results.
2. Optimal Control Model Based on the Effective Utilization Rate
2.1. The Basic Model
It is well known that renewable resources have their own life cycles, even if they are not exploited or consumed. In general, we expect that the natural growth rate is greater than the natural mortality rate for renewable resources. However, the maximum amount of resources cannot exceed the environmental carrying capacity. In general, it is considered that the growth of resources satisfies the following logistic equation :where denotes the resource biomass, is the carrying capacity of the ecosystem, and is the intrinsic growth rate of resources. If we let be the exploitation amount, then model (1) becomesLet the total revenue from the exploited resources in unit time be expressed aswhere is the revenue per unit, is the cost per unit, and is the instantaneous social rate of discount. In this case, the objective function of the total revenue can be stated as follows:
Thus, the optimal control problem involving the basic model of renewable resources can be expressed as follows:
2.2. The Proposed Model Based on the Effective Utilization Rate
Let be the effective utilization rate at time ; then should satisfy the following three assumptions.(A1)The effective utilization rate will increase gradually with respect to (i.e., with the development of technology); that is, .(A2)The increase in will become more difficult after it reaches a certain level; that is, must satisfy .(A3)The ideal or the best utilization of resources is achieved completely; that is, .By these assumptions, we take as our effective utilization rate, which satisfies the preceding assumptions. Furthermore, let the initial effective utilization rate be and the ultimate effective utilization rate be . Then, parameter , in can be obtained as follows:
Thus, our improved objective function for optimal control associated with the effective utilization rate can be written asand our proposed optimal control model can be established as follows:where represents the management objective, is the maximum amount of exploitation, and the initial and the terminal populations of renewable resources , are assumed to be known. In model (8), the meanings of , , , , are similar to those in model (5).
3. Solutions to the Model
In this section, we first analyze the singularity of the proposed model (8) and we then apply the Bang-Bang approach to obtain its solution.
3.1. Solution of the Hamiltonian Formulation
The Hamiltonian formulation to the optimal control problem (8) can be written as follows:Then, we have
By the maximal principle, the Hamiltonian will obtain the maximal value with respect to if the objective function obtains the maximal value. The linear relationship between the Hamiltonian function and control changes the computation of the maximization problem (10) into a Bang-Bang optimal control problem, as follows:
However, from the solution of (11), we cannot obtain any information about the optimal control by . Indeed, there are two possible solutions for . Firstly, if there is only a countable time set that satisfies , we can still use the Bang-Bang theory to solve this model. Secondly, if there is an interval that satisfies for , the problem will become more complicated and difficult to solve. Thus, to obtain a better solution to the model, we first need to analyze whether such an interval exists or not. Thus, we need to discuss the singularity [25, 26] of the optimal control model (8).
3.2. Singularity Analysis of the Model
Based on the preceding discussion, we can see that the singularity of model (8) will affect the choice of method used to solve this model. In this subsection, we discuss the singularity of model (8) by reduction to absurdity.
Assuming that the optimal control is a singularity, that is, there is an interval for all that satisfies the derivative of with respect to can be expressed as by the costate equation of model (8)we obtain the following equation:By substituting (12) into (15), we obtain Equation (16) admits the positive root given by which indicates that if we want (16) to hold at time , then must satisfy (17). By taking the derivative of (17) with respect to , we obtain
In fact, (18) does not equal (2). In other words, that satisfies the condition given above does not exist. Compared with the basic model (5), which is a singularity, our proposed model (8) involving the effective utilization rate is a normal model. Therefore, the optimal control of model (8) can be solved using the Bang-Bang approach.
3.3. Existence of the Switching Time
In this subsection, we discuss the existence of the switching time . First, we consider , which indicates that the inequalitywill hold. In this case, the optimal strategy is by the Bang-Bang method.
However, this is not a reasonable method for exploiting renewable resources because the renewable resources will be destroyed. Therefore, we make the following assumption: . Then, (14) is obtained if the condition that is satisfied and is monotonically decreasing for . In addition, we discuss the monotonicity of the following function:By taking the derivative of (20) with respect to Then, It can be seen that if and only if satisfies the following condition:
However, this condition cannot be obtained from (22). Thus, (20) is monotonically decreasing in time period , which indicates that there will be an appropriate intersection between and ; that is, the switching time exists. Then we can obtain that
3.4. Solution to Model (8)
After obtaining the switching time and proving the nonsingularity of model (8), the optimal control strategy to (8) can be expressed by the Bang-Bang approach as follows:Next, we find the expressions of and the optimal renewable resources function . First, if we let and , we obtain By solving the ordinary differential equation (26), we obtainwhere is expressed as follows: Second, if we let and , we can obtain Unlike (26), we cannot solve the ordinary differential equation (29) directly using the Bernoulli equation. However, if there are suitable parameters , , the following formula holds; that is,Then, we can solve (29) by the Bernoulli equation. In fact, by (29) and (30), we havethat is, the following quadratic equation holds:Irrespective of whether (32) is positive or 0, that is, the suitable parameters , will be obtained, which indicates that (29) can be solved by the Bernoulli equation. It is known that is the environmental carrying capacity, which is a very large number; thus the formula holds. Then, we can obtainTherefore, we can solve (31), where the parameters , are given as follows: Moreover, we can obtain the solution of (29) by the Bernoulli equation:That is, where is expressed asTo summarize, in the interval , the optimal renewable resources function can be expressed asThe total amount of resources is continuous with respect to time , even with ; that is,We can determine the time by (40). Thus, we obtain the optimal exploitation and the optimal level of total resources under the maximum revenue:
In summary, the optimal management strategy can be interpreted as follows. At the beginning of a time period , we should not exploit the renewable resources to reach a certain number. At time , we can exploit the renewable resources at the maximum capacity. Finally, we obtain the total benefit as follows:
4. Numerical Example
Based on the solution to the proposed model (8), we now present an example that verifies the model. The initial parameters for the fisheries resources are specified in Table 1, that is, the carrying capacity of the ecosystem , the intrinsic growth rate of the renewable resources , the revenue per unit , the cost per unit , the instantaneous social rate of discount , the maximum amount of exploitation , the time interval , the ultimate amount of renewable resources , the initial effective utilization rate , and the ultimate effective utilization rate . Thus, by (6), (28), (35), and (38), we can obtain the remaining parameters: , , , , , and , respectively. By substituting parameters into (38), the switching time is obtained (note that in another example for tons based on the same calculation used for ). Finally, we obtain the expressions for optimal exploitation and the optimal level of total resources under the maximum revenue as follows (we omit the case where ):which are shown in Figures 1 and 2, respectively. The total maximum revenue for the cases where and , respectively, can also be obtained from (42) as follows:
In Figure 2, for the case where , the switching time is , which indicates that we should allow the renewable resources to grow naturally for about 4 months of 1 year. Subsequently, the renewable resources can be exploited optimally where and the final total maximum revenue is 136,254 dollars. Figure 2 also shows that the renewable resources retain approximately the same growth rate as before. However, when the optimal exploitation is with the development of exploitation technology, we should delay the exploitation time; that is, , which is about 10 months of 1 year. In this case, the total maximum revenue is about 138,753 dollars, which is greater than that for and the renewable resources continue to grow, but at a slower growth rate.
In this study, we considered a model for the optimal control of renewable resources, which is affected by the effective utilization rate. We analyzed the singularity of the model. If the model is a singularity, the problem is difficult because the maximal principle cannot identify an optimal candidate explicitly. Indeed, we proved that our proposed optimal control model involving the effective utilization rate is normal and we solved the model with the aid of the maximal principle. Finally, we determined the optimal exploitation rate and the optimal level of the total renewable resources under the maximum revenue. This study may provide reference values to facilitate the development of renewable resources. However, it is clear that is not related directly to according to the present study. Thus, further research is required to address the case where and have a linear relationship.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors are grateful to the reviewers for many helpful suggestions which improved the presentation of the paper. The authors acknowledge supports by the National Science Foundation of China (no. 61174209) and the Oriented Award Foundation for Science and Technological Innovation, Inner Mongolia Autonomous Region, China (2012).
H. S. Gordon, “An economic approach to the optimum utilization of fishery resources,” Fisheries Research, vol. 10, pp. 442–447, 1953.View at: Google Scholar
H. S. Gordon, “The economic theory of a common-property resource: the fishery,” The Journal of Political Economy, vol. 62, no. 2, pp. 124–142, 1954.View at: Google Scholar
C. W. Clark, Mathematical Bio-Economics: The Optimal Management of Renewable Resource, John Wiley & Sons, New York, NY, USA, 1976.
C. W. Clark and G. P. Kirkwood, “Bio-economic model of the Gulf of Carpentaria prawn fishery,” Journal of the Fisheries Research Board of Canada, vol. 36, pp. 1303–1312, 1979.View at: Google Scholar
B. S. Goh, “The usefulness of optimal control theory to biological problem,” in Theoretical Systems Ecology; Advances and Case Studies, E. Halfon, Ed., pp. 385–399, Academic Press, New York, NY, USA, 1979.View at: Google Scholar
C. W. Clark, Bio-Economic Modeling and Fisheries Management, John Wiley & Sons, New York, NY, USA, 1985.
J. C. Seijo, O. Defeo, and S.Salas, “Fisheries bio-economics: theory, modeling and management,” FAO Fisheries Technical Papers 368, 1998.View at: Google Scholar
R. J. H. Beverton and S. J. Holt, “On the dynamics of exploited fish populations,” Fishery Investigations, vol. 2, no. 19, pp. 1–533, 1957.View at: Google Scholar
M. D. Smith, Bayesian Bio-Economics of Marine Reserves, American Agricultural Economics Association, 2005.
T. Das, R. N. Mukherjee, and K. S. Chaudhuri, “Harvesting of a prey-predator fishery in the presence of toxicity,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 33, no. 5, pp. 2282–2292, 2009.View at: Publisher Site | Google Scholar | MathSciNet
P. J. Sqarre, “Bio-economic analytical model no. 5,” in Proceedings of the 13th Eafe Annual Conference, Thessaloniki, Greece, March 2005.View at: Google Scholar
O. Flaaten, Fisheries Economics and Management, University of Tromsø, Tromsø, Norway, 2011.
C. W. Clark, Mathematical Bio-Economics: The Optimal Management of Renewal Resources, John Wiley & Sons, New York, NY, USA, 1976.