Mixed wettability shale pore network flow simulation method

The invention relates to a mixed wettability shale pore network flow simulation method, and belongs to the field of shale digital core and multiphase flow simulation. In order to overcome the defects in the prior art, the invention aims to provide a mixed wettability shale pore network flow simulation method, which comprises the following steps: on the basis of the characteristics of different lithofacies in a real shale scanning electron microscope image, combining a four-parameter random growth algorithm and a pore growth kernel classification method to construct five shale lithofacies and five shale pores in sequence; obtaining a mixed wettability shale pore network through a maximum ball algorithm; and carrying out flow simulation to obtain a capillary force curve and a relative permeability curve. The method is easy and convenient to operate, can be used for mixed wettability shale pore network flow simulation, can also be used for research such as stress sensitivity, mercury injection simulation, electr