Due to the highly developed natural fractures and horizontal beddings in shale reservoirs

it is possible to generate complex volume fractures during hydraulic fracture treatment. In order to investigate the propagation process of complex fracture network

a three dimensional finite element model for volume fracture propagation of hydraulic fracture in shale reservoirs is built using fluid-solid coupling basic equations in porous medium and basic theory of damage mechanics. The numerical simulation results are in good accordance with laboratory experiment of shale fracture propagation

which prove the reliability of the numerical model. After a series of numerical simulations

some conclusions are drawn as follws: (1) The horizontal beddings can open and form horizontal fractures in hydraulic fracturing. They can intersect with vertical fractures and finally generate complex volume fracture network. (2) As the horizontal stress difference increases

the length of stimulated reservoir volume (SRV) can increase and the width can decrease

which means the ratio of length to width of SRV can increase. The length is the distribution distance of volume fracture along horizontal maximum principal insitu stress. The width is the distribution distance of volume fracture along horizontal minimum principal insitu stress. (3) As the pump rate of fracturing increases

the SRV length can decrease

width can increase

and the ratio of length to width of SRV can decrease accordingly. (4) When the residual tensile strength of the natural fracture becomes stronger

the width of SRV can decrease

the length and the ratio of length to width of SRV can increase. Research results can offer some references for hydraulic fracture design and operation of shale gas in China.