An order parameter functional method is proposed which permits description of structural order-disorder phase transitions. A feature peculiar to the method is that the probabilities of all atomic groups (clusters) are represented in terms of a distribution function that is immediately related to the long-range order parameters. Consideration is given to applications of the proposed method to describe atomic ordering in non-stoichiometric niobium and tantalum carbides.