Recently, optimal investment problem for insurers has attracted more and more attention. For example, Browne  considered the optimal investment problem for an insurance company with diffusion risk model and the stock price was described by a geometric Brownian motion. Hipp and Plum  used the classical Cramér-Lundberg risk model and obtained the optimal investment strategy for an insurance company. For other literatures on this field, please refer to     . Although many researchers study optimal investment problem for general insurers, there are few literatures studying investment problem for concrete life insurance companies.
In this paper, we consider the optimal investment problem for a life insurance company. Life insurance is a contract between an insurance policy holder and an insurer. The insurers promise to pay a benefit in exchange for a premium, upon the death of a policyholder. Depending on the contract, other events such as terminal illness or critical illness can also trigger payment. At the same time as the rapid economic development, the problem of unbalanced development among different regions emerges. This disparity in development results in differences in the health level and life expectancy of residents between different regions. Thus, insurance companies should consider these differences between different regions when pricing life insurance products. We use the mortality of critical diseases to describe the differences between the different regions. Biagini et al.  studied the pricing and hedging of a life insurance portfolio with dependent mortality risk. They introduced the Gaussian random fields to describe the mortality intensities and obtained the analytically results of the hedging strategy. According to the analysis of the mortality of critical diseases in different regions, we introduce the concept of health-levels for different regions. We then analyze the impact of health-level on the amounts of claims and adjust premiums for life insurance in various regions.
Suppose the life insurance company choose the mean-variance criterion. There are two main approaches to solve the mean-variance problem. One is to obtain the precommitment strategy and the other is to study the time-consistent strategy. Bi and Guo  derived the optimal precommitment strategy for insurers under mean-variance criterion for various cases. Li et al.  considered the time-consistent investment and reinsurance strategies for an insurer under Hestons stochastic volatility model. For more detailed discussion, see    . In this paper, we consider the time-consistent strategy for this optimal investment problem.
This paper according to the regional differences of the health level studies the optimal investment strategy for a life insurance company. The surplus process of the company is described by the classical Cramér-Lundberg Model and the insurer can invest in a risk-free asset and a risky asset. To maximize the profits and minimize the risk, we take the mean-variance criterion into account. Furthermore, we establish the corresponding Hamilton-Jacobi-Bellmen(HJB) equations and obtain the time-consistent optimal investment strategy. Finally, we study the effect of the health-level on the value function and the extra premium by numerical stimulations.
This paper is organized as follows: Section 2 formulates the model. In Section 3, we derive the optimal investment strategy and the expectation of the terminal wealth explicitly. Section 4 gives the numerical stimulations. Section 5 concludes the paper.
2 Model and Optimization Problem
In this section, let be a complete probability space, where is a fixed time horizon; and is a probability measure; and is the information of the financial and insurance market until time t.
2.1. Surplus Process for a Life Insurance Company
In this section, we will model the heath-level by using the mortality of the critical diseases and describe the surplus process of the life insurance company.
Firstly, the health-level is a defined by the mortality of the critical diseases in different regions. Suppose that there are k regions and M critical diseases. So the health-level of the kth region is defined as:
where is the mortality of the ith critical diseases in the kth region; is the minimum mortality of the ith critical diseases in all regions.
And the surplus process of the life insurance company is described by the classical Cramér-Lundberg model:
where is the initial capital of the life insurance company, y is the health-level of the location of the company, and is the premium rate of the company, is a compound Poisson process which represents the cumulative claims up to time t, is a homogeneous Poisson process with intensity and are i.i.d. positive random variables with mean and finite variance . The extra claims caused by the different health-level is denoted by .
Assume that is described by
where is the safety loading of the insurer, and is the extra premium associated with health-level.
and we can find that is a quadratic function according to the practical data (see Appendix). To simplify the calculation, is assumed to be the following linear function
where a and c are positive constants satisfying . The relationship of constants guarantee that the extra claim is larger than the extra premium .
Substituting Equations (3) and (4) into (2) to get a new form of surplus process:
Differentiating (5) with respect to t, we get
According to the approximation of Cramér-Lundberg model,
where is a Brownian motion on the space .
2.2. Financial Model
This section will model the financial assets and the wealth process of the insurance company.
In this paper, the financial market consists of a risk-free asset and a risk asset. The price of the risk-free (i.e., cash) asset is the following:
where is the initial price of risk-free asset and r is the risk-free interest rate.
The second asset in the market is risk asset which is described by a standard Brownian motion
where is the expect return rate of the risk asset; is the volatility of this asset. is a standard Brownian motion on the space . And the correlation coefficient between and is . Suppose that there are no transaction costs and trading is continuous.
Moreover, donate as the money which investment in the risk asset , and , which means that a short sell of the bonds is not permitted.
Then the wealth process of the life insurance company is
Finally, we defined the optimization problem for the continuous-time model (9). We want to maximize the fund size and to minimize the volatility of the terminal wealth. So we choose the mean-variance utility as our main criterion. And the optimization problem under this criterion can be described as follow:
where is the investment strategy of the insurance company.
3. Solution of the Optimal Control Problem
In this section, we will find out the optimal solution of the problem (10). Using the methods in  , the mean-variance optimal control problem is equal to the following Markovian time inconsistent stochastic optimal control problem:
where is a coefficient representing the degree of risk aversion of the insurance company. And also helps establish the optimal strategy of mean-variance optimal control problem. And the optimal investment strategy satisfies . Denote
and from (11), the value function is
Theorem 1 [Verification Theorem] If there exist three real functions satisfying the following extended Hamilton-Jacobi-Bellman equtions:
then there exist is the optimal strategy of problem (11), and , , .
Proof The proof of this theorem is similar to the proof in He and Liang  , so we omit the details here. +
After giving the theorem, we will solve the HJB Equations ((14), (16), (17)). Firstly, we establish the optimal strategy of problem (11).
From (13), we get
substituting (18) into (15),
Taking (19) and (18) into (14) and differentiating (14) with respect to , we can obtain
from this equation, we have
Substituting (21) into (14) and (16),
Suppose and have the following form:
Differentiating (24) with respect to x and t, we get
and taking (25) into (22) and (23), then the equations become
Let the coefficient of x and constant term of (26) and (27) be 0, we obtain
And then we get the solution of these ordinary differential equations, the results are as follows:
Substituting (30) and (25) into (21), we have
And the optimal investment strategy is
And we also have
where is the unique solution of the following equation:
4. Numerical Simulation
In this section, we will study the effects of model parameters on the effect of value functions and the extra premium . Throughout the numerical simulation, the initial parameters are given in Table 1 unless otherwise stated.
4.1. Effects of the Health-Level on the Extra Premium
Figure 1 shows the trends of extra premiums and extra claims. The straight line shows the trend of the extra claim and the other line shows the trend of
Table 1. Value of parameters in this section.
Where r is the risk-free interest rate; is the mean of random claims and is the variance of ; is the expect return rate of the risk asset and is the volatility of this asset; is the risk aversion of the insurance company; is the safety loading of the insurer; is the correlation coefficient between and ; is the intensity of the Poisson process ; T is the terminal time of our investment; a, b and c are constants which are used to describe the expression of the extra claim and the extra premium. (Modify according to comment 1)
Figure 1. Effects of health-level y on extra premium k(y).
the extra premuim . We find that the extra claims increase with y. This is consistent with intuition. Larger y means lower health-level, thus the extra claims increase as the health-level declines. Moreover, the extra premium increases with y, too. This shows that when the health level of a region is high, the company can appropriately reduce its premium. While the insurance company will charge more for the policyholder from regions with low health-level. And because the insurance company should afford a part of the risk of the health-level y, the extra claim is a little more than the extra premium from beginning to end. (Modify according to comment 2)
4.2. Effects of Parameters on the Value Function
Figure 2 shows the sensitivity analyzes of the value function of the insurance company. We find that the value functions increase with y first and it will decrease when the value of y is large enough. As y increases, the health-level decreases which means the extra claim is increasing and thus the value function becomes smaller accordingly. We also see that the value function is a increasing function of time t. So the investment strategy is efficient. From this figure, we know that if the health-level y changes over time t, it may be difficult to estimate the trend of value function. So the investment strategy can be useless in this situation. Thus we should use different expressions of extra premium in different regions because of the difference of the health-level between the different regions. The different extra premium can make the value function keep nondecreasing with y and t. (Modify according to comment 3 and 4) Figure 2(b) plots the effect of the initial wealth x on the value function. We see that the value function of the insurance company increases with x. This is consistent with intuition. The more the initial wealth is, the larger the value function is. And the effect of y on the value function is similar to Figure 2(a). The initial wealth is a constant during the period of investment, so it can not mislead the trend of value function with the change of health-level y. (Modify according to comment 3)
Figure 2. Effects of parameters t, x and y on value function . (a) Parameters t and y; (b) Parameters x and y.
From the form of optimal strategy , we find that the health-level y doesn’t influence the optimal strategy. So we do not study the effects of parameters on optimal strategy.
In this paper, we study the optimal investment strategy for a life insurance company with considering the differences in the health-level between different regions. The health-level of a region is defined by the incidence rate of the critical diseases. This paper first defines the expression of the health-level y and fits the extra claims according to the actual data. (Modify according to comment 5). The surplus process of the insurance company is described by Cramér-Lundberg model and the insurer is allowed to invest in a risk-free asset and a risky asset. The price process of the risky asset follows the Brownian motion and the insurer considers the mean-variance criterion. By the dynamic programming approach, we establish the HJB equations and derive the optimal investment strategy explicitly. Finally, numerical simulation is provided to analyze the effects of health-level and other parameters on the extra premium and value function. We find that the value function decreases as the corresponding health-level declines. Thus the life insurance company should consider a better expression of the extra premium to reduce the loss from the difference of the health-level in different regions and make the company can operate in a long term. (Modify according to comment 4 and 5). From Figure 2(a), we can find that the value function may decrease with t and y. So the extra premium is not enough to hedge the health-level risk. Thus we need a better definition of extra premium. For future research, we will try to find another way to calculate the extra premium . (Modify according to comment 6).
This research was supported by grants from the National Natural Science Foundation of China (11771329, 11871052, 11301376).
Fitting of the Extra Claims f(y)
We use the data from the URL: www.cdc.gov, www.acli.com and www.census.gov to fit the extra claim . Some data are given in Tables A1-A3.
To fit the extra claims, we give the following assumptions.
・ The difference of claims between different regions is only caused by different health-levels.
・ The beginning time is 2000.
・ The cumulate claims of different regions are the same.
Define the total claims with considering the health-level as follow:
Table A1. The per capita insurance claims of four states from 2002 to 2006.
Table A2. The critical diseases’ incidence of four states from 2002 to 2006.
Table A3. The health-level of four states from 2002 to 2006.
Figure A1. The fit curve of the extra claims f(y).
According to the real data and (42), we provide the fitting of in Figure A1.
From Figure A1, we see that is a quadratic curve approximately. Thus, we assume in (4).
 Browne, S. (1995) Optimal Investment Policies for a Firm with a Random Risk Process: Exponential Utility and Minimizing the Probability of Ruin. Mathematics of Operations Research, 20, 937-958.
 Zhou, M., Yuen, K.C. and Yin, C. (2017) Optimal Investment and Premium Control in a Nonlinear Diffusion Model. Acta Mathematicae Applicatae Sinica, 33, 945-958.
 Wang, Y., Rong, X. and Zhao, H. (2017) Optimal Investment Strategies for an Insurer and a Reinsurer with a Jump Diffusion Risk Process under the CEV Model. Journal of Computational and Applied Mathematics, 328, 414-431.
 Biagini, F., Botero, C. and Schreiber, I. (2011) Risk-Minimization for Life Insurance Liabilities with Dependent Mortality Risk. Mathematics and Financial Economics, 27, 505-533.
 Li, Z., Zeng, Y. and Lai, Y. (2012) Optimal Time-Consistent Investment and Reinsurance Strategies for Insurers under Heston’s SV Model. Insurance Mathematics and Economics, 51, 191-203. ttps://doi.org/10.1016/j.insmatheco.2011.09.002
 Shen, Y. and Zeng, Y. (2015) Optimal Investment-Reinsurance Strategy for Mean-Variance Insurers with Square-Root Factor Process. Insurance Mathematics and Economics, 62, 118-137.
 Zeng, Y. and Li, Z. (2017) Optimal Time-Consistent Investment and Reinsurance Policies for Mean-Variance Insurers. Insurance Mathematics and Economics, 49, 145-154.
 Li, D., Shen, Y. and Zeng, Y. (2018) Dynamic Derivative-Based Investment Strategy for Mean-Variance Asset-Liability Management with Stochastic Volatility. Insurance Mathematics and Economics, 78, 72-86.
 He, L. and Liang, Z. (2013) Optimal Investment Strategy for the DC Plan with the Return of Premiums Clauses in a Mean-Variance Framework. Insurance Mathematics and Economics, 53, 643-649.