Dynamic simulation for tempo-spatial distribution of strain energy density in inhomogeneous stratified crust

Through simulating the research on dynamic variations of strain energy density (SED) in seismogenic model with hard inclusion, the authors have gained further knowledge to such problems as the process of earthquake preparation, initial rupture, conditions of the initial rupture and fracture propagation direction, etc. Results of the research show that SED (strain energy density) in soft inclusion is very high during the initial period of earthquake preparation. And the increment of SED in the soft area decreases at the later stage of the process. Meanwhile, the increment increases quickly in hard inclusion and in the intersection zone of the inclusion with an erecting fault, where the increment of SED is maximum. Thus, the intersection zone between hard inclusion with larger elastic modulus and erecting fault becomes the place where the initial rupture or earthquake occurs. The fracture in the end part of the hard inclusion spreads along a direction nearly vertical to the erecting fault, so the theoretical fracture direction is consistent with that calculated by digital simulation.