Formation interval determination method of MPD based on risk aversion and casing level optimization