Understanding the characteristics of hydraulic fracture,porous flow and heat transfer in fractured rock is critical for geothermal power generation applications,and numerical simulation can provide a powerful approach for systematically and thoroughly investigating these problems. In this paper,we present a fully coupled solid-fluid code using discrete element method(DEM) and lattice Boltzmann method(LBM). The DEM with bonded particles is used to model the deformation and fracture in solid,while the LBM is used to model the fluid flow. The two methods are two-way coupled,i.e.,the solid part provides a moving boundary condition and transfers momentum to fluid,while the fluid exerts a dragging force to the solid. Two widely used open source codes,the ESy S_Particle and the OpenL B,are integrated into one code and paralleled with Message Passing Interface(MPI) library. Some preliminary 2D simulations,including particles moving in a fluid and hydraulic fracturing induced by injection of fluid into a borehole,are carried out to validate the integrated code. The preliminary results indicate that the new code is capable of reproducing the basic features of hydraulic fracture and thus offers a promising tool for multiscale simulation of porous flow and heat transfer in fractured rock. Understanding the characteristics of hydraulic fracture, porous flow and heat transfer in fractured rock is critical for geothermal power generation applications, and numerical simulation can provide a powerful approach for systematically and thoroughly investigating these problems. In this paper, we present a fully coupled solid- fluid code using discrete element method (DEM) and lattice Boltzmann method (LBM). The DEM with bonded particles is used to model the deformation and fracture in solid, while the LBM is used to model the fluid flow. The two methods are two- way coupled, ie, the solid part provides a moving boundary condition and transfers momentum to fluid, while the fluid exerts a dragging force to the solid. Two widely used open source codes, the ESy S_Particle and the OpenL B, are integrated into one code and paralleled with Message Passing Interface (MPI) library. Some preliminary 2D simulations, including particles moving in a fluid and hydraulic fracturing induced by injection of fl The preliminary results that that new code is capable of reproducing the basic features of hydraulic fractures and therefore offers a promising tool for multiscale simulation of porous flow and heat transfer in fractured rock .